R paketleri nlme ve lme4'te çalışıyordum, çoklu rastgele etkilere sahip modelleri belirlemeye çalışıyordum. Sadece nlme'nin varyansın heterojen yapısını belirleyebildiğini buldum . Bu nedenle, sıcaklığın (Y) zamana bağlı (saat cinsinden), kesişimin tarihe ve yıla göre değiştiği ve varyansın da yıla göre değiştiği bir modelim var:
fit1 <- lme(Y ~ time, random=~1|year/date, data=X, weights=varIdent(form=~1|year))
Ancak, başka bir rastgele terim eklemem gerekirse (tarihe göre değişen zaman) ve modeli şöyle belirtin:
fit2 <- lme(Y ~ time, random=list(~1|year, ~time-1|date, ~1|date), data=X,
weights=varIdent(form=~1|year))
rastgele etkiler birbirlerine yuvalanmış hale gelir: yıl içindeki tarih; ve sonra tarih ve yılda tarih.
Ben de denedim
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))
ama bir hata veriyor:
Error in pdConstruct.pdBlocked(object, form = form, nam = nam, data = data, :
cannot have duplicated column names in a "pdMat" object
Benzer sorunla ilgili birçok sorunun zaten olduğunu anladım, ancak davamın cevabını gerçekten bulamadım. Modelin doğru tanımlamasına yardımcı olabilir misiniz?