Dear Hermann,
If you are actually using a GLMM, then I can believe the numbers
because the LogL value printed is not really what is being minimised.
i.e. it only holds the 'Conditional Normal' part of the model.
If you are fitting a Normal model, then there must be something
going wrong but you have not given enough information for me to detect it.
Arthur
------------- Begin Forwarded Message -------------
X-Authentication-Warning: lamb.chiswick.anprod.csiro.au: petidomo set sender to
asreml-owner@lamb.chiswick.anprod.csiro.au using -f
Date: Wed, 19 Jul 2000 15:04:47 -0400
From: "C. Balde" <cbalden@entelchile.net>
MIME-Version: 1.0
To: asreml <asreml@chiswick.anprod.csiro.au>
Subject: Local maximum of LogL
Content-Transfer-Encoding: 7bit
Dear Arthur, Dear All,
When fitting a univariate model I get the following problem: The LogL
function converge at a maximum value. Sometimes the value decreases in
every iteration and sometimes it first increases and then decreases, but
the final value is not the maximum. Changing the initial values or using
!GU option makes the LogL to increase in every itration, but the final
value is the same as in prior run.
So I think ASREML has found a local maximum, but not the global maximum.
Is this right ? Any sugestion on how to solve this problem ?
Hermann
PS: Following are the LogL values so you can understan my questions.
LogL=-1177.45 S2= 5.1233 802 df 0.1000 0.1000
1.000
LogL=-1177.60 S2= 4.9262 802 df 0.1277 0.1299
1.000
LogL=-1178.34 S2= 4.6714 802 df 0.1730 0.1699
1.000
LogL=-1179.18 S2= 4.4919 802 df 0.2145 0.1949
1.000
LogL=-1179.47 S2= 4.4425 802 df 0.2288 0.1993
1.000
LogL=-1179.50 S2= 4.4371 802 df 0.2305 0.1996
1.000
LogL=-1179.50 S2= 4.4366 802 df 0.2306 0.1996
1.000
----------
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml
------------- End Forwarded Message -------------
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD mailto:Arthur.Gilmour@agric.nsw.gov.au
Principal Research Scientist (Biometrics) fax: <61> 2 6391 3899
NSW Agriculture <61> 2 6391 3922
Orange Agricultural Institute telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA home: <61> 2 6362 0046
ASREML is still free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
To join the asreml discussion list, send the message
subscribe
mailto:asreml-request@chiswick.anprod.CSIRO.au
To send messages to the list, mailto:asreml@chiswick.anprod.CSIRO.au
Asreml list archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<> <> <> <> <> <> <>
While they watched, He was taken up,
and a cloud received Him out of their sight.
And while they looked steadfastly toward heaven as He went up,
behold two men stood by them in white apparel who also said
"Men of Galilee, why do you stand gazing up into heaven?
This same Jesus who was taken up from you into heaven,
will so come in like manner as you saw Him go into heaven."
Acts 1 9:11
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml