[R] Error message

[Bert Gunter]

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