如何在R中执行双尾双样本t检验

时间:2018-04-28 16:52:10

标签: r t-test

我在创建双尾双样本t检验时遇到了问题。给出的唯一值是下面的值。我们有两组S和L,分别有中期和期末考试。我的任务是在学期中期分数上进行一次双尾双样本t检验,然后是从学期中期到最后一学期考试的分数差异。

d <- read.table(text="
Cohort  N Midterm_Mean Midterm_SD Final_Mean Final_SD Diff_Mean Diff_SD
2016_L 38     77.4        3.0       73.7       4.2      -3.7        2.1
2017_S 37     81.9        2.1       70.0       4.6     -11.9        2.8
", header=TRUE)

结果应该是:

  • 对于中期:( t:-1.24,p:0.219> 0.01,科恩的d:-0.286)
  • 为决赛:( t:2.30,p:0.024> 0.01,科恩的d:0.532)

我没有得到这些结果。如何在R中正确实现测试?

编辑:我已经有了以下内容:

L1 <- rnorm(mean = 77.4, sd = 3,   n=38)
S1 <- rnorm(mean = 81.9, sd = 2.1, n=37)
L2 <- rnorm(mean = 73.7, sd = 4.2, n=38)
S2 <- rnorm(mean = 70,   sd = 4.6, n=37)
tab <- data.frame(result=c(L1,S1), 
                  group=c(rep(0, times=38), rep(1,times =37)) ,head=TRUE)
attach(tab)
t <- t.test (result ~ group, mu=0, alternative="two.sided", 
             conf=0.95, var.equal=T, paired=F)

上面的代码仅用于比较中间条款。首先,我认为我必须创建具有给定值的法向量。我将它们放在数据框中并用0或1标记它们是否来自(S)或(L)。这或多或少是解决任务的正确方法吗?因为结果得到的p值太小。

2 个答案:

答案 0 :(得分:1)

基本上,你根本没有真正走上正轨。我认为这是作业作业?在现实生活中,您几乎总是从原始数据开始,这是R所期望的。但是在这里你没有原始数据,你只有摘要统计数据。因此,您需要自己进行计算才能得到答案。 (在您的编辑中,您正在制作看起来有点像真实数据的新原始数据 - 但当然不完全相同。)

因为这就像家庭作业一样,我不能在这里提供这些计算;相反,在文本中查找公式并遵循它们。基本思想是获得差异,除以标准误差(您需要从这两个标准偏差和样本大小计算),然后使用普通表(或R函数)来获得p -值。尝试一下,然后回过头来处理任何后续问题 - 编辑原始问题。

答案 1 :(得分:0)

检查?t.test

你会发现你需要交给它两个向量(x和y)。默认是测试平均差为零。双面测试也是默认测试。您可以切换是否假设等效方差 - var.equal = TRUE。

您报告p值的方式是非标准的。我不太清楚结果应该是什么。