[R] nlme: problems with gls

Douglas Bates bates at stat.wisc.edu
Tue Dec 18 18:45:28 CET 2001


kjetil halvorsen <kjetilh at umsanet.edu.bo> writes:

> Hola!
> 
> > library(nlme)
> Loading required package: nls 
> > presion.gls <- gls( I(despues-antes) ~ 1 + I(antes-mean(antes)), 
> +                     correlation=corSymm(value=0.5, form= ~ 1 |
> subject) )
> > summary(presion.gls)
> Generalized least squares fit by REML
>   Model: I(despues - antes) ~ 1 + I(antes - mean(antes)) 
>   Data: NULL 
>        AIC      BIC    logLik
>   213.1707 218.4995 -102.5853
> 
> Correlation Structure: General
>  Formula: ~1 | subject 
>  Parameter estimate(s):
> Error in "[<-"(*tmp*, lower.tri(val), value = aux) : 
>         incompatible types
> > presion.gls
> Generalized least squares fit by REML
>   Model: I(despues - antes) ~ 1 + I(antes - mean(antes)) 
>   Data: NULL 
>   Log-restricted-likelihood: -102.5853
> 
> Coefficients:
>            (Intercept) I(antes - mean(antes)) 
>            -14.1000000             -0.1485486 
> 
> Correlation Structure: General
>  Formula: ~1 | subject 
>  Parameter estimate(s):
> Error in "[<-"(*tmp*, lower.tri(val), value = aux) : 
>         incompatible types
> 
> 
> What is the error?

I imagine it is in the print.summary.gls method.

Can you send me a copy of the data set that generated this problem so
I can check the problem?
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list