22.07.2013 Views

QMSim - CGIL - University of Guelph

QMSim - CGIL - University of Guelph

QMSim - CGIL - University of Guelph

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

INPUT PARAMETER FILE<br />

The program requires a parameter file, in which various parameters for the simulation should be<br />

specified. The parameter file consists <strong>of</strong> five main sections. The first part describes global<br />

parameters, the second part describes historical generations, the third part describes parameters for<br />

subpopulations and generations, the fourth part contains genome parameters and the fifth part is<br />

related to the output options. The order <strong>of</strong> commands within each section is not normally important.<br />

An example <strong>of</strong> parameter file is given below.<br />

/*******************************<br />

** Global parameters **<br />

*******************************/<br />

title = "Simulating two divergent lines - 50k SNP panel";<br />

nr = 1; //Number <strong>of</strong> replicates<br />

h2 = 0.2; //Overall heritability<br />

qtlh2 = 0.2; //QTL heritability<br />

phvar = 1.0; //Phenotypic variance<br />

/*******************************<br />

** Historical generations **<br />

*******************************/<br />

begin_hist; //Historical generations<br />

nhg = 200; //Number <strong>of</strong> historical generations<br />

sfhg = 420; //Size <strong>of</strong> the first historical generation<br />

sihg = 420; //Size <strong>of</strong> intermediate historical generation<br />

ihg = 200; //Intermediate historical generation number<br />

slhg = 420; //Size <strong>of</strong> the last historical generation<br />

nmlhg = 20; //Number <strong>of</strong> male in the last historical generation<br />

end_hist;<br />

/*******************************<br />

** Populations **<br />

*******************************/<br />

begin_pop = "Line 1";<br />

begin_founder;<br />

male [n = 20, pop = "hp"];<br />

female [n = 400, pop = "hp"];<br />

end_founder;<br />

ls = 1 2 [0.05]; //Litter size<br />

pmp = 0.5 /fix; //Proportion <strong>of</strong> male progeny (random, fix or fixwf)<br />

ng = 10; //Number <strong>of</strong> generations<br />

md = random; //Mating design<br />

sr = 0.4; //Replacement ratio for sires<br />

dr = 0.2; //Replacement ratio for dams<br />

sd = phenotypic /h; //Selection design<br />

cd = age; //Culling design<br />

begin_popoutput;<br />

ld /bin 10 /maf 0.1 /gen 0;<br />

data;<br />

genotype /snp_code;<br />

allele_freq /gen 10;<br />

end_popoutput;<br />

end_pop;<br />

begin_pop = "Line 2";<br />

Dairy Cattle Breeding and Genetics Committee Meeting, September 18, 2008. 4

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!