Acessibilidade / Reportar erro

ÍNDICE DE VALOR DE IMPORTÂNCIA DE ESPÉCIES ARBÓREAS DA FLORESTA NACIONAL DO TAPAJÓS VIA ANÁLISES DE COMPONENTES PRINCIPAIS E DE FATORES

IMPORTANCE VALUE INDEX OF TREE SPECIES OF TAPAJÓS NATIONAL FOREST BY ANALYSIS OF PRINCIPAL COMPONENTS AND FACTORS

RESUMO

Este trabalho avaliou, a partir das análises de componentes principais e de fatores, a aplicabilidade do modelo linear IVI = a1Ar + a2Dr + a3Fr para obter índices de valor de importância (IVI) para espécies arbóreas, em que Ar é a abundância relativa, Dr é a dominância relativa e Fr é a frequência relativa, comparando-os com o método univariado que considera a equação tradicional que considera a1 = a2 = a3 =1. Os dados foram analisados a partir de uma amostra de 104,04 ha, localizada na Floresta Nacional do Tapajós, município de Belterra, estado do Pará, abrangendo um conjunto de 36 conglomerados de 4 ha, formato cruz de malta, com quatro subparcelas de 1 ha. Para obter os IVIs por espécie dentro de cada conglomerado, as subparcelas foram subdivididas cada uma em 5 subsubparcelas de 0,2 ha. Os resultados mostraram que os valores obtidos pelos métodos multivariados são diferentes dos computados pela equação tradicional. Verificou-se, também, que os valores dos índices uni e multivariados variaram de forma significativa dentro da floresta, não sendo recomendável, como normalmente ocorre, inferir um valor único para a área total da floresta, pois espécies consideradas de menor importância para o total da floresta podem possuir grande importância em determinadas subáreas. O uso do modelo linear com a1 = a2 = a3 =1 não é adequado, pois, por espécie e de acordo com a sua distribuição na floresta, esses coeficientes apresentam grandes variações.

Palavras-chave:
fitossociologia; inventário florestal; amostra por conglomerado

ABSTRACT

Principal components and factor analysis were used to evaluate the applicability of the linear model IVI = a1Ar + a2Dr + a3Fr for obtaining indexes of value of importance (IVI) for tree species, where Ar is relative abundance, Dr is the relative dominance and Fr is the relative frequency, to compare with the standard univariate method that considers a1 = a2 = a3 =1. The data comprised a sample of 104.04 ha, located in Tapajós National Forest, Belterra municipality, state of Pará, covering of a set of 36 clusters of 4 ha each, Malta croos-shaped, with four subplots of 1ha each. To obtain the IVI's by species within each cluster, each subplot was subdivided each into 5 sub-subplots of 0.2 ha each. The analysis revealed that the values obtained by multivariable differ from those obtained by the standard equation. It was also found that the values of univariate and multivariate indices vary significantly within the forest, and therefore are not recommended, as they usually are, to infer a single value for the whole forest area, because species considered of minor importance for the overall forest may be of great importance in particular subareas. The use of linear model with a1 = a2 = a3 =1 is not suitable, because of the wide variation of the coefficients among species and their spatial distribution in the forest.

Keywords:
phytosociology; forest inventory; cluster sampling

INTRODUÇÃO

A fitossociologia possibilita obter análises da estrutura florestal a partir de diversos índices existentes na literatura. O índice de valor de importância (IVI) sugerido por Curtis e McIntosh (1950CURTIS, J. T.; MCINTOSH, R. P. The interrelations of certain analytic and synthetic phytosociological characters. Ecology, Washington, v. 31, p. 434-455, 1950.), definido como a soma aritmética dos valores relativos de abundância, dominância e frequência, tem sido muito utilizado para determinar a importância ecológica das espécies, através da hierarquização em termos do grau de ocupação de sua população dentro do espaço geométrico da floresta, o que é expresso pelo número (abundância), tamanho (dominância) e distribuição espacial (frequência) dos indivíduos da população. Quanto maiores os valores de abundância, frequência e dominância mais importância terá a espécie dentro do complexo florístico da área. Cain et al. (1956CAIN, S. A. et al. Application of some phytosociological techniques to Brazilian rain forest. American Journal of Botany, New York, v. 43, p. 911-941, 1956.) foram os pioneiros na utilização do IVI para caracterizar a estrutura de uma floresta tropical.

Lamprecht (1964) e Finol (1971FINOL, U. H. Nuevos parametros a considerarse en el analysis estructural de Las Selvas Virgenes Tropicales. Revista Forestal Venezolana, Mérida, v. 14, p. 29-42, 1971.) ressaltam que o conceito de estrutura florestal deve ser ampliado, devendo ser adicionado ao modelo matemático outros parâmetros como posição sociológica, regeneração natural, qualidade do fuste, etc.

Queiroz (1995QUEIROZ, W. T. Inventário Florestal: processo de pós-estratificação multidimensional. Boletim da Faculdade de Ciências Agrárias do Pará, Belém, v. 23, p. 7-21, 1995.) ressalta que uma espécie pode não ter importância ecológica na área florestal como um todo, mas ter uma relevância setorizada em determinados estratos. Este autor recomenda que o plano amostral contemple a possibilidade de obter o valor do IVI por unidade amostral, sugerindo para isso o uso de uma estrutura amostral por conglomerados com subdivisão das subparcelas e assim obter os índices de importância por subárea.

