[R-br] encontrar resposta no tempo anterior

Elias Teixeira Krainski eliaskrainski em yahoo.com.br
Quarta Fevereiro 27 12:45:03 BRT 2013


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





Mais detalhes sobre a lista de discussão R-br