Извлечение медианы и доверительных интервалов из нескольких переменных в отдельную таблицу в R - PullRequest
0 голосов
/ 07 ноября 2018

У меня есть data.frame в R, который выглядит так:

ID       A       B     C      D
1  31.63237      6     24    10.6
2       NA      1     62    10.4
3 18.04803      5      2     7.9
3 26.73194     24     19     7.6

Как я могу извлечь только срединные и 95% доверительные интервалы в отдельную таблицу? Это мой желаемый вывод

    A       B     C       D

Медиана 26,73201 16 35 8,35

5% 18,48319 19,39076 18,61792
95% 40,39744 39,89518 41,60440

1 Ответ

0 голосов
/ 07 ноября 2018

Вы можете сделать:

df <- read.table(header=TRUE, text=
"ID       A       B     C      D
1  31.63237      6     24    10.6
2       NA      1     62    10.4
3 18.04803      5      2     7.9
3 26.73194     24     19     7.6")

Q <- sapply(df[-1], quantile, probs=c(0.5, 0.05, 0.95), na.rm=TRUE)
Q
# > Q
#            A    B     C      D
# 50% 26.73194  5.5 21.50  9.150
# 5%  18.91642  1.6  4.55  7.645
# 95% 31.14233 21.3 56.30 10.570

Диапазон от 5% до 95% составляет 90% доверительного интервала.
Вы можете установить имена строк результата.

rownames(Q) <- c("Median", "05%", "95%")
Q
# > Q
#               A    B     C      D
# Median 26.73194  5.5 21.50  9.150
# 05%    18.91642  1.6  4.55  7.645
# 95%    31.14233 21.3 56.30 10.570
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...