A construção matemática do IVI, baseada na soma aritmética simples dos valores de abundância, dominância e frequência, que representa uma média aritmética, foi contestada por Queiroz (1995QUEIROZ, W. T. Inventário Florestal: processo de pós-estratificação multidimensional. Boletim da Faculdade de Ciências Agrárias do Pará, Belém, v. 23, p. 7-21, 1995.), visto que não considera as correlações entre as variáveis envolvidas. Este autor recomenda o uso de análises de componentes principais e de fatores para definir a combinação linear que servirá como índice, resultado que corresponde a obter uma média ponderada, tal que cada variável terá seu peso de importância na construção do índice. No caso do primeiro componente principal, ou do primeiro fator, explicar significativamente a estrutura de dispersão dos dados, então os mesmos podem ser considerados como funções lineares para expressar o índice. Por outro lado, no caso de haver necessidade de usar dois ou mais componentes ou fatores, então o índice pode ser definido, se usada à análise de fatores, como o módulo da resultante da soma vetorial dos escores fatoriais ponderados pela contribuição em porcentagem dos mesmos à comunalidade total. No caso de componentes principais, ter-se-á o módulo da resultante da soma vetorial dos escores dos componentes principais ponderados pela contribuição em porcentagem dos mesmos à variância total.

Este trabalho objetivou estudar e comparar o método tradicional univariado definido como a soma aritmética dos valores relativos de abundância (Ar ), frequência (Fr) e dominância (Dr ) para obtenção do índice de valor de importância de espécies arbóreas na análise da vegetação da Floresta Nacional do Tapajós versus índices calculados via análises de componentes principais e de fatores.

MATERIAL E MÉTODOS

Caracterização e localização da área de estudo

A área de estudo localiza-se na Floresta Nacional do Tapajós (2°40' - 4°10' S; 54°45' - 55°00' W), município de Santarém no estado do Pará, Brasil. A vegetação da área é constituída, basicamente, de uma floresta alta e relativamente aberta. A variação interna do conjunto florestal abrange áreas que fisionomicamente se aproximam estruturalmente às matas cipoálicas que provavelmente são representativas de estágios dinâmicos (QUEIROZ, 1984QUEIROZ, W. T. Análise de fatores ("Factor Analysis") pelo método da máxima verossimilhança: aplicação ao estudo da estrutura de florestas tropicais. 1984. 112 f. Tese (Doutorado em Estatística e Experimentação Agronômica) - Escola Superior de Agricultura "Luiz de Queiroz", Piracicaba, 1984.). De acordo com PROJETO RADAMBRASIL (1976) na sub-região dos Baixos Platôs da FLONA do Tapajós são encontradas em grande escala as espécies florestais Bertholletia excelsa Bonpl. - Lecythidaceae (castanheira-do-brasil), Bowdichia nitida Spruce ex Benth. - Fabaceae (sucupira), Ruizterania albiflora (Warm.) Marc.-Berti -Vochysiaceae (mandioqueira), Goupia glabra Aubl. - Goupiaceae (cupiúba) e Manilkara elata (Allemão ex Miq.) Monach - Sapotaceae (maçaranduba).

A área amostrada possui um tamanho de 3.745,44 ha, iniciando no km 62 da Rodovia BR-163, que liga Cuiabá a Santarém, estendendo-se no sentido norte-sul até o km 72, distanciando-se da margem direita da estrada um mínimo de 550 m, atingindo uma profundidade de 12,5 km em direção Leste - Oeste.

Planejamento da amostragem e coleta dos dados

A área de estudo foi dividida em 36 unidades primárias de 104,04 ha, correspondendo a um quadrado de 1.020 m de lado. Em cada unidade primária foi alocado um conglomerado, no qual a estrutura amostral possuía uma forma de cruz apresentando quatro subparcelas retangulares, orientadas segundo as direções Norte-Sul e Leste-Oeste (Figura 1), nas quais as subparcelas (unidades secundárias) foram tomadas com o tamanho de 500 m de comprimento por 20 m de largura, totalizando uma área de 1 ha. Cada subparcela foi subdividida em cinco subsubparcelas (unidades terciárias) de 20 m de largura e 100 m de comprimento. A área amostral total foi de 104,04 ha, que correspondeu à medição de 36 conglomerados de 4 ha cada ou 720 unidades terciárias de 20 m × 100 m.

FIGURA 1
Estrutura conglomerada cruz de malta utilizada no inventário florestal da FLONA Tapajós, Município de Belterra no estado do Pará, em 1977
FIGURE 1:
Cross of Malta shaped cluster structure used in the Tapajós National Forest inventory, Belterra municipality, state of Pará, in 1977.

