Question

我的数据集如下：

library(dplyr)
salary_old<-c(100,200,300,400,10000,100,10,20,30)
salary_new<-c(200,300,400,500,230,240,30,40,50)
d<-as.Date(c('2019-01-01','2019-01-02','2019-01-03'))
country<-c('USA','UK','IR')
id<-c('A','B','A')
data<-data.frame(id,country,d,salary_new,salary_old)
data<-data %>% arrange(id,country,d)

然后我想使用循环或套用函数来计算每个组的新旧工资的T.test，以检查每个组的p值是否小于0.001。

我写的代码如下：

    datlist <- split(data ,list(data$id , data$country) )
 datlist<-datlist[sapply(datlist, nrow)>0]

 results<-  lapply(
  1:length(datlist) ,
  FUN = function(x) {
   t.test(datlist[[c(x,4)]] , datlist[[c(x,5)]])
  })

x<-matrix(nrow=3,ncol=1)
for (i in 1:length(results)){
   x[i]<-results[[i]]$p.value
  x[i]<-(sum(x[i]<0.001))

}

我获得了每个组的p.value，但是我想将其连接到主数据框，以了解与id / country / combination关联的p_value是多少：

id   country  p-value         status
A     USA      0.5417366       0
A     IR       0.4321609       0
B     UK       0.7066187       0

Answer 1

如果我正确理解了您的问题，则希望在操作结束时将所有结果合并到一个data.frame中。首先，您为每个t检验构建迷你data.frame，然后再构建 rbind个。

您可以使用此：

results <- lapply(
    datlist,
    FUN = function(x) {
      tres <- t.test(x[["salary_new"]], x[["salary_old"]])
      lt05 <- sum(tres$p.value < 0.05)
      data.frame(x[1L, "id", drop = FALSE],
          x[1L, "country", drop = FALSE],
          pval = tres$p.value,
          status = lt05)
    }
)

do.call(rbind, unname(results))

注意：问题中的显着性水平为0.05。

如何计算R中数据帧中每个组的t检验

1 个答案: