Pertanyaan kemarin Menentukan akurasi model yang memperkirakan probabilitas acara membuat saya ingin tahu tentang penilaian probabilitas. The Brier skor adalah rata-rata kuadrat error ukuran. Apakah pengukuran kinerja kesalahan absolut rata-rata analog punya nama
Saya sangat baru di parsial kuadrat terkecil (PLS) dan saya mencoba memahami output dari fungsi R plsr()dalam plspaket. Mari kita mensimulasikan data dan menjalankan PLS: library(pls) n <- 50 x1 <- rnorm(n); xx1 <- scale(x1) x2 <- rnorm(n); xx2 <- scale(x2) y <- x1 + x2 +...
Saya menjalankan tes root unit berikut (Dickey-Fuller) pada serangkaian waktu menggunakan ur.df()fungsi dalam urcapaket. Perintahnya adalah: summary(ur.df(d.Aus, type = "drift", 6)) Outputnya adalah: ############################################### # Augmented Dickey-Fuller Test Unit Root Test...
Di R, jika saya set.seed (), dan kemudian menggunakan fungsi sampel untuk mengacak daftar, dapatkah saya menjamin saya tidak akan menghasilkan permutasi yang sama? yaitu... set.seed(25) limit <- 3 myindex <- seq(0,limit) for (x in seq(1,factorial(limit))) { permutations <-...
Saya memiliki sampel data yang dihasilkan Roleh rnorm(50,0,1), jadi data jelas mengambil distribusi normal. Namun, Rtidak "tahu" informasi distribusi ini tentang data. Apakah ada metode Ryang dapat memperkirakan jenis distribusi sampel saya berasal? Jika tidak, saya akan menggunakan...
Saya telah menetapkan pertanyaan berikut sebagai bagian dari pekerjaan rumah: Rancang dan laksanakan studi simulasi untuk menguji kinerja bootstrap untuk memperoleh interval kepercayaan 95% pada rata-rata sampel data yang univariat. Implementasi Anda bisa dalam R atau SAS. Aspek kinerja yang...
Saya memiliki banyak set data terkait. Korelasi pearson antara pasangan mereka biasanya lebih besar dari korelasi spearman. Itu menunjukkan korelasi apa pun yang linear, tetapi orang mungkin berharap bahwa bahkan jika pearson dan spearman itu sama. Apa artinya ketika ada kesenjangan yang pasti...
Saya mencoba menganalisis data tindakan berulang dan berusaha untuk membuatnya bekerja R. Data saya pada dasarnya adalah sebagai berikut, saya memiliki dua kelompok perlakuan. Setiap subjek dalam setiap kelompok diuji setiap hari dan diberi skor (persentase benar pada suatu tes). Data dalam format...
Mengapa bahkan memiliki prior non-informatif? Mereka tidak memberikan informasi tentang . Jadi mengapa menggunakannya? Mengapa tidak hanya menggunakan prior informatif? Misalnya, misalkan . Lalu apakah adalah non-informatif sebelum ?θ ∈ [ 0 , 1 ] θ ∼ U ( 0 , 1 ) θθθ\thetaθ ∈ [ 0 , 1 ]θ∈[0,1] \theta...
Dalam posting sebelumnya saya bertanya-tanya bagaimana cara menangani skor EQ-5D . Baru-baru ini saya menemukan regresi kuantil logistik yang disarankan oleh Bottai dan McKeown yang memperkenalkan cara yang elegan untuk menangani hasil yang terbatas. Rumusnya
Catatan: pertanyaan ini adalah repost, karena pertanyaan saya sebelumnya harus dihapus karena alasan hukum. Sambil membandingkan PROC CAMPURAN dari SAS dengan fungsi lmedari nlmepaket di R, saya menemukan beberapa perbedaan yang agak membingungkan. Lebih khusus lagi, derajat kebebasan dalam...
Saya telah menggunakan rjags untuk menjalankan MCMC pada model, yang ditentukan dalam bahasa JAGS. Apakah ada cara yang baik untuk mengekstrak model itu dan melakukan prediksi dengannya (menggunakan distribusi posterior parameter saya)? Saya dapat menentukan kembali model dalam R dan menyambungkan...
Diberikan dataset: x <-
Saya mencoba menguji berbagai pendekatan analisis data fungsional. Idealnya, saya ingin menguji panel pendekatan yang saya miliki pada data fungsional yang disimulasikan. Saya sudah mencoba untuk menghasilkan FD simulasi menggunakan pendekatan berdasarkan pada penjumlahan Gaussian noise (kode di...
Saya mengalami masalah kinerja menggunakan MCMCglmmpaket di R untuk menjalankan model efek campuran. Kode ini terlihat seperti ini: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn,...
Saya memiliki fungsi probabilitas berikut: Prob=11+e−zProb=11+e−z\text{Prob} = \frac{1}{1 + e^{-z}} dimana z=B0+B1X1+⋯+BnXn.z=B0+B1X1+⋯+BnXn.z = B_0 + B_1X_1 + \dots + B_nX_n. Model saya terlihat seperti Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])\Pr(Y=1) =...
Apakah lebih baik untuk membedakan rangkaian (dengan asumsi perlu) sebelum menggunakan Arima ATAU lebih baik menggunakan parameter d dalam Arima? Saya terkejut betapa berbedanya nilai-nilai yang dipasang tergantung pada rute mana yang diambil dengan model dan data yang sama. Atau apakah saya...
Para peneliti sering menggunakan dua ukuran yang memiliki item yang sangat mirip dan berpendapat bahwa mereka mengukur hal-hal yang berbeda (misalnya, "Saya selalu khawatir ketika saya berada di sekitar mobil"; "Saya takut mobil"). Mari kita memanggil langkah-langkah hipotetis Ketakutan Mobil...
Katakanlah saya memiliki data berikut dan saya menjalankan model regresi: df=data.frame(income=c(5,3,47,8,6,5), won=c(0,0,1,1,1,0), age=c(18,18,23,50,19,39), home=c(0,0,1,0,0,1)) Di satu sisi, saya menjalankan model linier untuk memprediksi pendapatan: md1 = lm(income ~ age + home + home,...
Saya menggunakan paket randomForest dalam R dan menggunakan data iris, hutan acak yang dihasilkan adalah klasifikasi tetapi ketika saya menggunakan dataset dengan sekitar 700 fitur (fitur-fiturnya adalah masing-masing piksel dalam gambar 28x28 piksel) dan kolom label dinamai label, yang...