par2 <- 'no' par1 <- '1' library(caret) mywarning <- '' par1 <- as.numeric(par1) if(is.na(par1)) { par1 <- 1 mywarning = 'Warning: you did not specify the column number of the target! The first column was selected by default.' } x <- na.omit(data.frame(t(x))) k <- length(x[1,]) n <- length(x[,1]) x <- as.data.frame(x) for(ii in 1:k) { x[,ii] <- as.numeric(x[,ii]) } myf <- formula(paste(colnames(x)[par1],' ~ .',sep='')) myf lm_grid <- expand.grid(intercept = c(TRUE, FALSE)) fitControl <- trainControl(method = 'repeatedcv', number = 10, repeats = 5) if(par2=='no') { lm_via_caret <- train(myf, data = x, method = 'lm', tuneGrid = lm_grid) } if(par2=='yes') { lm_via_caret <- train(myf, data = x, method = 'lm', tuneGrid = lm_grid, trControl = fitControl) } lm_via_caret load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Multiple Regression',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('',RC.texteval('lm_via_caret$results'),' ',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('',RC.texteval('lm_via_caret$finalModel$tuneValue'),' ',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('',RC.texteval('varImp(lm_via_caret )'),' ',sep='')) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('',RC.texteval('lm_via_caret$finalModel'),' ',sep='')) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|