diff --git a/README.md b/README.md index 3c31edf..315ce5e 100644 --- a/README.md +++ b/README.md @@ -9,11 +9,11 @@ If the user has a single genotype file, then user can make use of bcftools scrip from smcpp_pipeline import vcf2smc,estimate,plot 1) Function to convert vcf file to smc format - def vcf2smc(popfile,chromosome_info,gaps_bed,num_workers) + def vcf2smc(popfile,chromosome_info,gaps_bed,num_workers): 2) Function to estimate the demographic history of a single population - def estimate(popfile,mu = 1.25e-8,spline = "piecewise",start_time = 100, end_time=100000,knots = 8 + def estimate(popfile,mu = 1.25e-8,spline = "piecewise",start_time = 100, end_time=100000,knots = 8): 3) Function to generate plots from the model.json file obtained using estimate function - def plot(popfile,time_start=100,time_end=15000) + def plot(popfile,time_start=100,time_end=15000): popfile should be a tab separated file containing two columns