আমি lme4আমার পরীক্ষামূলক ডিজাইনের জন্য রৈখিক মিশ্র ইফেক্ট মডেলের (সহ ) সূত্র নির্দিষ্ট করার চেষ্টা করছি, তবে আমি নিশ্চিত নই যে আমি এটি সঠিকভাবে করছি।
নকশা: মূলত আমি উদ্ভিদের উপর একটি প্রতিক্রিয়া পরামিতি পরিমাপ করছি। আমি চিকিত্সার 4 স্তর, এবং 2 সেচ স্তর আছে। গাছগুলি 16 টি প্লটে বিভক্ত করা হয়েছে, প্রতিটি প্লটের মধ্যে আমি 4 টি সাব-প্লট নমুনা করি। প্রতিটি উপ-প্লটে আমি 15 থেকে 30 টি পর্যবেক্ষণ (পাওয়া উদ্ভিদের সংখ্যার উপর নির্ভর করে) নিই। অর্থাৎ মোট 1500 সারি রয়েছে।
প্রথমদিকে সাব্প্লট স্তরটি কেবলমাত্র নমুনার উদ্দেশ্যে ছিল এখানে, তবে আমি ভেবেছিলাম মডেলটিতে এটি বিবেচনায় নিতে চাই (64৪-স্তরের ভেরিয়েবল হিসাবে) কারণ আমি দেখেছি যে একটি উপ-প্লট থেকে অন্য উপ-প্লটে অনেকগুলি পরিবর্তনশীলতা রয়েছে এমনকি একই প্লটের অভ্যন্তরে (পুরো প্লটের মধ্যে পরিবর্তনের চেয়ে বড়)।
আমার প্রথম ধারণাটি ছিল:
library(lme4)
fit <- lmer(y ~ treatment*irrigation + (1|subplot/plot), data=mydata)
অথবা
fit <- lmer(y ~ treatment*irrigation + (1|subplot) + (1|plot), data=mydata)এটা কি ঠিক? আমি নিশ্চিত নই যে আমার সূত্রে আমাকে প্লট / সাবপ্লট উভয় স্তরই রাখতে হবে কিনা। কোনও স্থির প্রভাব উল্লেখযোগ্য নয় তবে এলোমেলো প্রভাবগুলি খুব তাৎপর্যপূর্ণ।