R/hcaPlot.R

In jasenfinch/OrbiFIEmisc: Miscellaneous functions for the analysis of FIE-HRMS data.

#' hcaPlot
#' @export

hcaPlot <- function(dn,dat.all,cls,Path,DF,HPC_mode=F){
  for (i in 1:length(dn)) {
    ldamod <- nlda(dat.all[[i]],cls)
    if (length(unique(cls))<3){
      bitmap(paste(Path,DF,paste(DF,"HCA",sep="_"),paste(paste(DF, paste("lda_hca",dn[i],"Tw1",as.character(round(ldamod$Tw[[1]]*10, digits=0)/10),sep="_"),sep="_"),".bmp",sep=""),sep="/"))
    } else {
      bitmap(paste(Path,DF,paste(DF,"HCA",sep="_"),paste(paste(DF, paste("lda_hca",dn[i],"Tw1",as.character(round(ldamod$Tw[[1]]*10, digits=0)/10),"Tw2",as.character(round(ldamod$Tw[[2]]*10, digits=0)/10),sep="_"),sep="_"),".bmp",sep=""),sep="/"))
    }
    hca.nlda(ldamod)
    paste("Tw1",as.character(round(ldamod$Tw[[1]]*10, digits=0)/10), sep="_")
    dev.off()
  }
}