[Biomod-commits] Biomod.RangeSize

[Ashley Brooks]

Hey!

So I would like to do a Biomod.RangeSize comparing an average of two GCM's
to the current climate data.

Here is what I've done so far:

> load("proj.cccm80a2a/Total_consensus_cccm80a2a")
> load("proj.had80a2a/Total_consensus_had80a2a")
> MyAverageA2A =
rowMeans(cbind(Total_consensus_cccm80a2a[,1,"prob.mean.weighted"],
Total_consensus_had80a2a[,1,"prob.mean.weighted"]))
> MyThresA2A =
mean(consensus_cccm80a2a_results$Cflorida$thresholds["prob.mean.weighted",
c(2:4, 6:8)])
>  MyAverageBinA2A=BinaryTransformation(MyAverageA2A, MyThresA2A)

For the current climate data all that I could think to do was put it in a
similar 'format' as MyAverageBinA2A.

Here is what I tried:

> load("proj.present/Total_consensus_present_Bin")
> Present<-Total_consensus_present_Bin[,1,"prob.mean.weighted"]
> Biomod.RangeSize(CurrentPred=Present, FutureProj=MyAverageBinA2A,
SpChange="Cflorida.2080a2a")
Error in `dimnames<-.data.frame`(`*tmp*`, value = list(NULL, c("Loss",  :
  invalid 'dimnames' given for data frame


I've tried to use this function a couple of different ways (each
incorporating MyAverageBinA2A) and keep getting the same error message.
 Please let me know if you have any suggestions or let me know if I took a
wrong step somewhere.

Thank you so much for all your help!  It is greatly appreciated!
Ashley
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/biomod-commits/attachments/20120227/6c8838d1/attachment.html>