【问题标题】:t.test from summary statistics in R?t.test 来自 R 中的汇总统计数据?
【发布时间】:2014-02-02 19:10:59
【问题描述】:

在 R 中,prop.test 可以采用汇总统计表(交叉表中的单元格计数)而不是所有数据(参见参数 x here):

x - 成功计数的向量,具有两个条目的一维表,或具有 2 列的二维表(或矩阵),给出成功和失败的计数,分别。

有没有办法给t.test(或一些等效的函数)两个组的平均值、标准差和组大小的汇总统计数据,并让它执行适当的计算?

【问题讨论】:

  • 为什么这个问题被否决了?我对答案很感兴趣。
  • 像什么?如果您有组的均值/方差,只需使用适当的公式计算统计量。并且t.test 将参数x 作为数据,就像prop.test 一样。包含成功和失败的表格不是汇总统计。
  • 您能解释一下您的问题背后的原因吗?你想达到什么目的?
  • 我假设您有每个组的均值、标准差和样本量。我不知道这是否已经写在某个地方,但是通过实现计算合并方差等的基本规则很容易编写;您将需要 pt 来计算 p 值。 stats:::t.test.default 拥有所有必要的机制,但很复杂,因为它处理许多不同的情况并进行大量输入检查。
  • @Glen_b 如果您将其作为答案提交,我会接受。

标签: r


【解决方案1】:

CRAN 已经至少有一个包执行 t 检验的汇总统计版本,请参阅包 BSDA 中的函数 tsum.test

This question on stats.SE 包含几组 R 代码也应该足够了。

注意事项:

那里的答案似乎没有利用htest 类。我可能会保持一致性,尽管这不是一个大问题。

我实际上喜欢看到的是t.test 能够混合使用(两个样本都有观察或两个样本都有总结,或者每个样本一个,我不认为tsum.test 也不会),我希望看到 t.test 正确处理只有一次观察的一个样本的情况(当然只有var.equal=TRUE) ;就目前而言,所有这些添加都很简单,而且并不经常需要,令我惊讶的是,少量的附加行并不是简单地插入的。

(有趣的是,这个问题从你的问题转移到另一个方向,我认为它们非常相似;在临界情况下,我想会有一些问题被迁移到迁移另一端的人think 应该保持原状。可能是我对您的问题以及您的回答实际上是触发这种情况的原因。)

【讨论】:

猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2020-12-03
  • 2012-06-04
相关资源
最近更新 更多