안녕하세요,, 항상 많은 도움을 받고 있습니다. 

약제의 각 그룹에서 non-response rate을 우선 cumulate incidence로 산출하였으며,

이에 대한 incidence rate per 100 PY 및 reference 그룹대비 rate ratio를 산출하고자 합니다.

incidence rate per 100 PY에 대해서는 아래와 같이 수치를 산출하였습니다. 

데이터 %>% group_by(drug) %>% summarise(N = n(), event = sum(NP), rate = event/ N*100, person_year =sum(yr), incidence = event /person_year*1000)

다만 제가 궁금한 사항은 95% CI를 산출하는 것인데요. 

혹시 어떤 방법으로 R에서 이를 산출할 수 있을지 문의드립니다.