R/dynRB_Vn.R

In dynRB: Dynamic Range Boxes

dynRB_Vn <-
function(A=A,steps=201,correlogram=FALSE, row_col = c(2, 2)){
  steps0<-steps
  names(A)[1]<-"Species"
  dims<-ncol(A)-1
  
  if(dims==1){
    
    insects<-levels(factor(A$Species))
    S2<-A[2:(dims+1)]
    insects<-as.character(levels(factor(A$Species)))
    G<-data.frame(V1=insects)
    for(i in 1:dims){
      G<-cbind(G,data.frame(x=0))
      names(G)[i+1]<-names(A)[i+1]
    }
    plot_data_volume<-list()
    names_data_S1<-c()
    for(i in 1:nrow(G)){
      S1<-subset(A,A$Species==G$V1[i])[,2:(dims+1)]
      
      EE<-.volumeA_coordinates_full_dim1(S1,S2,steps=steps0 )
      G[i,2:(dims+1)]<-EE$integral_coord
      plot_data_volume[[i]]<-EE$plot_volume  
      names_data_S1<-c(names_data_S1,names(EE$integral_coord))
      print(i)
    }
    plot_alpha_grid<-EE$alpha_grid
    print(G)
    r<-list(plot_data_volume=plot_data_volume, names_data_S1=names_data_S1, plot_alpha_grid=plot_alpha_grid, result=G)
    invisible(r) 
    
  }else{
    
    insects<-levels(factor(A$Species))
    if(correlogram==TRUE){
      par(mfrow=row_col)
      G<-expand.grid(insects)
      h<-1
      k<-1
      for(i in 1:length(insects)){
        S1<-subset(A,A$Species==G$Var1[i])[,2:(dims+1)]
        M<-cor(S1,use="pairwise.complete.obs")
        corrplot(M, method="color",   
                 type="upper", order="hclust", addCoef.col = "black", 
                 tl.col="black", tl.srt=45,
                 sig.level = 0.01, insig = "blank", diag=FALSE, mai = c("",as.character(G$Var1[i])),mar = c(0,2,0,0))
        if(k/h==row_col[1]*row_col[2]){
          title("Correlogram for each species", outer = T)
          h<-h+1
          readline(prompt = "Press <Enter> to continue...")}
        k<-k+1
      }
      title("Correlogram for each species", outer = T)
    }
    S2<-A[2:(dims+1)]
    insects<-as.character(levels(factor(A$Species)))
    G<-data.frame(V1=insects)
    for(i in 1:dims){
      G<-cbind(G,data.frame(x=0))
      names(G)[i+1]<-names(A)[i+1]
    }
    plot_data_volume<-list()
    names_data_S1<-c()
    for(i in 1:nrow(G)){
      S1<-subset(A,A$Species==G$V1[i])[,2:(dims+1)]
      
      EE<-.volumeA_coordinates_full(S1,S2,steps=steps0 )
      G[i,2:(dims+1)]<-EE$integral_coord
      plot_data_volume[[i]]<-EE$plot_volume  
      names_data_S1<-c(names_data_S1,names(EE$integral_coord))
      print(i)
    }
    plot_alpha_grid<-EE$alpha_grid
    print(G)
    r<-list(plot_data_volume=plot_data_volume, names_data_S1=names_data_S1, plot_alpha_grid=plot_alpha_grid, result=G)
    invisible(r) 
  
  }
  
}