Solved
Contributor
Posts: 23

# 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
Posts: 5,529

## Re: How to compute incidence rate and confidence interval

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

All Replies
Solution
‎10-08-2013 09:51 PM
Posts: 5,529

## Re: How to compute incidence rate and confidence interval

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.