Hi Kristen

Just type
AICc.lme

at the command prompt. That shows the code that is used to build up the AIC
formula. In this case, the relevant bit looks like
        if (identical(nobs, NULL)) {
            n <- nrow(mod$fitted)
        }
        else {
            n <- nobs
        }
        LL <- logLik(mod)[1]
        K <- attr(logLik(mod), "df")

So it uses the total number of observations, unless you specify the
argument nobs in the call to AICc.

hth


On Mon, Feb 4, 2013 at 7:38 PM, Kristen Gorman <kgorman@sfu.ca> wrote:

> Dear all,
> I am curious, how would I find out the exact equation used by the
> AICcmodavg package in R for calculating AICc with lme models? I understand
> what the AICc equation is from Burnham and Anderson (2002), but I am unsure
> as to how the package deals with the number of observations within lme
> class models given random effects.
>
> Thanks for any guidance,
>
> Kristen Gorman
>
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
>
>


-- 
Drew Tyre

School of Natural Resources
University of Nebraska-Lincoln
416 Hardin Hall, East Campus
3310 Holdrege Street
Lincoln, NE 68583-0974

phone: +1 402 472 4054
fax: +1 402 472 2946
email: atyre2@unl.edu
http://snr.unl.edu/tyre
http://aminpractice.blogspot.com
http://www.flickr.com/photos/atiretoo

	[[alternative HTML version deleted]]

