এটি কিছুটা স্ব-প্রচারের মতো দেখতে (এবং আমি মনে করি এটি)। তবে আমি আর ( সিআরএএন- তে উপলব্ধ) জন্য একটি lsmeans প্যাকেজ তৈরি করেছি যা এই ধরণের পরিস্থিতিটি হ্যান্ডেল করার জন্য ডিজাইন করা হয়েছে। এটি আপনার উদাহরণের জন্য কীভাবে কাজ করে তা এখানে:
> sample.data <- data.frame(IV=rep(1:4,each=20),DV=rep(c(-3,-3,1,3),each=20)+rnorm(80))
> sample.aov <- aov(DV ~ factor(IV), data = sample.data)
> library("lsmeans")
> (sample.lsm <- lsmeans(sample.aov, "IV"))
IV lsmean SE df lower.CL upper.CL
1 -3.009669 0.2237448 76 -3.4552957 -2.564043
2 -3.046072 0.2237448 76 -3.4916980 -2.600445
3 1.147080 0.2237448 76 0.7014539 1.592707
4 3.049153 0.2237448 76 2.6035264 3.494779
> contrast(sample.lsm, list(mycon = c(-3,-1,1,3)))
contrast estimate SE df t.ratio p.value
mycon 22.36962 1.000617 76 22.356 <.0001
আপনি চাইলে তালিকায় অতিরিক্ত বিপরীতে উল্লেখ করতে পারেন। এই উদাহরণস্বরূপ, আপনি অন্তর্নির্মিত রৈখিক বহুপদী বিপরীতে একই ফলাফল পাবেন:
> con <- contrast(sample.lsm, "poly")
> con
contrast estimate SE df t.ratio p.value
linear 22.369618 1.0006172 76 22.356 <.0001
quadratic 1.938475 0.4474896 76 4.332 <.0001
cubic -6.520633 1.0006172 76 -6.517 <.0001
এটি নিশ্চিত করার জন্য, নোট করুন যে "poly"
স্পেসিফিকেশনটি এটি কল করার নির্দেশ দেয় poly.lsmc
, যা এই ফলাফলগুলি উত্পাদন করে:
> poly.lsmc(1:4)
linear quadratic cubic
1 -3 1 -1
2 -1 -1 3
3 1 -1 -3
4 3 1 1
আপনি যদি কয়েকটি বিপরীতে একটি যৌথ পরীক্ষা করতে চান তবে এর সাথে test
ফাংশনটি ব্যবহার করুন joint = TRUE
। উদাহরণ স্বরূপ,
> test(con, joint = TRUE)
এটি একটি "ধরণের III" পরীক্ষা উত্পন্ন করবে। বিপরীতে car::Anova()
, এটি মডেল-ফিটিং পর্যায়ে ব্যবহৃত বিপরীতে কোডিং নির্বিশেষে এটি সঠিকভাবে করবে। এর কারণ এটি হচ্ছে যে লিনিয়ার ফাংশনগুলি পরীক্ষা করা হচ্ছে তা মডেল হ্রাসের মাধ্যমে স্পষ্টত না বরং সরাসরি নির্দিষ্ট করা হয়। একটি অতিরিক্ত বৈশিষ্ট্য হ'ল এমন একটি ক্ষেত্রে যেখানে পরীক্ষা করা হচ্ছে বৈসাদৃশ্যগুলি রৈখিকভাবে নির্ভরশীল তা সনাক্ত করা হয় এবং সঠিক পরীক্ষার পরিসংখ্যান এবং স্বাধীনতার ডিগ্রি উত্পন্ন হয়।