Dear R-List,

I'm using the package biomod to identify areas of potential distribution of a butterfly. I don't have independent data to test, so I'm using a partition of the data 80/20. Maybe I should say as well that I chose to run the Models with only 1 repetitions but 2 runs of Pseudo-absences. 

I have 3 questions.
1. In the tutorial for this package is described the function level.plot. It is useful to plot the results of every model but in the tutorial they ran the Models with "independent data", Is there any possibility to use this function without having to run the models with independent data?

load("pred/Pred_Ag_cing")
level.plot(Pred_Ag_cing[,"RF",1,1],coorxy,title='Ag_cingulata RF')

Error in level.plot(Pred_Ag_cing[, "RF", 1, 1], coorxy, title = "Ag_cingulata RF") : 
  
 data and coordinates should be of the same length 

length(coorxy$Longitude)
[1] 2697182

length(Pred_Ag_cing[,"RF",1,1])
[1] 1350

It gives me an error saying that the coordinates and the data must be the same length, (the truth is I don't get why without running the models with independent data it doesn't get the same length but with independent it does.)


2. Is there any possibility to export the data to make
 the plots as an ascii files so can be imported into GIS?

3. Is the only form to plot the outputs of the models in GIS by extracting the model and using Raster Calculator? If it is the case, how do I get the final models of RandomForest, MARS and GAM?

Any help would be very much appreciated, 

Thanks a lot in advanced,

Liliana Ballesteros
PhD Student
University of Basel





      
	[[alternative HTML version deleted]]

