Resumos
Para avaliar a influência do gene halotano sobre a curva de crescimento de suínos, bem como sua interação com o sexo do animal, foi proposta uma modelagem hierárquica Bayesiana. Nesta abordagem, os parâmetros dos modelos não-lineares de crescimento (Logístico, Gompertz e von Bertalanffy) foram estimados conjuntamente com os efeitos de sexo e genótipos do gene halotano. Foram utilizados 344 animais F2(Comercial x Piau) pesados ao nascer, aos 21, 42, 63, 77, 105 e 150 dias. O modelo Logístico foi aquele que apresentou melhor qualidade de ajuste por apresentar menor DIC (Deviance Information Criterion) que os demais. As amostras das distribuições marginais a posteriori para as diferenças entre as estimativas dos parâmetros do modelo Logístico indicaram que o peso dos machos à idade adulta com genótipo heterozigoto (HalNn) foi superior ao dos homozigotos (HalNN). A título de comparação, também foi considerada a abordagem frequentista tradicional, baseada em dois passos distintos, a qual, por apresentar um menor poder de discernimento estatístico, não mostrou diferenças significativas
A hierarchical Bayesian modeling was used to evaluate the influence of halothane gene and its interaction with sex on pig´s growth curves. Under this approach, the parameters from growth models (Logistic, Gompertz and von Bertalanffy) were estimated jointly with the effects of halothane gene and sex. A total of 344 F2 (Commercial x Piau) animals were weighted at birth, 21, 42, 63, 77, 105 and 150 days in life. The Logistic model has presented the best fit based on DIC (Deviance Information Criterion). Thus, the samples from marginal posterior distributions for the differences between the parameters estimates of Logistic model have indicated that the maturity weight of males with heterozygous genotypes (HalNn) was superior to males with homozygous genotypes (HalNN). In order to realize a comparison with the traditional methodology, the frequentist approach based on two distinct steps also was used, but there was not identified significant differences between growth curve parameter estimates from each group (combinations of halothane genotypes and sex)
improvement; nonlinear models; pig stress syndrome
INTRODUÇÃO:
Uma das estratégias adotadas na suinocultura para aumentar a eficiência da produção e
melhorar a qualidade dos produtos foi o uso de genótipos especializados para ganho de
peso. Devido ao aumento da frequência do alelo mutado do gene do estresse suíno
(halotano) (FUJII et al., 1991FUJII, J. et al. Identification of a mutation in porcine ryanodine
receptor associated with malignant hyperthermia. Science, v.253, n.2, p.448-451,
1991. Disponível em: <http://dx.doi.org/ 10.1126/science.1862346>. Acesso em:
15 nov. 2012. doi: 10.1126/science.1862346.
http://dx.doi.org/
...
) nas linhagens
terminais selecionadas para produção, surgiram hipóteses de que este gene esteja
diretamente relacionado com crescimento (SATHER et al.,
1991SATHER, A.P. et al. The effect of the halothane genotype on pork
production and meat quality of pigs reared under commercial conditions. Canadian
Journal of Animal Science, v.71, n.4, p.959-967, 1991. Disponível em:
<http://dx.doi.org/10.4141/cjas91-117>. Acesso em: 29 set. 2012. doi:
10.4141/cjas91-117.
http://dx.doi.org/10.4141/cjas91-117...
; BRIDI et al., 2008BRIDI, A.M. et al. Efeito da ractopamina e do gênero no desempenho e na
carcaça de suínos de diferentes genótipos halotano. Seminário de Ciências Agrárias,
Londrina, v.29, n.3, p.713-722, 2008.). Animais com
genótipos heterozigotos (HalNn) para o gene halotano são extremamente
sensíveis ao estresse causado, por exemplo, pela mistura com outros animais, transporte,
desmame e manejo (CULAU et al., 2002CULAU, P.O.V. et al. Influência do gene halotano sobre a qualidade da
carne suína. Revista Brasileira de Zootecnia, v.31, p.954-961, 2002. Disponível em:
<http://dx.doi.org/10.1590/S1516-35982002000400019>. Acesso em: 10 dez. 2012.
doi:10.1590/S1516-35982002000400019.
http://dx.doi.org/10.1590/S1516-35982002...
).
Alguns estudos, como o de BAND et al. (2005)BAND, G.O. et al. Relationship between the Porcine Stress Syndrome gene
and carcass and performance traits in F2 pigs resulting from divergent crosses.
Genetics and Molecular Biology, v.28, n.1, p.92-96, 2005. Disponível em:
<http://dx.doi.org/10.1590/S1415-47572005000100016>. Acesso em: 12 jan. 2013.
doi:10.1590/S1415-47572005000100016
http://dx.doi.org/10.1590/S1415-47572005...
,
foram realizados com o intuito de avaliar a influência do gene halotano sobre o peso dos
animais em diferentes períodos de desenvolvimento, do desmame ao abate, porém nenhuma
diferença significativa foi encontrada. KUSEC et al.
(2007)KUSEC, G. et al. Influence of MHS genotype and feeding regime on
allometric and temporal growth of pigs assessed by magnetic resonance imaging.
Livestock Science,v.110, p.89-100, 2007. Disponível em:
<http://dx.doi.org/10.1016/j.livsci.2006.10.007>. Acesso em: 7 out. 2012. doi:
10.1016/j.livsci.2006.10.007.
http://dx.doi.org/10.1016/j.livsci.2006....
realizaram estudo pioneiro ao avaliar o efeito do gene halotano sobre
as estimativas dos parâmetros do modelo Logístico ajustado a dados de peso-idade de
animais cruzados (Pietrain×Hampshire e Large White×German). Estes autores utilizaram uma
metodologia contemplando dois passos independentes, no primeiro, ajustou-se o modelo
não-linear para dados individuais de peso-idade e, no segundo, as estimativas dos
parâmetros provenientes desse ajuste foram consideradas como variáveis dependentes em um
modelo linear (ANOVA), com efeito de gene halotano (com dois níveis: genótipos
HalNN e HalNn). De acordo com BLASCO et al. (2003)BLASCO, A et al. A Bayesian analysis of the effect of selection for
growth rate on growth curves in rabbits. Genetic Selection Evolution, v.35,n p.21-41,
2003. Disponível em: <http://dx.doi.org/10.1051/gse:2002034>. Acesso em: 12
nov. 2012. doi: 10.1186/1297-9686-35-1-21
http://dx.doi.org/10.1051/gse:2002034...
, algumas críticas a esta metodologia são: a estimação
dos parâmetros e a estimação do efeito do gene são realizadas em duas etapas distintas,
portanto a estimação dos parâmetros é realizada sem levar em consideração a diferença
entre os animais quanto à presença do gene; na segunda etapa, as variáveis dependentes
são estimativas de parâmetros de modelos não-lineares, portanto sujeitas à falta de
normalidade requerida pela utilização de testes tradicionais como F e t-Student.
Segundo SILVA et al. (2005)SILVA, F.F. et al. Abordagem Bayesiana da curva de lactação de cabras
Saanen de primeira e segunda ordem de parto. Pesquisa agropecuária brasileira, v.40,
n.1, p.27-33, 2005. Disponível em:
<http://dx.doi.org/10.1590/S0100-204X2005000100004>. Acesso em: 15 nov. 2012.
doi: 10.1590/S0100-204X2005000100004.
http://dx.doi.org/10.1590/S0100-204X2005...
, em estudos
envolvendo ajuste de modelos de regressão não-linear a dados zootécnicos (curvas de
crescimento e de lactação), a inferência Bayesiana pode ser utilizada com sucesso, pois
reduz, substancialmente, o número de estimativas viesadas na presença de poucas
observações, e fornece estimação por intervalo mais precisas em relação à obtida pela
metodologia Frequentista. Além disso, sua principal vantagem é permitir que as
estimativas dos parâmetros do modelo não-linear sejam estimadas simultaneamente com os
efeitos de possíveis fatores que as afetam, isto é, sob o enfoque de modelagem
hierárquica Bayesiana (SILVA et al., 2007SILVA, F.F. et al. Abordagem Bayesiana da curva de lactação de cabras
Saanen de primeira e segunda ordem de parto. Pesquisa agropecuária brasileira, v.40,
n.1, p.27-33, 2005. Disponível em:
<http://dx.doi.org/10.1590/S0100-204X2005000100004>. Acesso em: 15 nov. 2012.
doi: 10.1590/S0100-204X2005000100004.
http://dx.doi.org/10.1590/S0100-204X2005...
). Esta
modelagem permite ainda a comparação estatística direta e simplificada entre curvas
provenientes de diferentes tratamentos por meio da obtenção de uma distribuição para as
diferenças entre as estimativas dos parâmetros.
Diante do exposto, o principal objetivo do presente estudo foi investigar a influência do gene halotano sobre as estimativas dos parâmetros das curvas de crescimento de 344 animais (machos castrados e fêmeas) de uma população F2 (Comercial x Piau), considerando dados de pesagens do nascimento aos 150 dias de idade. Para tanto, foi proposta uma abordagem hierárquica Bayesiana contemplando três diferentes modelos de regressão não-linear (Logístico, Gompertz e von Bertalanffy). A título de comparação com a metodologia proposta, também foi empregada a análise frequentista tradicional baseada em dois passos distintos.
MATERIAL E MÉTODOS:
Os dados utilizados são provenientes de 344 animais F2 (Comercial x Piau), machos castrados e fêmeas, da Granja de Melhoramento de Suínos do Departamento de Zootecnia da Universidade Federal de Viçosa (UFV). Os animais foram pesados nas seguintes fases de desenvolvimento: ao nascimento, aos 21, 42, 63, 77, 105 e 150 (abate) dias. As genotipagens para o gene halotano foram realizadas no Laboratório de Biotecnologia Animal (Labtec) do Departamento de Zootecnia (DZO) - UFV. Observou-se que 149 machos e 141 fêmeas não eram portadores do alelo mutado do gene halotano (HalNN), enquanto que 20 machos e 34 fêmeas eram portadores deste alelo (HalNn).
Independentemente da abordagem utilizada (Bayesiana ou frequentista), foram considerados os modelos de regressão não-linear Gompertz , Logístico e von Bertalanffy
. Em todos estes modelos, yij é o peso do animal i na idade j, Ai é a estimativa de peso na maturidade (peso adulto) do animal i, bi é o parâmetro de integração (não tem interpretação biológica, é uma constante matemática cuja função é ajustar o peso inicial em relação à origem) do animal i e Ki é a taxa de maturidade (representa uma medida de precocidade, ou de velocidade de crescimento) do animal i, ɛij é o termo de erro aleatório, considerado independente e identicamente distribuído, com distribuição normal com média zero e variância σe 2.
Na abordagem hierárquica Bayesiana, o efeito dos quatro grupos considerados (macho e fêmea com genótipos HalNN e macho e fêmea com genótipos HalNn) sobre as estimativas dos parâmetros das curvas de crescimento foram estimados separadamente para cada um dos três modelos considerados (Gompertz, Logístico e von Bertalanffy). No teorema de Bayes, a distribuição a posteriori dos parâmetros é proporcional ao produto da distribuição conjunta dos dados amostrais (função de verossimilhança) e das distribuições a priori para os parâmetros de interesse. Assumindo normalidade dos dados de crescimento, as funções de verossimilhança para os modelos Gompertz, Logístico e von Bertalanffy são dadas, respectivamente, por:
Nas expressões (1), (2) e (3), o vetor θ representa os parâmetros A, b e K para todos os indivíduos, além da variância do erro. No primeiro estágio da hierarquia, foram assumidas as seguintes distribuições a priori para os parâmetros dos modelos de crescimento , e . No segundo e último estágio, foram especificadas as distribuições a priori para os parâmetros das distribuições assumidas no primeiro estágio. Tais distribuições foram as seguintes:
Além disso, especificou-se também uma distribuição a priori para a variância do erro:
.Os valores dos hiperparâmetros
foram escolhidos com base em análises prévias e em revisões de literatura.
A modelagem Bayesiana apresentada foi realizada no software R (R DEVELOPMENT CORE TEAM, 2012R DEVELOPMENT CORE TEAM. R: a language and environment for statistical
computing. Vienna, Austria, 2012. ISBN 3-900051-07-0. Disponível em:
<http://www.R-project.org/>. Online. Acesso em: 5 dez. 2012.
http://www.R-project.org/...
) através do pacote
R2OpenBUGS, no qual foram implementados os algoritmos Gibbs sampler e o
Metropolis-Hastings. Para avaliar a convergência, foram utilizados
os testes de GEWEKE (1992)GEWEKE, J. Evaluating the accuracy of sampling-based approaches to
calculating posterior moments. New York: Oxford University, 1992.
631p. e RAFTERY & LEWIS (1992)RAFTERY, A.E.; LEWIS, S. How many iterations in the gibbs sampler? In:
BERNARDO, J.M. et al. (Eds.). Bayesian statistics. Oxford: University, 1992.
p.763-773. mediante o pacote
Bayesian Output Analysis (boa) do R.
As amostras das distribuições marginais a posteriori para as diferenças
entre as estimativas dos parâmetros A e K (aqueles com interpretação biológica direta)
para cada um dos quatro grupos considerados foram obtidas com o intuito de verificar
significância estatística por meio do intervalo HPD (Highest Posterior
Density) de 95%. Se o valor zero estiver contido neste intervalo, conclui-se
que os parâmetros das duas populações envolvidas no contraste são estatisticamente
iguais. Essa metodologia foi apresentada por SILVA et
al. (2005)SILVA, F.F. et al. Abordagem Bayesiana da curva de lactação de cabras
Saanen de primeira e segunda ordem de parto. Pesquisa agropecuária brasileira, v.40,
n.1, p.27-33, 2005. Disponível em:
<http://dx.doi.org/10.1590/S0100-204X2005000100004>. Acesso em: 15 nov. 2012.
doi: 10.1590/S0100-204X2005000100004.
http://dx.doi.org/10.1590/S0100-204X2005...
para comparar parâmetros de curvas de lactação de cabras referentes
a duas diferentes populações, e posteriormente utilizada para comparar curvas de
crescimento de gado Nelore de diferentes grupos genéticos (SILVA et al., 2007SILVA, N.A.M. et al. Aplicação do método bayesiano na estimação de curva
de Crescimento em animais da raça Nelore. Ceres, v.54, n.312, p.191-198, 2007.
<http://www.redalyc.org/articulo.oa?id=305226740005 >. Acesso em: 3 jan. 2013.
doi: 305226740005.
http://www.redalyc.org/articulo.oa?id=30...
).
Na abordagem frequentista tradicional em dois passos, no primeiro, os modelos Gompertz, Logístico e von Bertalanffy foram ajustados para cada animal por meio do método dos quadrados mínimos ordinários, com soluções obtidas por meio do processo iterativo de Gauss-Newton, utilizando o pacote nls do R. No segundo passo, as estimativas dos parâmetros obtidas no primeiro passo foram consideradas como variáveis dependentes em um modelo ANOVA, sendo os tratamentos as combinações de genótipos e sexo.
Os códigos das rotinas implementadas no software R são de livre acesso em http://www.det.ufv.br/~moyses/links.php.
RESULTADOS E DISCUSSÕES:
Na tabela 1, estão apresentadas, para cada modelo de crescimento considerado, as estimativas (média a posteriori) para os parâmetros populacionais µA (peso adulto) e µK (taxa de maturidade) em cada combinação dos níveis de sexo e genótipos do halotano. Também são apresentadas as estimativas das diferenças entre estes parâmetros para cada um destes grupos (combinações), bem como os intervalos HPDs obtidos. Nesta tabela, observa-se que, para os modelos Gompertz e Logístico, há uma diferença significativa entre a estimativa da média do peso adulto dos machos HalNn e machos HalNN, uma vez que o intervalo HPD para a diferença não contém o valor zero. Além disso, como essas diferenças são maiores que zero, conclui-se que o peso adulto dos machos HalNn é maior do que o peso adulto dos machos HalNN, uma vez que foram consideradas distribuições para a diferença HalNn - HalNN. Para o modelo von Bertalanffy, não foram observadas diferenças significativas.
Por meio do teste de RAFTERY & LEWIS (1992)RAFTERY, A.E.; LEWIS, S. How many iterations in the gibbs sampler? In: BERNARDO, J.M. et al. (Eds.). Bayesian statistics. Oxford: University, 1992. p.763-773., pode-se dizer que a convergência foi atingida, pois o período de aquecimento (burn in), o número total de iterações e o intervalo entre as amostragens (thin) usados foram superiores ao mínimo recomendado por este teste. Além disso, todos os p-valores obtidos para o critério de GEWEKE (1992)GEWEKE, J. Evaluating the accuracy of sampling-based approaches to calculating posterior moments. New York: Oxford University, 1992. 631p. indicaram convergência das cadeias, visto que H0 postula convergência, e esses p-valores variaram de 0,19 a 0,99.
Sob o enfoque Bayesiano, para escolher o modelo que melhor descreveu o crescimento dos
animais, foi utilizado o Deviance Information Criterion (DIC), o qual
indicou o modelo Logístico, uma vez que este apresentou um menor DIC em todos os grupos
considerados. Esse resultado está de acordo com o encontrado por FREITAS (2005)FREITAS, A.R. de. Curvas de crescimento na produção animal. Revista
Brasileira de Zootecnia, v.34, n.3, p.786-795, 2005. Disponível em:
<http://dx.doi.org/10.1590/S1516-35982005000300010>. Acesso em: 2 nov. 2012.
doi: 10.1590/S1516-35982005000300010.
http://dx.doi.org/10.1590/S1516-35982005...
, que, ao estudar sete modelos não-lineares,
considerando-se o ajuste de curvas de crescimento na produção animal, verificou que o
modelo Logístico apresentou melhor ajuste para a maioria das espécies estudadas, dentre
elas a suína. KUSEC et al. (2007)KUSEC, G. et al. Influence of MHS genotype and feeding regime on
allometric and temporal growth of pigs assessed by magnetic resonance imaging.
Livestock Science,v.110, p.89-100, 2007. Disponível em:
<http://dx.doi.org/10.1016/j.livsci.2006.10.007>. Acesso em: 7 out. 2012. doi:
10.1016/j.livsci.2006.10.007.
http://dx.doi.org/10.1016/j.livsci.2006....
, ao estudar
suínos provenientes de diferentes cruzamentos, também apontaram para o modelo Logístico
como sendo o mais adequado para descrever as curvas de crescimento dos animais.
Como o modelo Logístico foi o que melhor se ajustou aos dados, as estimativas dos
parâmetros deste modelo foram utilizadas com a finalidade de avaliar a influência do
gene halotano sobre a curva de crescimento dos animais. Desse modo, conforme mostrado na
tabela 1, conclui-se que o peso adulto dos suínos machos HalNn é maior do que
o peso adulto dos suínos machos HalNN. Em contraste com este trabalho, BAND et al. (2005)BAND, G.O. et al. Relationship between the Porcine Stress Syndrome gene
and carcass and performance traits in F2 pigs resulting from divergent crosses.
Genetics and Molecular Biology, v.28, n.1, p.92-96, 2005. Disponível em:
<http://dx.doi.org/10.1590/S1415-47572005000100016>. Acesso em: 12 jan. 2013.
doi:10.1590/S1415-47572005000100016
http://dx.doi.org/10.1590/S1415-47572005...
, ao avaliarem a influência do
gene halotano sobre características de desempenho (pesos a diferentes idades, inclusive
ao abate) em uma população de suínos F2, resultantes de cruzamentos divergentes
(comercial x Piau), utilizando um modelo linear, concluíram que, para a maioria das
características, não se observaram diferenças significativas entre os genótipos do
halotano. Do mesmo modo, KUSEC et al. (2007)KUSEC, G. et al. Influence of MHS genotype and feeding regime on
allometric and temporal growth of pigs assessed by magnetic resonance imaging.
Livestock Science,v.110, p.89-100, 2007. Disponível em:
<http://dx.doi.org/10.1016/j.livsci.2006.10.007>. Acesso em: 7 out. 2012. doi:
10.1016/j.livsci.2006.10.007.
http://dx.doi.org/10.1016/j.livsci.2006....
não
encontraram diferenças significativas entre os genótipos do halotano sobre as
estimativas dos parâmetros do modelo Logístico, considerando também dois regimes de
alimentação (intensivo e restritivo). Na figura 1, está apresentada a densidade
a posteriori para a diferença entre a estimativa da média do peso
adulto de suínos machos com genótipo heterozigoto e homozigoto dominante (μAmNn
- μAmNN). Nota-se que o valor zero realmente
encontra-se fora (à esquerda) do intervalo HPD de 95% para a diferença entre as médias
do peso adulto, indicando que genótipos HALNn apresentam um maior peso à
idade adulta.
: Densidade a posteriori para a diferença entre a estimativa do peso adulto dos suínos machos com genótipo homozigoto dominante e a estimativa do peso adulto dos suínos machos com genótipo heterozigoto (μAmNn- μAmNN).
Na abordagem Frequentista, após obter as estimativas individuais dos parâmetros dos modelos Gompertz, Logístico e von Bertalanffy, estas foram consideradas como variáveis dependentes em um modelo linear (ANOVA), cujos níveis dos tratamentos foram as combinações entre os genótipos do gene halotano (com dois níveis, genótipos HalNN e HalNn) e sexo (com dois níveis, macho e fêmea). Ao contrário do resultado encontrado pela abordagem Bayesiana, na abordagem Frequentista, não foi possível detectar diferenças significativas entre as estimativas dos parâmetros em relação aos quatro grupos considerados (machos e fêmeas com diferentes genótipos), uma vez todos os p-valores para o teste F foram maiores que 0,05. Este resultado está de acordo com os resultados encontrados pelos autores citados anteriormente, os quais também se basearam em ANOVA para tirar as conclusões a respeito da influência do halotano sobre o crescimento de suínos. Nessa abordagem Frequentista, para escolher o modelo que melhor descreveu o crescimento dos animais, utilizou-se o coeficiente de determinação (R2) e o Critério de informação de Akaike (AIC). Nesta abordagem, também foi verificado que o modelo Logístico foi o que melhor se ajustou aos dados, pois este apresentou, respectivamente, maiores e menores valores para o R2 e AIC.
As médias das estimativas dos parâmetros do modelo Logístico, Gompertz e von Bertalanffy, as quais não diferiram estatisticamente entre grupos, estão apresentadas na tabela 2. Nesta tabela, nota-se ainda que a abordagem Frequentista superestimou as estimativas dos parâmetros das curvas de crescimento, principalmente para o peso adulto (A).
: Média das estimativas dos parâmetros A (peso adulto) e K (taxa de maturidade), para os modelos Logístico, Gompertz e von Bertalanffy (abordagem Frequentista).
CONCLUSÃO:
O modelo que melhor se ajustou aos dados apresentados neste trabalho foi o modelo Logístico, pois apresentou um melhor ajuste, tanto na abordagem Frequentista (via AIC e R2) quanto na Bayesiana (via DIC). Esta última abordagem mostrou-se mais apropriada para avaliar a influência do gene halotano sobre as curvas de crescimento dos animais, permitindo identificar que o genótipo HALNn exerce influência significativa apenas no peso adulto dos machos, sendo que estes apresentam peso superior àqueles de genótipo homozigoto.
- BAND, G.O. et al. Relationship between the Porcine Stress Syndrome gene and carcass and performance traits in F2 pigs resulting from divergent crosses. Genetics and Molecular Biology, v.28, n.1, p.92-96, 2005. Disponível em: <http://dx.doi.org/10.1590/S1415-47572005000100016>. Acesso em: 12 jan. 2013. doi:10.1590/S1415-47572005000100016
» https://doi.org/10.1590/S1415-47572005000100016» http://dx.doi.org/10.1590/S1415-47572005000100016 - BLASCO, A et al. A Bayesian analysis of the effect of selection for growth rate on growth curves in rabbits. Genetic Selection Evolution, v.35,n p.21-41, 2003. Disponível em: <http://dx.doi.org/10.1051/gse:2002034>. Acesso em: 12 nov. 2012. doi: 10.1186/1297-9686-35-1-21
» https://doi.org/10.1186/1297-9686-35-1-21» http://dx.doi.org/10.1051/gse:2002034 - BRIDI, A.M. et al. Efeito da ractopamina e do gênero no desempenho e na carcaça de suínos de diferentes genótipos halotano. Seminário de Ciências Agrárias, Londrina, v.29, n.3, p.713-722, 2008.
- CULAU, P.O.V. et al. Influência do gene halotano sobre a qualidade da carne suína. Revista Brasileira de Zootecnia, v.31, p.954-961, 2002. Disponível em: <http://dx.doi.org/10.1590/S1516-35982002000400019>. Acesso em: 10 dez. 2012. doi:10.1590/S1516-35982002000400019.
» https://doi.org/10.1590/S1516-35982002000400019» http://dx.doi.org/10.1590/S1516-35982002000400019 - FREITAS, A.R. de. Curvas de crescimento na produção animal. Revista Brasileira de Zootecnia, v.34, n.3, p.786-795, 2005. Disponível em: <http://dx.doi.org/10.1590/S1516-35982005000300010>. Acesso em: 2 nov. 2012. doi: 10.1590/S1516-35982005000300010.
» https://doi.org/10.1590/S1516-35982005000300010» http://dx.doi.org/10.1590/S1516-35982005000300010 - FUJII, J. et al. Identification of a mutation in porcine ryanodine receptor associated with malignant hyperthermia. Science, v.253, n.2, p.448-451, 1991. Disponível em: <http://dx.doi.org/ 10.1126/science.1862346>. Acesso em: 15 nov. 2012. doi: 10.1126/science.1862346.
» https://doi.org/10.1126/science.1862346» http://dx.doi.org/ 10.1126/science.1862346 - GEWEKE, J. Evaluating the accuracy of sampling-based approaches to calculating posterior moments. New York: Oxford University, 1992. 631p.
- KUSEC, G. et al. Influence of MHS genotype and feeding regime on allometric and temporal growth of pigs assessed by magnetic resonance imaging. Livestock Science,v.110, p.89-100, 2007. Disponível em: <http://dx.doi.org/10.1016/j.livsci.2006.10.007>. Acesso em: 7 out. 2012. doi: 10.1016/j.livsci.2006.10.007.
» https://doi.org/10.1016/j.livsci.2006.10.007» http://dx.doi.org/10.1016/j.livsci.2006.10.007 - R DEVELOPMENT CORE TEAM. R: a language and environment for statistical computing. Vienna, Austria, 2012. ISBN 3-900051-07-0. Disponível em: <http://www.R-project.org/>. Online. Acesso em: 5 dez. 2012.
» http://www.R-project.org/ - RAFTERY, A.E.; LEWIS, S. How many iterations in the gibbs sampler? In: BERNARDO, J.M. et al. (Eds.). Bayesian statistics. Oxford: University, 1992. p.763-773.
- SATHER, A.P. et al. The effect of the halothane genotype on pork production and meat quality of pigs reared under commercial conditions. Canadian Journal of Animal Science, v.71, n.4, p.959-967, 1991. Disponível em: <http://dx.doi.org/10.4141/cjas91-117>. Acesso em: 29 set. 2012. doi: 10.4141/cjas91-117.
» https://doi.org/10.4141/cjas91-117» http://dx.doi.org/10.4141/cjas91-117 - SILVA, F.F. et al. Abordagem Bayesiana da curva de lactação de cabras Saanen de primeira e segunda ordem de parto. Pesquisa agropecuária brasileira, v.40, n.1, p.27-33, 2005. Disponível em: <http://dx.doi.org/10.1590/S0100-204X2005000100004>. Acesso em: 15 nov. 2012. doi: 10.1590/S0100-204X2005000100004.
» https://doi.org/10.1590/S0100-204X2005000100004» http://dx.doi.org/10.1590/S0100-204X2005000100004 - SILVA, N.A.M. et al. Aplicação do método bayesiano na estimação de curva de Crescimento em animais da raça Nelore. Ceres, v.54, n.312, p.191-198, 2007. <http://www.redalyc.org/articulo.oa?id=305226740005 >. Acesso em: 3 jan. 2013. doi: 305226740005.
» https://doi.org/305226740005» http://www.redalyc.org/articulo.oa?id=305226740005
Datas de Publicação
-
Publicação nesta coleção
Out 2014
Histórico
-
Recebido
26 Set 2013 -
Aceito
24 Jan 2014