R/5adboosting.R

In 13308204545/Chord: Chord

##V2:设置了可选参数
DBboostTrain<-function(mattest,mfinal){
  require(adabag)
  str(mattest)
  mattest$true<-as.factor(mattest$true)
  if (ncol(mattest)>5) {
    name.method<-""
    for (i in colnames(mattest)[6:ncol(mattest)]) {
      name.method<-paste0(name.method,"+",i)
    }
    DBboost<- boosting(formula(paste0("true ~bcds_s+cxds_s+dbf_s+scran_s",name.method)),data = mattest,boos=TRUE, mfinal=mfinal)
  }else{
    DBboost<- boosting(true~bcds_s+cxds_s+dbf_s+scran_s,data = mattest,boos=TRUE, mfinal=mfinal)
  }
# DBboost<- boosting(true~bcds_s+cxds_s+dbf_s+scran_s,data = mattest,boos=TRUE, mfinal=mfinal)
  return(DBboost)
}
DBboostPre<-function(DBboost,seu=seu,sce=sce,mattest=mattrain,outname="out"){
  require(adabag)

  bcds_s<-sce$bcds_score
  cxds_s<-sce$cxds_score
  dbf_s<-seu@meta.data[,ncol(seu@meta.data)-1]
  scran_s<-seu@meta.data[,"scran"]
  mattest<-as.data.frame(cbind(bcds_s,cxds_s,dbf_s,scran_s))
  mattest$bcds_s<-as.numeric(mattest$bcds_s)
  mattest$cxds_s<-as.numeric(mattest$cxds_s)
  mattest$dbf_s<-as.numeric(mattest$dbf_s)
  mattest$scran_s<-as.numeric(mattest$scran_s)

  str(mattest)
  pre_DB <- predict(DBboost,newdata = mattest)
  mattest$chord<-pre_DB$prob[,1]-pre_DB$prob[,2]

  write.csv(mattest,file = paste0("finalScore.",outname,".csv"))

  pdf(paste0(outname,"hist.pdf"))
  p<-hist(x = mattest$chord ,breaks = 50, xlim = c(-1,1))
  dev.off()
  return(mattest)
}

DBboostPre2<-function(DBboost,seu=seu,sce=sce,mattest=mattrain,outname="out"){
  require(adabag)
  pre_DB <- predict(DBboost,newdata = mattest)
  mattest$chord<-pre_DB$prob[,1]-pre_DB$prob[,2]

  write.csv(mattest,file = paste0("finalScore.",outname,".csv"))

  pdf(paste0(outname,"hist.pdf"))
  p<-hist(x = mattest$chord ,breaks = 50, xlim = c(-1,1))
  dev.off()
  return(mattest)
}