400论坛币求如何用R实现Stata 配对t检验输出结果

这个其实没什么复杂，主要是你从t.test里面提取相应的数据合成一个数据框就好了，另外，我同意楼上的看法，你那个one sample t test看起来是实际上才是paired t test的结果。

1. library(dplyr)
   
2. library(purrr)
   
3. 
   
4. t_test_df <- function(dat){
   
5.   
   
6.   t_test <- apply(dat, 2, t.test)
   
7.   
   
8.   t_pair <- data.frame(obs = map_dbl(t_test, 2) + 1, # df + 1 is num of obs
   
9.                      avg = map_dbl(t_test, 5),
   
10.                      std_err = map_dbl(t_test, 7),
    
11.                      std_dev = apply(dat, 2, sd),
    
12.                      ci_low = map(t_test, 4) %>% map_dbl(., 1),
    
13.                      ci_up = map(t_test, 4) %>% map_dbl(., 2)
    
14.                      )
    
15.   
    
16.   t_pair_info <- data.frame(t = as.numeric(t_test$d[1]),
    
17.                             df = as.numeric(t_test$d[2]),
    
18.                             less = as.numeric(t_test$d[3])/2,
    
19.                             two_sided = as.numeric(t_test$d[3]),
    
20.                             greater = 1 - as.numeric(t_test$d[3])/2)
    
21.   
    
22.   # should be result on t_test$x1, t_test$x2
    
23.   t_one <- t_pair[3,]
    
24.   t_one_info <- t_pair_info
    
25.   
    
26. }

复制代码