Esse esquema amostral possibilitou obter o índice de valor de importância setorizado por conglomerado e para o total da área de estudo. No primeiro caso, foram consideradas como área mínima para a avaliação dos parâmetros estruturais as unidades terciárias de 0,2 ha (20 m × 100 m), o que totalizou 20 unidades por conglomerado. No segundo caso, considerou-se área mínima o próprio conglomerado com área de 4 ha, perfazendo 36 unidades, possibilitando obter os valores dos índices para a área total.

Os dados foram coletados por Queiroz (1977QUEIROZ, W. T. Efeitos da variação estrutural em unidades amostrais na aplicação de processos de amostragem em conglomerados nas Florestas do Planalto do Tapajós. 1977. 109 f. Dissertação (Mestrado em Engenharia Florestal) - Universidade Federal do Paraná, Curitiba, 1977.) e consistiram na medição do volume comercial de todas as árvores que apresentaram DAP acima de 30 cm. As árvores foram identificadas por parabotânico e em seguida armazenadas no herbário da EMBRAPA, Belém - PA.

Índice de valor de importância (IVI) segundo Curtis e McIntosh (1950CURTIS, J. T.; MCINTOSH, R. P. The interrelations of certain analytic and synthetic phytosociological characters. Ecology, Washington, v. 31, p. 434-455, 1950.)

Curtis e McIntosh (1950CURTIS, J. T.; MCINTOSH, R. P. The interrelations of certain analytic and synthetic phytosociological characters. Ecology, Washington, v. 31, p. 434-455, 1950.) propuseram o índice de valor de importância como a soma aritmética dos valores relativos de abundância (Ar ), frequência (Fr) e dominância (Dr ) para i-ésima espécie, ou seja, IVIi= Ari + Dri + Fri . Esse índice, considerado absoluto, leva em consideração apenas os parâmetros da estrutura horizontal das espécies, os quais foram calculados conforme Salomão et al. (2012SALOMÃO, R. P. et al. Análise fitossociológica de floresta ombrófila densa e determinação de espécies-chave para recuperação de área degradada através da adequação do índice de valor de importância. Boletim do Museu Paraense Emílio Goeldi. Ciências Naturais, Belém, v. 7, n. 1, p. 57-102, 2012.).

A abundância é definida como o número de indivíduos de uma espécie por unidade de área, sendo expresso em termos de relativa (Ar), conforme a expressão: Ari = (ni/N) x 100, em que ni = número de indivíduos da i-ésima espécie na amostragem; N é a abundância total ou número total de indivíduos por área amostral (soma das abundâncias absolutas de todas as espécies amostradas).

A frequência é um indicador do padrão de distribuição espacial dos indivíduos de uma população na floresta, sendo definida como a razão entre o número de parcelas em que ocorre pelo menos um indivíduo de uma espécie e o número total de parcelas da amostragem. Fri é a frequência relativa da espécie i definida como a porcentagem de ocorrência da espécie em relação ao somatório das porcentagens de ocorrências considerando o número total de espécies estudadas.

A dominância é definida como o grau de ocupação que uma espécie exerce sobre o espaço volumétrico do ecossistema. Em virtude da dificuldade de avaliação desse volume, a dominância é expressa em termos de área basal. A dominância relativa Dri é obtida pelo quociente entre a área basal da i-ésima espécie, em m2/ha, em relação à área amostrada em ha.

Para melhor explicar o comportamento da espécie serão considerados dois procedimentos, o primeiro corresponde ao IVI obtido dentro de cada um dos 36 conglomerados, enquanto que no segundo, o índice é computado acumulando os conglomerados.

Índice de valor de importância pelas análises de componente principal e de fatores

A análise de componentes principais consiste em, dado um conjunto de variáveis aleatórias, obter, por transformação linear, um número igual de novas variáveis não correlacionadas (componentes principais), tal que a soma das variâncias seja invariante. A análise de fatores estuda as relações internas de um conjunto de variáveis, substituindo-as por um número menor de variáveis hipotéticas denominadas fatores que podem ser correlacionados (fatores oblíquos) ou não correlacionados (fatores ortogonais). Esses fatores devem explicar parte significativa da variância do conjunto das variáveis originais. Valente et al. (2011VALENTE, M. D. R. et al. Modelo de predição para o volume total de Quaruba (Volchysia inundata ducke) via análise de fatores e regressão. Revista Árvore, Viçosa, v. 35, n. 2, p. 307-317, 2011.) utilizaram a análise de fatores para predizer o volume de madeira para ajuste de regressão.

No caso dos primeiros componentes principais ou dos fatores explicarem parte significativa da variância total, pode-se, então, avaliá-los para cada unidade de amostra e considerá-los como os valores do índice de valor de importância.

O IVI obtido pelas análises de componente principal e de fatores foi calculado somente em relação aos valores por conglomerado. Nesse caso, os dados foram arranjados em uma matriz de dados Xn x p , por espécie, em que n = 36 conglomerados e p são as três variáveis respostas Ar, Fr e Dr. Os componentes principais das observações da matriz Xn x p são funções lineares IVIi= a1i Ar+a2iDr+a3iFr com variância igual a , tal que , sendo que S = [sij] é a matriz amostral de covariâncias.

Para obter o primeiro componente principal calcula-se av1 , tal que , dada a condição , atinja o valor máximo.

