Error message
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
<jdnewmil at 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:<jdnewmil at 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 <15264831 at sun.ac.za>" <15264831 at 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 at 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 at 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