search for: glsestimate

Hi I have a huge array with series of data. For each cell in the array I fit a linear model, either using lm() or gls() with lm() there is no problem, but with gls() I get an error: Error in glsEstimate(glsSt, control = glsEstControl) : computed gls fit is singular, rank 2 as soon as there are data like this: > y1 <- c(0,0,0,0) > x1 <- c(0,1,1.3,0) > gls(y1~x1) Error in glsEstimate(glsSt, control = glsEstControl) : computed gls fit is singular, rank 2 of cour...

...yze 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]]