我在创建双尾双样本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)
结果应该是:
我没有得到这些结果。如何在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值太小。
答案 0 :(得分:1)
基本上,你根本没有真正走上正轨。我认为这是作业作业?在现实生活中,您几乎总是从原始数据开始,这是R所期望的。但是在这里你没有原始数据,你只有摘要统计数据。因此,您需要自己进行计算才能得到答案。 (在您的编辑中,您正在制作看起来有点像真实数据的新原始数据 - 但当然不完全相同。)
因为这就像家庭作业一样,我不能在这里提供这些计算;相反,在文本中查找公式并遵循它们。基本思想是获得差异,除以标准误差(您需要从这两个标准偏差和样本大小计算),然后使用普通表(或R函数)来获得p -值。尝试一下,然后回过头来处理任何后续问题 - 编辑原始问题。
答案 1 :(得分:0)
检查?t.test
你会发现你需要交给它两个向量(x和y)。默认是测试平均差为零。双面测试也是默认测试。您可以切换是否假设等效方差 - var.equal = TRUE。
您报告p值的方式是非标准的。我不太清楚结果应该是什么。