Saya memiliki kumpulan data deret waktu yang saya coba paskan dengan Hidden Markov Model (HMM) untuk memperkirakan jumlah status laten dalam data. Kode pseudo saya untuk melakukan ini adalah sebagai berikut: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
Saya menyesuaikan model efek acak dengan glmerbeberapa data bisnis. Tujuannya adalah untuk menganalisis kinerja penjualan oleh distributor, dengan mempertimbangkan variasi regional. Saya memiliki variabel berikut: distcode: ID distributor, dengan sekitar 800 level region: ID geografis tingkat...
Saya mendapat kesan bahwa fungsi lmer()dalam lme4paket tidak menghasilkan nilai-p (lihat lmer, nilai-p dan semua itu ). Saya telah menggunakan MCMC menghasilkan nilai p sebagai per pertanyaan ini: Efek signifikan dalam lme4model campuran dan pertanyaan ini: Tidak dapat menemukan nilai-p dalam...
Saya memiliki satu set data yang sangat kecil tentang kelimpahan lebah soliter yang saya kesulitan analisis. Ini menghitung data, dan hampir semua penghitungan berada dalam satu perawatan dengan sebagian besar nol dalam perawatan lainnya. Ada juga beberapa nilai yang sangat tinggi (masing-masing di...
Data saya dijelaskan di sini Apa yang dapat menyebabkan "Error () model singular error" di aov saat memasang ANOVA tindakan berulang? Saya mencoba untuk melihat efek dari interaksi menggunakan lmersehingga kasus dasar saya adalah: my_null.model <- lmer(value ~ Condition+Scenario+...
Saya adalah pengguna yang lebih terbiasa dengan R, dan telah mencoba memperkirakan lereng acak (koefisien seleksi) untuk sekitar 35 individu selama 5 tahun untuk empat variabel habitat. Variabel respons adalah apakah suatu lokasi merupakan habitat "bekas" (1) atau "tersedia" (0) ("digunakan" di...
Baru-baru ini saya diberitahu bahwa tidak mungkin untuk menggabungkan kovariat yang bervariasi waktu dalam model campuran longitudinal tanpa memperkenalkan jeda waktu untuk kovariat ini. Bisakah Anda mengkonfirmasi / menolak ini? Apakah Anda punya referensi tentang situasi ini? Saya mengusulkan...
Bagaimana seseorang dapat memperoleh bobot regresi standar (efek tetap) dari regresi multilevel? Dan, sebagai "add-on": Apa cara termudah untuk mendapatkan bobot standar ini dari mer-object (dari lmerfungsi lme4paket di
Saya sedang mengerjakan kumpulan data untuk mengevaluasi dampak pengeringan pada aktivitas mikroba sedimen. Tujuannya adalah untuk menentukan apakah dampak pengeringan bervariasi di antara jenis sedimen dan / atau kedalaman sedimen. Desain eksperimental adalah sebagai berikut: Faktor pertama...
Saya sudah cocok dengan model lmer dengan yang berikut ini (walaupun dibuat keluaran): Random effects: Groups Name Std.Dev. day:sample (Intercept) 0.09 sample (Intercept) 0.42 Residual 0.023 Saya benar-benar ingin membangun interval kepercayaan untuk setiap efek menggunakan rumus...
Saya berharap seseorang dapat membantu dengan apa yang saya pikir merupakan pertanyaan yang relatif sederhana, dan saya pikir saya tahu jawabannya tetapi tanpa konfirmasi itu sudah menjadi sesuatu yang saya tidak bisa yakin. Saya memiliki beberapa data hitungan sebagai variabel respon dan saya...
The mgcvpaket untuk Rmemiliki dua fungsi untuk pas interaksi produk tensor: te()dan ti(). Saya memahami pembagian kerja dasar antara keduanya (menyesuaikan interaksi non-linear vs menguraikan interaksi ini menjadi efek utama dan interaksi). Yang tidak saya mengerti adalah mengapa te(x1, x2)dan...
Saya menjalankan glmm dengan variabel respons binomial dan prediktor kategoris. Efek acak diberikan oleh desain bersarang yang digunakan untuk pengumpulan data. Data terlihat seperti ini: m.gen1$treatment [1] sucrose control protein control no_injection ..... Levels: no_injection control sucrose...
Pertanyaan saya didasarkan pada respons ini yang menunjukkan lme4::lmermodel mana yang sesuai dengan langkah-langkah berulang ANOVA dua arah: require(lme4) set.seed(1234) d <- data.frame( y = rnorm(96), subject = factor(rep(1:12, 4)), a = factor(rep(1:2, each=24)), b = factor(rep(rep(1:2,...
Saya sedang menyelesaikan beberapa analisis pada sejumlah besar data. Saya ingin mengambil model linier yang digunakan di bagian pertama pekerjaan dan menyesuaikannya menggunakan model campuran linier (LME). LME akan sangat mirip dengan pengecualian bahwa salah satu variabel yang digunakan dalam...
Diberikan tiga variabel, ydan x, yang positif kontinu, dan z, yang kategorikal, saya memiliki dua model kandidat yang diberikan oleh: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) dan fit.fe <- lm( y ~ 1 + x ) Saya berharap untuk membandingkan model ini untuk menentukan model mana...
Saya ingin memodelkan dua variabel waktu yang berbeda, beberapa di antaranya sangat collinear dalam data saya (usia + kohort = periode). Melakukan ini saya mengalami beberapa masalah dengan lmerdan dan interaksi poly(), tetapi mungkin tidak terbatas lmer, saya mendapat hasil yang sama dengan...
Saya melakukan tes post-hoc pada model efek campuran linear dalam R( lme4paket). Saya menggunakan multcomppaket ( glht()fungsi) untuk melakukan tes post-hoc. Desain eksperimental saya adalah tindakan berulang, dengan efek blok acak. Model ditentukan sebagai: mymod <- lmer(variable ~ treatment...
Demi contoh sederhana berasumsi bahwa ada dua model regresi linier Model 1 memiliki tiga prediktor, x1a, x2b, danx2c Model 2 memiliki tiga prediktor dari model 1 dan dua prediktor tambahan x2adanx2b Ada persamaan regresi populasi di mana varians populasi yang dijelaskan adalah untuk Model 1 dan...
Saya mencoba mereproduksi beberapa uji interaksi antara keduanya lmdan lmerpada tindakan berulang (2x2x2). Alasan saya ingin membandingkan kedua metode ini adalah karena GLM SPSS untuk tindakan berulang menghasilkan hasil yang sama persis dengan lmpendekatan yang disajikan di sini, jadi pada...