Aplicação do resíduo específico na análise de grupos de experimentos

Application of specific error in the analysis of groups of experiments

Resumos

A análise de grupos de experimentos tem como objetivo principal produzir informações que possam ser aplicadas a uma população, seja no espaço e/ou no tempo. O objetivo deste trabalho é aplicar a técnica do resíduo específico para contornar o problema da heterocedasticidade presente nas interações dos experimentos com tratamentos, e obter as esperanças dos quadrados médios, para os efeitos considerados na análise de variância. Os dados utilizados referem-se a grupo de 11 experimentos, de adubação, no esquema fatorial 3³ com cana de açúcar, desenvolvidos em diferentes localidades do Estado de São Paulo. Os resultados revelaram que a técnica do resíduo específico é adequada quando a heterocedasticidade está presente nas interações.

interação; heterogeneidade de variâncias; contraste ortogonal; grupos de experimentos

The main aim of the analysis of groups of experiments is to produce information that can be applied to a population either in space or in time and in both. The aim of this paper is to apply the specific error technique in order to solve the heterogeneity problem of the interactions of experiments with treatments, and to obtain the expectations of the mean squares for the considered effects in the analysis of variance. The data that were used refer to a group of 11 experiments, in a 3³ factorial design with sugar cane, conducted at different places of Sao Paulo State, Brazil. The results showed that the specific error technique is adequate when the heterogeneity is present in the interactions.

interaction; heterogeneity of variances; orthogonal contrast; groups of experiments

METODOLOGIA E TÉCNICAS EXPERIMENTAIS

Aplicação do resíduo específico na análise de grupos de experimentos

Application of specific error in the analysis of groups of experiments

Wilson Alves de OliveiraI, ^*; Maria Cristina Stolf Nogueira^II

^ICentro de Ciências Exatas e Tecnológicas da Universidade Estadual do Oeste do Paraná – UNIOESTE, Caixa Postal 801, 85814-110 Cascavel(PR). E-mail: [e-mail]

^IIDepartamento de Ciências Exatas – ESALQ/USP, Caixa Postal 09, 13418-900 Piracicaba (SP). E-mail: [e-mail]

RESUMO

A análise de grupos de experimentos tem como objetivo principal produzir informações que possam ser aplicadas a uma população, seja no espaço e/ou no tempo. O objetivo deste trabalho é aplicar a técnica do resíduo específico para contornar o problema da heterocedasticidade presente nas interações dos experimentos com tratamentos, e obter as esperanças dos quadrados médios, para os efeitos considerados na análise de variância. Os dados utilizados referem-se a grupo de 11 experimentos, de adubação, no esquema fatorial 3³ com cana de açúcar, desenvolvidos em diferentes localidades do Estado de São Paulo. Os resultados revelaram que a técnica do resíduo específico é adequada quando a heterocedasticidade está presente nas interações.

Palavras-chave: interação, heterogeneidade de variâncias, contraste ortogonal, grupos de experimentos.

ABSTRACT

The main aim of the analysis of groups of experiments is to produce information that can be applied to a population either in space or in time and in both. The aim of this paper is to apply the specific error technique in order to solve the heterogeneity problem of the interactions of experiments with treatments, and to obtain the expectations of the mean squares for the considered effects in the analysis of variance. The data that were used refer to a group of 11 experiments, in a 3³ factorial design with sugar cane, conducted at different places of Sao Paulo State, Brazil. The results showed that the specific error technique is adequate when the heterogeneity is present in the interactions.

Key words: interaction, heterogeneity of variances, orthogonal contrast, groups of experiments.

1. INTRODUÇÃO

Experimentos agrícolas envolvendo o mesmo fator ou fatores podem ser desenvolvidos em vários locais distintos e/ou em várias ocasiões diferentes. Algumas vezes, o objetivo da pesquisa é produzir informações para uma população que se estende, no espaço e/ou no tempo.

A análise estatística apropriada para dados de grupos de experimentos vai variar com o objetivo da pesquisa. Todavia, os passos preliminares da análise tendem a ser os mesmos em todos os casos e uma aplicação, não crítica, dos métodos apropriados para análise conjunta, pode levar a conclusões errôneas.

YATES e COCHRAN (1938) sugerem um exame cuidadoso dos resultados experimentais antes de iniciar qualquer análise combinada, e atestam que o procedimento de análise de variância usual, apropriado para analisar os resultados de um só experimento, pode requerer modificação, devido à heterocedasticidade presente entre os diferentes experimentos e, podendo também, ocorrer a heterocedasticidade entre os componentes da interação de tratamentos com locais. Assim, sugerem a partição das somas de quadrados devido a tratamentos e da interação de tratamentos com locais em componentes ortogonais apropriados.

