Hello again...,


I have a question regarding the habitat suitability maps in R.

I performed MADIFA on my data, but I can't execute the command:

pred <- predict(mad, map)

_my problem is:
_
MADIFA wants an object of class SpatialPixelsDataFrame (map).

The problem is that I have 3 different areas with 3 different extents , 
so, I have 3 different SpatialPixelsDataFrame with different coordinates 
and different extents.

I tried to put together these 3 SPDF to have 1, but I didn't find the 
right command.... I also tried to put togheter all the rasters (with the 
stack command) at the beginning of the process before to do the SPDFs, 
but, again it was not possible...

Below I reported my MADIFA results:

MADIFA

$call: madifa(dudi = pc_mad, pr = pr_mad, scannf = FALSE, nf = 7)

  

eigen values: 15.58 3.255 0.9557 0.6351 0.2338 ...

$nf: 7 axes saved

  

   vector   length mode     content                       

1 $pr      418353 numeric vector of presence            

2 $mahasu 418353 numeric squared Mahalanobis distances

3 $lw      418353 numeric row weights                   

4 $cw      418353 numeric column weights                

5 $eig     8       numeric eigen values                  

  

   data.frame nrow    ncol content                                          

1 $tab        418353 8     modified array                                   

2 $li         418353 7     row coordinates                                   

3 $l1         418353 7     row normed scores (variance weighted by $pr = 1)

4 $co         8       7     column coordinates                               

5 $cor        8       7     cor(habitat var., scores) for available points

  

  

  

Can anybody help me?

  

Thank you very much again.

  

Michela

-- 
*****************************************************************
Dr. Michela Giusti
ISPRA
Istituto Superiore per la Protezione e la Ricerca Ambientale
Via Brancati, 60
00144 Roma
Tel. +390650074746
michela.giusti@isprambiente.it
*****************************************************************


	[[alternative HTML version deleted]]

