[R-sig-eco] Help: Error running Projection.raster in BIOMOD
Marco
wul1984 at gmail.com
Sun Feb 20 04:32:59 CET 2011
Dear Community:
I am running Projection.raster on my species, and got one simple error
that I had no idea why.
My running code is simply like this:
>library(BIOMOD)
>setwd("F:/QTP_Layer/model_run/allpresence")
>load(paste(getwd(),"/BIOMOD_run.RData",sep=""))
>qtp_all<-brick("F:/QTP_Layer/data/qtp_all")
>qtp_sel<-stack(qtp_all[[1]],qtp_all[[4]],qtp_all[[15]],qtp_all[[18]],qtp_all[[20]]) ### selected variables
>layerNames(qtp_sel)<-Biomod.material$VarNames
>Sp.env.sel<-qtp_sel
>Projection.raster(RasterProj=Sp.env.sel,Proj.name='qtp.proj',
>GLM = T, GBM = T, GAM = T, CTA = F, ANN = F, SRE = F, quant=0.025, FDA =F, MARS = F, RF = T, BinRoc=F, BinKappa=F, BinTSS=T, FiltRoc=F, FiltKappa=F, FiltTSS=T, repetition.models=T,stack.out=TRUE)
# Here is my raster stack input information:
class : RasterStack
dimensions : 1800, 4200, 5 (nrow, ncol, nlayers)
resolution : 0.008333333, 0.008333333 (x, y)
extent : 70, 105, 25, 40 (xmin, xmax, ymin, ymax)
projection : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
min values : -227 4125 33 11 83
max values : 242 11533 156 2510 8469
#and the error is simply:
Error in x at layers[[i]] : subscript out of bounds
When traced the error back, I found that I am able to generate the
projection for each "g" (algorithm), but then came the error:
pile.proj <- stack(pile.proj, run.name2 = g)
#where pile.proj is a NULL stack, and g is the projection result.
Could anybody explain why?
Thanks in advance~
Sincerely
Marco
--
-------------------------------------------------------
Ph.D candidate,
Chinese Ecological Research Network,
Institute of Geographical Sciences and Natural Resources Research,
Chinese Academy of Sciences,
100101, Beijing
Tel: (86)-10-6488 9809
More information about the R-sig-ecology
mailing list