### 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)

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.