[R-br] encontrar resposta no tempo anterior

[Elias Teixeira Krainski]

n <- 1000
x <- (arima.sim(list(ar=.9), n)>0) + 1
ant1 <- c(NA, x[-n])
ant2 <- c(NA, NA, x[-(-1:0+n)]) 
cor(cbind(x, ant1, ant2), u='p')


On Wed, 2013-02-27 at 11:43 -0300, Maurício Lordêlo wrote:
> Saudações a todos desta lista.
> Estou precisando de uma ajuda no sentido de definir uma nova variável
> em função de outra já existente.
> No exemplo ilustrativo, tenho uma variável resposta binária
> (denominada de "resp") observada em quatro indivíduos (id) e medida em
> quatro momentos distintos (time).
> Preciso desta resposta no tempo anterior (inicialmente de primeira e
> segunda ordem). 
> A resposta prévia de primeira ordem (resp_ant1) indica o valor da
> resposta daquele indivíduo no tempo imediatamente anterior.
> De forma análoga, tem-se a resposta prévia de segunda ordem
> (resp_ant2). 
> Neste exemplo com poucas observações fica "fácil" visualizar estas
> respostas prévias e criar estas novas variáveis.
> Alguém teria uma sugestão de uma forma mais simples para encontrar as
> variáveis "resp_ant1" e "resp_ant2"?
> Desde já agradeço,
> Maurício 
> 
> 
> id=rep(1:4, each = 4) 
> time=rep(1:4, 4)
> resp=c(1,1,0,1,0,0,0,1,0,1,1,1,0,1,0,1)
> resp_ant1=c(NA,1,1,0,NA,0,0,0,NA,0,1,1,NA,0,1,0)
> resp_ant2=c(NA,NA,1,1,NA,NA,0,0,NA,NA,0,1,NA,NA,0,1)
> 
> 
> _______________________________________________
> R-br mailing list
> R-br em 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.