Основная идея состоит в том, чтобы создать линейную модель в R таким образом, чтобы FinH объяснялся SoilNkh, dDDSP, dDDSP2, Provenance, Site, Genotype и Block, где SoilNkh, dDDSP и dDDSP2 являются непрерывными ковариатами, Provenance, Site, Genotype иБлоки - это факторы, Сайт и Провенанс фиксированы, а Генотип и Блок случайны.Кроме того, Genotype вложен в Provenance и Block на сайте.
Поскольку порядок, в котором расположены переменные, имеет важное значение, это должен быть тип I типа Anova с параметрами в следующем порядке: FinH ~ SoilNkh, Site, dDDSP, dDDSP2, Provenance, Site: Provenance, Provenance /Генотип, Сайт / Блок
Для фиксированной части у меня все в порядке:
test31 <-lm(FinH~SoilNkh + Site + dDDSP + dDDSP2 + Provenance + Site:Provenance ,data=d1)
test32 <-aov(FinH~SoilNkh + Site + dDDSP + dDDSP2 + Provenance + Site:Provenance ,data=d1
При попытке указать случайную часть, взяв вышеуказанный текст в качестве отправной точки, начинаются проблемы:)
Мне кажется, это должно быть в форме:
test64 <- lme(FinH~SoilNkh + Site + dDDSP + dDDSP2 + Provenance + Site:Provenance,
random = ~1|Provenance/Genotype + ~1|Site/Block,data=d1)
, но я не могу избежать ошибки
"Ошибка в getGroups.data.frame (dataMix,groups): неверная формула для групп "
Я потерян для подсказок, правда, поэтому любой совет будет отличным!Если будут предоставлены какие-либо данные, я, конечно, с радостью предоставлю, но пока не могу понять, как ...
Заранее спасибо за ваше время.