set.seed(2018)
data=data.frame(
patient=1:10000,
primdiag=sample(390:398,10000,replace=TRUE),
onset=sample(seq.Date(as.Date("2015/2/1"),
as.Date("2016/2/1"),"1 day"),10000,replace=TRUE),
sex=sample(c("M","F"),10000,replace=TRUE),
county=sample(c("Albany","New York"),10000,replace=TRUE)
)
output.series=case.series(
data,ICD=392:396,diagnosis="primdiag",
date="onset",start="2015/1/1",end="2016/12/31",by1="sex")
head(output.series)
Run the code above in your browser using DataLab