Se o primeiro componente principal explicar parte significativa da variância total, ou seja, mais de 80%, de acordo com Hair Junior et al. (2005), então o mesmo deverá ser considerado como um bom índice de valor de importância.

O segundo componente principal (CP2) é obtido calculando-se av2, tal que atinja o valor máximo atendendo às condições , enquanto o terceiro componente principal (CP3) consiste na obtenção de av3 tal que atinja o valor máximo atendendo às condições

Por outro lado, no caso de haver necessidade de usar dois ou mais componentes, então o índice pode ser definido como o módulo da resultante da soma vetorial dos escores dos componentes principais ponderados pela contribuição em porcentagem (Pi) dos mesmos à variância total (CPsv).

Seja o modelo m-dimensional linear escrito na forma matricial: , sendo o vetor correspondente a p = 3 variáveis respostas.

A é a matriz de ordem 3 × m dos pesos correspondentes às três variáveis respostas, sendo m o número de fatores fatores fatores fatores do modelo, tal que m ≤ p

é o vetor correspondente aos fatores comuns (variáveis latentes) aleatórios, independentes e de variâncias iguais a um; é o vetor inerente aos resíduos específicos aleatórios de média zero.

Seja a matriz populacional de covariâncias de Ar, Dr e Fr. Considerando pressuposições: Resulta: que é a matriz de covariâncias das varáveis Ar, Dr e Fr sob o modelo fatorial.

Os elementos da matriz são denominados de comunalidades das variáveis respostas Ar, Dr e Fr, ou seja, é a porção da variância que a variável contribui para o fator.

No uso da análise de fatores, após verificar a normalidade dos dados, é importante identificar na matriz de correlações, o número de correlações que são estatisticamente diferentes de zero. Hair Junior et al. (2005) citam que para aplicar a análise de fatores o número de correlações com valores menores que 0,30 não deve ser grande.

Mingoti (2005MINGOTI, S. A. Análise de dados através de Métodos de Estatística Multivariada: uma abordagem aplicada. Belo Horizonte: UFMG, 2005. 300 p.) cita que alguns autores recomendam que para obter um ajustamento adequado do modelo de análise de fatores é importante que a matriz de correlações inversa seja próxima da matriz diagonal. Com base nessa condição, pode-se usar o coeficiente de KMO (Kaiser-Meyer-Olkin) que mede o grau de correlação parcial entre as variáveis. De acordo com Maroco (2007MAROCO, J. Análise estatística com utilização do SPSS. 3. ed. Lisboa: Edições Sílabo, 2007. 824 p.), a estatística de KMO (1.1) compara as correlações simples com as correlações parciais observadas entre as variáveis, sendo dada pela fórmula abaixo, para i ≠ j, tal que rij é o coeficiente de correlação entre as variáveis i e j e aij é o coeficiente de correlação parcial, em que os aij deverão assumir valores próximos de zero, uma vez que se pressupõe que os fatores são ortogonais entre si.

1.1

Os valores desse coeficiente variam de 0 a 1 e, segundo Fávero et al. (2009FÁVERO, L. P. et al. Análise de dados: modelagem para análise de dados. Rio de janeiro: Elsevier, 2009. 646 p.), quanto mais próximo da unidade mais recomendável é a aplicação da análise de fatores, pois quando as correlações parciais são próximas de zero, o coeficiente KMO tende a 1, significando que a matriz de correlações inversa tende à diagonal. Conforme estes autores os valores abaixo de 0,50 indicam não ser recomendável à aplicação da análise de fatores.

Além do uso do coeficiente de KMO, é recomendável também aplicar o teste de esfericidade de Bartlett para testar a hipótese nula de que as variáveis não sejam correlacionadas, ou seja, que a matriz de correlação entre as variáveis Ar, Fr e Dr é uma identidade, pois se essa hipótese não for rejeitada, então não tem sentido usar a técnica de análise de fatores, pois as variáveis originais são consideradas independentes. A expressão (1.2) indicada por Bartlett (1954) com 1/2×p(p+1) graus de liberdade foi usada para testar a hipótese .

1.2

Segundo Fávero et al. (2009FÁVERO, L. P. et al. Análise de dados: modelagem para análise de dados. Rio de janeiro: Elsevier, 2009. 646 p.), calcula-se para cada variável, similarmente à estatística KMO, uma medida de adequação da amostra (Measure of Sampling Adequacy - MSA). Para decidir sobre a necessidade de eliminação de variáveis no modelo de análise de fatores, pode-se obter a matriz Anti-Imagem O output do software SPSS apresenta os valores de MSA na diagonal principal dessa matriz. Hair Junior et al. (2005) definem a matriz de correlações anti-imagem como a matriz das correlações parciais entre as variáveis após a análise de fatores. A diagonal contém as medidas de adequação da amostra para cada variável e os demais valores são correlações parciais entre as variáveis.

A consistência interna ou confiabilidade entre as variáveis considerando as escalas dos dados pode ser avaliada pelo coeficiente de confiabilidade Alfa de Cronbach, citado em Cronbach (1951CRONBACH, L. J. Coefficient alpha and the internal structure of tests. Psychometrika, London, v. 16, n. 3, p. 297-334, 1951. ), dado pela expressão 1.3, tal que, p é o número de variáveis, Cov é a média das covariâncias e V é a média das variâncias.

