Menggunakan paket ramalan yang sangat baik oleh Rob Hyndman, saya menemukan perlunya tidak hanya memiliki interval prediksi, tetapi untuk mensimulasikan sejumlah jalur masa depan, mengingat pengamatan masa lalu dari serangkaian waktu dengan musiman kompleks. Ada sesuatu untuk deret waktu yang...
Saya memiliki serangkaian tanggapan yang dinyatakan sebagai interval seperti contoh di bawah ini. > head(left) [1] 860 516 430 1118 860 602 > head(right) [1] 946 602 516 1204 946 688 di mana kiri adalah batas bawah dan kanan adalah batas atas dari respons. Saya ingin memperkirakan...
Saya mencoba menghitung distribusi posterior
Dalam peramalan deret waktu menggunakan berbagai model seperti AR, MA, ARMA, dll, kami biasanya fokus pada pemodelan data dalam perubahan waktu. Tetapi ketika kita memiliki 2 seri waktu yang menunjukkan koefisien korelasi Pearson mereka sangat berkorelasi, apakah mungkin untuk memodelkan nilai...
Saya perlu menghitung CDF Dirichlet , tetapi saya hanya dapat menemukan implementasi dari PDF . Apakah Anda tahu ada perpustakaan (terutama di R) yang
Saya memiliki beberapa set data EEG yang saya uji terhadap dua kelas. Saya bisa mendapatkan tingkat kesalahan yang layak dari LDA (distribusi kelas-kondisional bukan Gaussian, tetapi memiliki ekor yang sama dan pemisahan yang cukup baik), jadi saya ingin memplot ROC dari prediktor LDA terhadap set...
Ini adalah pertanyaan pemula tentang latihan dalam “Bayesian Computation with R” karya Jim Albert. Perhatikan bahwa sementara ini mungkin pekerjaan rumah, dalam kasus saya tidak, karena saya belajar metode Bayesian di R karena saya pikir saya mungkin menggunakannya dalam analisis masa depan...
Saya menulis pertanyaan ini dengan merujuk pada contoh pada hal 138-142 dari dokumen berikut: ftp://ftp.software.ibm.com/software/analytics/spss/documentation/amos/20.0/id/Manuals/IBM_SPSS_Amos_User_Guide.pdf . Berikut adalah gambar ilustrasi dan tabel: Saya mengerti bahwa variabel laten tidak...
Ketika saya menjalankan kode ini: require(nlme) a <- matrix(c(1,3,5,7,4,5,6,4,7,8,9)) b <- matrix(c(3,5,6,2,4,6,7,8,7,8,9)) res <- lm(a ~ b) print(summary(res)) res_gls <- gls(a ~ b) print(summary(res_gls)) Saya mendapatkan koefisien yang sama dan signifikansi statistik yang...
Saya ingin melakukan analisis kekuatan untuk uji hipotesis (non-) kesetaraan proporsi di mana proporsinya sangat kecil. Saya ingin melakukannya tanpa menggunakan perkiraan normal (atau Poisson) dari distribusi binomial. Ada beberapa jenis pertanyaan kekuatan yang ingin saya atasi. Pasca-hoc :...
Tolong izinkan saya untuk mengajukan pertanyaan dasar. Saya memahami mekanisme Naive Bayes untuk variabel diskrit, dan dapat mengulang perhitungan "dengan tangan". (kode HouseVotes84 semua jalan per di bawah ini). Namun - Saya berjuang untuk melihat bagaimana mekanisme bekerja untuk variabel...
Saya punya parameter θθ\theta yang terletak di antara [ 0 , 1 ][0,1][0,1]. Katakanlah saya dapat menjalankan percobaan dan memperolehnya θ^= θ + wθ^=θ+w\hat{\theta} = \theta + wdimana wwwadalah Gaussian standar. Yang saya butuhkan adalah perkiraanθθ\thetayang 1) tidak bias 2) hampir pasti dibatasi....
Saya mencoba menjalankan regresi Cox pada dataset 2.000.000 sampel baris sebagai berikut hanya menggunakan R. Ini adalah terjemahan langsung dari PHREG di SAS. Sampel tersebut mewakili struktur dataset asli. ## library(survival) ### Replace 100000 by 2,000,000 test <-
Saya sedang mengerjakan model yang mengandalkan fungsi parametrized jelek yang berfungsi sebagai fungsi kalibrasi pada bagian model. Menggunakan pengaturan Bayesian, saya perlu mendapatkan prior non-informatif untuk parameter yang menggambarkan fungsi saya. Saya tahu bahwa idealnya, saya harus...
Saya memiliki 2 penelitian yang mengamati respons pasien terhadap obat yang sama. Studi 1 menemukan 10.000 gen yang diekspresikan di atas latar belakang dan 500 di antaranya diekspresikan secara berbeda dan disebut sebagai tanda respon obat. Studi 2 menemukan 1.000 gen yang mewakili tanda tangan...
Saya memerlukan satu set data imputed (misalnya untuk membuat dummy grup negara dari data pendapatan per kapita negara imputed). R menawarkan paket paket untuk membuat beberapa data imputasi (misalnya Amelia) dan menggabungkan hasil dari beberapa dataset (seperti dalam MItools). Kekhawatiran saya...
Saya ingin memasukkan bobot sampel ke model regresi kuantil saya, tetapi saya tidak yakin bagaimana melakukan ini. Saya sudah menentukan berat badan saya, yang merupakan replikasi bobot yang sudah diberikan dalam dataset survei (dihitung dalam paket survei):
Saya sendiri akan selalu menggunakan mean geometrik untuk memperkirakan median lognormal. Namun, di dunia industri, kadang-kadang menggunakan median sampel memberikan hasil yang lebih baik. Dengan demikian pertanyaannya adalah, adakah rentang cutoff / titik mulai dari mana median sampel dapat...
Adakah yang tahu metode untuk membandingkan dua variabel dengan uji chi square jika variabelnya dari survei yang berbeda dengan svydesign()pernyataan yang berbeda ? Saya mencari untuk menguji perbedaan dalam distribusi variabel di dua gelombang survei, tetapi svychisq()pernyataan tersebut terbatas...
Saya telah melakukan pengelompokan titik koordinat (bujur, lintang) dan menemukan hasil yang mengejutkan dan merugikan dari kriteria pengelompokan untuk jumlah cluster yang optimal. Kriteria diambil dari clusterCrit()paket. Poin-poin yang saya coba klaster pada plot (karakteristik geografis dari...