Как выполнить парный t-тест без всей выборки данных?

Question

У меня есть несколько сводок данных для парного t-теста из трех разностных источников (больниц), что означает, что у меня нет полных данных выборки.То, что я должен выполнить для парного t-теста, - это среднее значение (разность данных каждой парной выборки).
Данные выглядят как таблица ниже (каждая строка - один пациент), и у меня есть mean(diff) и sd(diff).
Есть ли какая-нибудь функция, с помощью которой я мог бы легко выполнить парный t-тест и получить 95% CI?Из-за проблем с конфиденциальностью пациента просто невозможно получить весь набор данных.

3 months 6 months  diff
  1         3       -2
  2         1        1
  5         9       -4

Answer 1

Парный t-критерий в точности эквивалентен t-критерию, который проверяет среднее значение разностей против нулевой гипотезы нуля.

diffs <- c(-2,1,-4)
m <- mean(diffs)
s <- sd(diffs)
n <- length(diffs)

Теперь предположим, что вы не знаете diffs,просто m, s и n (все три необходимы)

pval <- 2*pt(abs(m/s),lower.tail=FALSE,df=n)
## 0.555
qq <- qt(0.975,df=n)
ci <- m + c(-1,1)*qq*s
## [1] -9.675648  6.342314