Hi,
I want to get predictions for the subrace effects for four different
sites (traits), using the following code:
!PART 21 #Full multivariate Four sites
!maxit 12
massflow bawdflow shalflow sso3flow ~ Trait Trait.subrace at(Trait,1).rep,
at(Trait,2).bigrep at(Trait,3).rep,
!r at(Trait,2).plot 38 -at(Trait,3).rep.block,
Trait.csiroID !f mv
1 2 1 !STEP .1
0 0 0
Trait 0 US !+10 !GUFUFFUFFFP
734
0 346
0 0 192
0 0 0 913
Trait.csiroID 2
Trait 0 CORR !+10 !GUUUFFUUUUU
1.08
0.89 0.68
0 0 1.11
267 335 174 245
csiroID
predict subrace
The run converges but it does not produce a .pvs file with the
predictions for subrace, although there are solutions for Trait and
Trait.subrace in the .sln file.
The .as file says:
Forming 6505 equations: 218 dense
Initial updates will be shrunk by factor 0.100
Prediction of at(Trait,3) 1 0 0 -13 0 -1 3 not yet enabled
Prediction of at(Trait,2) 1 0 0 -13 0 -1 2 not yet enabled
Prediction of at(Trait,1) 1 0 0 -13 0 -1 1 not yet enabled
Prediction of at(Trait,2) 1 0 0 -13 0 -1 2 not yet enabled
Prediction of mv_estimates 3922 0 0 -4 0 0 0 not yet enabled
NOTICE: 4540 (more) singularities,
LogL=-4165.24 S2= 1.0000 5063 df
LogL=-4163.48 S2= 1.0000 5063 df
LogL=-4159.68 S2= 1.0000 5063 df
LogL=-4157.94 S2= 1.0000 5063 df
LogL=-4157.71 S2= 1.0000 5063 df
LogL=-4157.71 S2= 1.0000 5063 df
Unable to predict at(Trait,3) 4 0 0 1 -13 0 -1 3
LogL=-4157.71 S2= 1.0000 5063 df
Not all subraces are in all sites. There are 11, 11, 18 and 8 in each
site respectively. What is the problem?
Luis
--
Dr Luis A. Apiolaza
Quantitative Geneticist
CRC for Sustainable Production Forestry
School of Plant Science
University of Tasmania
GPO Box 252-55
Hobart TAS 7001
Australia
phone: +61-3-6226 2213
fax: +61-3-9229 2698
email: mailto:luis.apiolaza@utas.edu.au
'There is no intellectual exercise which is not
ultimately useless'
Jorge Luis Borges
http://www.utas.edu.au/docs/plant_science/
;-)
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml