Saya telah menjelajahi Google dan situs ini dan saya masih bingung tentang fungsi lmer di perpustakaan lme4. Saya memiliki beberapa data yang dikumpulkan dari bangsal psikiatri yang berbeda, yang memiliki struktur bertingkat. Untuk menyederhanakan, saya akan memilih dua level 2 dan dua variabel...
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...
Pertimbangkan percobaan dengan beberapa partisipan manusia, masing-masing diukur beberapa kali dalam dua kondisi. Model efek campuran dapat dirumuskan (menggunakan sintaks lme4 ) sebagai: fit = lmer( formula = measure ~ (1|participant) + condition ) Sekarang, katakan saya ingin membuat interval...
Misalkan, dalam penelitian terhadap 15 subjek, variabel respons (res) dimodelkan dengan dua variabel penjelas, satu (level) dikategorikan dengan 5 level dan lainnya (waktu respons: RT) kontinu. Dengan lmer dalam paket lme4 dari R, saya punya: fm1 <- lmer(res ~ level * RT + (level-1 | subject),...
Pertanyaan ini mungkin terdengar bodoh, tetapi ... apakah benar bahwa efek acak hanya dapat berlaku untuk variabel kategori (seperti id individu, id populasi, ...), misalkan katakanlah adalah variabel kategori:xsayaxix_i ysayayiy_i ~βxsayaβxi\beta_{x_i} βxsayaβxi\beta_{x_i} ~No r m ( μ ,...
Saya memiliki tipe data berikut. Saya telah mengevaluasi 10 orang yang masing-masing diulang 10 kali. Saya memiliki 10x10 matriks relasi (hubungan antara semua kombinasi individu). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld =...
Rekan saya dan saya sedang menyesuaikan serangkaian model efek campuran linier dan nonlinear dalam R. Kami diminta untuk melakukan validasi silang pada model yang dipasang sehingga orang dapat memverifikasi bahwa efek yang diamati relatif dapat digeneralisasikan. Ini biasanya tugas yang sepele,...
Dalam R, bagaimana cara menentukan model lmer tanpa efek tetap global? Misalnya, jika saya mengatakan sesuatu seperti lmer(y ~ (1 | group) + (0 + x | group), data = my_df) model yang pas akan ysaya j= a +αsaya+βsayaxsaya jyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Bagaimana saya...
Saya mencoba untuk berurusan dengan analisis waktu ke acara menggunakan hasil biner berulang. Misalkan waktu-ke-peristiwa diukur dalam beberapa hari tetapi untuk saat ini kami memutuskan waktu ke minggu. Saya ingin memperkirakan estimator Kaplan-Meier (tetapi memungkinkan untuk kovariat)...
Cangkok berikut diambil dari artikel ini . Saya pemula untuk bootstrap dan mencoba mengimplementasikan bootstrap parametrik, semiparametrik, dan nonparametrik untuk model campuran linier dengan R bootpaket. Kode R Ini Rkode saya : library(SASmixed) library(lme4) library(boot) fm1Cult...
Saya mencoba untuk menganalisis efek Tahun pada variabel logInd untuk kelompok individu tertentu (saya punya 3 kelompok). Model paling sederhana: > fix1 = lm(logInd ~ 0 + Group + Year:Group, data = mydata) > summary(fix1) Call: lm(formula = logInd ~ 0 + Group + Year:Group, data =...
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...
Dalam studi longitudinal, hasil dari unit berulang kali diukur pada titik waktu dengan total pengukuran tetap kali (tetap = pengukuran pada unit diambil pada waktu yang sama). i t mYi tYitY_{it}sayaiitttmmm Unit ditugaskan secara acak baik untuk perawatan, , atau ke kelompok kontrol, . Saya ingin...
Saya punya pertanyaan dasar tentang pendekatan untuk pemodelan rata-rata menggunakan kriteria IT untuk model berat dalam satu set kandidat. Sebagian besar sumber yang saya baca tentang model rata-rata menganjurkan rata-rata estimasi koefisien parameter berdasarkan bobot model (baik menggunakan...
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),...
Haruskah saya memasukkan efek acak dalam model bahkan jika mereka tidak signifikan secara statistik? Saya memiliki desain eksperimental tindakan berulang, di mana setiap individu mengalami tiga perawatan yang berbeda dalam urutan acak. Saya ingin mengontrol efek individu dan ketertiban, tetapi...
Contoh: Saya memiliki kalimat dalam deskripsi pekerjaan: "Java senior engineer in UK". Saya ingin menggunakan model pembelajaran yang mendalam untuk memperkirakannya sebagai 2 kategori: English dan IT jobs. Jika saya menggunakan model klasifikasi tradisional, hanya dapat memprediksi 1 label dengan...
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...
Perbedaan antara variabel yang benar-benar acak (yang harus dimodelkan seperti itu) dan variabel non-acak yang beberapa orang berpendapat dapat dimodelkan sebagai acak jika itu adalah model hierarkis / bertingkat, menjadi buram bagi saya. Bates dan Bolker mencontohkan efek acak dengan kasus...
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...