1.3

O valor α pode variar de 0 a 1, significando que o valor 1 demonstra presença de consistência interna de 100 % e zero significa total ausência. A literatura recomenda como limite de aceitação valores acima de 0,70 (CORRAR; PAULO; DIAS FILHO, 2007CORRAR. L. J.; PAULO, E.; DIAS FILHO, J. M. Análise multivariada para os cursos de administração, ciências contábeis e economia. São Paulo: Atlas, 2007. 344 p.).

Neste trabalho será aplicado o método de análise de fatores denominado de máxima verossimilhança. Lawley (1940LAWLEY, D. N. The estimation of factor looding by the method of maximum likelihood. Proceedings of the Royal Society of Edinburgh, Edinburgh, v. 60, p. 64-82, 1940.) aborda originalmente o problema da estimação na análise de fatores, tornando-a um método estatístico propriamente dito. Este autor emprega o método da máxima verossimilhança para estimar os parâmetros do modelo fatorial , mostrando que as estimativas dos pesos fatoriais e as estimativas das variâncias específicas , tal que , são soluções das seguintes equações matriciais implícitas:

1.4

1.5

A principal vantagem do Método de Máxima Verossimilhança é verificar se um determinado número de fatores m, considerando uma amostra de tamanho n, está adequado para reproduzir o modelo. Segundo Bartlett (1954BARTLETT, M. S. A note on multiplying factors for various chi-squared approximations. Jornal of the Royal Statistical Society, London, v. 16, p. 296-298, 1954.), H0 é verdadeira quando a variável Um, dada em 1.6, tem aproximadamente uma distribuição qui-quadrado com v = 1/2[(p-m)2-p-m] graus de liberdade.

1.6

Os escores fatoriais não são estimados no sentido estatístico usual, pois não são parâmetros do modelo de fatores, mas valores atribuídos às variáveis hipotéticas. Atribuindo o princípio de quadrados mínimos para predizer os escores dos fatores, Thonsom (1951THONSOM, G. H. The factorial analysis of human ability. London: University of London, 1951.) estabeleceu o método para obter m regressões lineares dos fatores sobre as variáveis respostas, resultando a fórmula para obter os escores fatoriais que é uma matriz de ordem N ×n.

No caso do primeiro fator (F1 ) explicar significativamente a estrutura de dispersão dos dados, então o mesmo pode ser considerado como função linear para expressar o índice. Por outro lado, no caso de haver necessidade de usar dois ou mais fatores, então o índice pode ser definido como o módulo da resultante da soma vetorial dos escores fatoriais ponderados pela contribuição em porcentagem dos mesmos à comunalidade total (Fsv).

RESULTADOS E DISCUSSÃO

Análise segundo Curtis e McIntosh (1950CURTIS, J. T.; MCINTOSH, R. P. The interrelations of certain analytic and synthetic phytosociological characters. Ecology, Washington, v. 31, p. 434-455, 1950.)

A espécie Lecythis lúrida (Miers) S. A. Mori foi uma das mais importantes na estrutura da floresta estudada, pois ocorreu em todos os 36 conglomerados. Observando-se a Figura 2, verifica-se que a espécie Lecythis lúrida (Miers) S.A. Mori apresentou valores de IVI de 34,49, 36,07 e 42,16, respectivamente, para os conglomerados 5º, 20º e 30º, enquanto para as áreas acumuladas de 20 ha (5 conglomerados), 80 ha (20 conglomerados) e 120 ha (30 conglomerados) os valores de IVI foram de 23,45, 21,64 e 21,6, respectivamente. Estes resultados denotam uma grande variação no valor IVI quando avaliado por conglomerado, o que é uma evidência da existência de grande variação na ocorrência da espécie na área, ou seja, essa espécie não possui uma distribuição uniforme na floresta. Se o objetivo da amostragem é identificar e caracterizar essas comunidades, então o IVI por conglomerados é o que deve ser usado. Todavia, se o objetivo da amostragem é obter uma estimativa do IVI das espécies em relação à área total, o IVI acumulado deve ser considerado, pois ele indicará a importância da espécie em toda a área.

FIGURA 2
Distribuição do IVI absoluto para a espécie Lecythis lurida considerando os conglomerados isoladamente e na forma acumulada na área inventariada FLONA Tapajós, Município de Belterra PA, em 1977.
FIGURE 2:
Distribution of absolute IVI for Lecythis lurida considering clusters individually and cumulatively in Tapajós National Forest inventory, Belterra municipality, in 1977.

A espécie Cordia bicolor A.DC. foi a espécie que apresentou a menor importância em termos dos valores de IVI na área total, inclusive com vários conglomerados apresentando IVI iguais a zero (Figura 3).

FIGURA 3
Distribuição do IVI absoluto para a espécie Cordia bicolor considerando os conglomerados isoladamente e na forma acumulada na área inventariada da FLONA Tapajós, Município de Belterra - PA, em 1977
FIGURE 3
Distribution of absolute IVI for Cordia bicolor considering clusters individually. And cumulatively in Tapajós National Forest inventory, Belterra municipality-PA state, in 1977.

