I run SampleBranchLengths.sh in Relate software.
Command is as follow:
${PATH_TO_RELATE}/scripts/SampleBranchLengths/SampleBranchLengths.sh \ -i ${INPUT_DIR}/n_ENe_chr${CHR} \ -o n_ENe_SBL_chr${CHR} \ -m 1.25e-8 \ --coal ${INPUT_DIR}/n_ENe.coal \ --format b \ --num_samples 5
And got the timeb suffix files for every chromosome.
Then how can I get the single SNP likelihood by palm/lik.py ?
Where can I set the pos of the snp, and calculate the snp ?
I run SampleBranchLengths.sh in Relate software.
Command is as follow:
${PATH_TO_RELATE}/scripts/SampleBranchLengths/SampleBranchLengths.sh \ -i ${INPUT_DIR}/n_ENe_chr${CHR} \ -o n_ENe_SBL_chr${CHR} \ -m 1.25e-8 \ --coal ${INPUT_DIR}/n_ENe.coal \ --format b \ --num_samples 5And got the timeb suffix files for every chromosome.
Then how can I get the single SNP likelihood by palm/lik.py ?
Where can I set the pos of the snp, and calculate the snp ?