Displaying 1 result from an estimated 1 matches for "tpmen".
Did you mean:
tpmex
2011 Feb 26
0
Help: Error en model.frame.defaul
...t 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=&quo...