আমি ব্যবহার করার চেষ্টা করছি lmeথেকে nlmeথেকে প্রতিলিপি ফলাফলে প্যাকেজ aovপুনরাবৃত্তি পরিমাপ করে ANOVAs জন্য। আমি এটি একটি একক-ফ্যাক্টর পুনরাবৃত্তি ব্যবস্থা পরীক্ষার জন্য এবং একটি-বিষয় ফ্যাক্টর এবং একটি-বিষয় ফ্যাক্টরের মধ্যে একটি দ্বি-ফ্যাক্টর পরীক্ষার জন্য করেছি, তবে এর মধ্যে দুটি সহ একটি দ্বি-ফ্যাক্টর পরীক্ষার জন্য এটি করতে আমার সমস্যা হচ্ছে - সাবজেক্টস ফ্যাক্টর।
একটি উদাহরণ নীচে প্রদর্শিত হয়। Aএবং Bস্থির-প্রভাব কারণ এবং subjectএলোমেলো-প্রভাব ফ্যাক্টর।
set.seed(1)
d <- data.frame(
Y = rnorm(48),
subject = factor(rep(1:12, 4)),
A = factor(rep(1:2, each=24)),
B = factor(rep(rep(1:2, each=12), 2)))
summary(aov(Y ~ A*B + Error(subject/(A*B)), data=d)) # Standard repeated measures ANOVA
library(nlme)
# Attempts:
anova(lme(Y ~ A*B, data=d, random = ~ 1 | subject)) # not same as above
anova(lme(Y ~ A*B, data=d, random = ~ 1 | subject/(A+B))) # gives error
আমি পিনহিরো এবং বেটস বইতে এর ব্যাখ্যা দেখতে পেলাম না তবে আমি এটি উপেক্ষা করেছি।
aovকলটি কেবল একটি স্ট্যান্ডার্ড পুনরাবৃত্তি-পরিমাপ নকশা, যা এক হিসাবে হালকা বিশ্লেষণ করবেlmer(Y~A*B+(1|subject))। (যদিও আরও জটিল মডেলগুলির জন্য এই উত্তরটি দেখুন যেগুলি পার্শ্ব -এস-এফেক্টের বৈচিত্র এবং পারস্পরিক সম্পর্কের প্রাক্কলন অনুমান করে: stats.stackexchange.com/questions/13166/rs-lmer-cheat-sheet/… )