[R-sig-ME] binary trait
Mohamed Salem
mmisalem at gmail.com
Thu Dec 1 12:00:12 CET 2016
Dears,
I am trying to use MCMCglmm to estimate heritability for binary trait.
I used this model
"
prior <- list(R = list(V = 1, fix = 1), G = list(G1 =
list(V = 1, nu = 1000, alpha.mu = 0, alpha.V = 1)))
model1 <- MCMCglmm(SB ~ 1 + Farm +year, random = ~animal, family =
"ordinal",
prior = prior, pedigree = Ped, data = Data, nitt = 1e+06,burnin = 10000,
thin = 100)
and when I diagnosed the MCMC work by autocorr.diag(model1$VCV)
I found this results
animal units
Lag 0 1.0000000 NaN
Lag 100 0.9786790 NaN
Lag 500 0.9092071 NaN
Lag 1000 0.8360430 NaN
Lag 5000 0.4860764 NaN
how can I avoid this problem?
--
*Mohamed M. I. Salem*
Lecturer
Department of Animal Production
Faculty of Agriculture, Alexandria University
Aflaton St., El-Shatby, P. O. Box 21545
Alexandria, Egypt
Tel: (+203) 5915427, (+203) 5924391
Fax (+203)5922780
e-mail: mmisalem at gmail.com
mohamed.salim at alexu.edu.eg
skype: mmisalem
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list