Another alternative might be to post the the r-sig-mixed-models list gls() is part of nlme.
-- Bert On Tue, Apr 30, 2013 at 11:12 AM, Bert Gunter <bgun...@gene.com> wrote: > Actually, Jeff, it may not be. > > The OP has likely overfit -- as a guess, Depth is essentially > continuous and as.factor() produces too many levels. I would recommend > that the OP find a local statistical consultant to help him, as he > seems out of his depth, statistically. Failing that, he should > probably post to a statistical list like stats.stackexchange.com and > ask for statistical counsel there. > > -- Bert > > On Tue, Apr 30, 2013 at 8:59 AM, Jeff Newmiller > <jdnew...@dcn.davis.ca.us> wrote: >> That is unlikely in the absence of a reproducible example [1]. >> >> [1] >> http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example >> --------------------------------------------------------------------------- >> Jeff Newmiller The ..... ..... Go Live... >> DCN:<jdnew...@dcn.davis.ca.us> Basics: ##.#. ##.#. Live Go... >> Live: OO#.. Dead: OO#.. Playing >> Research Engineer (Solar/Batteries O.O#. #.O#. with >> /Software/Embedded Controllers) .OO#. .OO#. rocks...1k >> --------------------------------------------------------------------------- >> Sent from my phone. Please excuse my brevity. >> >> "Havenga, BS, Mnr <15264...@sun.ac.za>" <15264...@sun.ac.za> wrote: >> >>>Hi there >>> >>>I am a Masters student at the University of Stellenbosch. I have been >>>using R to analyze the data, using the GLS model, of one of my >>>experiments. >>>The problem that I am having is that whenever I run my model using: >>>fit.glsmodel1<-gls(Number~as.factor(Season)+as.factor(Depth)+as.factor(Orientation), >>>data=Number, weights=varPower(), method="ML") >>> >>>I get the error: >>>Error in glsEstimate(object, control = control) : computed "gls" fit >>>is singular, rank 8 >>> >>>I would like to know if anyone can tell me what is causing this error. >>>Any advice will be much appreciated. >>> >>>Kind regards >>>Brendan >>> >>> [[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. >> >> ______________________________________________ >> 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. > > > > -- > > Bert Gunter > Genentech Nonclinical Biostatistics > > Internal Contact Info: > Phone: 467-7374 > Website: > http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm ______________________________________________ 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.