0

我有两个单独的列表，每个列表包含 4 个 data.frame。我需要在两个列表中的每个 data.frames 之间执行学生 t 检验（ t.test） 。rainfall

这里的列表：

lst1 = list(data.frame(rnorm(20), rnorm(20)), data.frame(rnorm(25), rnorm(25)), data.frame(rnorm(16), rnorm(16)), data.frame(rnorm(34), rnorm(34)))    
lst1 = lapply(lst1, setNames, c('rainfall', 'snow'))

lst2 = list(data.frame(rnorm(19), rnorm(19)), data.frame(rnorm(38), rnorm(38)), data.frame(rnorm(22), rnorm(22)), data.frame(rnorm(59), rnorm(59)))
lst2 = lapply(lst2, setNames, c('rainfall', 'snow'))

我需要做的是：

t.test(lst1[[1]]$rainfall, lst2[[1]]$rainfall)
t.test(lst1[[2]]$rainfall, lst2[[2]]$rainfall)
t.test(lst1[[3]]$rainfall, lst2[[3]]$rainfall)
t.test(lst1[[4]]$rainfall, lst2[[4]]$rainfall)

我可以通过编写 4 个 data.frame 中的每一个来做到这一点（我实际上有 40 个带有我的真实数据），但我想知道是否存在更智能、更快捷的方法来做到这一点。

下面是我尝试过的（没有成功）：

myfunction = function(x,y) {
  test = t.test(x, y)
  return(test)
}

result = mapply(myfunction, x=lst1, y=lst2)
4

2 回答 2

1 
x <- NULL
for (i in seq_along(lst1)){
  x[[i]] <- t.test(lst1[[i]]$rainfall, lst2[[i]]$rainfall)
}
x
于 2017-02-11T23:56:18.480 回答
0

为我工作。不过，我会用它simplify = FALSE来更好地格式化结果。

lst1 <- list()
lst1[[1]] <- data.frame(rainfall = rnorm(10))
lst1[[2]] <- data.frame(rainfall = rnorm(10))
lst2 <- list()
lst2[[1]] <- data.frame(rainfall = rnorm(10))
lst2[[2]] <- data.frame(rainfall = rnorm(10))

myfunction = function(x,y) {
  test = t.test(x$rainfall, y$rainfall)
  return(test)
}

mapply(myfunction, x = lst1, y = lst2, SIMPLIFY = FALSE)
于 2017-02-11T22:54:30.460 回答