A grande variação nos valores do IVI dentro dos conglomerados foi verificada para todas as espécies encontradas na área de estudo, principalmente para as espécies pouco abundantes como C. bicolor. Os resultados apresentados mostram que espécies consideradas pouco importantes na estrutura de uma área florestal podem apresentar grande relevância quando analisadas dentro de determinadas subáreas. Esse fato deveria ser levado em consideração no planejamento da produção no manejo florestal (PEREIRA JUNIOR, 2004PEREIRA JÚNIOR, R. A. A utilização da fitossociologia como uma ferramenta na seleção de espécies para exploração e definição de tratamentos silviculturais em floresta de terra-firme sob manejo florestal, na Amazônia brasileira. 2004. 115 f. Dissertação (Mestrado em Ciências Florestais) - Universidade Federal Rural da Amazônia, Belém, 2004.).

Índices segundo as análises de componentes principais e de fatores

Os resultados da análise para o primeiro componente principal, considerando as dez espécies que apresentaram os maiores valores, segundo a expressão IVIi= Ari + Dri + Fri , assim como os autovetores e suas porcentagens de explicação (Tabela 1) demonstram que o uso da função linear IVIi= Ari + Dri + Fri ou escrita na forma IVIi= 0,577350Ari + 0,577350Dri + 0,577350Fri , segundo a condição , ou seja, norma do vetor dos coeficientes igual a unidade, não é recomendável, visto que os coeficientes do componente principal mostram valores diferentes entre si, resultantes das inter-relações entre variáveis aleatórias Ar, Dr e Fr.

TABELA 1:
Resultados da análise do primeiro componente principal para as 10 espécies de maior IVI na FLONA Tapajós, Município de Belterra - PA, em 1977.
TABLE 1
Results of the analysis of the first main component for the 10 species of highest IVI in Tapajós National Forest, Belterra, PA state, in 1977.

As análises de componentes principais para as 157 espécies mostraram que apenas Swartzia arborescens (Aubl.) Pittier, Pachira aquatica Aubl., Tachigali alba Ducke, Glycydendron amazonicum Dicke, Bowdichia nitida Spruce ex Berth., Mezilaurus itauba (Meisn.) Taub. ex Mez, Cedrela odorata L. e Virola michelii Heckel tiveram o 1º componente principal explicando menos de 90 % da variância total, ou seja, apenas 8 em um universo de 157 espécies, o que demonstra que o primeiro componente principal pode ser considerado um excelente Índice de Valor de Importância.

Para ilustrar a metodologia aplicada neste trabalho, utilizou-se a espécie Lecythis lurida, devido a sua ocorrência em todos os 36 conglomerados.

Na matriz de correlações da espécie Lecythis lurida (Tabela 2), verifica-se que não ocorreram valores menores que 0,30. A aplicação do teste F de Snedecor mostrou que todas as correlações são diferentes de zero (p-valor < 0,00001), demonstrando serem significativas as inter-relações entre as variáveis Ar, Fr e Dr . O teste de esfericidade de Bartlett apresentou um valor calculado igual a 83,4322 (p-valor < 0,00001), o qual permite concluir pela rejeição da hipótese de que a matriz de correlações corresponde a uma identidade. Nesse caso, podem-se aplicar as técnicas de análise de componentes principais e de fatores, pois as variáveis originais são consideradas dependentes.

TABELA 2:
Matrizes de covariâncias e de correlações de Lecythis lurida.
TABLE 2
Covariance and correlation matrices of Lecythis lurida

Os valores do coeficiente de confiabilidade de Cronbach foram iguais a 0,841155 para os dados da matriz de covariâncias e de 0,926466 considerando a matriz de correlações, demonstrando que as variáveis apresentam consistência interna ou confiabilidade com relação às escalas utilizadas para medição dos dados.

O valor da estatística de Kaiser-Meyer-Olkin (KMO) dada pela fórmula 1.3 apresentou um valor igual a 0,71678183, o que indica a existência de correlações entre as variáveis, recomendando, então, a aplicação da análise de fatores, pois KMO apresentou valor superior a 0,60.

O teste que mede o grau de correlação parcial entre as variáveis (Measure of Sampling Adequacy-MSA) apresentou os seguintes resultados: 0,64908; 0,77881050 e 0,74669263, respectivamente para as variáveis Ar, Dr e Fr. Como os valores das MSA foram superiores a 0,6, então deve-se usar todas as variáveis (Ar, Dr e Fr) nas análises, ou seja, não é necessário eliminar nenhuma delas.

