Hi,

I was trying to do cross-validation using the crossval function (bootstrap package), with the following code:

---------------------------------------------------------------------------------------------------------

 theta.fit <- function(x,y){     
      model <- svm(x,y,kernel = "linear") 
 }
 
 theta.predict <- function(fit,x){
      prediction <- predict(fit,x)            
      return(prediction)
 }
 
 x <- matrix(rnorm(5100),102,50)
 rownames(x) <- paste('a',1:102,sep='')
 colnames(x) <- paste('b',1:50,sep='')
 y <- factor(sample(1:2,102,replace=T))
 
 results <- crossval(x,y,theta.fit,theta.predict)   # LOOCV 

-----------------------------------------------------------------------------------------------------------

I get the following error:

Error in scale(newdata[, object$scaled, drop = FALSE], center = object$x.scale$"scaled:center",  : 
  (subscript) logical subscript too long

It seems to work alright if I use 10 fold cross validation (e.g. results <- crossval(x,y,theta.fit,theta.predict,ngroup = 10), but gives the error for LOOCV.

What am I doing wrong?

thanks!

My session info is:
> sessionInfo()
R version 2.7.1 (2008-06-23) 
i386-pc-mingw32 

locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] rpart_3.1-41     lattice_0.17-8   ROCR_1.0-2       gplots_2.6.0    
[5] gdata_2.4.2      gtools_2.5.0     e1071_1.5-18     class_7.2-42    
[9] bootstrap_1.0-21

loaded via a namespace (and not attached):
[1] grid_2.7.1  tools_2.7.1


---------------------------------------------------------------------------------------------------------



      
	[[alternative HTML version deleted]]

