自定义R汇总函数中q1/q3分位数结果格式修正问题
解决q1、q3列带标签向量的问题
方案一:修改自定义汇总函数
直接在计算q1、q3时将结果转为纯数值,修改funcio_resum函数中的对应行即可:
funcio_resum <-function(x) { n<-NROW(x) mean<-mean(x) median <- median (x) min<-min (x) max<-max(x) # 将quantile结果转为纯数值 q1 <- as.numeric(quantile(x, .25)) q3 <- as.numeric(quantile(x, .75)) SD<- sd(x) IQR <- IQR(x) quantile <- quantile(x, probs=seq(.1,.9, by =.1)) hist<-hist(x) summary<-list(n=n,mean=mean, median = median, max=max,min=min, q1 = q1, q3 = q3, SD=SD, IQR = IQR) return(summary) }
修改后重新运行tapply和do.call(rbind, resum),q1、q3列就会显示纯数值。
方案二:对已生成的结果做后续处理
如果不想修改原函数,可以在得到resum后,单独转换q1、q3列的数值:
# 转换q1列 resum[,"q1"] <- sapply(resum[,"q1"], as.numeric) # 转换q3列 resum[,"q3"] <- sapply(resum[,"q3"], as.numeric)
这样处理后,原本带标签的向量就会被提取出纯数值,替换原列内容。
原理说明
quantile()函数默认返回带分位数命名的向量(比如25%对应q1),as.numeric()会自动去除标签只保留数值部分;用sapply遍历每一行的q1/q3元素,可批量完成分组结果的转换。
内容的提问来源于stack exchange,提问作者FrancoBattiato
相关产品推荐
相关产品推荐

