Displaying 4 results from an estimated 4 matches for "glsst".
Did you mean:
glist
2003 Sep 25
1
Error from gls call (package nlme)
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 course, th...
2017 May 10
2
bug report: nlme model-fitting crashes with R 3.4.0
...# object 'inner_perc_table' not found
gls(distance ~ age + Sex, data = Orthodont,
correlation = corCompSymm( form = ~ 1 | Subject))
# Error in corMatrix.corCompSymm(object) :
# object 'compSymm_matList' not found
#
#Upon debugging, the error is thrown by line 60 of gls code
#
#glsSt <- Initialize(glsSt, dataMod, glsEstControl)
R.version
# _
# platform x86_64-w64-mingw32
# arch x86_64
# os mingw32
# system x86_64, mingw32
# status
# major 3
# minor 4.0
# year 2017
# month 04
# day 21
# svn...
2004 Oct 28
0
Auxilliary args in gls
...iary variable "knots" in my model formula.
It should not be part of the data frame, as it is a vector of a
different length. How can this be done?
The gls() code suggests it may be difficult to keep the auxillary arg
out of the model.frame but in the evaluation frame of the formula:
glsSt <- glsStruct(corStruct = correlation, varStruct =
varFunc(weights))
mfArgs <- list(formula = asOneFormula(formula(glsSt), model,
groups), data = data, na.action = na.action)
if (!missing(subset)) {
mfArgs[["subset"]] <- asOneSidedFormula(Call[["subse...
2017 May 11
0
bug report: nlme model-fitting crashes with R 3.4.0
...;
> gls(distance ~ age + Sex, data = Orthodont,
> correlation = corCompSymm( form = ~ 1 | Subject))
>
> # Error in corMatrix.corCompSymm(object) :
> # object 'compSymm_matList' not found
> #
> #Upon debugging, the error is thrown by line 60 of gls code
> #
> #glsSt <- Initialize(glsSt, dataMod, glsEstControl)
>
> R.version
>
> # _
> # platform x86_64-w64-mingw32
> # arch x86_64
> # os mingw32
> # system x86_64, mingw32
> # status
> # major 3
> # minor 4.0
> # year...