[Biomod-commits] PA strategy: disk

[Damien Georges]

Dear Semra,

Your can work with rasters in meters or whatever units you want if your 
resp.xy are expressed same units.

I made some test and disk PA selection seems to work well on my computer..

Are you sure to have correctly extract your species occurances?
Did you keep only the presences?
Could you check your /myResp object summary?/

If you want, you can send me your raster objects (in private mail), I 
will test if I get the same issue.

All the best,

Damien.


On 09/10/2012 09:13, semra yalcin wrote:
> 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 <http://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/
> /
> /
> ------------------------------------------------------
>
>
>
> _______________________________________________
> Biomod-commits mailing list
> Biomod-commits at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits

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