Bayesian Approach For Animal Breeding Data Analysis
Bayesian Approach For Animal Breeding Data Analysis
G. R. Gowane
ICAR-CSWRI Avikanagar
Choice of estimating VC
REML
DFREML, MTDFREML, VCE, Wombat
Why Bayes?
Combine information from prior to obtain the
posterior distribution.
MORE ACCURACY!!
It is an alternate approach
In Bayesian theory we
know that all problems
are reduced to a single
pathway: we should
look for a posterior
distribution, given the
distribution of the data
and the prior
distribution.
Laplace (1800)
Normal distribution
Many applications, including census
[sampling models]
Gauss (1800)
Least squares
Applications to astronomy [measurement error models]
Prior???
Litter Size of Landrace Pig ~ 10
Spanish Landrace?
Mean = 5
Prior mean = 10
What to do?
Prior information is the
information about the
parameters we want to
estimate that
exists before we perform our
experiment.
MTGSAM
Van Tassell and Van Vleck (1995)
Model Assumptions
y = X + Zu + e
where is the vector of fixed effects associated with
records in y by X, and u is the vector of random effects
associated with records in y by Z, and e is the vector of
random residual effects.
Variance Components:
The MTGSAM posterior mean estimate for
(co)variance components is based on the
expected value of the IW RV
The mean of a (co)variance component is
calculated as the average of expected values
over the length of the post burn-in chain.
Actual
results
are
given in
the
MTGS8
3
MTGS72
Animal effect
Second
animal effect
Pulldat.exe
You need to have pulldat.exe
in the same folder where
MTGS61, MTGS62 and
MTGS63 are present.
Only after files are
subjected to pulldat they
can be further used for
Gibanal or other
analyses.
Open the CMD window and
change path directory to
pulldat
Give the name for output
fille.
Follow the options as given
except at the point Enter
the number of variables
in each record to be read
from MTGS61. Here in
our case, 5 variables
needs to be extracted,
however, the input will
change according to the
data in consideration and
Gibanal.exe
the observed values of the fixed and random
effects are written to unit 61
We can use Gibanal to see
Serial correlation
Burn in
Chain length
Bloco Freqncia
2.18
1
2.46
3
2.73
5
3.01
10
3.29
13
3.57
22
3.85
23
4.13
9
4.41
2
Thank You