ฉันพยายามวิเคราะห์ข้อมูลบางอย่างโดยใช้แบบจำลองเอฟเฟกต์ผสม ข้อมูลที่ฉันรวบรวมนั้นเป็นตัวแทนของน้ำหนักของสัตว์เล็กที่มีจีโนไทป์ต่างกันตามเวลา
ฉันใช้แนวทางที่เสนอที่นี่: https://gribblelab.wordpress.com/2009/03/09/repeated-measures-anova-using-r/
โดยเฉพาะฉันใช้โซลูชัน # 2
ดังนั้นฉันมีสิ่งที่ชอบ
require(nlme)
model <- lme(weight ~ time * Genotype, random = ~1|Animal/time,
data=weights)
av <- anova(model)
ตอนนี้ฉันต้องการเปรียบเทียบหลายอย่าง ใช้multcomp
ฉันสามารถทำได้:
require(multcomp)
comp.geno <- glht(model, linfct=mcp(Genotype="Tukey"))
print(summary(comp.geno))
และแน่นอนฉันสามารถทำเช่นเดียวกันกับเวลา
ฉันมีสองคำถาม:
- ฉันจะใช้
mcp
เพื่อดูปฏิสัมพันธ์ระหว่างเวลาและจีโนไทป์ได้อย่างไร เมื่อฉันวิ่ง
glht
ฉันได้รับคำเตือนนี้:covariate interactions found -- default contrast might be inappropriate
มันหมายความว่าอะไร? ฉันสามารถเพิกเฉยได้อย่างปลอดภัยหรือไม่? หรือฉันควรทำอย่างไรเพื่อหลีกเลี่ยง
แก้ไข: ฉันพบPDF นี้ที่ระบุว่า:
เพราะมันเป็นไปไม่ได้ในการกำหนดค่าพารามิเตอร์ของดอกเบี้ยโดยอัตโนมัติในกรณีนี้ MCP () ใน multcomp จะเริ่มต้นด้วยการสร้างการเปรียบเทียบผลกระทบหลักเพียงอย่างเดียวโดยไม่สนใจตัวแปรและการมีปฏิสัมพันธ์ ตั้งแต่เวอร์ชัน 1.1-2 เราสามารถระบุค่าเฉลี่ยของเงื่อนไขการใช้งานร่วมกันและ covariates โดยใช้อาร์กิวเมนต์ interaction_average = TRUE และ covariate_average = TRUE ตามลำดับในขณะที่รุ่นที่เก่ากว่า 1.0-0 โดยเฉลี่ยโดยอัตโนมัติผ่านข้อกำหนดของการโต้ตอบ เราแนะนำให้ผู้ใช้ทราบว่าพวกเขาเขียนชุดความแตกต่างที่พวกเขาต้องการด้วยตนเองเราควรทำสิ่งนี้เมื่อใดก็ตามที่มีข้อสงสัยเกี่ยวกับสิ่งที่การวัดความคมชัดเริ่มต้นซึ่งมักเกิดขึ้นในรุ่นที่มีเงื่อนไขการโต้ตอบเพื่อที่สูงขึ้น เราอ้างถึง Hsu (1996), บทที่ 7 และ Searle (1971), ตอนที่ 7.3 สำหรับการสนทนาและตัวอย่างเพิ่มเติมเกี่ยวกับปัญหานี้
ฉันไม่มีสิทธิ์เข้าถึงหนังสือเหล่านั้น แต่อาจมีบางคนที่นี่หรือ