1
Karışık efekt modellerinde tahmin: rastgele efektlerle ne yapmalı?
Bu varsayımsal veri kümesini ele alalım: set.seed(12345) num.subjects <- 10 dose <- rep(c(1,10,50,100), num.subjects) subject <- rep(1:num.subjects, each=4) group <- rep(1:2, each=num.subjects/2*4) response <- dose*dose/10 * group + rnorm(length(dose), 50, 30) df <- data.frame(dose=dose, response=response, subject=subject, group=group) Kullanabileceğimiz lmerastgele bir etki modeli ile tepkisini modellemek için: require(nlme) model <- lme(response …