O objetivo deste projeto é desenvolver e reunir metodologias estatísticas ligadas a problemas em Engenharia Florestal.
Modelos Lineares Generalizados, Testes Paramétricos e Não-Paramétricos, Análise de Sobrevivência, Análise de Séries Temporais, Modelos de Crescimento, Análise de Dados Longitudinais, Regressão por Quantis
Comparação entre duas famílias/tratamento (há dez séries temporais para cada família) . O objetivo é comparar o crescimento das famílias, distribuídas em 2 blocos (dois ambientes no local). Foram selecionadas 5 famílias de cada planta por bloco.
A luciana quer comparar estas duas famílias (seleção genética, pegou a melhor e a pior em 32 famílias). O que está em questão é medir a diferença entre o crescimento.
Variável dependente - Crescimento pelo dendômetro Variáveis independentes - Copa da árvore (não é medida frequentemente, mas pode ser medida de tempos em tempos),tensão de agua no solo a 30cm de profundidade e a 60cm de profundidade, temperatura, umidade do ar, precipitação(diária), insolação, velocidade do vento (talvez não interfira), dia de geada, estação do ano.
Dar uma olhada em modelos de função de transferência. delineamento original -32 famílias, 5 repetições
Dados sobre clima no Estado de Santa Catarina, locais : Ponte Serrada, Cacçador e Chapecó (não há média histórica). As variáveis analisadas são Temperatura, Número de Geadas, …
Primeiro Trabalho - Comparar observações climáticos dos locais acima, feita a partir de 2002, com a média histórica dos ultimos 30 anos.
Joel, estou enviando o arquivo com os dados para determinar as correlações existentes entre as inclinações diárias e os dados meteorológicos do ano de 2007. Neste arquivo tem duas planilhas, uma somente com os dados do ano de 2007 e a outra com os dados de 2007 e algumas informações de dezembro de 2006 para determinar as correlações com os dados meteorológicos dos meses anteriores. Qualquer dúvida me escreva.
Dados sobre crescimento de Pinus Taeda de acordo com diferentes tipos de tratamentos.
Este é um problema em que a determinação do plano amostral e do tamanho de amostra estão restritos aos custos monetários e, principalmente, operacionais ligados à coleta de dados. As variáveis sob as quais desejamos retirar informações são características das raízes finas (fine roots) e grossas (coarse roots) de Pinus Taeda plantadas na área experimental da empresa COMFLORESTA. A principal questão operacional vinculada a retirada das amostras é explicada pela dificuldade de locomoção na área com o "trator", equipamento que arranca a árvore do solo e a coloca em suspensão para que seja analisada a arquitetura da raiz grossa.
A importância no levantamento destas informações está no fato de que a absorção de nutrientes é feita pelas raízes finas e a sustentabilidade no crescimento está vinculada à raiz grossa. As características das raízes finas são expressas pelo peso e forma, enquanto a arquitetura da raiz grossa é mais dificil de ser definida ( ver Danjon e Reubens, 2008 ).
O que chamamos de amostragem aqui, na verdade é uma sub-amostragem feita em cima de árvores localizadas numa área experimental.
Planos amostrais sugeridos :
Divisão da dissertação do Valdeci em capítulos.
O objetivo é verificar a influência dos fatores no desenvolvimento geral da árvore, caracterizado pelo diâmetro e altura.
Os tratamentos foram aleatorizados em 4 blocos com o objetivo de controlar a variação espacial pois há um leve declive no terreno. No total foram criadas 32 parcelas para receber os 8 tratamentos resultante das possíveis combinações entre os fatores controlados no experimento( Resp.plantio e Sist.prod).
Metodologia Estatística
Será aplicado o modelo ANOVA para dois fatores em blocos inteiramente casualizados, adotando-se algum procedimento para controlar o número de mortos nas parcelas e seu efeito espacial nos resultados, caso este efeito exista. A hipótese existente para este efeito espacial é a de que parcelas com menor número de plantas oferecem mais espaço para o desenvolvimento da planta.
Uma primeira abordagem para minimizar o efeito da mortalidade é considerar o número de mortos na vizinhança (no máximo podem haver 8 mortos na vizinhança) como uma covariável (sugestão PJ). Devemos ter cuidado para declaração do modelo no R.
As análises de variância serão feitas para diâmetro, altura e volume de forma UNIVARIADA. Como volume é função de diâmetro e altura, devemos discutir se há necessidade realmente de repetir a análise para esta variável.
Definir alguns contrastes de interesse (Por exemplo: média dos tubetes x raiz nua)
Para o capítulo 2, serão analisadas as características das raízes finas. Para esta finalidade, foi selecionada uma árvore por parcela, totalizando 4 por tratamento. Foi selecionada a árvore no centro da parcela e, quando esta estava morta, selecionou-se a próxima na ordem de numeração crescente de acordo com o croqui do experimento.
Algumas hipóteses levantadas neste caso. Fatores ambientais podem levar a diferença entre os pesos das raízes finas nas diferentes posições(p1,p2,p3 e p4). Pode haver diferença nos pesos das raízes finas em função do sistema de produção pois espera-se que a raiz nua>tubete/126>tubete/55 10 meses > tubete 55 6 meses. Por hipótese, a raiz nua tem a tendência em ter melhor desenvolvimento pelas suas características de desenvolvimento no viveiro. Devemos verificar a interação entre a posição e os tratamentos (Será que existe ?????)
Matrizes de correlação para verificar se os pesos.
Definir alguns contrastes de interesse (Por exemplo: média dos tubetes x raiz nua)
Algumas dúvidas: (As observações nas posições são teoricamente correlacionadas espacialmente !!! Como superar isto ???)
No capítulo 3 foram amostradas, de maneira destrutiva, todas as árvores do primeiro bloco e cada ávore foi fotografada de 3 ângulos diferentes. Em cada ângulo, serão avaliadas as características da raiz grossa de acordo com escores pré-determinados. 64 árvores no total.
Verificar junto ao Cesar Taconeli a possibilidade de aplicar árvores de classificação multivariadas com o objetivo de identificar quais fatores levam aos escores atribuídos as raízes. Com esta metodologia, talvez consigamos encontrar a conjunção dos fatores que levem ao melhor arquitetura de raiz grossa.
Uma possibilidade (mais pobre) é aplicar testes não paramétricos para a nota em cada ângulo fotografado.
Neste capítulo, a primeira abordagem estatística consiste de utilizar uma técnica não paramétrica para encontrar alguma evidência de diferença entre os métodos e entre as empresas na constituição da raiz grossa. A justificativa para utilizar técnicas não-paramétricas é a característica da variável resposta que é qualitativa.
# Leitura dos dados
avalia<-read.csv2('http://www.leg.ufpr.br/~joel/dados/grossaclassif.csv')
# Resumo das variaveis
summary(avalia)
# Attachando os dados
attach(avalia)
# Nomes das variaveis
names(avalia)
# Carrega pacote para comparacoes multiplas nao parametricas
require(pgirmess)
# O gráfico de interacao é fundamental pois vai ser um instrumento para verificar
# a sua possivel existencia
interaction.plot(metodo,plantio,soma)
interaction.plot(plantio,metodo,soma)
# Alguns graficos exploratorios para entender melhor as interacoes
boxplot(soma[plantio=="empresa"]~metodo[plantio=="empresa"])
boxplot(soma[plantio=="terceiro"]~metodo[plantio=="terceiro"])
boxplot(soma[metodo=="M1"]~plantio[metodo=="M1"])
boxplot(soma[metodo=="M2"]~plantio[metodo=="M2"])
boxplot(soma[metodo=="M3"]~plantio[metodo=="M3"])
boxplot(soma[metodo=="M4"]~plantio[metodo=="M4"])
# Testes para verificar diferenças nos plantios dentro dos métodos
# Observação : como são dois niveis de plantio, o teste U de Mann-Whitney é um caso
# particular do kruskall-wallis e , portanto, não precisa fazer comparação multipla aqui
# Nomes dos metodos
m<-levels(metodo)
# loop que troca de metodos e testa diferencas entre os plantios
for (i in 1:4)
{
s<-soma[metodo==m[i]]
p<-plantio[metodo==m[i]]
print(paste("metodo",m[i]))
print(kruskal.test(s~p))
}
# aqui vamos trocar de plantios e verificar as diferenças entre os métodos
n<-levels(plantio)
# loop que troca de plantios e testa as diferenças entre os métodos...aqui
# já aproveito o embalo e faço as comparações multiplas com nível de significância
# de 10%
for (i in 1:2)
{
s<-soma[plantio==n[i]]
p<-metodo[plantio==n[i]]
print(paste("plantio",n[i]))
print(kruskal.test(s~p))
print(kruskalmc(s,p,prob=0.1))
}
Artigo sobre Modelos Lineares Generalizados Mistos
Artigo sobre Comparação Múltipla Não-paramétrica
Avaliação de Comparações Múltiplas
Análise de Sobrevivência Aplicada à Germinação
Sobre a Arquitetura de Raízes de Árvores
Características Morfológicas de Raízes de Pinus Taeda