Рассчитать стандартное отклонение по группе, исключая текущее наблюдение в R

Question

Легко рассчитать среднее значение каждой группы в кадре данных R.Если вы хотите исключить текущее наблюдение, это почти так же просто .

Есть ли простой способ исключить текущее наблюдение при расчете стандартного отклонения?

ДляНапример, когда у меня есть эта таблица

data.frame(country = c(rep("A",3), rep("B",3)), weight = c(10,11,12,20,25,30))

, мне нужна следующая таблица:

data.frame(country = c(rep("A",3), rep("B",3)), weight = c(10,11,12,20,25,30), standarddeviation = c(sd(c(11,12)), sd(c(10,12)), sd(c(10,11)), sd(c(25,30)), sd(c(20,30)), sd(c(20,25))))

Answer 1

Не очень красивое решение, но оно должно работать

library(dplyr)

data = data.frame(country = c(rep("A",3), rep("B",3)), weight = c(10,11,12,20,25,30))

cdata = list()

for(k in 1:length(unique(data$country))){
cdata[[k]] = filter(data,country==unique(country)[k])
}

for(i in 1:length(unique(data$country))){
  for(j in 1:nrow(cdata[[1]])){
    aux=cdata[[i]][-j,]
    cdata[[i]][j,"StandardDeviation"] = sd(aux$weight)
  }
}

rbind(cdata[[1]],cdata[[2]])

Answer 2

Можно использовать dplyr и mapply.mapply выполняется для каждой строки (группы), а sd вычисление исключает текущую строку.

library(dplyr)

df %>% group_by(country) %>%
  mutate(Sp_SD = mapply(function(x)sd(weight[-x]), 1:n()))


# # A tibble: 6 x 3
# # Groups: country [2]
# country weight Sp_SD
# <fctr>   <dbl> <dbl>
# 1 A         10.0 0.707
# 2 A         11.0 1.41 
# 3 A         12.0 0.707
# 4 B         20.0 3.54 
# 5 B         25.0 7.07 
# 6 B         30.0 3.54