TAMANHO AMOSTRAL ENVOLVENDO MÚLTIPLOS CARACTERES: UMA APLICAÇÃO À DIGESTIBILIDADE EM BOVINOS

SAMPLE SIZE INVOLVING MULTIPLE TRAITS: AN APPLICATION TO DIGESTIBILITY IN CATTLE

Resumos

O objetivo deste trabalho foi estimar o tamanho amostral mínimo (n) para comparar tratamentos em experimentos de consumo e de digestibilidade com bovinos, envolvendo múltiplos caracteres. Foram utilizados dados de digestibilidade de 72 novilhas com média de 18 meses de idade e 250 kg de peso. O experimento foi realizado na Embrapa-Centro de Pesquisa de Pecuária do Sudeste, São Carlos, SP, de 1988 a 1989, em delineamento inteiramente casualizado, com nove tratamentos organizados em esquema fatorial 3 x 3 (três grupos genéticos: Canchim, ½ Canchim + ½ Nelore e Nelore, e três níveis de proteína bruta: 6, 10 e 13%, com oito repetições cada, sendo a unidade experimental a novilha). Foram analisados consumo de ração (g/kg0,75) por quilograma de peso metabólico, energia digestível, nitrogênio retido (NR), NR (mg/kg0,75) e digestibilidades da matéria seca, proteína bruta, fibra em detergente neutro e fibra em detergente ácido. O valor mínimo de n, que permite detectar diferenças significativas (delta) entre vetores de médias de tratamentos, foi obtido por meio de um programa SAS (Statistical Analysis System), considerando modelo de distribuição normal t-variada, média zero e matriz de covariância sigma, estatística T² de Hotelling, distribuição F com parâmetro de não centralidade (<FONT FACE=Symbol>d²D</FONT>), erros do tipo I (alfa), poder do teste (1 beta) e delta. O valor de n variou de 6 a 47, sendo mais influenciado por alteração nos valores de delta, do que nos valores de alfa e poder do teste.

distribuição F; diferença mínima significativa; Hotelling

The objective of this work was to estimate minimum sample size (n) for comparison of treatments in experiments of consumption and digestibility of cattle in a multitrait analysis. Data from digestibility of 72 heifers aging approximately 18 months and with 250 kg of live weight were used. A completely randomized experiment was conducted at Embrapa-Centro de Pesquisa de Pecuária do Sudeste, São Carlos, SP, Brazil, from 1988 to 1989, considering the combination of nine treatments organized in a factorial 3 x 3 (three genetic groups: Canchim, Nellore and ½ Canchim + ½ Nellore, and three dietary levels of crude protein: 6, 10 and 13%, and eight heifers each). The end points analyzed were feed consumption (g/kg0.75), digestible energy, retained nitrogen (RN), RN (mg/kg0.75), and digestibilities of dry matter, crude protein, neutral detergent fiber and acid detergent fiber. The minimum sample size required to detect significant differences (delta) between mean vectors was calculated by a SAS program considering a t-variate normal distribution with mean vector 0 and covariance matrix sigma, Hotelling's T² test, F-distribution with noncentrality parameter (<FONT FACE=Symbol>d²D)</FONT> and combinations of type I error (alpha) and power of the test (1 beta) and delta. The calculated n ranged from 6 to 47, and was affected more by variation in delta than variation in alpha or power of the test.

F-distribution; minimum significative difference; Hotelling's

TAMANHO AMOSTRAL ENVOLVENDO MÚLTIPLOS CARACTERES:

UMA APLICAÇÃO À DIGESTIBILIDADE EM BOVINOS¹

ALFREDO RIBEIRO DE FREITAS², AIRTON MANZANO² E SÉRGIO NOVITA ESTEVES³

