[R] nlme: problems with gls

[kjetil halvorsen]

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?

Kjetil Halvorsen
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._