微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

R函数查找大于或等于特定数字的均值差

如何解决R函数查找大于或等于特定数字的均值差

我刚刚使用 R 开始了我的基础统计课程,我们正在学习使用 R 进行配对 t 检验。我遇到过这样的问题,我们得到两组数据,我们被要求找出均值的差异是等于 0 还是大于 0 等等。我们对方差未知的两个样本 x 和 y 使用的函数类似于下面的函数

t.test(x,y,var.equal=TRUE,alternative="greater")

我的问题是,如果我们想测试均值的差异是否大于或等于指定数字与小于特定数字而不是 0 的替代方案,我们将如何执行此操作。

例如,假设我们有两个数据,分别是 10 人前后的权重。我们如何测试体重的平均差异是否大于或等于 3 公斤,而不是体重平均差异小于 3 公斤的替代方案。有没有办法做到这一点?非常感谢有关此事的任何指导。

解决方法

如果您需要更多理论证明,也可能值得在 https://stats.stackexchange.com/ 上发帖。是否可以从 xy 中添加/减去 3kg,然后使用 t 检验来检查相似性?我认为这至少会告诉你哪种结果更有可能,如果那是最终目标。获得有关此方面的反馈会很好

# number of obs,and rnorm dist for simulating
N <- 10
mu <- 70
sd <- 10
set.seed(1)
x <- round(rnorm(N,mu,sd),1)

# three outcomes
# (1) no change
y_same <- x + round(rnorm(N,5),1)
# (2) average increase of 3
y_imp <- x + rnorm(N,3,5)
# (3) average decrease of 3
y_dec <- x + rnorm(N,-3,5)

# say y_imp is true
y_act <- y_imp
# can we test whether we're closer to the output by altering
# the original data? or conversely,altering y_imp
t_inc <- t.test(x+3,y_act,var.equal=TRUE,alternative="two.sided")
t_dec <- t.test(x-3,alternative="two.sided")

t_inc$p.value
[1] 0.8279801
t_dec$p.value
[1] 0.0956033

# one with the highest p.value has the closest distribution,so
# +3 kg more likely than -3kg
,

您可以设置 mu=3 以将原假设从 0 更改为 3,假设您的 x 变量采用上述单位。

t.test(x,y,mu=3,alternative="greater",paired=TRUE)

有关 Stack Exchange 的更多(一般)信息 [此处]。(https://stats.stackexchange.com/questions/206316/can-a-paired-or-two-group-t-test-test-if-the-difference-between-two-means-is-l/206317#206317)

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。