RESUMO - O objetivo deste trabalho foi estimar o tamanho amostral mínimo (n) para comparar tratamentos em experimentos de consumo e de digestibilidade com bovinos, envolvendo múltiplos caracteres. Foram utilizados dados de digestibilidade de 72 novilhas com média de 18 meses de idade e 250 kg de peso. O experimento foi realizado na Embrapa-Centro de Pesquisa de Pecuária do Sudeste, São Carlos, SP, de 1988 a 1989, em delineamento inteiramente casualizado, com nove tratamentos organizados em esquema fatorial 3 x 3 (três grupos genéticos: Canchim, ½ Canchim + ½ Nelore e Nelore, e três níveis de proteína bruta: 6, 10 e 13%, com oito repetições cada, sendo a unidade experimental a novilha). Foram analisados consumo de ração (g/kg^0,75) por quilograma de peso metabólico, energia digestível, nitrogênio retido (NR), NR (mg/kg^0,75) e digestibilidades da matéria seca, proteína bruta, fibra em detergente neutro e fibra em detergente ácido. O valor mínimo de n, que permite detectar diferenças significativas (D) entre vetores de médias de tratamentos, foi obtido por meio de um programa SAS (Statistical Analysis System), considerando modelo de distribuição normal t-variada, média zero e matriz de covariância S, estatística T² de Hotelling, distribuição F com parâmetro de não centralidade (d²_D), erros do tipo I (a), poder do teste (1 b) e D. O valor de n variou de 6 a 47, sendo mais influenciado por alteração nos valores de D, do que nos valores de a e poder do teste.

Termos para indexação: distribuição F, diferença mínima significativa, Hotelling.

SAMPLE SIZE INVOLVING MULTIPLE TRAITS:

AN APPLICATION TO DIGESTIBILITY IN CATTLE

ABSTRACT - The objective of this work was to estimate minimum sample size (n) for comparison of treatments in experiments of consumption and digestibility of cattle in a multitrait analysis. Data from digestibility of 72 heifers aging approximately 18 months and with 250 kg of live weight were used. A completely randomized experiment was conducted at Embrapa-Centro de Pesquisa de Pecuária do Sudeste, São Carlos, SP, Brazil, from 1988 to 1989, considering the combination of nine treatments organized in a factorial 3 x 3 (three genetic groups: Canchim, Nellore and ½ Canchim + ½ Nellore, and three dietary levels of crude protein: 6, 10 and 13%, and eight heifers each). The end points analyzed were feed consumption (g/kg^0.75), digestible energy, retained nitrogen (RN), RN (mg/kg^0.75), and digestibilities of dry matter, crude protein, neutral detergent fiber and acid detergent fiber. The minimum sample size required to detect significant differences (D) between mean vectors was calculated by a SAS program considering a t-variate normal distribution with mean vector 0 and covariance matrix S, Hotelling's T² test, F-distribution with noncentrality parameter (d²_D) and combinations of type I error (a) and power of the test (1 b) and D. The calculated n ranged from 6 to 47, and was affected more by variation in D than variation in a or power of the test.

Index terms: F-distribution, minimum significative difference, Hotelling's.

INTRODUÇÃO

Em ensaios de digestibilidade em bovinos, fatores como o uso de animais fistulados, instalações, disponibilidade de animais, entre outros, limitam ou dificultam a condução de pesquisas na área. Deste modo, deve-se determinar o número de animais a ser utilizado nos experimentos, de modo a fazer inferências precisas.

É fundamental determinar o número de indivíduos necessários (n) para testar os efeitos entre tratamentos, quando um conjunto de t caracteres é analisado no animal. Para a determinação de n, admitiu-se que a amostra era proveniente de uma população com distribuição normal multivariada, sendo adotado o modelo Y=Xb + e_i de dimensões n x p, n x k, k x p e n x p, respectivamente, sendo p o número de variáveis dependentes, k o número de parâmetros e n o número de observações. Definem-se ainda y_ij =(y_ij⁽¹⁾, y_ij⁽²⁾,..., y_ij^(p)) como sendo o vetor de resposta do j-ésimo indivíduo do i-ésimo tratamento e avaliado em p variáveis e, m_i (i = 1, ..., t), o vetor de média p x 1 associado ao tratamento i.

O objetivo deste trabalho foi estimar o número de indivíduos necessários para detectar significância entre vetores de médias de tratamentos em experimentos de consumo e de digestibilidade em bovinos.

MATERIAL E MÉTODOS

Obtenção dos dados

