calculation of AIC
On 2/02/2009, at 8:41 AM, Adam D. I. Kramer wrote:
On Sun, 1 Feb 2009, orzack wrote:
Speaking of this, does anybody know how to change the default rounding for glm (and lmer) OR for an R session in general (e.g., so that a regular call to glm would generate AIC values with more digits)?
1/7
[1] 0.1428571
options(digits=22) 1/7
[1] 0.1428571428571428
options(digits=23)
Error in options(digits = 23) : invalid 'digits' parameter, allowed 1...22
options(digits=22)
...but this of course won't help if there is explicit rounding programmed into glm/lmer. I also do not understand what would motivate this code, instead of a more straightforward round(aic,0).
R is open source. Look at the code (for print.glm()). Or better,
look at
args(print.glm)
and you will see that these are ``x'' and ``digits''.
This suggests that
print(fit,digits=7)
or
print(fit,digits=17)
(where ``fit'' is the object returned by glm()) may get you somewhere.
You could also simply type ``fit$aic'' to see the AIC printed to
whatever
number of significant digits your (global) options are set for.
cheers,
Rolf Turner
######################################################################
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}