আমি R
( lme4
প্যাকেজ) রৈখিক মিশ্র-প্রভাব মডেলটিতে পোস্ট-হক পরীক্ষা করছি । পোস্ট-হক পরীক্ষা করার জন্য আমি multcomp
প্যাকেজ ( glht()
ফাংশন) ব্যবহার করছি ।
আমার পরীক্ষামূলক নকশাটি এলোমেলো ব্লক প্রভাব সহ পুনরাবৃত্তি-ব্যবস্থাগুলি। মডেলগুলি হিসাবে নির্দিষ্ট করা হয়:
mymod <- lmer(variable ~ treatment * time + (1|block), data = mydata, REML = TRUE)
এখানে আমার ডেটা সংযুক্ত করার পরিবর্তে, আমি প্যাকেজের warpbreaks
মধ্যে ডেকে আনা ডেটা থেকে কাজ করছি multcomp
।
data <- warpbreaks
warpbreaks$rand <- NA
আমার "ব্লক" এফেক্টটি নকল করতে আমি একটি অতিরিক্ত র্যান্ডম ভেরিয়েবল যুক্ত করেছি:
warpbreaks$rand <- rep(c("foo", "bar", "bee"), nrow(warpbreaks)/3)
এটি আমার মডেলটিকে নকল করে:
mod <- lmer(breaks ~ tension * wool + (1|rand), data = warpbreaks)
আমি "উদাহরণ সচেতন অতিরিক্ত Multcomp Examples- 2 ওয়ে Anova" আপনি মাত্রা মধ্যে টান মাত্রা একটি তুলনামূলক এই উদাহরণ বিশালাকার wool
।
আমি যদি এর বিপরীতটি করতে চাই - এর স্তরের wool
মধ্যে স্তরগুলির তুলনা করব tension
? (আমার ক্ষেত্রে, এটি সময়ের স্তরের (তিন - জুন, জুলাই, আগস্ট) এর মধ্যে চিকিত্সার স্তরের (দুটি - 0, 1) তুলনা করা হবে।
আমি এটি করার জন্য নীচের কোডটি নিয়ে এসেছি, তবে এটি কাজ করছে বলে মনে হচ্ছে না (নীচে ত্রুটির বার্তা দেখুন)।
প্রথমত, উদাহরণ থেকে ( wool
এবং tension
স্থানগুলি বদলে নেওয়া):
tmp <- expand.grid(wool = unique(warpbreaks$wool), tension = unique(warpbreaks$tension))
X <- model.matrix(~ tension * wool, data = tmp)
glht(mod, linfct = X)
Tukey <- contrMat(table(warpbreaks$wool), "Tukey")
K1 <- cbind(Tukey, matrix(0, nrow = nrow(Tukey), ncol = ncol(Tukey)))
rownames(K1) <- paste(levels(warpbreaks$tension)[1], rownames(K1), sep = ":")
K2 <- cbind(matrix(0, nrow = nrow(Tukey), ncol = ncol(Tukey)), Tukey)
rownames(K2) <- paste(levels(warpbreaks$tension)[2], rownames(K2), sep = ":")
এখান থেকে নীচে আমার নিজস্ব কোড:
K3 <- cbind(matrix(0, nrow = nrow(Tukey), ncol = ncol(Tukey)), Tukey)
rownames(K2) <- paste(levels(warpbreaks$tension)[3], rownames(K3), sep = ":")
K <- rbind(K1, K2, K3)
colnames(K) <- c(colnames(Tukey), colnames(Tukey))
> summary(glht(mod, linfct = K %*% X))
Error in summary(glht(mod, linfct = K %*% X)) :
error in evaluating the argument 'object' in selecting a method for function 'summary': Error in K %*% X : non-conformable arguments