Foram utilizados dados de digestibilidade de 72 novilhas com média de idade de 18 meses e peso de 250 kg. O experimento, em que a unidade experimental foi representada pela novilha, foi realizado na Embrapa-Centro de Pesquisa de Pecuária do Sudeste (CPPSE), São Carlos, SP, de 1988 a 1989, e constou de três períodos. Em cada período, com 24 animais, foram avaliados três grupos genéticos (Canchim, ½ Canchim + ½ Nelore e Nelore) e uma dieta. O primeiro período (6% de PB) iniciou em fevereiro de 1988, o segundo (10% de PB), em agosto de 1988, e o terceiro (13% de PB), em maio de 1989. A relação concentrado: volumoso das rações de cada período foi de 25:75, 40:60 e 60:40, respectivamente. Cada período constou de 121 dias, sendo 21 para adaptação dos animais à dieta e ao confinamento, 90 para avaliação do consumo voluntário, e dez para ganho de peso e conversão alimentar, sendo três para adaptação e sete para colheita de fezes e urina, para determinação dos coeficientes de digestibilidade e do balanço de N. Os caracteres analisados foram o consumo de ração por quilograma de peso metabólico (kg^0,75) (CPM); diges-tibilidade da energia bruta (%) (DEB); digestibilidade da fibra em detergente ácido (%) (DFDA); digestibilidade da fibra em detergente neutro (%) (DFDN); digestibilidade da matéria seca (%) (DMS); digestibilidade da proteína bruta (%) (DPB); N retido (g/dia) (NR) e N retido por quilograma de peso metabólico (mg/kg^0,75) (NPM).

Método de análise

Os dados de digestibilidade, considerando confundimento entre efeitos de períodos e dietas, foram analisados por meio do procedimento GLM do SAS (SAS, 1993), conforme modelo abaixo:

em que

y_ijk = valor observado no k-ésimo animal do j-ésimo grupo genético da i-ésima dieta;

m = efeito médio;

D_i = efeito principal da dieta;

G(D)_j(i) = efeito do j-ésimo grupo genético dentro da i-ésima dieta;

b(X_ijk ^¾) = coeficiente de regressão linear (b) entre y_ijk e X_ijk (peso do k-ésimo animal da j-ésima dieta do i-ésimo grupo genético); é a média dos valores de X;

e_jk(i) = efeito aleatório, suposto ter distribuição normal, média zero e variância ,e_i ~ N (0, s²).

Para propósitos de estimação do tamanho amostral, interessa obter a matriz de covariância amostral (S) e sua inversa (S^-1), sem a influência de tratamentos específicos; assim, por meio do procedimento GLM eliminaram-se dos dados observados os efeitos devidos a D_i e G(D)_j(i) e b(X_ijk ^¾).

Para testar a hipótese da igualdade de dois vetores m₁ e m₂, H₀: m₁= m₂ versus H₁: m₁ ¹ m₂, tem-se:

A estatística F possui graus de liberdade p e (n₁ + n₂¾ p ¾ 1) e parâmetro de não-centralidade d² = n₁n₂(m₁ - m₂)¢ å^-1(m₁ - m ₂) /(n₁ + n₂). Para um dado a, rejeita-se H₀ se F > F(p, n₁ + n₂¾ p ¾ 1; d²).

Uma vez que não se conhecem todos os valores de m₁, ... , m_t e S para os quais se deseja aceitar ou rejeitar H₀ na comparação de dois vetores i e j, especificou-se um vetor de diferença mínima significativa (D), em valor absoluto, dada por ½m_i ^¾m_j½ > D, para algum j ¹ k (Scheffé, 1959), que se deseja detectar a um nível a de probabilidade e poder do teste (1 ^¾b). O valor mínimo de d²sujeito à restrição ½m_j ^¾m_k½ = D, definido por d²_D^,é igual a nD²/max _j<k {s²_j+s²_k^¾ 2s_jk }, em que s²_j e s²_k (j < k) são as variâncias e s_jk, a covariância, associadas às características j e k, respectivamente. Vonesh & Schork (1986) demonstram que essa expressão é apropriada para estimar o tamanho amostral n.

Estimativas de n, para t > 2 tratamentos , em função de v₁, v₂ , a, poder do teste (1 ^¾b) e de d²_D, foram obtidas iterativamente por meio de um programa SAS que utiliza as funções FPROB e FINV (Hardison et al., 1983; Vonesh & Schork, 1986 ). O valor de n foi obtido em função de a (0,01 e 0,05), poder do teste (0,95 e 0,90) e d²_D. Para o cálculo da D, considerou-se os valores 0,5s, 1,0s, 1,5s, 2,0s e 2,5s, sendo s vetor de desvio obtido da matriz de covariância (Tabela 1).

RESULTADOS E DISCUSSÃO

