Hey!<div><br></div><div>So I would like to do a Biomod.RangeSize comparing an average of two GCM's to the current climate data.</div><div><br></div><div>Here is what I've done so far:</div><div><br></div><div><div>
> load("proj.cccm80a2a/Total_consensus_cccm80a2a")</div><div>> load("proj.had80a2a/Total_consensus_had80a2a")</div><div>> MyAverageA2A = rowMeans(cbind(Total_consensus_cccm80a2a[,1,"prob.mean.weighted"], Total_consensus_had80a2a[,1,"prob.mean.weighted"]))</div>
<div>> MyThresA2A = mean(consensus_cccm80a2a_results$Cflorida$thresholds["prob.mean.weighted", c(2:4, 6:8)])</div></div><div>> MyAverageBinA2A=BinaryTransformation(MyAverageA2A, MyThresA2A)</div><div><br>
</div><div>For the current climate data all that I could think to do was put it in a similar 'format' as MyAverageBinA2A.</div><div><br></div><div>Here is what I tried:</div><div><br></div><div><div>> load("proj.present/Total_consensus_present_Bin")</div>
<div>> Present<-Total_consensus_present_Bin[,1,"prob.mean.weighted"]</div><div>> Biomod.RangeSize(CurrentPred=Present, FutureProj=MyAverageBinA2A, SpChange="Cflorida.2080a2a")</div></div><div>
Error in `dimnames<-.data.frame`(`*tmp*`, value = list(NULL, c("Loss", : </div><div><div> invalid 'dimnames' given for data frame</div></div><div><br></div><div><br></div><div>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.</div>
<div><br></div><div>Thank you so much for all your help! It is greatly appreciated!</div><div>Ashley</div>