KEMPTHORNE (1952) considera que na análise conjunta, existem duas possíveis dificuldades para a interpretação dos resultados na análise de variância. A primeira dificuldade é que a variância residual(s²)não é constante entre os experimentos, e a segunda, é que o componente de variância referente à interação de tratamentos com locais() depende da combinação de tratamentos e locais.

COCHRAN e COX (1957) fazem duas críticas sobre a análise preliminar. A primeira crítica, é que alguns componentes da soma de quadrados da interação de tratamentos com experimentos, podem ser maiores que outros, ou seja, a variância da interação não é homogênea. Nesse caso, o teste F para tratamentos tendo como denominador a interação de tratamentos com experimentos é viciado, no sentido de que o valor de F, obtido nas tabelas é muito baixo, isto é, obtêm-se muitos resultados significativos. Segundo os autores, um método para contornar essa dificuldade é decompor a soma de quadrados de tratamentos em um conjunto de componentes ortogonais, que fornecerá toda ou a maioria das informações de interesse. A soma de quadrados das interações é decomposta do mesmo modo, obtendo-se um resíduo específico, para testar cada um dos componentes de tratamentos. Uma segunda crítica é que as variâncias dos erros experimentais não são as mesmas em todos os experimentos, podendo essa variação invalidar o teste F para interações tendo como denominador o Resíduo Médio.

Em GOMES (1990), constata-se aplicação do método da decomposição da interação de tratamentos com locais, em análise de grupos de experimentos de adubação de algodão. O autor, ao fazer a análise conjunta dos diversos experimentos, considera individualmente os graus de liberdade relativos aos contrastes de tratamentos e estima os contrastes em cada local, obtendo assim, o quadrado médio apropriado para testar cada um dos contrastes. Alerta, porém, para o fato que esse método tem a desvantagem de reduzir excessivamente o número de graus de liberdade para os testes, recomendando seu uso somente quando for grande o número de locais.

NOGUEIRA (2004) realizou revisão sobre definições e conceitos envolvendo contrastes, buscando justificar que a técnica de contrastes ortogonais, com um número de graus de liberdade, é simples e bastante eficiente na análise de dados experimentais, tais como: obtenção dos efeitos principais, de efeito de interação, de efeitos aninhados, nas comparações entre grupos de médias e na obtenção do resíduo específico.

Neste trabalho, é aplicada a técnica do resíduo específico para contornar o problema da heterocedasticidade das interações e, são demonstradas as esperanças dos quadrados médios para os efeitos considerados na análise de variância.

2. MATERIAL E MÉTODOS

2.1 Material

Os dados analisados referem-se a rendimentos de cana-de-açúcar, em t ha^-1, obtidos de grupo de 11 experimentos, de adubação N x P x K, instalados em terra roxa legítima, em diferentes localidades do Estado de São Paulo, retirado de Malavolta et al (1963).

O nitrogênio foi utilizado nas doses de 0, 60 e 120 kg ha^-1, sendo 1/3 no plantio, 1/3 em cobertura, aos 3 meses, e 1/3 aos 6 meses, quando ainda era usado como adubo nitrogenado o salitre-do-chile. O fósforo, usado nas doses de 0, 75 e 150 kg ha^-1 de P₂O₅, foi fornecido todo no plantio, metade como superfosfato simples e metade como hiperfosfato (atualmente não utilizado). O potássio, aplicado sob a forma de cloreto, todo no plantio, teve as doses de 0, 75 e 150 kg.ha^-1 de K₂O.

Os experimentos foram instalados no esquema fatorial 3³, com confundimento de 2 graus de liberdade da interação tripla N x P x K, de acordo com o grupo W (YATES, 1937), gerando três blocos envolvendo nove tratamentos para cada bloco, com uma repetição. Neste caso, segundo BANZATTO e KRONKA (1989), entre outros, a parte não confundida da interação tripla, com seis números de graus de liberdade, funciona como resíduo. Visando obter resultados mais amplos e detalhados, desdobraram-se as interações duplas, isolando-se apenas a interação entre os componentes lineares e juntando-se os outros três componentes (ou números de graus de liberdade) ao resíduo, que passou a ser 15 graus de liberdade.

2.2 Métodos

2.2.1 Modelo Matemático

O modelo matemático adotado para a análise conjunta foi o seguinte:

no qual:

