grafico_dados longitudinais

Boa Noite Pessoal, Estou precisando muito de uma ajuda com graficos... Tenho um conjunto de dados longitudinais desbalanceados, onde cada amostra foi medida duas ou tres vezes (dependendo da amostra) com intervalos de 1 a 2 anos entre as medições. Eu gostaria de plotar em um gráfico X1~X2, de modo que os pontos que correspondem a mesma unidade amostral, porem em tempos diferentes, possam ser ligados, unidos um ao outro, mostrando o que aconteceu com a unidade amostral de um ano pro outro (aumentou X1 em função do aumento de X2, por exemplo). Lembrando que nem X1 nem X2 correspondem a variavel tempo. Quero fazer algo semelhante a figura que esta disponivel no seguinte link: http://www.datafilehost.com/download-e7370438.html Alguem sabe como fazer isso no R? Abraços Ayuni

Embora existam funcoes especificas em pacotes especificos aqui vao uma possivel solucao usanbdo apensa funções básucas do R foo <- data.frame(ind = c(1,1, 2,2,2, 3,3, 4, 5,5,5,5), tempo = c(1,3, 2,4,5, 1,4, 3, 1,2,3,5), dados = c(21, 23, 12, 14, 18, 23, 27, 13, 10, 11, 13, 19) ) with(foo, plot(dados ~ tempo, type="n")) with(foo, by(foo, ind, function(x) lines(x$dados ~ x$tempo, type="b"))) On Thu, 28 Apr 2011, Ayuni Sena wrote:
Boa Noite Pessoal, Estou precisando muito de uma ajuda com graficos... Tenho um conjunto de dados longitudinais desbalanceados, onde cada amostra foi medida duas ou tres vezes (dependendo da amostra) com intervalos de 1 a 2 anos entre as medições. Eu gostaria de plotar em um gráfico X1~X2, de modo que os pontos que correspondem a mesma unidade amostral, porem em tempos diferentes, possam ser ligados, unidos um ao outro, mostrando o que aconteceu com a unidade amostral de um ano pro outro (aumentou X1 em função do aumento de X2, por exemplo). Lembrando que nem X1 nem X2 correspondem a variavel tempo. Quero fazer algo semelhante a figura que esta disponivel no seguinte link: http://www.datafilehost.com/download-e7370438.html Alguem sabe como fazer isso no R? Abraços Ayuni

Ayuni, Ainda não estou seguro sobre o que você quer, mas aí um CMR da <- expand.grid(tempo=1:4, x=sort(runif(15,0,20))) da$y <- rnorm(da$x) require(lattice) xyplot(y~x, data=da, groups=tempo, type=c("p","l")) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (3)
-
Ayuni Sena
-
Paulo Justiniano
-
Walmes Zeviani