Displaying 1 result from an estimated 1 matches for "pibtrimestr".
Did you mean:
pibtrimestral
2004 Apr 25
1
ts's in lm()
...rep(3, length(xx))), frequency=12, start=
c(start(xx)[1], (start(xx)[2]-1)*3+1) )
common <- ts.intersect(xx, yy)
colnames(common) <- c(name.x, name.y)
xx <- common[,1]
yy <- common[,2]
sea <- seasonal.dummies(xx)
return(lm(yy ~ 0+sea+xx))
}
> mod1 <- modelo1(PIBtrimestral, gasolina.eq)
Error in "storage.mode<-"(`*tmp*`, value = "double") :
invalid time series parameters specified
lm() works with time series objects made directly with ts(), and the
arguments to modelo1() is of this type.
Kjetil Halvorsen