y_ijklm é a produção de cana-planta, em t.ha-1, observada no experimento i, bloco j, para a klm-ésima combinação dos fatores N, P e K respectivamente;

i = 1, ... , I; j = 1, ... , J; k = 1, ... , K; l = 1, ... , L e m = 1, ... , M;

µ é a média geral;

E_i é o efeito do experimento i, considerado aleatório;

B_j (E_i) é o efeito do bloco j dentro do experimento i, considerado aleatório;

N_k é o efeito do fator nitrogênio, considerado fixo;

P_l é o efeito do fator fósforo, considerado fixo;

K_m é o efeito do fator potássio, considerado fixo;

(NP)k_l é o efeito da interação nitrogênio com fósforo, considerado fixo;

(NP)k_m é o efeito da interação nitrogênio com potássio, considerado fixo;

(PK)_lm é o efeito da interação fósforo com potássio, considerado fixo;

(EN)_ik é o efeito da interação experimento com nitrogênio, considerado aleatório;

(EP)_il é o efeito da interação experimento com fósforo, considerado aleatório;

(EK)_im é o efeito da interação experimento com potássio, considerado aleatório;

(ENP)_ikl é o efeito da interação experimento com N e com P, considerado aleatório;

(ENK)_ikm é o efeito da interação experimento com N e com K, considerado aleatório;

(EPK)_ilm é o efeito da interação experimento com P e com K, considerado aleatório;

e_ijklm é o erro aleatório, normal e independentemente distribuído, com média zero variância s2.

2.2.2. Resíduo Específico

O método proposto para contornar a heterocedasticidade das variâncias das interações é decompor as somas de quadrados de tratamentos em um conjunto de componentes ortogonais. As somas de quadrados das interações, também, são decompostas do mesmo modo, para se obter um resíduo específico, o qual servirá para testar cada um dos componentes de tratamentos com o componente obtido das interações dos tratamentos com experimentos.

Como os níveis dos fatores de tratamentos, para o conjunto de dados considerado são equidistantes, essa decomposição pode ser feita de modo simples, utilizando os coeficientes dos polinômios ortogonais.

O cálculo das somas de quadrados para a regressão linear e quadrática, para os efeitos principais, é obtido da seguinte forma:

e ;

em que e

c_1n e c_2n são os coeficientes dos polinômios ortogonais de 1º e 2º graus respectivamente;

T_n é o total das observações do n-ésimo nível do efeito considerado;

r é o número de parcelas somadas para obter cada total (T_n);

a é o número de níveis do efeito considerado;

I é o número de experimentos.

As somas de quadrados para as regressões linear e quadrática, para as interações são obtidas por:

SQRLinear ´ Exper.

e SQRQuadratica ´ Exper.

Um programa SAS foi elaborado para realizar as análises estatísticas propostas, utilizando o comando CONTRAST/ PROC GLM.

3. RESULTADOS E DISCUSSÃO

3.1 Esperança dos quadrados médios

Para se testar a hipótese de interesse, através do teste F, é necessário estudar a natureza dos quadrados médios na análise de variância. Usando o modelo matemático adotado no item 2.2.1, foram consideradas as seguintes restrições:

(i) Para os efeitos fixos:

(ii) Para os efeitos aleatórios:

E_i, B_j(E_i), (EN)_ik, (EP)_il,, (EK)_im, (ENP)_ikl, (ENK)_ikm, (EPK)_ilm e e_ijklm são normais e independentemente distribuídos, com médias zero e variâncias , , , , , , , , , respectivamente.

Além dessas restrições, considera-se que os efeitos aleatórios são não correlacionados.

Na tabela 1, estão apresentadas as esperanças matemáticas para os quadrados médios dos efeitos considerados no modelo, obtidas através do método dos momentos (Hicks, 1973).

A composição para aplicação do teste F, considerando os resultados obtidos para as E(QM) foi realizada da seguinte forma:

(a) Para o efeito de Experimentos (E):

com (I-1) e v₁ graus de liberdade, sendo v₁, obtido pela fórmula de SATTERTHWAITE (1946).

(b) Para o efeito de doses de N:

com (K-1) e (I-1)(K-1) graus de liberdade.

(c) Para o efeito de doses de P:

com (L-1) e (I-1)(L-1) graus de liberdade.

(d) Para o efeito de doses de K:

com (M-1) e (I-1)(M-1) graus de liberdade.

(e) Para o efeito da interação NxP:

com (K-1)(L-1) e (I-1)(K-1)(L-1) graus de liberdade.

(f) Para o efeito da interação NxK:

