Saya mencoba untuk menghasilkan model campuran linier kode R adalah sebagai berikut. lme (Average.payoff ~ Game + Type + Others.Type + Game: Type + Game: Others.Type + Type: Others.Type, random = ~ 1 | Subjek, metode = "REML", data = Subjectsm1) -> lme1 Istilah respons Average.payoff...
Misalkan saya punya satu sampel frekuensi dari 4 peristiwa yang mungkin: Event1 - 5 E2 - 1 E3 - 0 E4 - 12 dan saya memiliki probabilitas yang diharapkan dari peristiwa saya terjadi: p1 - 0.2 p2 - 0.1 p3 - 0.1 p4 - 0.6 Dengan jumlah frekuensi yang diamati dari empat acara saya (18) saya dapat...
Intro Saya memiliki peserta yang berulang kali menyentuh permukaan yang terkontaminasi dengan E. coli dalam dua kondisi ( A = memakai sarung tangan, B = tanpa sarung tangan). Saya ingin tahu apakah ada perbedaan antara jumlah bakteri di ujung jari mereka dengan dan tanpa sarung tangan tetapi juga...
Saya memasang glmermodel dalam lme4paket R. Saya mencari tabel anova dengan nilai p yang ditunjukkan di sini, tetapi saya tidak dapat menemukan paket yang cocok. Apakah mungkin melakukannya di R? Model yang saya pas adalah dalam bentuk: model1<-glmer(dmn~period*teethTreated+(1|fullName),...
library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1...
Saya memiliki beberapa kesulitan untuk mendapatkan hasil yang setara antara aovmodel tindakan berulang antara-dalam dan lmermodel campuran. Data dan skrip saya terlihat sebagai berikut data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id...
Dalam Model Generalized Linear Mixed logistik (family = binomial), saya tidak tahu bagaimana menafsirkan varians efek acak: Random effects: Groups Name Variance Std.Dev. HOSPITAL (Intercept) 0.4295 0.6554 Number of obs: 2275, groups: HOSPITAL, 14 Bagaimana cara menginterpretasikan hasil...
Saya mensimulasikan uji coba Bernoulli dengan acak logitθ∼N(logitθ0,12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)antara kelompok dan kemudian saya cocokkan dengan model yang sesuai dengan lme4paket: library(lme4) library(data.table) I <- 30 # number of...
Saya menggunakan lmer di R untuk memeriksa efek kondisi ( cond) pada beberapa hasil. Berikut adalah beberapa data yang dibuat, di mana s adalah pengidentifikasi subjek dan a, bdan ckondisi. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30,...
Saya menyesuaikan model efek campuran dengan istilah spline dalam aplikasi di mana tren dari waktu ke waktu dikenal sebagai kurva-linier. Namun, apa yang ingin saya nilai adalah apakah tren kurva-linier terjadi karena penyimpangan individu dari linieritas, atau apakah itu efek pada tingkat kelompok...
Kejadian yang tidak biasa ketika berhadapan dengan model campuran maksimal yang kompleks (memperkirakan semua efek acak yang mungkin untuk data dan model yang diberikan) adalah sempurna (+1 atau -1) atau korelasi yang hampir sempurna di antara beberapa efek acak. Untuk tujuan diskusi, mari kita...
Saya punya beberapa pertanyaan tentang spesifikasi dan interpretasi GLMM. 3 pertanyaan pasti statistik dan 2 lebih spesifik tentang R. Saya posting di sini karena pada akhirnya saya pikir masalahnya adalah interpretasi hasil GLMM. Saat ini saya mencoba untuk menyesuaikan GLMM. Saya menggunakan...
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 Sedimen...
Saya telah membaca bahwa distribusi untuk estimator efek acak di lme4 sangat miring dan untuk alasan ini kesalahan standar tidak dilaporkan. Saya ingin tahu apakah ada yang bisa memberikan referensi untuk ini? Saya memiliki akses ke buku karya Bates dan Pinherio, tetapi tidak untuk Raudenbush dan...
Misalkan A adalah faktor tetap dan B adalah faktor acak. A dan B dilintasi. Faktor acak C lainnya bersarang di B. Mempertimbangkan interaksi AB dan AC, bagaimana saya harus menyesuaikan model dengan menggunakan fungsi lme dalam paket nlme? Terima
Saya mencoba mendapatkan prediksi untuk pengamatan dari objek lme. Ini seharusnya cukup mudah. Namun, karena saya mendapatkan berbagai jenis kesalahan untuk uji coba yang berbeda, menurut saya saya kehilangan sesuatu. Model saya adalah sebagai berikut: model <- lme(log(child_mortality) ~...
Saya memiliki 2 faktor Adan B(5 × 3) dan satu kovariat Xdalam desain subjek. Inilah cara saya menentukan keseluruhan model saya: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Interpretasi saya adalah bahwa saya melihat grafik y~x, di mana kemiringan berubah karena kovariat, dan garis...
Saya memiliki masalah dalam menemukan solusi tentang cara menjalankan tes post-hoc (Tukey HSD) setelah 2-faktor (keduanya dalam-mata pelajaran) mengukur ANOVA berulang dalam R. Untuk ANOVA, saya telah menggunakan fungsi aov: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1)) Setelah...
Saya punya masalah dengan menafsirkan interaksi 2-arah dan 3-arah di lmer. DV saya adalah tinggi yang merupakan variabel kontinu. Semua IV adalah variabel kategori. Faktor pertama adalah hewan, baik tikus atau singa. Faktor kedua adalah jenis kelamin, baik pria maupun wanita. Faktor ketiga adalah...
Saya ingin tahu tentang bagaimana paket lmerTest di R, khususnya fungsi "rand", menangani tes efek acak. Pertimbangkan contoh dari pdf lmerTest tentang CRAN yang menggunakan kumpulan data "wortel" bawaan : #import lme4 package and lmerTest package library(lmerTest) #lmer model with correlation...