如何解决R中的卡方结构
我将平均光谱分为两组,使用卡方检验测试两组之间差异的显着性。我在R中构建了测试,得出的值约为468,这似乎太大了。每行代表所有光谱在波长下的测量值。这是代码:group1 <- dt[,c("p2","p8","p9","p10")]
group2 <- dt[,c("p1","p3","p4","p5","p6","p7","p11","p12","p13","p14")]
group1$Mean <- rowMeans(group1)
group2$Mean <- rowMeans(group2)
for (i in 1:nrow(dt)){
group1$sd[i] <- sd(dt[i,2:5])}
for (i in 1:nrow(dt)){
group2$sd[i] <- sd(dt[i,2:10])}
chi <- c()
for (i in 1:nrow(dt)){
chi[i] <- ((group1$Mean[i] - group2$Mean[i])^2) / sqrt(group1$sd[i]^2 + group2$sd[i]^2)}
sum(chi)
问题出在哪里?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。