Displaying 2 results from an estimated 2 matches for "l10area".
2009 Sep 01
1
understanding the output from gls
...ng "Correlation:"? At first, I thought that this was
describing possible correlations among the predictor variables, but then I
saw that it also included the model intercept. What do these correlation
value mean?
##More detailed information
##function calls:
sppl.i.xx = gls(all.all.rch~l10area+newx, data = gtemp, method="ML")
sppl.i.ex = gls(all.all.rch~l10area+newx, data = gtemp, method="ML",
correlation = corExp(c(20,.8), form=~x+y|area, nugget=TRUE))
##model summaries
> summary(sppl.i.xx)
Generalized least squares fit by maximum likelihood
Mo...
2009 Sep 22
1
odd (erroneous?) results from gls
...se I have made an error, in the
hope that one of you folks might be able to correct me. Thanks in advance
for your time.
The issue is in 2 parts.
(A) I've used gls to fit a model with two fixed effects and a corExp
object. By my count, this fitting process estimates 5 parameters:
(Intercept), l10area, newx, range, and nugget. With 118 total df, there
should be 118 - 5 = 113 residual df. However, the output from summary.gls
reports 115 residual degrees of freedom. Is this an error in summary or
gls, or is there an error in my count?
(B) Summary.gls reports logLik=-273.6. Using my count of 5 est...