Tengo una tabla que va de la siguiente forma
dx<-data.frame("clinica"<-as.factor(c("A","A","A","A","A","A","A","A","B","B","B","B","B","B","B","B")),
"Diabetes"<-c(0,1,0,1,1,1,0,0,1,1,1,1,1,0,0,1),
"Hipertension"<-c(1,1,1,0,1,0,0,0,1,1,1,1,1,0,1,0))
Necesito calcular el total de casos con diabetes y hipertensión para cada clinica y el porcentaje por clinica. Lo he intentado de la siguiente forma
df5<-df %>%
group_by(clinica)%>%
summarize(Diabetes=sum(Diabetes), Hipertension=sum(Hipertension))
df_5<-df5 %>%
pivot_longer(c("Diabetes","Hipertension"), names_to = "Comorbilidades", values_to = "valor")
df_5<-df_5 %>%
group_by(sede_dx)%>%
mutate(per=paste0(round(valor/sum(valor)*100, 2), "%")) %>%
ungroup
El problema es que calcula el porcentaje de la suma de las dos enfermedades y no del total de casos, y como una persona puede tener las dos enfermedades, pues me da un resultado erróneo, les agradezco mucho la ayuda