load("/home/yduan/soft/bio/seq/RNA_seq/CPAT-1.2.2/dat/Fly_logitModel.RData")
test <- read.table(file="/home/yduan/work_file/lncRNAdb/methods_assement/prediction_results/random5k/refseq_mouse/cpat_fly.out.dat",sep="\t",col.names=c("ID","mRNA","ORF","Fickett","Hexamer"))
test$prob <- predict(mylogit,newdata=test,type="response")
attach(test)
output <- cbind("mRNA_size"=mRNA,"ORF_size"=ORF,"Fickett_score"=Fickett,"Hexamer_score"=Hexamer,"coding_prob"=test$prob)
write.table(output,file="/home/yduan/work_file/lncRNAdb/methods_assement/prediction_results/random5k/refseq_mouse/cpat_fly.out",quote=F,sep="\t",row.names=ID)