Hi dinhdiep,
I am very interested in your great work and want to use your method to generate MHL for liver cancer. When running your script I found there are some missing in the examples. For example, in this command:
./scripts/bam2cghap.sh allcpg/cpg.small.txt.gz BAMfiles/Colon_primary_tumor_sept9_promoter.bam test
I can not find a allcpg directory in your repository. With out the file, I can't figure out the format and generate my own file that is compatible with the command.
Hi dinhdiep, I am very interested in your great work and want to use your method to generate MHL for liver cancer. When running your script I found there are some missing in the examples. For example, in this command:
./scripts/bam2cghap.sh allcpg/cpg.small.txt.gz BAMfiles/Colon_primary_tumor_sept9_promoter.bam test
I can not find a allcpg directory in your repository. With out the file, I can't figure out the format and generate my own file that is compatible with the command.
And in this command,
./scripts/cghap2mhbs.sh results/HaploInfo/chr22.sub.hapInfo.txt example/N37_10_tissue_pooled.autosomes.RD10_80up.genomecov.bed 0.3 chr22
The results/HaploInfo/chr22.sub.hapInfo.txt is missing. I assume that this file is generated from one step of your script.
Can you help me address this problem? Thanks!