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

数据如下：
x1        x2        d
10        12        -2
11        13        -2
12        15        -3
13        14        -1
14        15        -1
15        17        -2
16        19        -3


数据也可见附件。


输出结果如下：

最好在R中输出dataframe格式个结果，然后以网页的格式呈现结果。

请将code以附件出售的形式上传，谢谢！

aaa.xls (9.5 KB)

2021-1-31 12:44:08 上传

溜了溜了，营养跟不上。。你给T检验杠上了

乐天天12300 发表于 2021-1-31 16:23
溜了溜了，营养跟不上。。你给T检验杠上了

###set work directory
setwd ("文档路径") ###路径用("C:/Users/.....") ---是“/"不是”\"
getwd()
###install.packages("readxl")
library(readxl)
###install.packages("tidyverse")
library(tidyverse)
###nstall.packages("ggpubr")
library(ggpubr)
#install.packages("rstatix")
library(rstatix)

###get data from the file
mydata<-read_excel('aaa.xls')
str(mydata)

## Transform wide data frame to long data frame
mydata.long <- mydata%>%
  gather(key = "group", value = "value", x1, x2)

## Only get all rows for x1 and x2
transform_data<-mydata.long[,2:3]
## Read data
transform_data

# Compute t-test
results <-t.test(value ~ group, data = transform_data, paired = TRUE)
##Show results
results

#### One-sample t-test
mean(transform_data$value)##calculate mu value
one_sample<-t.test(transform_data$value,mu=14)
one_sample

试一下这个吧。
不过你的one-sample t-test应该不是那样算的。

joven123 发表于 2021-2-1 09:07
###set work directory
setwd ("文档路径") ###路径用("C:/Users/.....") ---是“/"不是”\"
getwd()

能否将以上代码写入一个函数中，函数参数为输入数据集，变量1，变量2，输出。

继续请教！