QMSim - CGIL - University of Guelph
QMSim - CGIL - University of Guelph
QMSim - CGIL - University of Guelph
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