[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