瀏覽代碼

上传文件至 '.'

master
yingyu 6 年之前
父節點
當前提交
c0505b1b4f
共有 1 個檔案被更改,包括 5 行新增5 行删除
  1. +5
    -5
      RNAseq_1_expr_ballgown.R

+ 5
- 5
RNAseq_1_expr_ballgown.R 查看文件

logexpr_out<-cbind(rownames(logexpr),round(logexpr,3)) logexpr_out<-cbind(rownames(logexpr),round(logexpr,3))
colnames(logexpr_out)[1]<-"Gene" colnames(logexpr_out)[1]<-"Gene"


message("output log2 expression file\n")
message("Write log2 expression file\n")


write.table(logexpr_out,file = paste(out_dir,opt$project_code,"_geneexp_log2fpkm_floor0p01_c",ncol(logexpr),"r",nrow(logexpr),"_",Sys.Date(),".txt",sep=""),sep="\t",row.names=F,quote=F)
write.table(logexpr_out,file =file = paste(out_dir,opt$project_code,"_geneexp_log2FPKM.txt",sep=""),sep="\t",row.names=F,quote=F)
}else{ }else{


#output expression file with fpkm #output expression file with fpkm
expr<-cbind(rownames(expr),round(expr,3)) expr<-cbind(rownames(expr),round(expr,3))
colnames(expr)[1]<-"Gene" colnames(expr)[1]<-"Gene"


message("output fpkm expression file\n")
message("Write fpkm expression file\n")


write.table(expr,file = paste(out_dir,opt$project_code,"_geneexp_fpkm_c",ncol(expr),"r",nrow(expr),"_",Sys.Date(),".txt",sep=""),sep="\t",row.names=F,quote=F)
write.table(expr,file = paste(out_dir,opt$project_code,"_geneexp_FPKM.txt",sep=""),sep="\t",row.names=F,quote=F)
} }


write.table(expr_fpkm_out,file = paste(out_dir,opt$project_code,"_geneexp_FPKM.txt",sep=""),sep="\t",row.names=F,quote=F)





Loading…
取消
儲存