Question

这是df：

    # A tibble: 6 x 5
      t      a      b       c       d
  <dbl>   <dbl>   <dbl>   <dbl>   <dbl>
1 3999. 0.00586 0.00986 0.00728 0.00856
2 3998. 0.0057  0.00958 0.00702 0.00827
3 3997. 0.00580 0.00962 0.00711 0.00839
4 3996. 0.00602 0.00993 0.00726 0.00875

除了不包括第一列之外，我想获得所有行的均值。我写的代码：

df$Mean <- rowMeans(df[select(df, -"t")])

我得到的错误：

    Error: Must subset columns with a valid subscript vector.
x Subscript `select(group1, -"t")` has the wrong type `tbl_df<
  p2 : double
  p8 : double
  p10: double
  p9 : double
>`.
ℹ It must be logical, numeric, or character.

我尝试将df转换为矩阵，但随后出现另一个错误。我该怎么解决？

现在，我正在尝试使用以下代码计算标准误：

se <- function(x){sd(df[,x])/sqrt(length(df[,x]))}
sapply(group1[,2:5],se)

我尝试指出应该使用哪些列来计算错误，但又会弹出错误：

 Error: Must subset columns with a valid subscript vector.
x Can't convert from `x` <double> to <integer> due to loss of precision.

我使用了有效的列下标，所以我不知道为什么会出错。

Answer 1

类似的base R解决方案是：

df$Mean <- rowMeans(df[,-1],na.rm=T)

输出：

     t       a       b       c       d      Mean
1 3999 0.00586 0.00986 0.00728 0.00856 0.0078900
2 3998 0.00570 0.00958 0.00702 0.00827 0.0076425
3 3997 0.00580 0.00962 0.00711 0.00839 0.0077300
4 3996 0.00602 0.00993 0.00726 0.00875 0.0079900

Answer 2

我们可以使用setdiff返回不是't'的列，然后获取rowMeans。假设列“ t”可以在任何地方，而不是基于列的位置

df$Mean <- rowMeans(df[setdiff(names(df), "t")], na.rm = TRUE)
df
#     t       a       b       c       d      Mean
#1 3999 0.00586 0.00986 0.00728 0.00856 0.0078900
#2 3998 0.00570 0.00958 0.00702 0.00827 0.0076425
#3 3997 0.00580 0.00962 0.00711 0.00839 0.0077300
#4 3996 0.00602 0.00993 0.00726 0.00875 0.0079900

select中的

dplyr返回data.frame的子集，而不返回列名或索引。因此，我们可以直接应用rowMeans

library(dplyr)
rowMeans(select(df, -t), na.rm = TRUE)

或在管道中

df <- df %>%
         mutate(Mean = rowMeans(select(., -t), na.rm = TRUE))

更新

如果需要获取每行的标准错误，可以将apply与MARGIN一起使用

apply(df[setdiff(names(df), 't')], 1, 
            function(x) sd(x)/sqrt(length(x)))

或者使用rowSds中的matrixStats

library(matrixStats)
rowSds(as.matrix(df[setdiff(names(df), 't')]))/sqrt(ncol(df)-1)

数据

df <- structure(list(t = c(3999, 3998, 3997, 3996), a = c(0.00586, 
0.0057, 0.0058, 0.00602), b = c(0.00986, 0.00958, 0.00962, 0.00993
), c = c(0.00728, 0.00702, 0.00711, 0.00726), d = c(0.00856, 
0.00827, 0.00839, 0.00875)), class = "data.frame", row.names = c("1", 
"2", "3", "4"))

如何计算数据框中每一行的均值？ [R]

2 个答案:

更新

数据