r - 如何将 na.rm 作为参数传递给 tapply？

Question

我想从数据框中计算均值和标准差，其中一列用于参数，一列用于组标识符。使用 tapply 时如何计算它们?我可以用 sd(v1, group, na.rm=TRUE) ，但不适合 na.rm=TRUE使用 tapply 时进入语句.omit.na是没有选择。我有一大堆参数，并且在排除具有一个缺失值的所有行时，必须逐步检查它们而不会丢失一半的数据帧。

data("weightgain", package = "HSAUR")
tapply(weightgain$weightgain, list(weightgain$source, weightgain$type), mean)

这同样适用于 by陈述。

x<-c(1,2,3,4,5,6,7,8,9,NA)
y<-c(2,3,NA,3,4,NA,2,3,NA,2)
group<-rep((factor(LETTERS[1:2])),5)
df<-data.frame(x,y,group)
df

by(df$x,df$group,summary)
by(df$x,df$group,mean)

sd(df$x) #result: NA
sd(df$x, na.rm=TRUE) #result: 2.738613

任何想法如何完成这项工作？

Answer 1

我认为这应该做你想做的。

选择你想要的列:

v = c("x", "y")#or
v = colnames(df)[1:2]

使用 sapply迭代 v并将值传递给 tapply :

sapply(v, function(i) tapply(df[[i]], df$group, sd, na.rm=TRUE))