আমি বর্তমানে আর প্যাকেজ lme4 ব্যবহার করছি ।
আমি এলোমেলো প্রভাব সহ লিনিয়ার মিশ্র প্রভাবগুলির মডেলগুলি ব্যবহার করছি:
library(lme4)
mod1 <- lmer(r1 ~ (1 | site), data = sample_set) #Only random effects
mod2 <- lmer(r1 ~ p1 + (1 | site), data = sample_set) #One fixed effect +
# random effects
mod3 <- lmer(r1 ~ p1 + p2 + (1 | site), data = sample_set) #Two fixed effects +
# random effects
মডেলগুলির তুলনা করতে, আমি anova
ফাংশনটি ব্যবহার করছি এবং সর্বনিম্ন এআইসির মডেলের তুলনায় এআইসির পার্থক্যগুলি দেখছি:
anova(mod1, mod2, mod3)
মডেলগুলির তুলনা করার জন্য উপরেরটি ঠিক আছে।
যাইহোক, প্রতিটি মডেলের জন্য উপযুক্ত ব্যবস্থার সদ্ব্যবহারের ব্যাখ্যা করার জন্য আমার কিছু সহজ উপায়ও প্রয়োজন। এমন ব্যবস্থা নিয়ে কারও কি অভিজ্ঞতা আছে? আমি কিছু গবেষণা করেছি এবং মিশ্রিত প্রভাবগুলির মডেলগুলির স্থির প্রভাবগুলির জন্য আর স্কোয়ারে জার্নাল পত্রিকা রয়েছে:
- চেং, জে।, এডওয়ার্ডস, এলজে, মালদোনাদো-মোলিনা, এমএম, কোমরো, কেএ, এবং মুলার, কেই (২০১০)। বাস্তব লোকের জন্য বাস্তব অনুদায়ী ডেটা বিশ্লেষণ: একটি ভাল যথেষ্ট মিশ্র মডেল তৈরি করা। মেডিসিনে পরিসংখ্যান, 29 (4), 504-520। doi: 10.1002 / sim.3775
- এডওয়ার্ডস, এলজে, মুলার, কেই, ওল্ফঞ্জার, আরডি, কাকিশ, বিএফ, এবং স্ক্যাবেনবার্গার, ও। (২০০৮)। রৈখিক মিশ্র মডেলের স্থির প্রভাবগুলির জন্য একটি আর 2 পরিসংখ্যান। মেডিসিনে পরিসংখ্যান, 27 (29), 6137-6157। doi: 10.1002 / sim.3429
তবে মনে হয়, উপরের কাগজপত্রে প্রস্তাবিত পদক্ষেপের মতো পদক্ষেপের ব্যবহারকে ঘিরে কিছু সমালোচনা রয়েছে।
কেউ দয়া করে আমার মডেলগুলিতে প্রযোজ্য এমন উপযুক্ত পদক্ষেপের ধার্মিকতা ব্যাখ্যা করার জন্য কয়েকটি সহজ পরামর্শ দিতে পারেন?
mixed()
আমার এএফেক্স প্যাকেজে ফাংশনটি দেখুন ( বিকাশ সংস্করণটিতে প্যারামেট্রিক বুটস্ট্র্যাপও রয়েছে )। কিছু রেফারেন্স জন্য এখানে দেখুন ।
KRmodcomp
প্যাকেজ থেকে ব্যবহার করে pbkrtest
। আপনি KRmodcomp
মডেলদের তুলনা করতে সরাসরি ব্যবহার করতে পারেন ।