com (K-1)(M-1) e (I-1)(K-1)(M-1) graus de liberdade.

(g) Para o efeito da interação PxK:

com (L-1)(M-1) e (I-1)(L-1)(M-1) graus de liberdade.

(h) Para o efeito da interação ExN:

, com (I-1)(K-1) e v₂ graus de liberdade, sendo v₂, obtido pela fórmula de Satterthwaite ( 1946).

(i) Para o efeito da interação ExP:

, com (I-1)(L-1) e v₃ graus de liberdade, sendo v₃, obtido pela fórmula de Satterthwaite ( 1946).

(j) Para o efeito da interação ExK:

, com (I-1)(M-1) e v₄ graus de liberdade, sendo v₄, obtido pela fórmula de Satterthwaite ( 1946).

(k) Para o efeito da interação ExNxP:

, com (I-1)(K-1)(L-1) e I(15) graus de liberdade.

(l) Para o efeito da interação ExNxK:

, com (I-1)(K-1)(M-1) e I(15) graus de liberdade.

(m) Para o efeito da interação ExPxK:

, com (I-1)(L-1)(M-1) e I(15) graus de liberdade.

3.1.1 Esperança dos quadrados médios para os componentes lineares e quadráticos de N, P e K

No caso dos experimentos de adubação, existe interesse em determinar a curva de produção da cultura em estudo, em função das doses crescentes aplicadas de N, P e K. Para isso, é usual decompor os efeitos de N, P e K em componentes lineares e quadráticos. A técnica utilizada é a dos polinômios ortogonais, e designa-se como c_1k e c_2k os coeficientes dos polinômios linear e quadrático para o nitrogênio respectivamente.

A soma de quadrados para o efeito linear do nitrogênio é dada pela seguinte expressão:

, tal que .

Assim,

em que,.

Portanto,

, (1) pois, a SO (N') está associada a 1 número de grau de liberdade e é o valor teórico do contraste que corresponde ao efeito linear de N.

Para o conjunto de dados em questão, o fator N tem três níveis, quantitativos, e os valores de c_1k são referentes aos coeficientes do polinômio ortogonal de grau 1, isto é, c₁₁ = -1, c₁₂ = 0 e c₁₃ = 1. Logo,

Analogamente, utilizando d_il como coeficientes do polinômio ortogonal para o efeito linear do fósforo e e_lm para o efeito linear do potássio, obtém-se:

e

,

Sendo:

e os valores teóricos dos contrastes que correspondem aos efeitos lineares de, respectivamente, P e K.

Para obter as esperanças dos quadrados médios para os efeitos quadráticos de N, P e K, basta substituir os coeficientes lineares c_1k, d_1l e e_1m pelos coeficientes quadráticos c_2k , d_2l e e_2m respectivamente.

3.1.2 Esperanças dos quadrados médios para N' x P', N' x K' e P' x K'.

Nos experimentos de adubação, com a finalidade de verificar se o efeito de um nutriente tem influência sobre o efeito do outro, existe interesse em estudar as interações de primeira ordem entre os nutrientes utilizados. As interações de interesse prático, segundo GOMES (1990), são somente aquelas entre os efeitos lineares, isto é, N' x P', N'x K' e P'x K'.

A soma de quadrados para N'x´P' é dada pela seguinte expressão:

,

tal que .

Seguindo procedimento análogo à N', P' e K', obtém-se que:

Sendo:

, e os valores teóricos dos contrastes que correspondem aos efeitos de N' x P' x , N' x K' e P' x K' respectivamente.

3.1.3 Esperança dos quadrados médios para as interações E x N', E x P' e E x K'

A soma de quadrados para E x N' é dada pela seguinte expressão:

, tal que,.

Aplicando-se esperança matemática, obtém-se que, e assim, (2)

Portanto, observando as equações (1) e (2), verifica-se que QM(EN') é o resíduo específico para testar o efeito linear do nitrogênio (N').

Analogamente tem-se,

e ,

que correspondem, respectivamente, aos resíduos específicos para testar os efeitos lineares do fósforo (P') e do potássio (K').

3.1.4 Esperança dos quadrados médios para E x N' x P', E x N' x K' e E x P' x K'

A soma de quadrados para EN'P' é dada pela seguinte expressão:

, tal que, .

Aplicando-se esperança matemática, obtém-se que,

.

e assim,

.

Analogamente tem-se,

e

.

que correspondem, respectivamente, aos resíduos específicos para testar as interações

3.2 Verificação da homocedasticidade entre experimentos

O teste utilizado para verificar a homogeneidade das variâncias residuais, entre experimentos, foi o teste de Bartlett, sendo essa hipótese rejeitada, a de 5% de probabilidade. Na tabela 2, observam-se os quadrados médios dos resíduos para os experimentos.

3.3 Análise de variância conjunta

O teste F, aplicado na análise de variância, tem os valores obtidos de acordo com os resultados das respectivas esperanças dos quadrados médios, apresentados em 3.1, cujos resultados são mostrados na tabela 3.

Pela tabela 3, verifica-se que todas as interações envolvendo o fator experimentos não foram significativas, isto é, os efeitos de N, P e K foram semelhantes nos diferentes experimentos (locais). As interações N'P', N'K' e P'K' não foram significativas, posto que o efeito de um fertilizante não interfere no efeito do outro, em média. Observou-se, também, que os efeitos médios para N e K foram significativos. Assim, realizou-se a decomposição dos dois graus de liberdade associados à soma de quadrados de cada um deles, ou seja, um referente ao efeito linear e outro, ao efeito quadrático do fertilizante. O teste F, para os efeitos linear e quadrático de N e K, tem como denominadores EN e EK respectivamente. Neste trabalho, o interesse está em decompor as somas de quadrados das interações EN e EK para a obtenção dos resíduos específicos, para se testar cada um desses efeitos, devido à heterocedasticidade presente nas variâncias residuais entre experimentos. Os resultados são apresentados na tabela 4.

Na tabela 5, constata-se a aplicação do teste F, sem utilizar o resíduo específico. Confrontando as tabelas 4 e 5, observa-se que o denominador na composição do teste F, para N' e N'' é o QM(EN), que representa o valor médio entre QM(EN') e QM(EN''). Se estes valores são relativamente diferentes, o denominador com o QM(EN) não é o mais adequado.

