Doğrusal karışık efekt modelinin sonuçlarını gösteren grafikler


13

R'deki doğrusal karışık efekt modellemesini kullanarak bazı verileri analiz ediyorum. Sonuçlarla bir poster yapmayı planlıyorum ve sadece karışık efekt modelleri ile deneyimli herhangi bir kişinin, sonuçları göstermek için hangi grafikleri kullanacağını önerebileceğini merak ediyorum modeli. Kalan arsalar, orijinal değerlere karşı uygun değerlerin arsalarını vb. Düşünüyordum.

Bunun verilerime çok bağlı olacağını biliyorum ama sadece doğrusal karma efekt modellerinin sonuçlarını göstermenin en iyi yolunu hissetmeye çalışıyordum. R'de nlme paketini kullanıyorum.

Teşekkürler


3
Tarafından kitap Pinheiro ve Bates çeşitli örnekler içerir. §4.3, "Gömme Modelin İncelenmesi" bölümüne bakabilirsiniz.
Sergio

Yanıtlar:


22

Modelinize bağlıdır, ancak, tecrübelerime göre, karışık efekt modelleri hakkında iyi bir anlayışa sahip olmayan iş arkadaşları bile, tahminleri farklı gruplama düzeyleriyle çizdiğinizde olduğu gibi:

library(nlme)
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1|Subject)

newdat <- expand.grid(Sex=unique(Orthodont$Sex),
                  age=c(min(Orthodont$age),
                            max(Orthodont$age)))

library(ggplot2)
p <- ggplot(Orthodont, aes(x=age, y=distance, colour=Sex)) +
  geom_point(size=3) +
  geom_line(aes(y=predict(fm2), group=Subject, size="Subjects")) +
  geom_line(data=newdat, aes(y=predict(fm2, level=0, newdata=newdat), size="Population")) +
  scale_size_manual(name="Predictions", values=c("Subjects"=0.5, "Population"=3)) +
  theme_bw(base_size=22) 
print(p)

resim açıklamasını buraya girin


@ Roland, cevabın için teşekkürler. Modelim, bazı değişkenlerle bağımsız ve bağımlı değişkenler içeren doğrusal karma efekt modelidir.
John_dydx

@ Roland, sadece fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1|Subject)aynı olup olmadığını sorabilir miyim fm3 <- lme(distance ~ age*Sex, data = Orthodont, random = ~ 1|Subject) ? Sex'i model için bir eş değişken olarak kullanmaya çalışıyorum.
John_dydx

No age * Sexaynıdır age + Sex + age:Sexbu etkileşimi içerir, diğer bir deyişle.
Roland

1
Evet tabi ki. Kaldırmanız gerekecek colour=Sex.
Roland

1
Evet, ancak bu temel ggplot2 işlevselliğidir. Belgeleri ve eğiticileri inceleyin. Kullanmak isteyebilirsiniz scale_colour_manual.
Roland
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.