DATA Step, Macro, Functions and more

How to compute incidence rate and confidence interval

Accepted Solution Solved
Reply
Contributor
Posts: 23
Accepted Solution

How to compute incidence rate and confidence interval

I have a dataset from a longitudinal study to compute incidence rate and the confidence interval.I need to compute the number of new cases(the number of the variable 'case'=1) divided by the sum of the total person-years.Need your kindly help on the coding.Thanks!

The example dataset is like

case    pyr

1          3

0          5

0         5

1         3.4

0         5

0        5

0        5

1       2.1

0       5


Accepted Solutions
Solution
‎10-08-2013 09:51 PM
Respected Advisor
Posts: 4,919

Re: How to compute incidence rate and confidence interval

Posted in reply to jing2000yr

Assuming you want to sum the pyr only for the new cases :

data have;
input ase    pyr;
datalines;
1          3
0          5
0         5
1         3.4
0         5
0        5
0        5
1       2.1
0       5
;

data want(keep=newCases sumPyr ratio);
do until(done);
     set have end=done;
     if ase then do;
          newCases + 1;
          sumPyr + pyr;
          end;
     end;

ratio = newCases / sumPyr;
run;

proc print data=want noobs; run;

PG

PG

View solution in original post


All Replies
Solution
‎10-08-2013 09:51 PM
Respected Advisor
Posts: 4,919

Re: How to compute incidence rate and confidence interval

Posted in reply to jing2000yr

Assuming you want to sum the pyr only for the new cases :

data have;
input ase    pyr;
datalines;
1          3
0          5
0         5
1         3.4
0         5
0        5
0        5
1       2.1
0       5
;

data want(keep=newCases sumPyr ratio);
do until(done);
     set have end=done;
     if ase then do;
          newCases + 1;
          sumPyr + pyr;
          end;
     end;

ratio = newCases / sumPyr;
run;

proc print data=want noobs; run;

PG

PG
Contributor
Posts: 23

Re: How to compute incidence rate and confidence interval

Thanks!

🔒 This topic is solved and locked.

Need further help from the community? Please ask a new question.

Discussion stats
  • 2 replies
  • 268 views
  • 0 likes
  • 2 in conversation