Hi and thanks to all R developers and helpers,

I wanrt to take  the 1 and 3 dimensions of this object to create a new one without the years and with the next functions to do some specific calculations:

> str(PMpa)
 num [1:499105, 1:60, 1:12] 29.8 55.8 29.7 25.1 25 ...
 - attr(*, "dimnames")=List of 3
  ..$ punt: NULL
  ..$ any : chr [1:60] "1950" "1951" "1952" "1953" ...
  ..$ mes : chr [1:12] "Gener" "Febrer" "Març" "Abril" ...


I'm writing a function where appears the next error message:

Error en model.frame.default(formula = x ~ c("1950:2009"), drop.unused.levels = TRUE) : 
variable lengths differ (found for 'c (1950:2009)')


Script:

library(robustbase)
load("PMpa.Rdat")

R2 <- function(x) var(predict(x))/var(predict(x)+resid(x))

TPMEN<-t(apply(PMpa,c(1,3),function(x){
kk<-lmrob(x~c("1950:2009")); 
c(kk$coefficients[2],R2(kk), (1-summary(kk)$coefficients[2,4])*100, (confint(kk,level=.95)*10)[2,])}))

dimnames(TPMEN)[[3]][1:3]<-c("Pendent","R2","Significancia")

save(TPMEN,file="TPMEN.Rdat")

working:

> TPMEN<-t(apply(PMpa,c(1,3),function(x){
+ kk<-lmrob(x~c("1950:2009")); 
+ c(kk$coefficients[2],R2(kk),(1-summary(kk)$coefficients[2,4])*100,(confint(kk,level=.95)*10)[2,])}))
Error en model.frame.default(formula = x ~ c("1950:2009"), drop.unused.levels = TRUE) : 
  las longitudes de las variables son diferentes (encontrada para 'c("1950:2009")')


Thanks to all for your time and help.


	[[alternative HTML version deleted]]

