Sobrevivência 3 meses
Análises de Sobrevivência (16_abr_2011)
No caso dos dados de sobrevivência tem que usar análises específicas, que tem como variável resposta o “failure-time”. Referências: Scheiner & Gurevitch e R-Book.
Primeiro contato com as análises de sobrevivência:
Primeiros gráficos (linha cheia:RA; linha tracejada: RB)
27 março 2011
A minha dúvida é sobre qual análise usar para dados coletados em diferentes tempos. (sobrevivência, altura, número de folhas…)
Relembrando… Pergunta: Desempenho das espécies durante a fase de plântula difere entre RA e RB?
H0: Sobrevivência não difere entre RA e RB
HA: espécies dominantes como adulto em um ambiente terão menor mortalidade (maior sobrevivência) de plântulas neste ambiente em relação ao outro. P.e. plântulas de espécies dominantes na RB como adultos sobreviveram melhor na RB em relação à RA.
Já tenho os dados de 3 meses de coleta semanal de sobrevivência média de cada parcela para 2 sps (N inicial= 15 plântulas de cada sp em cada parcela), sendo 10 parcelas na RA e 10 na RB.
Se considerarmos apenas a sobrevivência final a análise é uma ANOVA simples (variável resposta contínua: sobrevivência; variável preditora categórica: ambiente, com dois níveis: RA e RB). Mas como os dados foram coletados ao longo do tempo não devemos jogar fora esses dados e usar só a sobrev. final, certo? (ver figuras)
Pelo que eu estive lendo tem 2 soluções:
1) fazer um modelo linear para cada ambiente e comparar as inclinações (linhas tracejadas nas figuras). Mas isso está mais relacionado com a “velocidade” com que as plântulas morrem, certo?
2) Usar ANOVA com medidas repetidas. Mas como faço isso no R?? A idéia é semelhante à uma ANOVA split-plot, certo? (ver script) Não sei bem como construir o modelo nesse caso…
Além disso, no caso da sobrevivência a minha variável resposta é a sobrevivência média por parcela, que se baseia no N inicial, não sei se essa é a melhor forma de representar a variável resposta. E o fato de ser %? Precisa fazer transformações?
Arquivos: