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 ______________________________________________ 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.
