У меня есть кадр данных, который выглядит так:
#df
ID DRUG FED AUC0t Tmax Cmax
1 1 0 100 5 20
2 1 1 200 6 25
3 0 1 NA 2 30
4 0 0 150 6 65
Ans и так далее. Я хочу обобщить некоторые статистические данные по AUC, Tmax и Cmax препаратом DRUG
и FED STATUS FED
. Я использую dplyr. Например: для AUC:
CI90lo <- function(x) quantile(x, probs=0.05, na.rm=TRUE)
CI90hi <- function(x) quantile(x, probs=0.95, na.rm=TRUE)
summary <- df %>%
group_by(DRUG,FED) %>%
summarize(mean=mean(AUC0t, na.rm=TRUE),
low = CI90lo(AUC0t),
high= CI90hi(AUC0t),
min=min(AUC0t, na.rm=TRUE),
max=max(AUC0t,na.rm=TRUE),
sd= sd(AUC0t, na.rm=TRUE))
Однако вывод не сгруппирован с помощью DRUG и FED. Он дает только одну строку, содержащую статистику всех, не разбитую на DRUG и FED.
Любая идея, почему? и как я могу заставить его поступать правильно?