พล็อตเพื่อแสดงผลลัพธ์ของโมเดลเอฟเฟกต์เชิงเส้นผสม


13

ฉันกำลังวิเคราะห์ข้อมูลบางอย่างโดยใช้การสร้างแบบจำลองเอฟเฟกต์แบบผสมเชิงเส้นในอาร์ฉันวางแผนที่จะสร้างโปสเตอร์ที่มีผลลัพธ์และฉันก็สงสัยว่าถ้าใครมีประสบการณ์กับโมเดลเอฟเฟกต์ผสมสามารถแนะนำแผนการที่จะใช้ แบบ ฉันกำลังคิดเกี่ยวกับพล็อตที่เหลือพล็อตของค่าติดตั้งกับค่าเดิม ฯลฯ

ฉันรู้ว่าสิ่งนี้จะขึ้นอยู่กับข้อมูลของฉันเป็นอย่างมาก แต่ฉันแค่พยายามทำความเข้าใจถึงวิธีที่ดีที่สุดในการอธิบายผลลัพธ์ของตัวแบบผสมเอฟเฟกต์เชิงเส้น ฉันใช้แพ็คเกจ nlme ใน R

ขอบคุณ


3
หนังสือโดยPinheiro และ Batesมีหลายตัวอย่าง คุณสามารถดู§4.3 "การตรวจสอบรุ่นติดตั้ง"
Sergio

คำตอบ:


22

มันขึ้นอยู่กับแบบจำลองของคุณ แต่จากประสบการณ์ของฉันแม้กระทั่งเพื่อนร่วมงานที่ไม่มีความเข้าใจที่ดีเกี่ยวกับแบบจำลองเอฟเฟกต์แบบผสมเช่นถ้าคุณวางแผนการทำนายด้วยระดับการจัดกลุ่มที่แตกต่างกัน:

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)

ป้อนคำอธิบายรูปภาพที่นี่


@ Roland ขอบคุณสำหรับคำตอบของคุณ โมเดลของฉันเป็นแบบจำลองเอฟเฟกต์แบบเชิงเส้นที่ประกอบด้วยตัวแปรอิสระและตัวแปรตามที่มีตัวแปรร่วมบางตัว
John_dydx

@ Roland ผมสามารถเพียงแค่ถามว่าเป็นเช่นเดียวกับfm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1|Subject) fm3 <- lme(distance ~ age*Sex, data = Orthodont, random = ~ 1|Subject) ฉันกำลังพยายามใช้เซ็กส์เป็นผู้ร่วมทำงานให้กับโมเดล
John_dydx

ไม่age * Sexเหมือนกันage + Sex + age:Sexคือมันรวมถึงการมีปฏิสัมพันธ์
Roland

1
ใช่แน่นอน colour=Sexคุณจะต้องลบ
Roland

1
ใช่ แต่นั่นคือฟังก์ชัน ggplot2 พื้นฐาน ศึกษาเอกสารและแบบฝึกหัด scale_colour_manualคุณอาจต้องการที่จะใช้
Roland
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.