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.glsGeneralized 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._