Da matriz de variância-covariância dos dados (Tabela 1), calculou-se max _j<k (s²_j+s²_k^¾ 2s_jk ) = (217,4826 + 50,8182 ^¾ 2 x 15,8558) = 236,6392. Utilizando-se este valor e ainda = 1,0 para = 14,7473, obteve-se o valor de d²_D(função de n, valor a ser estimado), como sendo d²_D= 0,9204n.

As Tabelas 2 a 4 apresentam as estimativas de n para a variando de 0,01 a 0,05, poder do teste (1^¾ b) de 0,80 a 0,95 e D de 1,0s a 2,5s. Uma vez que estas estimativas foram obtidas considerando a variabilidade de cada caráter, sem a influência de tratamentos específicos, e ainda os dois maiores valores da matriz de variância-covariância, inferências sobre os valores de n podem ser feitas para o conjunto de caracteres estudado (Vonesh & Schork, 1986).

Como exemplo, seja o planejamento de um experimento de digestibilidade com bovinos em que se deseja comparar 10 tratamentos, considerando ainda a avaliação dos caracteres citados neste trabalho, que são normalmente utilizados neste tipo de pesquisa. Se pretender detectar diferença significativa entre dois vetores de médias de tratamentos a 5% de probabilidade (a = 0,05), considerando que a diferença entre os dois vetores seja no máximo igual a um desvio-padrão (D = 1,0s) e poder do teste F (1-b) igual a 0,80, há necessidade de 28 indivíduos (Tabela 2). Para esta mesma situação, ao adotar maior rigor na comparação de vetores de médias de tratamentos, ou seja, diferença significativa a 1% de probabilidade (a = 0,01) e poder do teste igual a 0,95, seriam necessários 47 indivíduos, conforme observa-se na Tabela 4.

Para número de tratamentos superior a 10, uma aproximação razoável do número de indivíduos a ser utilizado, independentemente de qual seja a combinação de a, b e D, seria acrescentar 1,0 (um) ao valor de n das Tabelas de 2 a 4, para cada tratamento adicional. Por exemplo, para detectar diferença mínima significativa entre vetores de médias de tratamentos para D =1,0s, a = 0,02 e poder do teste igual a 0,95 (Tabela 3), considerando-se 12 tratamentos, deve-se utilizar 45 indivíduos.

A Tabela 5 apresenta as estimativas dos vetores da média dos caracteres e da diferença mínima significativa de 1,0s a 2,5s, obtidas dos dados observados. Com base nesta situação, quando se deseja detectar diferença significativa entre duas médias de um experimento com oito tratamentos, mantendo-se o poder do teste a 0,90 e níveis de probabilidade (a) variando de 0,05 a 0,01, verifica-se a necessidade de 13 a 16 indivíduos, caso a diferença entre os dois vetores seja bastante evidente (D = 2,5s, por exemplo), ou 30 a 38 indivíduos, quando a diferença for pequena (D = 1,0s). Assim, verifica-se que o número de indivíduos a ser considerado é bastante influenciado pela magnitude da D entre os vetores a serem comparados. Vale salientar que as estimativas obtidas nesse trabalho são válidas para experimentos de digestiblidade envolvendo os oito caracteres aqui estudados, admitindo-se ainda a pressuposição de homogeneidade da matriz de variância-covariância entre os tratamentos.

CONCLUSÕES

1. O número de indivíduos necessários para detectar significância entre vetores de médias de tratamentos de dados de digestibilidade com bovinos é influenciado pelo grau adotado da diferença mínima significativa do erro do tipo I e do poder do teste.

2. O número de indivíduos é mais influenciado por alteração no grau da diferença mínima significativa do que no grau do erro tipo I e do poder do teste.

REFERÊNCIAS

• HARDISON, C.D.; QUADE, D.; LANGSTON, R.E. Nine functions for probability distributions. In: SAS INSTITUTE. SUGI supplemental library user´s guide Cary, NC, 1983. p.229-236.
• SAS INSTITUTE. SAS/STAT user´s guide: statistics, versão 6. 4.ed. Cary, NC, 1993. v.2, 1686p.
• SCHEFFÉ, A. The analysis of variance Toronto: John Wiley, 1959. 477p.
• VONESH, E.F.; SCHORK, M.A. Sample size analysis of repeated measurement. Biometrics, Baltimore, v.42, n.3, p.601-610, Sept. 1986.

Datas de Publicação

• Publicação nesta coleção
  21 Nov 2005
• Data do Fascículo
  Nov 1999