Синтаксис для УНИАНОВА в R

Question

Я пытаюсь реплицировать следующий синтаксис SPSS для unianova типа I в R

UNIANOVA FinH BY Site Provenance Block Genotype WITH SoilNkh dDDSP dDDSP2
/RANDOM=Block Genotype
/METHOD=SSTYPE(1)
/INTERCEPT=INCLUDE
/POSTHOC=Site (SNK) 
/SAVE=PRED RESID
/CRITERIA=ALPHA(0.05)
/DESIGN=SoilNkh Site dDDSP dDDSP2 Provenance Provenance*Site Genotype(Provenance)
Block(Site).

Если сайт и провенанс фиксированы, генотип, вложенный в провенанс, и блок, вложенный в сайт, являются случайными, а dDDSP, dDDSP2 иSoilNkh.Это должен быть тип I ANOVA.

Лучшее, что я могу придумать в отношении R-синтаксиса, это

m11 <- lme(FinH ~ 1 + SoilNkh + Site + dDDSP + dDDSP2 + Provenance + Site:Provenance, 
    random = (1|Provenance/Genotype) + (1|Site/Block), data=spss)
Anova(m11,type="I")

и довольно много вариантов (что не должно быть возможным:)), но безрезультатно.Честно говоря, я заблудился во всем, что Google может предложить по этому вопросу.

Я был бы очень признателен за любую помощь в переводе синтаксиса.

Заранее спасибо!

Answer 1

Это, кажется, приближает меня:

FinH <- lmer(FinH~SoilN + Site + dDDSP + dDDSP2 + Provenance +
                  Site:Provenance +(1|Provenance/Genotype) + (1|Site/Block),
              REML = FALSE,
              data=subset(spss,Year=="2011"))

Но, как вы видите, не совсем идентично.Также мне не удалось рассчитать статистику для случайных эффектов.Я, конечно, продолжу работать над копированием анализа, но всегда приветствуются дальнейшие советы (и приемы).