Also, gls1 was fitted with an REML method:you should use an ML method if you
want to compute AIC and friends.

R

2011/9/4 Roland Sookias <r.sookias@gmail.com>

> Just realized was because I had an object called AICc.
>
> On Sun, Sep 4, 2011 at 2:08 AM, Roland Sookias <r.sookias@gmail.com>
> wrote:
>
> > Hi
> >
> > I get the following error when I try and get the AICc for a gls
> regression
> > using qpcR:
> >
> >
> > > AICc(gls1)
> > Loading required package: nlme
> > Error in n/(n - p - 1) : 'n' is missing
> >
> > My gls is like this:
> >
> > > gls1
> > Generalized least squares fit by REML
> >   Model: thercarnmax ~ therherbmax
> >   Data: NULL
> >   Log-restricted-likelihood: 2.328125
> >
> > Coefficients:
> > (Intercept) therherbmax
> >   1.6441405   0.2998528
> >
> > Degrees of freedom: 6 total; 4 residual
> > Residual standard error: 0.1506473
> >
> >
> > Cheers
> >
> > Roland
> >
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



-- 
Renaud Lancelot
EDENext Project, coordinator: http://www.edenext.eu/

CIRAD, UMR15, Campus International de Baillarguet TA A-DIR / B
F34398 Montpellier

Tel.  +33 4 67 59 37 17  -  Fax  +33 4 67 59 37 95
Secr. +33 4 67 59 37 37  - Cell. +33 6 77 52 08 69

	[[alternative HTML version deleted]]

