====== Ocotea pulchella ===== ===== Planilhas ===== {{rgr_ocotea_final.xls|}} {{rgr_ocotea.xls|Dados tempo zero para todos indiv}} {{dados_regressao.xls|Peso seco para regressão inicial}} {{dados_regressao.txt| Peso seco inicial formato txt}} ===== Análises ===== ==== Regressão inicial ==== **Objetivo** -> estimar massa inicial das plantas que serão secas com 6 meses de idade [[ Dúvidas/respostas do Alê]] ==== Análises iniciais ==== {{regressao_ocotea.txt|script R}} {{dados_regressao.txt|}} - Distribuição dos dados de peso seco : Não parece normal {{massa_ocotea_hist.jpg?300|}} - Plotando contra distribuição normal: parece ter mesmo desvio da normalidade {{qq_ocotea.jpg? 300|}} - Se transformar dados para **escala log**: parece ter menos desvio! {{qq_ocotea_log.jpg? 300|}} ==== Regressão simples massa seca x altura ==== primeiro sem transformar os dados: {{massa_altura.jpg? 400|}} **Modelo:** m1<-lm(tab$massa_tot_mg~tab$altura) summary(lm(tab$massa_tot_mg~tab$altura)) Call: lm(formula = tab$massa_tot_mg ~ tab$altura) Residuals: Min 1Q Median 3Q Max -41.161 -14.571 -3.247 8.264 72.796 Coefficients: ||Estimate Std.| Error| t value |Pr(>t)| |(Intercept)| 17.503| 15.610| 1.121| 0.27109| |tab$altura| 17.915| 4.487| 3.992| 0.00039 | Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 25.7 on 30 degrees of freedom Multiple R-squared: 0.3469, **Adjusted R-squared: 0.3252 ** F-statistic: 15.94 on 1 and 30 DF, p-value: 0.0003901 **Gráficos diagnósticos:** {{diagnostico_m1.jpg? 500|}} Pelo grafico 1 (residuos x valores observados) -> relaçao nao linear! Pelo grafico 2(grafico QQ) -> desvio da normalidade nos resíduos ==== Regressão simples transformando para log ==== {{massa_altura_log.jpg? 400|}} **Modelo** : m2<-lm(log(tab$massa_tot_mg)~log(tab$altura)) summary(m2) Call: lm(formula = log(tab$massa_tot_mg) ~ log(tab$altura)) Residuals: Min 1Q Median 3Q Max -0.586234 -0.215341 -0.004106 0.154675 0.681365 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 3.3324 0.2007 16.607 < 2e-16 *** log(tab$altura) 0.8115 0.1676 4.841 3.65e-05 *** Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 0.3026 on 30 degrees of freedom Multiple R-squared: 0.4385, **Adjusted R-squared: 0.4198** F-statistic: 23.43 on 1 and 30 DF, p-value: 3.652e-05 **Gráficos diagnósticos**: {{diagnostico_m2.jpg ? 500|}} ==== Numero de folhas x massa seca ==== {{massa_nfolha.jpg|}} numero de folhas é uma variavel categorica, nao da para usar regressao, certo?