「R」如何计算几何平均数
刚遇到一个有意思的问题，如何用R计算几何平均数。如果数字少，简单，计算很容易，直观上，先用prod函数连乘，然后开方即可。
但我的数值大，连乘几十个之后R结果就是INF了，然后开方就还是INF，算不出来！
聪明人就会动脑筋了，转个弯，先取对数，再指数化！Stackoverflow上的解答让我大开眼界，下面给一个通用的计算函数：
gm_mean = function(x, na.rm=TRUE, zero.propagate = FALSE){
if(any(x < 0, na.rm = TRUE)){
return(NaN)
}
if(zero.propagate){
if(any(x == 0, na.rm = TRUE)){
return(0)
}
exp(mean(log(x), na.rm = na.rm))
} else {
exp(sum(log(x[x > 0]), na.rm=na.rm) / length(x))
}
}最后一个参数指定是否容忍0的存在。
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2020-06-30，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读
腾讯云开发者
