refseq_human_5000.result: refseq_human_5000.svm refseq_human_5000.predict
	perl /home/yduan/usr/bin/bin/svmresult.pl refseq_human_5000.svm refseq_human_5000.predict > refseq_human_5000.result

refseq_human_5000.predict: refseq_human_5000.scaled /home/yduan/soft/bio/seq/RNA_seq/iSeeRNA-1.2.2/model/hg19.model
	/home/yduan/usr/bin/bin/svm-predict -b 1 refseq_human_5000.scaled /home/yduan/soft/bio/seq/RNA_seq/iSeeRNA-1.2.2/model/hg19.model refseq_human_5000.predict
refseq_human_5000.scaled: refseq_human_5000.svm /home/yduan/soft/bio/seq/RNA_seq/iSeeRNA-1.2.2/model/hg19.scale.param
	/home/yduan/usr/bin/bin/svm-scale -r /home/yduan/soft/bio/seq/RNA_seq/iSeeRNA-1.2.2/model/hg19.scale.param refseq_human_5000.svm > refseq_human_5000.scaled

refseq_human_5000.svm: refseq_human_5000.list refseq_human_5000.orf refseq_human_5000.feature all_refseq_human_5000.consv
	perl /home/yduan/usr/bin/bin/merge.pl refseq_human_5000.list all_refseq_human_5000.consv refseq_human_5000.orf refseq_human_5000.feature > refseq_human_5000.svm

refseq_human_5000.gtf: /home/yduan/work_file/lncRNAdb/methods_assement/prediction_results/random5k/refseq_human/wd/refseq_human_5000.gtf
	/home/yduan/usr/bin/bin/gffread -o /dev/stdout /home/yduan/work_file/lncRNAdb/methods_assement/prediction_results/random5k/refseq_human/wd/refseq_human_5000.gtf | perl -ne 'print unless /^#/' | perl -ne '$$_="chr$$_" unless /^chr/; $$_=~s/Parent=/transcript_id /; print ' > refseq_human_5000.gtf

refseq_human_5000.fa: refseq_human_5000.gtf
	/home/yduan/usr/bin/bin/gffread -g /home/yduan/soft/bio/seq/RNA_seq/iSeeRNA-1.2.2/genome/hg38 -w refseq_human_5000.fa refseq_human_5000.gtf

refseq_human_5000.list: refseq_human_5000.fa
	perl /home/yduan/usr/bin/bin/fa2list.pl refseq_human_5000.fa refseq_human_5000.list

refseq_human_5000.feature: refseq_human_5000.fa
	perl /home/yduan/usr/bin/bin/fa2feature.pl refseq_human_5000.fa refseq_human_5000.feature

refseq_human_5000.orf: refseq_human_5000.fa
	/home/yduan/usr/bin/bin/txCdsPredict -anyStart refseq_human_5000.fa refseq_human_5000.orf

consv/Makefile: refseq_human_5000.gtf
	cd consv && perl /home/yduan/usr/bin/bin/prepare_consv.pl ../refseq_human_5000.gtf /home/yduan/soft/bio/seq/RNA_seq/iSeeRNA-1.2.2/consv/hg38 /home/yduan/usr/bin/bin/calcConsv all_refseq_human_5000.consv > Makefile && cd ..

all_refseq_human_5000.consv: consv/Makefile
	cd consv && make -j 1 && cd ..

clean:
	rm -rf consv/* *refseq_human_5000.*