文章/答案/技术大牛

发布

社区首页 >问答首页 >从r中manova分析中提取误差矩阵

问从r中manova分析中提取误差矩阵
EN

Stack Overflow用户

提问于 2018-06-22 16:02:01

回答 1查看 212关注 0票数 1

我使用car软件包进行了anova分析。这里，我以矩阵的形式得到了误差和因子(=Site)的平方和。密码在这里

library(car)
## One-Way MANOVA
anv<-summary(Anova(lm(cbind(Al, Fe, Mg, Ca, Na) ~ Site, data=Pottery)))
anv

现在，我想从这个结果中提取这两个矩阵，以便进行下一次分析。我想要的输出是这样的-

 >err_ss

       Al          Fe          Mg          Ca         Na
Al 48.2881429  7.08007143  0.60801429  0.10647143 0.58895714
Fe  7.0800714 10.95084571  0.52705714 -0.15519429 0.06675857
Mg  0.6080143  0.52705714 15.42961143  0.43537714 0.02761571
Ca  0.1064714 -0.15519429  0.43537714  0.05148571 0.01007857
Na  0.5889571  0.06675857  0.02761571  0.01007857 0.19929286



>Site_ss 

        Al          Fe          Mg         Ca         Na
Al  175.610319 -149.295533 -130.809707 -5.8891637 -5.3722648
Fe -149.295533  134.221616  117.745035  4.8217866  5.3259491
Mg -130.809707  117.745035  103.350527  4.2091613  4.7105458
Ca   -5.889164    4.821787    4.209161  0.2047027  0.1547830
Na   -5.372265    5.325949    4.710546  0.1547830  0.2582456

使用下面的代码后，我得到了错误矩阵的结果，但未能提取最后一个错误矩阵。

str(anv)
err_ss<-anv$"SSPE"
> err_ss
       Al          Fe          Mg          Ca         Na
Al 48.2881429  7.08007143  0.60801429  0.10647143 0.58895714
Fe  7.0800714 10.95084571  0.52705714 -0.15519429 0.06675857
Mg  0.6080143  0.52705714 15.42961143  0.43537714 0.02761571
Ca  0.1064714 -0.15519429  0.43537714  0.05148571 0.01007857
Na  0.5889571  0.06675857  0.02761571  0.01007857 0.19929286

有人能帮我做一下吗？提前谢谢。

matrix

manova

回答 1

Stack Overflow用户

回答已采纳

发布于 2018-06-22 16:12:04

我们需要根据结构提取。如果我们看一下str(anv)，它是一个由7个元素组成的list。"SSPE“、"SSPH”元素位于"multivariate.tests“的”站点“元素之下。因此，用$ of [[提取这两个矩阵，然后得到两个矩阵。

out <- anv$multivariate.tests$Site[c("SSPE", "SSPH")]
out
#$SSPE
#           Al          Fe          Mg          Ca         Na
#Al 48.2881429  7.08007143  0.60801429  0.10647143 0.58895714
#Fe  7.0800714 10.95084571  0.52705714 -0.15519429 0.06675857
#Mg  0.6080143  0.52705714 15.42961143  0.43537714 0.02761571
#Ca  0.1064714 -0.15519429  0.43537714  0.05148571 0.01007857
#Na  0.5889571  0.06675857  0.02761571  0.01007857 0.19929286

#$SSPH
#            Al          Fe          Mg         Ca         Na
#Al  175.610319 -149.295533 -130.809707 -5.8891637 -5.3722648
#Fe -149.295533  134.221616  117.745035  4.8217866  5.3259491
#Mg -130.809707  117.745035  103.350527  4.2091613  4.7105458
#Ca   -5.889164    4.821787    4.209161  0.2047027  0.1547830
#Na   -5.372265    5.325949    4.710546  0.1547830  0.2582456

与其拥有多个对象，不如将其保存在list中。

票数 1

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/50991759

复制

相似问题

问从r中manova分析中提取误差矩阵
EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问从r中manova分析中提取误差矩阵EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问从r中manova分析中提取误差矩阵
EN