Saya memiliki data berikut yang terletak di sini . Saya mencoba menghitung interval kepercayaan 95% pada kemurnian rata-rata ketika persentase hidrokarbon adalah 1,0. Di R, saya masukkan yang berikut ini. > predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95) fit lwr...
Saya menganalisis set data tentang komunitas intertidal. Data berupa persen tutupan (dari rumput laut, teritip, kerang, dll) dalam kuadrat. Saya terbiasa memikirkan analisis korespondensi (CA) dalam hal jumlah spesies , dan analisis komponen utama (PCA) sebagai sesuatu yang lebih berguna untuk tren...
Saya telah melihat dalam literatur yang diterbitkan (dan diposting di sini) bahwa efisiensi relatif asimptotik dari uji peringkat Wilcoxon menandatangani setidaknya 0,864 bila dibandingkan dengan uji t. Saya juga pernah mendengar bahwa ini hanya berlaku untuk sampel besar, meskipun beberapa buku...
Apakah mungkin untuk memplot plot dependensi parsial untuk menampilkan probabilitas kelas dan memperkirakan efek prediktor untuk model GBM ? Sesuatu yang mirip dengan partialPlotdari randomForestpaket. Menurut artikel ini , sebagian plot dapat dilakukan dengan gbm. Terima kasih sebelumnya atas...
Dalam Dixon, Coles ( 1997 ), mereka telah menggunakan estimasi kemungkinan maksimum untuk dua model Poisson independen yang dimodifikasi dalam (4.3) untuk memodelkan skor dalam sepak bola. Saya mencoba menggunakan R untuk "mereproduksi" alfa dan beta serta parameter efek rumah (hal 274, Tabel 4)...
Saya tertarik menemukan prosedur untuk mensimulasikan data yang konsisten dengan model mediasi yang ditentukan. Menurut kerangka model persamaan struktural linear umum untuk menguji model mediasi pertama kali diuraikan oleh Barron dan Kenny (1986) dan dijelaskan di tempat lain seperti Judd,...
Saya sedang mengerjakan estimasi vektor auto-regression (VAR) dan impulse response function (IRFs) berdasarkan data panel dengan 33 orang lebih dari 77 kuartal. Bagaimana seharusnya jenis situasi ini dianalisis? Algoritma apa yang ada untuk tujuan ini? Saya lebih suka melakukan analisis ini dalam...
Saya mencoba memperkirakan efek 2 obat ( drug1, drug2) pada kemungkinan pasien jatuh ( event). Para pasien dapat terjatuh lebih dari satu kali dan dapat memakai atau melepas obat kapan saja. Pertanyaan saya adalah bagaimana data harus disusun berkaitan dengan periode waktu (hari), khususnya...
Saya memiliki model regresi logistik (cocok melalui glmnet di R dengan regularisasi jaring elastis), dan saya ingin memaksimalkan perbedaan antara positif benar dan positif palsu. Untuk melakukan ini, prosedur berikut muncul dalam pikiran: Sesuai dengan model regresi logistik standar Dengan...
Saya ingin menilai keakuratan tes normalitas pada ukuran sampel yang berbeda dalam R (Saya menyadari bahwa tes normalitas dapat menyesatkan ). Sebagai contoh, untuk melihat tes Shapiro-Wilk, saya sedang melakukan simulasi berikut (dan juga merencanakan hasilnya) dan berharap bahwa ketika ukuran...
Saya baru dalam bidang statistik dan saya mencoba memahami perbedaan antara ANOVA dan regresi linier. Saya menggunakan R untuk mengeksplorasi ini. Saya membaca berbagai artikel tentang mengapa ANOVA dan regresi berbeda tetapi masih sama dan bagaimana dapat divisualisasikan dll. Saya pikir saya...
Indeks stabilitas populasi mengukur perubahan distribusi variabel dengan membandingkan sampel data dalam dua periode waktu. Ini sangat umum digunakan untuk mengukur pergeseran skor. Ini dihitung sebagai berikut: 1) Sampel dari periode dasar diskritisasi. Biasanya dipartisi ke dalam desil 2)...
Seorang kolega saya mengirimi saya masalah ini yang tampaknya membuat putaran di internet: If $3 = 18, 4 = 32, 5 = 50, 6 = 72, 7 = 98$, Then, $10 =$ ? Jawabannya sepertinya 200. 3*6 4*8 5*10 6*12 7*14 8*16 9*18 10*20=200 Ketika saya melakukan regresi linier di R: data <-...
Saya melakukan regresi linier menggunakan fungsi Rm: x = log(errors) plot(x,y) lm.result = lm(formula = y ~ x) abline(lm.result, col="blue") # showing the "fit" in blue tapi itu tidak pas. Sayangnya saya tidak bisa memahami manualnya. Bisakah seseorang mengarahkan saya ke arah yang benar agar...
Saya melihat lembar excel yang mengklaim menghitung χ2χ2\chi^2 , tetapi saya tidak mengenali cara melakukannya, dan saya bertanya-tanya apakah saya kehilangan sesuatu. Berikut adalah data yang dianalisis: +------------------+----------+----------+ | Total Population | Observed | Expected...
Saya ingin menyesuaikan DLM dengan koefisien yang bervariasi waktu, yaitu perpanjangan untuk regresi linier biasa, .yt=θ1+θ2x2yt=θ1+θ2x2y_t = \theta_1 + \theta_2x_2 Saya memiliki prediktor ( ) dan variabel respon ( y t ), marine & tangkapan ikan pedalaman tahunan masing-masing dari 1950 -...
Saya memiliki banyak rangkaian waktu dengan periode: hari, minggu atau bulan. Dengan stl()fungsi atau dengan loess(x ~ y)saya dapat melihat bagaimana tren seri waktu tertentu terlihat. Saya perlu mendeteksi apakah tren deret waktu meningkat atau menurun. Bagaimana saya bisa mengaturnya? Saya...
Buku Elemen Pembelajaran Statistik (tersedia dalam PDF online) membahas bias optimisim (7.21, halaman 229). Ini menyatakan bahwa bias optimisme adalah perbedaan antara kesalahan pelatihan dan kesalahan dalam sampel (kesalahan diamati jika kita sampel nilai-nilai hasil baru di masing-masing poin...
Saya sedang mengerjakan contoh-contoh dalam Analisis Data Does Bayesian Kruschke , khususnya ANOVA eksponensial Poisson di bab. 22, yang ia sajikan sebagai alternatif untuk uji chi-square independensi untuk tabel kontingensi. Saya dapat melihat bagaimana kita mendapatkan informasi tentang...
R in Action (Kabacoff, 2011) menyarankan rutinitas berikut untuk menguji penyebaran berlebihan dalam regresi logistik: Fit regresi logistik menggunakan distribusi binomial: model_binom <- glm(Species=="versicolor" ~ Sepal.Width, family=binomial(), data=iris) Fit regresi logistik menggunakan...