Silakan pertimbangkan data ini: dt.m <- structure(list(id = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), occasion = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("g1", "g2"), class =...
Saya mencoba menggunakan lmedari nlmepaket untuk mereplikasi hasil dari aovuntuk ANOVA tindakan berulang. Saya telah melakukan ini untuk percobaan pengukuran berulang faktor tunggal dan untuk eksperimen dua faktor dengan satu faktor antar-subjek dan satu faktor dalam-subyek, tetapi saya mengalami...
Saya memiliki GLMM dengan distribusi binomial dan fungsi tautan logit dan saya merasa bahwa aspek penting dari data tidak terwakili dengan baik dalam model. Untuk menguji ini, saya ingin tahu apakah data dijelaskan dengan baik oleh fungsi linear pada skala logit. Oleh karena itu, saya ingin tahu...
Ada perbedaan yang membuat saya bingung dengan model campuran, dan saya bertanya-tanya apakah saya bisa mendapatkan kejelasan tentang itu. Mari kita asumsikan Anda punya model campuran data hitungan. Ada variabel yang Anda tahu Anda inginkan sebagai efek tetap (A) dan variabel lain untuk waktu (T),...
Katakanlah kami sedang mengerjakan model efek acak dari beberapa data jumlah dari waktu ke waktu, dan kami ingin mengendalikan beberapa tren. Biasanya, Anda akan melakukan sesuatu seperti: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") untuk memasukkan bentuk kuadrat untuk...
Pertimbangkan data berikut dari dua arah dalam desain subjek: df <- "http://personality-project.org/r/datasets/R.appendix4.data" df <- read.table(df,header=T) head(df) Observation Subject Task Valence Recall 1 1 Jim Free Neg 8 2 2 Jim Free Neu 9 3 3 Jim Free Pos 5 4 4 Jim Cued Neg 7 5 5 Jim...
Katakanlah kita harus GLMM mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Model-model ini tidak bersarang dalam arti biasa: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) jadi kita tidak bisa...
Ini adalah bingkai data saya: Group <- c("G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3") Subject <-...
Latar Belakang: Catatan: Kumpulan data dan kode-r saya termasuk di bawah teks Saya ingin menggunakan AIC untuk membandingkan dua model efek campuran yang dihasilkan menggunakan paket lme4 di R. Setiap model memiliki satu efek tetap dan satu efek acak. Efek tetap berbeda di antara model, tetapi...
Saya tahu bahwa salah satu keuntungan dari model campuran adalah bahwa mereka memungkinkan untuk menentukan matriks varians-kovarians untuk data (simetri gabungan, autoregresif, tidak terstruktur, dll.) Namun, lmerfungsi dalam R tidak memungkinkan spesifikasi mudah dari matriks ini. Adakah yang...
Saya mengerti bahwa kami menggunakan model efek acak (atau efek campuran) ketika kami percaya bahwa beberapa parameter model bervariasi secara acak di beberapa faktor pengelompokan. Saya memiliki keinginan untuk menyesuaikan model di mana respons telah dinormalisasi dan terpusat (tidak sempurna,...
Saya memiliki data dari pasien yang dirawat dengan 2 jenis perawatan selama operasi. Saya perlu menganalisis efeknya pada detak jantung. Pengukuran detak jantung dilakukan setiap 15 menit. Mengingat bahwa lama operasi dapat berbeda untuk setiap pasien, setiap pasien dapat memiliki antara 7 dan 10...
Saya telah menggunakan MCMCglmmpaket baru-baru ini. Saya bingung dengan apa yang disebut dalam dokumentasi sebagai struktur-R dan struktur-G. Ini tampaknya berhubungan dengan efek acak - khususnya menentukan parameter untuk distribusi sebelumnya, tetapi diskusi dalam dokumentasi tampaknya...
Rumus model dalam R seperti y ~ x + a*b + c:d didasarkan pada apa yang disebut notasi Wilkinson : Wilkinson dan Rogers 1973, Deskripsi Simbolik Model Faktorial untuk Analisis Varians . Makalah ini tidak membahas notasi untuk model campuran (yang mungkin belum ada saat itu). Jadi di mana formula...
Saya mengalami kesulitan memahami hasil lmer()model saya . Ini adalah model sederhana dari variabel hasil (Dukungan) dengan berbagai penyadapan Negara / efek acak Negara: mlm1 <- lmer(Support ~ (1 | State)) Hasilnya summary(mlm1)adalah: Linear mixed model fit by REML Formula: Support ~ (1 |...
Dalam studi saya, saya akan mengukur beban kerja dengan beberapa metrik. Dengan variabilitas denyut jantung (HRV), aktivitas electrodermal (EDA) dan dengan skala subjektif (IWS). Setelah normalisasi, IWS memiliki tiga nilai: Beban kerja lebih rendah dari biasanya Beban kerja rata-rata Beban kerja...
Katakanlah kita punya model mod <- Y ~ X*Condition + (X*Condition|subject) # Y = logit variable # X = continuous variable # Condition = values A and B, dummy coded; the design is repeated # so all participants go through both Conditions # subject = random effects for different subjects...
Saya menggunakan glmerfungsi dari lme4paket dalam R, dan saya menggunakan bobyqapengoptimal (yaitu default dalam kasus saya). Saya mendapat peringatan, dan saya ingin tahu apa artinya. Warning message: In optwrap(optimizer, devfun, start, rho$lower, control = control, : convergence code 3 from...
Saya mencoba menggunakan lme4::glmer()agar sesuai dengan model campuran binomial umum (GLMM) dengan variabel dependen yang bukan biner, tetapi variabel kontinu antara nol dan satu. Orang dapat menganggap variabel ini sebagai probabilitas; sebenarnya itu adalah probabilitas sebagaimana dilaporkan...
Saya bereksperimen dengan algoritma mesin peningkat gradien melalui caretpaket di R. Menggunakan dataset penerimaan perguruan tinggi kecil, saya menjalankan kode berikut: library(caret) ### Load admissions dataset. ### mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") ###...