R软件的t检验，无原始数据怎么办？

已知均数、标准差、样本量，未知原始数据。

## 建立公式

T.test <- function(n, mean, sd) {

s <- sum((n – 1) * sd^2) / (sum(n) – 2) # weighted variance

t <- sqrt(prod(n) / sum(n)) * (diff(mean) / sqrt(s)) # t statistic

df <- sum(n) – 2 # degrees of freedom

p <- (1 – pt(abs(t), df)) * 2 # p value

c(t = t, p = p)

}

## 建立数据框，输入均数、标准差、样本量

dat <- data.frame(mean = c(14.32, 15.67), sd = c(1.22, 1.39), n = c(20, 21))

dat

## 分析结果

T.test(dat$n, dat$mean, dat$sd)