Rodolpho edited subsubsection_M_todos_textit_rea__.tex  over 8 years ago

Commit id: 9cae7aa4368f5fad49166e6f1ce270fe6d3b5971

deletions | additions      

       

\textit{Análises estatísticas}  O modelo de ocupação para comunidades (multi-species occupancy model) foi elaborado de forma a representar três processos distintos, mas relacionados: presença ou ausência das espécies na metacomunidade, presença ou ausência das espécies na localidade amostrada e detecção ou não das espécies em cada período amostral \cite{Dorazio_&_2005}. \cite{Dorazio_&_Royle_2005}.  Estes três processos são considerados de forma hierárquica no modelo, de modo que a detecção das espécies depende da ocupação das localidades por estas espécies, que por sua vez depende das espécies estarem presente no pool regional de espécies considerado como a metacomunidade \cite{Dorazio_et_al._2006}. Como as ocasiões temporais de amostragem podem ser modeladas por covariáveis de detecção e também a ocupação de cada espécie também pode ser relacionada a covariáveis de ocupação, os modelos podem se tornar bastante complexos e com muitos parâmetros a serem estimados. Nestes casos, tem sido indicado o uso da abordagem bayesiana, onde os parâmetros destes modelos são mais facilmente estimados usando o algoritmo computacional Markov Chain Monte Carlo method (ou método MCMC). Entretanto, como no presente momento nós estamos construindo os modelos e buscando entender melhor esta abordagem, nós optamos por estimar o número de espécies por sítio a partir de um modelo básico, que não inclui nenhuma covariável. Desta forma, para estimar e avaliar a relação entre as variáveis dependentes e independentes, nós utilizamos a abordagem de verossimilhança, a qual é mais conhecida, relativamente mais simples e facilmente implementável em softwares estatísticos como o R \cite{team2014r}. As variáveis dependentes (riqueza observada e estimada) foram relacionadas com as independentes (dois primeiros eixos da PCA da vegetação) por meio de um modelo linear generalizado (GLM) com distribuição dos erros Poisson, a qual é indicada para modelar variáveis dependentes que são contagens, como no nosso caso. Os modelos construídos para representar a relação entre as variáveis incluem relações lineares entre as variáveis dependentes e cada uma das variáveis independentes, além de um modelo onde a variável dependente possa ser relacionada com a soma dos efeitos das duas variáveis independentes. Além da relação linear entre riqueza de espécies e a estrutura da vegetação, é esperado também que a riqueza possa estar relacionada com a estrutura vegetacional por uma relação quadrática. Neste caso, a hipótese ecológica representada seria a de uma maior riqueza de espécies em estruturas de vegetação intermediárias e as áreas com estruturas menores e maiores com riquezas mais baixas. Assim, nós também construímos modelos em que cada uma das variáveis independentes se relacionam de forma quadrática com as variáveis dependentes, além de um modelo onde ambas as variáveis independentes possam se relacionar de forma quadrática com a variável dependente.  Para a seleção dos modelos que melhor representem a variação da riqueza observada e estimada em nossas amostras, nós utilizamos uma adaptação do critério de informação de Akaike (AIC), indicada para lidar com pequenas amostras (AICc). Este critério permite escolhermos o modelo que mais se aproxima do modelo “original” que possa ter gerado os dados, mesmo sendo este modelo desconhecido. Esta escolha é baseada na diferença dos AICc's entre os modelos concorrentes (deltaAICc), sendo o melhor modelo aquele que apresenta menor valor de AICc (deltaAICc=0). Outros modelos que apresentem deltaAICc igual ou menor que dois em comparação com o melhor, são considerados igualmente plausíveis e não devem ser descartados como possível explicação para os dados \cite{burnham2002model}.