首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >R函数对样本均值、SD和n进行方差分析和TukeyHSD

R函数对样本均值、SD和n进行方差分析和TukeyHSD
EN

Stack Overflow用户
提问于 2015-03-25 15:39:24
回答 1查看 1.9K关注 0票数 0

我想知道是否有从样本均值,SD和样本大小(n)进行方差分析和TukeyHSD的R函数。在样本统计中,我从ind.oneway.second中找到了rpsychi函数来进行方差分析,但在样本统计上却找不到执行TukeyHSD的函数。

代码语言:javascript
复制
library("rpsychi")

Mean <- c(90,85,92,100,102,106)
SD <- c(9.035613,11.479667,9.760268,7.662572,9.830258,9.111457)
SampleSize <- c(9,9,9,9,9,9)

fm1 <- ind.oneway.second(Mean, SD, SampleSize)
fm1
names(fm1)
TukeyHSD(fm1$anova.table)
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2015-03-25 16:12:01

这应该能行。基本上,它只是模拟与给定平均值、标准差和样本大小完全匹配的数据。在此之后,您可以对感兴趣的模型进行拟合,并进行您感兴趣的任何特殊测试。

代码语言:javascript
复制
Mean <- c(90,85,92,100,102,106)
SD <- c(9.035613,11.479667,9.760268,7.662572,9.830258,9.111457)
SampleSize <- c(9,9,9,9,9,9)


gen_data <- function(means, sds, samplesizes){
  n.grp <- length(means)
  grps <- factor(rep(1:n.grp, samplesizes))
  dat <- lapply(1:n.grp, function(i) {scale(rnorm(samplesizes[i]))*sds[i] + means[i]})
  y <- do.call(rbind, dat)
  out <- data.frame(group = grps, y = y)
  out
}

simulated_data <- gen_data(Mean, SD,SampleSize)
av <- aov(y ~ group, data = simulated_data)
summary(av)
TukeyHSD(av)
票数 4
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/29260139

复制
相关文章

相似问题

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