[R] The profile log-likelihood problem
Duncan Mackay
dulcalma at bigpond.com
Fri Nov 22 00:18:41 CET 2013
Hi
Without knowing anything about the data and the warning messages, however
the p.vec value seems a little strange.
try
p.vec= seq(1.35, 1.75, by = 0.1)
It could be that you have run out of memory to do 50 sets of calculations
If ok with the result then repeat with a narrower focus. Most probably
quicker
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mackay at northnet.com.au
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of smart hendsome
Sent: Thursday, 21 November 2013 12:43
To: r-help at r-project.org
Cc: roslinazairimah at ump.edu.my
Subject: [R] The profile log-likelihood problem
Hi, Im still new to R and I have a problem regarding the code below. I don't
know how to solve that. I want to compute profile log-likelihood using the
rainfall data. I use the code below from Dunn:
"
tweedie.profile(Amount~Year*Month,p.vec=seq(1.35,1.75,length=50),method="ser
ies",do.ci=TRUE) "
I got the message below when I run the code above. Anyone know why the p.vec
taken from 1.35 to 1.75? Can someone help me?
"1.35 1.358163 1.366327 1.374490 1.382653 1.390816 1.398980 1.407143
1.415306 1.423469 1.431633 1.439796 1.447959 1.456122 1.464286 1.472449
1.480612 1.488776 1.496939 1.505102 1.513265 1.521429 1.529592 1.537755
1.545918 1.554082 1.562245 1.570408 1.578571 1.586735 1.594898 1.603061
1.611224 1.619388 1.627551 1.635714 1.643878 1.652041 1.660204 1.668367
1.676531 1.684694 1.692857 1.701020 1.709184 1.717347 1.725510 1.733673
1.741837 1.75 ..................................................Done.
No valid values of the likelihood computed: smooth aborted
Consider trying another value for the input method.
Error in if ((xi.max > 0) & (xi.max < 1)) { :
missing value where TRUE/FALSE needed
In addition: There were 50 or more warnings (use warnings() to see the first
50)"
Thanks for yourattention and help.
Regard,
Zuhri
[[alternative HTML version deleted]]
More information about the R-help
mailing list