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 > [email protected] > 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. _______________________________________________ R-br mailing list [email protected] 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.
