27 Fev
2013
27 Fev
'13
15:45
### complementando porque vc tem individuos... n.times <- 5; n.ind <- 3 dat <- data.frame(expand.grid(time=1:n.times, indiv=1:n.ind)) dat$resp <- as.vector(replicate(n.ind, arima.sim(list(ar=.9),n.times)>0)+0) str(dat) lagFun <- function(x) cbind(ant1 = c(NA, x[-length(x)]), ant2 = c(NA, NA, x[-(-1:0+length(x))])) dat <- data.frame(dat, do.call('rbind', tapply(dat$res, dat$ind, lagFun))) str(dat)