আমি একটি মিশ্র প্রভাব মডেল ব্যবহার করে কিছু তথ্য বিশ্লেষণ করার চেষ্টা করছি। আমি যে তথ্য সংগ্রহ করেছি তা সময়ের সাথে সাথে বিভিন্ন জিনোটাইপের কিছু তরুণ প্রাণীর ওজন উপস্থাপন করে।
আমি এখানে প্রস্তাবিত পদ্ধতির ব্যবহার করছি: https://gribblelab.wordpress.com/2009/03/09/repected-measures-anova-used-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
এর মানে কী? আমি কি এটিকে নিরাপদে উপেক্ষা করতে পারি? বা এড়াতে আমার কী করা উচিত?
সম্পাদনা: আমি এই পিডিএফটি পেয়েছি যা বলছে:
যেহেতু এই ক্ষেত্রে স্বয়ংক্রিয়ভাবে আগ্রহের পরামিতিগুলি নির্ধারণ করা অসম্ভব, মাল্টকম্পে এমসিপি () ডিফল্টরূপে কোভারিয়েট এবং ইন্টারঅ্যাকশন উপেক্ষা করে কেবলমাত্র মূল প্রভাবগুলির জন্য তুলনা তৈরি করে । সংস্করণ ১.১-২ থেকে, যথাক্রমে ইন্টারঅ্যাকশন_এভারেজ = ট্রু এবং কোভারিয়েট_ভেরেজ = ট্রু ব্যবহার করে ইন্টারঅ্যাকশন শর্তাবলী ও কোভেরিয়েটগুলি যথাক্রমে গড় নির্দিষ্ট করতে পারে, যেখানে 1.0-0 এর চেয়ে পুরানো সংস্করণগুলি ইন্টারঅ্যাকশন শর্তগুলির চেয়ে স্বয়ংক্রিয়ভাবে গড় হয়। আমরা ব্যবহারকারীদের পরামর্শ দিই, তারা ম্যানুয়ালি, তারা চান তার বিপরীতে সেট লেখেন।যখনই ডিফল্ট কি পরিমাপ করে তার তুলনায় সন্দেহ থাকে যখনই এটি করা উচিত, যা সাধারণত উচ্চতর অর্ডার ইন্টারঅ্যাকশন শর্তাদি মডেলগুলিতে ঘটে। আমরা এই বিষয়ে আরও আলোচনা এবং উদাহরণের জন্য হু (1996), অধ্যায় 7 ~ এবং সেরেল (1971), অধ্যায় 7.3 দেখুন refer
আমার কাছে এই বইগুলিতে অ্যাক্সেস নেই, তবে সম্ভবত এখানে কেউ আছে?