A Tabela 3 mostra os valores das cargas ou pesos fatoriais, comunalidades e os coeficientes dos escores fatoriais obtidos pela análise de fatores pelo método da máxima verossimilhança a partir da matriz de correlações. Verificou-se, a partir do critério baseado no percentual de variância explicada, a necessidade de apenas um fator, que apresentou uma comunalidade igual a 2,44068178, correspondendo a uma explicação de 81,36 % da variação total. Ressalta-se que, nesta pesquisa, não foi possível aplicar o teste da razão da máxima verossimilhança, segundo Bartlett (1954BARTLETT, M. S. A note on multiplying factors for various chi-squared approximations. Jornal of the Royal Statistical Society, London, v. 16, p. 296-298, 1954.), expressão 1.6, para conhecer qual o número de fatores m está adequado para reproduzir o modelo, haja vista que o valor do número graus de liberdade de , dado por v =1\2[(p-m)2-p-m] , ser menor ou igual a zero para m. Então, a expressão 1.8 pode ser considerada como uma boa equação para construir índices.

1.8

TABELA 3:
Valores das cargas fatoriais, comunalidades e dos coeficientes fatoriais.
TABLE 3
Values of factor loadings, communalities and factor coefficients

A espécie Lecythis lurida apresentou os seguintes resultados para os três componentes principais:

CP1 = 0,506444Ar + 0,829354Dr + 0,235979Fr (autovalor igual a 23,3962)

CP2 = 0,724313Ar - 0,557663Dr + 0,405441Fr (autovalor igual a 1,8772)

CP3 = -0,467851Ar + 0,034410Dr + 0,883137Fr (autovalor igual a 0,4214)

O primeiro componente principal (CP1 ) apresentou um autovalor igual a 23,3962 que correspondeu a uma contribuição de 91,05 % para a variância total. O valor do autovalor para o segundo componente principal (CP2 ) foi de 1,8772, com participação de 7,31 %, enquanto o autovalor do terceiro componente principal (CP3 ) foi de 0,4214 que consistiu em uma contribuição de apenas 1,64 %. Os valores de CP1 e CP2 explicaram conjuntamente 98,36 % da dispersão total. Os valores dos índices de valor de importância baseados em CP1 = 0,506444Ar + 0,829354Dr + 0,235979Fr, considerando os 36 conglomerados, estão apresentados na Tabela 4. Na coluna CCP1 está a classificação hierárquica em função da magnitude dos resultados. Portanto, o primeiro componente principal pode ser considerado como uma boa equação para expressar o índice de valor de importância, pois explicou 91,05 % da dispersão total.

TABELA 4:
Comparação da importância de Lecythis lurida em termos de IVI (CIVI), escore do CP1 (CCP1), soma vetorial CPsv dos escores dos CP1, CP2 e CP3 (CCpv), escore F1 (CF1).
TABLE 4
Comparison of the importance of Lecythis lurida according to IVI (CIVI), score CP1 (CCP1),vetorial sum CPsv of the scores of CP1,CP2 e CP3 (CCpv), score F1 (CF1).

Por outro lado, optou-se pela obtenção do índice dado pela expressão 1.7 e definido como o módulo da resultante da soma vetorial dos escores dos componentes principais, ponderados pela contribuição em porcentagem dos respectivos componentes principais à variância total. Os valores deste índice (CPsv) para os 36 conglomerados são mostrados na Tabela 4. Na coluna CCPsv está a classificação hierárquica em função da magnitude dos resultados.

1.7

Para mostrar as metodologias descritas e efetivar uma comparação entre os índices IVI, escore do CP1, a soma vetorial considerando os escores dos três componentes principais (CP1, CP2, CP3) e o escore F1 , em termos de classificação pela magnitude para os 36 conglomerados, considerou-se a espécie Lecythis lurida. As observações concernentes às variáveis Ar, Fr e Dr estão inseridas na Tabela 4.

As porcentagens de coincidências considerando o índice com relação aos índices CP1 (36,1%), CPsv (36,1) e F1 (25,0%) (Tabela 5) são inexpressivas para afirmar que podem produzir resultados similares. Os índices CP1 e CPsv podem ser considerados similares, pois apresentaram uma coincidência nas classificações de 88,9 %. Este resultado deve-se ao fato do primeiro componente ter contribuído com 91,05 % da variância total, demonstrando que o CP1 é uma boa escolha como índice de valor de importância. Comparando-se CP1 com F1 obteve-se apenas 19,4 % de coincidências na classificação, concluindo-se que, apesar de serem técnicas multivariadas, apresentaram resultados diferentes. Apenas seis conglomerados tiveram a mesma classificação, considerando-se todos os índices, ou seja, uma porcentagem de coincidências de 16,7 %.

TABELA 5:
Número e porcentagens de classificações coincidentes dos diversos índices.
TABLE 5
Number and percentage of coincident classifications of various rates

Os valores das cargas fatoriais expressam que as correlações simples entre as variáveis originais Ar, Dr e Fr e a latente F1 correspondem, respectivamente, a 98,74 %, 84,5 % e 86,7 %.

CONCLUSÕES

Conclui-se que o planejamento da estrutura amostral por conglomerado com subdivisões das subparcelas propiciou uma estratificação estrutural-florística do povoamento, principalmente caracterizando que espécies consideradas como de menor importância para a estrutura florestal na área como um todo, apresentam grande relevância em determinadas subáreas da população;

As equações lineares obtidas através da análise de componentes principais explicaram mais de 90 % da variância dos dados originais, podendo-se concluir que o primeiro componente principal é uma boa escolha como um índice de valor de importância para caracterizar a estrutura florestal, pois apresenta a vantagem de absorver a maior parte da variabilidade apresentada pelas variáveis estudadas;

