मैं 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)क्या वो सही है? मुझे यकीन नहीं है कि मुझे अपने फॉर्मूले में प्लॉट / सबप्लॉट का स्तर रखना चाहिए। कोई निश्चित प्रभाव महत्वपूर्ण नहीं है लेकिन यादृच्छिक प्रभाव बहुत महत्वपूर्ण हैं।