r语言 - 计算组均值,同时排除每个个案的单个值



我有一个包含70个案例(研究参与者)的数据集。是否有一个函数可以计算这 70 个个案的平均值,以便每个个案不包括在分析中。这看起来像:

"mean for case x = (value(1) + ... value(n) - value(x))/n"

任何信息都会有所帮助。

您可以按照您的建议进行操作,并从总数中删除每个案例:

x <- c(1:10)
(sum(x) - x) / (length(x) - 1)
#[1] 6.000000 5.888889 5.777778 5.666667 5.555556 5.444444 5.333333 5.222222 5.111111 5.000000
mean(2:10)
#[1] 6
mean(1:9)
#[1] 5

编辑:更新以尝试解决评论中的后续问题:

set.seed(123)
df <- data.frame(group = rep(letters[1:3], each = 3), 
                 value = rnorm(9), stringsAsFactors = F)
df
#group       value
#1     a -0.56047565
#2     a -0.23017749
#3     a  1.55870831
#4     b  0.07050839
#5     b  0.12928774
#6     b  1.71506499
#7     c  0.46091621
#8     c -1.26506123
#9     c -0.68685285
df$loo_mean <- unlist(tapply(df$value, df$group, 
                      function(x) (sum(x) - x) / (length(x) - 1)))
df
  #group       value    loo_mean
#1     a -0.56047565  0.66426541
#2     a -0.23017749  0.49911633
#3     a  1.55870831 -0.39532657
#4     b  0.07050839  0.92217636
#5     b  0.12928774  0.89278669
#6     b  1.71506499  0.09989806
#7     c  0.46091621 -0.97595704
#8     c -1.26506123 -0.11296832
#9     c -0.68685285 -0.40207251
mean(df$value[2:3])
#[1] 0.6642654
mean(df$value[c(7,9)])
#[1] -0.1129683

这是一种矢量化方法,以避免一次平均每个子集:

x <- runif(70)
sapply(seq_along(x), function(i) mean(x[-i]))

相关内容

  • 没有找到相关文章

最新更新