लमे [बंद] में कई (अलग) यादृच्छिक प्रभाव निर्दिष्ट करना

Question

मैं आर संकुल nlme और lme4 में काम कर रहा था , कई यादृच्छिक प्रभावों के साथ मॉडल निर्दिष्ट करने की कोशिश कर रहा था । मैंने पाया, कि केवल विक्म ही विचरण की विषम संरचना को निर्दिष्ट करने की अनुमति देता है। इसलिए, मुझे एक मॉडल मिला, जहां तापमान (वाई) समय (घंटों में) पर निर्भर करता है, अवरोधन तारीख और वर्ष के अनुसार बदलता रहता है, और विचरण भी वर्ष के अनुसार बदलता रहता है:

fit1 <- lme(Y ~ time, random=~1|year/date, data=X, weights=varIdent(form=~1|year))

हालाँकि, अगर मुझे एक और यादृच्छिक शब्द जोड़ने की जरूरत है (समय तिथि के अनुसार बदलती है), और इस तरह मॉडल निर्दिष्ट करें:

fit2 <- lme(Y ~ time, random=list(~1|year, ~time-1|date,  ~1|date), data=X, 
            weights=varIdent(form=~1|year))

यादृच्छिक प्रभाव एक दूसरे में निहित हो जाते हैं: वर्ष में तारीख; और फिर तारीख में तारीख और वर्ष में।

मैंने भी कोशिश की

one  <- rep(1, length(Y))
fit3 <- lme(Y ~ time, random=list(one=pdBlocked(list(pdSymm(~1|year/date), 
            pdSymm(~time-1|year)))), data=X, weights=varIdent(form=~1|year))

लेकिन यह एक त्रुटि देता है:

Error in pdConstruct.pdBlocked(object, form = form, nam = nam, data = data,  :
  cannot have duplicated column names in a "pdMat" object

मैं समझता हूं कि इसी तरह की समस्या से संबंधित पहले से ही कई सवाल हैं, लेकिन मुझे वास्तव में मेरे मामले का जवाब नहीं मिला। क्या आप मुझे मॉडल के सही विनिर्देश के साथ मदद कर सकते हैं?

Answer 1

कई संघर्षों के बाद मुझे अपनी समस्या का हल मिला, जिसे मैं यहाँ पोस्ट कर रहा हूँ अगर किसी के पास ऐसे ही सवाल होंगे:

fit <- lme(Y ~ time, random=list(year=~1, date=~time), data=X, weights=varIdent(form=~1|year))