首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >R中UNIANOVA的语法

R中UNIANOVA的语法
EN

Stack Overflow用户
提问于 2018-11-20 04:35:18
回答 1查看 120关注 0票数 1

我正在尝试在R中为类型I unianova复制以下SPSS语法

代码语言:javascript
运行
复制
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语法而言,我能想到的最好的方法是

代码语言:javascript
运行
复制
m11 <- lme(FinH ~ 1 + SoilNkh + Site + dDDSP + dDDSP2 + Provenance + Site:Provenance, 
    random = (1|Provenance/Genotype) + (1|Site/Block), data=spss)
Anova(m11,type="I")

和相当多的变种(这是不可能的:) ),但无济于事。坦率地说,我也迷失在谷歌提供的关于这个主题的所有东西中。

我非常感谢在翻译语法方面的任何帮助。

提前感谢!

EN

回答 1

Stack Overflow用户

发布于 2019-02-02 03:37:56

这似乎让我更接近了:

代码语言:javascript
运行
复制
FinH <- lmer(FinH~SoilN + Site + dDDSP + dDDSP2 + Provenance +
                  Site:Provenance +(1|Provenance/Genotype) + (1|Site/Block),
              REML = FALSE,
              data=subset(spss,Year=="2011"))

但正如你所看到的,它们并不完全相同。此外,我还没有成功地计算随机效果的统计数据。当然,我将继续工作,以复制分析,但欢迎进一步的提示(和技巧)。

票数 -1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/53382250

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档