[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

*To*: <asreml@chiswick.anprod.csiro.au>*Subject*: rescaling design matrix*From*: richardm@calm.wa.gov.au*Date*: Thu, 9 Mar 2000 15:54:40 +0800*Reply-To*: <richardm@calm.wa.gov.au>*Sender*: asreml-owner@lamb.chiswick.anprod.csiro.au

Hi all, I am attempting to run a simple univariate "tree" model. Asreml has put out a warning message : WARNING: 1 singularities in AI matrix. The suggestion is then made that the design matrix be rescaled. I am a little uncertain as to how this may be achieved. Can this be done from within asreml? Any pointers in the right direction would be greatly appreciated. .as file is as follows Merridin P95 tree analysis tree !p fem 25 !I male rep 20 prov 3 !a oil ellmer95.prn ellmer95.prn oil~ mu rep prov !r tree relevent part of .asr file is below: ASREML [ 1 Jul 1999] Merredin P95 tree analysis 09 Mar 2000 11:49:27.360 8.00 Mbyte C:\asreml\ellmer95.as A-inverse retrieved from ainverse.bin PEDIGREE [ellmer95.prn ] has 1303 identities, 2581 Non zero elements Reading ellmer95.prn FREE FORMAT skipping 0 lines Univariate analysis of oil Using 1278 records [of 1278 read from 1278 lines of ellmer95.prn ] Model term Size Type COL Minimum Mean Maximum #zero #miss 1 tree 1303 Direct 1 2.000 663.9 1303. 0 0 2 fem 25 Factor 2 1 13.2574 25 0 0 3 male 1 Covariat 3 0.000 0.000 0.000 1278 0 4 rep 20 Factor 4 1 10.6729 20 0 0 5 prov 3 Factor 5 1 2.2926 3 0 0 6 oil 1 Variate 6 0.1000 2.653 5.900 0 0 Forming 1328 equations: 25 dense Initial updates will be shrunk by factor 0.548 NOTICE: 2 (more) singularities, LogL=-506.043 S2= 0.19031E-04 1256 df 0.5090E+05 1.000 LogL=-506.043 S2= 0.19019E-04 1256 df 0.5093E+05 1.000 Final parameter values 50933. 1.0000 WARNING: 1 singularities in AI matrix. Source Model terms Gamma Component Comp/SE % C tree 1303 1303 50933.5 0.968694 25.06 0 P Variance 1278 1256 1.00000 0.190188E-04 0.00 0 S WARNING: Code B - fixed at a boundary (!GP) C - Constrained by user (!CON) S - Singular Information matrix S means there is no information in the data for this parameter. Very small components with Comp/SE ratios of zero sometimes indicate poor scaling. Consider rescaling the design matrix in such cases Analysis of Variance DF F-incr F-adj StndErrDiff 4 rep 20 35.83 8.03 0.1515 5 prov 2 6.99 6.99 0.2748 Solution Standard Error T-value T-prev 5 prov NB 0.721356E-01 0.299888 0.24 WE -0.709678 0.295945 -2.40 -3.53 4 rep thanks in advance Richard ======================================= Richard Mazanec Dept of Conservation and Land Management CALMScience Division Brain street, Manjimup, 6258. Phone +61 8 9771 7995 Fax +61 8 9777 1183 Email: richardm@calm.wa.gov.au -- Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml

- Prev by Date:
**genetic groups** - Next by Date:
**Running ASReml in MS-DOS mode** - Prev by thread:
**genetic groups** - Next by thread:
**Re: rescaling design matrix** - Index(es):