I'm trying to simulate trend data over a five year period.  I want
different trend profiles...the simplest being a linear trend.  I've been
using the following code:

 

patBdta1 <- NULL

for(i in 1:100)

patBdta1 <- rbind(patBdta1,c(yr1= mean(rbinom(50,1,.50)),

                             yr2 =mean(rbinom(50,1,.51)),

                             yr3 =mean(rbinom(50,1,.52)),

                             yr4 =mean(rbinom(50,1,.53)),

                             yr5 =mean(rbinom(50,1,.54))))

 

This code creates 100 data sets each with a 5 yr binomial trend profile
with a slope of approximately .01.

 

Now, what I want to do is pass this code (or some code) in such a way
that I can simulate various trend slopes (i.e., pass in a loop or
vectorize vs. copying/repeating code for each slope).

 

Any guidance is much appreciated.  Thanks!

 

Marc

 

Marc W. Zodet, MS

Senior Health Statistician

Agency for Healthcare Research and Quality


	[[alternative HTML version deleted]]