4. CONCLUSÕES

1. As interações E x N', E x P', E x K', E x N'', E x P'' e E x K'' são os resíduos específicos para se testar os efeitos de N', P', K', N'', P'' e K'' respectivamente. Da mesma forma que as interações E x N'P', E x N'K' e E x P'K' são os resíduos específicos para se testar os efeitos de N' x P', N' x K' e P' x K' respectivamente.

2. A técnica do resíduo específico, aplicada nos dados experimentais apresentados, é alternativa adequada, de análise de dados de grupo de experimentos, na presença de heterocedasticidade entre experimentos.

Recebido para publicação em 8 de julho de 2005 e aceito em 9 de maio de 2007.

• BANZATTO, D.A.; KRONKA, S.N. Experimentação Agrícola. Jaboticabal: FUNEP-UNESP, 1989. 247p.
• COCHRAN, W.G. The combination of estimates from different experiments. Biometrics, Baltimore, v.10, n.1, p.101-129, 1954.
• COCHRAN, W.G.; COX, G.M. Experimental Designs 2.ed. New York: John Wiley and Sons, 1957. 611p.
• GOMES, F.P. Curso de Estatística Experimental. 13.ed. Piracicaba: Nobel, 1990. 467p.
• HICKS, C.R. Fundamental Concepts in the Design of Experiments. 2.ed. New York: Holt, Rinehart and Winston, 1973. 349p.
• KEMPTHORNE, O. The design and analysis of experiments. New York: John Wiley and Sons, 1952. 631p.
• MALAVOLTA, E. et al. A diagnose foliar na cana-de-açucar IV. Resultados de 40 ensaios fatoriais NPK 333, primeiro corte no estado de São Paulo. ESALQ-USP, Piracicaba, 1963. 47p.
• NOGUEIRA, M.C.S. Orthogonal Contrasts: Definitions and Concepts. Scientia Agricola, Piracicaba, v. 61, n. 1, p. 118-124, 2004.
• SAS Institute Inc., SAS/STAT^TM User's Guide, Release 9.1.ed. Cary, NC: SAS Institute Inc., 2006.
• SATTERTHWAITE, F.E. An approximate distribution of estimates of variance components. Biometrics Bulletin, Washington, v.2, p.110-114, 1946.
• YATES, F. The design and analysis of factorial experiments. Imperial Bureau of Soil Science, Harpenden, v.35, 1937.
• YATES, F.; COCHRAN, W.G. The analysis of groups of experiments. Journal of Agricultural Science, Harpenden, v.28, p.556-580, 1938.

*

Autor para correspondência.

Datas de Publicação

• Publicação nesta coleção
  24 Jan 2008
• Data do Fascículo
  2007

Histórico

• Aceito
  09 Maio 2007
• Recebido
  08 Jul 2005