[Biomod-commits] Projection() Error

[Brenna Forester]

Hello all,

I have successfully run Projection() on 10 different climate data sets.  I just encountered the following error, however, while projecting using an eleventh data set:

Error in if (sum(ProbData) != 0) { :
  missing value where TRUE/FALSE needed

I see from previous posts that this error can sometimes be avoided by changing FDA=F.  I did so and Projections() ran with no problems.  Ensemble.Forecasting() also ran fine, however ProjectionBestModel() gave the same error as above.  I am using all models except SRE.

Any feedback?  The code and version info are below.

Thank you!
Brenna

> load("Rhodiola_run.Rdata")

> Projection(Proj=hadc21[,3:10],Proj.name='HADC.21k',GLM = T, GBM = T, GAM = T, CTA = T, ANN = T, SRE = F, FDA =T, MARS = T, RF = T, BinRoc=T, BinKappa=T, BinTSS=T, FiltRoc=T, FiltKappa=T, FiltTSS=T, repetition.models=T)

Error in if (sum(ProbData) != 0) { :
  missing value where TRUE/FALSE needed

## Rerun using FDA=F

> Projection(Proj=hadc21[,3:10],Proj.name='HADC.21k',GLM = T, GBM = T, GAM = T, CTA = T, ANN = T, SRE = F, FDA =F, MARS = T, RF = T, BinRoc=T, BinKappa=T, BinTSS=T, FiltRoc=T, FiltKappa=T, FiltTSS=T, repetition.models=T)

> ProjectionBestModel(Proj.name='HADC.21k', Bin.trans=T, Filt.trans=T,method='all')

Error in if (sum(ProbData) != 0) { :
  missing value where TRUE/FALSE needed

> Ensemble.Forecasting(Proj.name= "HADC.21k", weight.method='Kappa', decay='proportional', PCA.median=F, binary=T, bin.method='Kappa', Test=T, repetition.models=T)

Rhodiola

 consensus_HADC.21k_results
$Rhodiola
$Rhodiola$weights
             ANN    CTA    GAM    GBM    GLM   MARS FDA     RF SRE
PA1       0.1373 0.1266 0.1515 0.1585 0.1191 0.1499   0 0.1569   0
PA1_rep1  0.1359 0.1457 0.1457 0.1552 0.1171 0.1361   0 0.1644   0

## Results truncated...

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
I am running BIOMOD Version 1.1-5 on the following system:

> version
               _
platform       i386-pc-mingw32
arch           i386
os             mingw32
system         i386, mingw32
status
major          2
minor          12.1
year           2010
month          12
day            16
svn rev        53855
language       R
version.string R version 2.12.1 (2010-12-16)











-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/biomod-commits/attachments/20110323/d5015d78/attachment.htm>