我正在尝试在R中为类型I unianova复制以下SPSS语法
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型方差分析。
就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")
和相当多的变种(这是不可能的:) ),但无济于事。坦率地说,我也迷失在谷歌提供的关于这个主题的所有东西中。
我非常感谢在翻译语法方面的任何帮助。
提前感谢!
发布于 2019-02-02 03:37:56
这似乎让我更接近了:
FinH <- lmer(FinH~SoilN + Site + dDDSP + dDDSP2 + Provenance +
Site:Provenance +(1|Provenance/Genotype) + (1|Site/Block),
REML = FALSE,
data=subset(spss,Year=="2011"))
但正如你所看到的,它们并不完全相同。此外,我还没有成功地计算随机效果的统计数据。当然,我将继续工作,以复制分析,但欢迎进一步的提示(和技巧)。
https://stackoverflow.com/questions/53382250
复制相似问题