[R] Error message

[Bert Gunter]

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 <bgunter at 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
> <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



-- 

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