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.sook...@gmail.com> > Just realized was because I had an object called AICc. > > On Sun, Sep 4, 2011 at 2:08 AM, Roland Sookias <r.sook...@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]] ______________________________________________ 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.
