Если вы хотите, чтобы процент пациентов с высоким кровяным давлением в каждой группе ИМТ был приведен ниже:
library(dplyr)
bmi <- data.frame(bmi = c(rep('uw',5),rep('norm',5), rep('over',5),rep('obese',5)),
hpb = rbinom(20,1,.5) , stringsAsFactors = FALSE)
bmi2 <- bmi %>%
group_by(bmi, hpb) %>%
tally %>%
group_by(bmi) %>%
mutate(pct= n/ sum(n))
Если вы хотите, чтобы общий процент пациентов с повышенным кровяным давлением здесь представлял собой код:
bmi2 <- bmi %>%
group_by(bmi, hpb) %>%
tally %>%
ungroup %>%
mutate(pct= n/ sum(n))
В приведенном выше коде вы пытаетесь разделить на n()
, но если вы хотите использовать свой код, возможно, замените n()
только на переменную n
.