A equação linear correspondente ao primeiro fator obtido pelo método da máxima verossimilhança para a espécie Lecythis lurida é uma boa escolha como um índice de valor de importância, visto que apresentou uma comunalidade igual a 2,44068178, correspondendo a uma explicação de 81,36 % da variação total

REFERÊNCIAS

  • BARTLETT, M. S. A note on multiplying factors for various chi-squared approximations. Jornal of the Royal Statistical Society, London, v. 16, p. 296-298, 1954.
  • CAIN, S. A. et al. Application of some phytosociological techniques to Brazilian rain forest. American Journal of Botany, New York, v. 43, p. 911-941, 1956.
  • CORRAR. L. J.; PAULO, E.; DIAS FILHO, J. M. Análise multivariada para os cursos de administração, ciências contábeis e economia. São Paulo: Atlas, 2007. 344 p.
  • CRONBACH, L. J. Coefficient alpha and the internal structure of tests. Psychometrika, London, v. 16, n. 3, p. 297-334, 1951.
  • CURTIS, J. T.; MCINTOSH, R. P. The interrelations of certain analytic and synthetic phytosociological characters. Ecology, Washington, v. 31, p. 434-455, 1950.
  • FÁVERO, L. P. et al. Análise de dados: modelagem para análise de dados. Rio de janeiro: Elsevier, 2009. 646 p.
  • FINOL, U. H. Nuevos parametros a considerarse en el analysis estructural de Las Selvas Virgenes Tropicales. Revista Forestal Venezolana, Mérida, v. 14, p. 29-42, 1971.
  • HAIR JUNIOR, J. F. et al. Ensaio sobre la Estructura Florística de la Parte Sub-Oriental del Bosque Universitário "El Caimital" Estado Barinas. Revista Forestal Venezolana , Mérida, v. 7, n. 10-11, p. 77-119, 1964.
  • LAWLEY, D. N. The estimation of factor looding by the method of maximum likelihood. Proceedings of the Royal Society of Edinburgh, Edinburgh, v. 60, p. 64-82, 1940.
  • MAROCO, J. Análise estatística com utilização do SPSS. 3. ed. Lisboa: Edições Sílabo, 2007. 824 p.
  • MINGOTI, S. A. Análise de dados através de Métodos de Estatística Multivariada: uma abordagem aplicada. Belo Horizonte: UFMG, 2005. 300 p.
  • PEREIRA JÚNIOR, R. A. A utilização da fitossociologia como uma ferramenta na seleção de espécies para exploração e definição de tratamentos silviculturais em floresta de terra-firme sob manejo florestal, na Amazônia brasileira. 2004. 115 f. Dissertação (Mestrado em Ciências Florestais) - Universidade Federal Rural da Amazônia, Belém, 2004.
  • QUEIROZ, W. T. Efeitos da variação estrutural em unidades amostrais na aplicação de processos de amostragem em conglomerados nas Florestas do Planalto do Tapajós. 1977. 109 f. Dissertação (Mestrado em Engenharia Florestal) - Universidade Federal do Paraná, Curitiba, 1977.
  • QUEIROZ, W. T. Análise de fatores ("Factor Analysis") pelo método da máxima verossimilhança: aplicação ao estudo da estrutura de florestas tropicais. 1984. 112 f. Tese (Doutorado em Estatística e Experimentação Agronômica) - Escola Superior de Agricultura "Luiz de Queiroz", Piracicaba, 1984.
  • QUEIROZ, W. T. Inventário Florestal: processo de pós-estratificação multidimensional. Boletim da Faculdade de Ciências Agrárias do Pará, Belém, v. 23, p. 7-21, 1995.
  • PROJETO RADAMBRASIL. Folha S/A 21-Santarém. Rio de Janeiro: Departamento Nacional de Produção Mineral, 1976. (Levantamento de Recursos Naturais, 10).
  • SALOMÃO, R. P. et al. Análise fitossociológica de floresta ombrófila densa e determinação de espécies-chave para recuperação de área degradada através da adequação do índice de valor de importância. Boletim do Museu Paraense Emílio Goeldi. Ciências Naturais, Belém, v. 7, n. 1, p. 57-102, 2012.
  • THONSOM, G. H. The factorial analysis of human ability. London: University of London, 1951.
  • VALENTE, M. D. R. et al. Modelo de predição para o volume total de Quaruba (Volchysia inundata ducke) via análise de fatores e regressão. Revista Árvore, Viçosa, v. 35, n. 2, p. 307-317, 2011.

Datas de Publicação

  • Publicação nesta coleção
    Mar 2017

Histórico

  • Recebido
    09 Nov 2012
  • Aceito
    16 Jun 2015
Universidade Federal de Santa Maria Av. Roraima, 1.000, 97105-900 Santa Maria RS Brasil, Tel. : (55 55)3220-8444 r.37, Fax: (55 55)3220-8444 r.22 - Santa Maria - RS - Brazil
E-mail: cienciaflorestal@ufsm.br