[Biomod-commits] PA strategy: disk

[semra yalcin]

Dear Damien,

I would like to generate PA by using disk strategy in Biomod2. But I got
error while doing this. When I select PA.dist.min = 0, it can select PAs.
On the other hand, when selecting a value bigger than zero for PA.dist.min,
it shows an error *below*.  My layers' linear units are meter, can it be a
problem?

Thank you in advance,

Semra

------------------------------------------------------
*Used codes (similar to default code) :*

*myBiomodData <- BIOMOD_FormatingData(resp.var = myResp,
expl.var = datacurr,
resp.xy = myRespXY,
resp.name = myRespName,
PA.nb.rep = 5,
PA.nb.absences = 10000,*
*PA.strategy = 'disk',*
* PA.dist.min = 20000)  **# 20 km; there are enough pixels for selecting
10000 PA with this rule*
*
*
------------------------------------------------------

Error message:

* |
|                                        |   0%Error in setValues(r, values
= callGeneric(getValues(e1), e2)) :
  length(values) is not equal to ncell(x), or to 1
Calls: BIOMOD_FormatingData ... disk.pseudo.abs.selection -> .local -> > ->
> -> setValues -> setValues*

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

*Layers:*

*> datacurr
class       : RasterStack
dimensions  : 710, 1675, 1189250, 7  (nrow, ncol, ncell, nlayers)
resolution  : 1000, 1000  (x, y)
extent      : -125022.9, 1549977, 3967607, 4677607  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=36 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0
+units=m +no_defs
names       : aspect, bio04, bio05, bio06, bio12, bio14, wetness
min values  : -1.00, 5240.00, 93.00, -250.00, 252.00, 0.00, 0.52
max values  : 360, 10293, 415, 85, 2184, 84, 22

> myResp.ras
class       : RasterLayer
dimensions  : 710, 1675, 1189250  (nrow, ncol, ncell)
resolution  : 1000, 1000  (x, y)
extent      : -125022.9, 1549977, 3967607, 4677607  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=36 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0
+units=m +no_defs
data source : /home/semra/biomod/speciesdata/kn500_.img
names       : kn500_
values      : 1, 1  (min, max)
Raster Attribute Table
     fields : ID COUNT BinValues Value
        min :  0   500         1     1
        max :  0   500         1     1*
*
*
------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/biomod-commits/attachments/20121009/41abdb8b/attachment.html>