|
|
|
|
|
|
|
|
|
|
|
|
|
|
message("Write log2 expression file\n") |
|
|
message("Write log2 expression file\n") |
|
|
|
|
|
|
|
|
write.table(logexpr_out,file =file = paste(out_dir,opt$project_code,"_geneexp_log2FPKM.txt",sep=""),sep="\t",row.names=F,quote=F) |
|
|
|
|
|
|
|
|
write.table(logexpr_out,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 |
|
|
|
|
|
|
|
|
write.table(expr,file = paste(out_dir,opt$project_code,"_geneexp_FPKM.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) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|