Muito obrigado! Problema resolvido. Abraço 2013/2/27 Elias Teixeira Krainski <eliaskrainski@yahoo.com.br>
### 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.