# Phenotypic mean = 0.000 (SD= 8.602) # QTL allele prop = 0.250 # QTL genot means = -5.888 0.000 5.888 # VQ = 13.000 ( 17.6%) # VA = 13.000 ( 17.6%) # VC = 0.000 ( 0.0%) # VE = 48.000 ( 64.9%) set locus mztwin qua set locus dBP qua set loc marker001 marker 10.00 set loc marker002 marker 15.00 set loc marker003 marker 20.00 set loc marker004 marker 25.00 set loc marker005 marker 30.00 set loc marker006 marker 35.00 set loc marker007 marker 40.00 set loc marker008 marker 45.00 set loc marker009 marker 50.00 set locus qtl marker 50.00 set loc marker010 marker 55.00 set loc marker011 marker 60.00 set loc marker012 marker 65.00 set loc marker013 marker 70.00 set loc marker014 marker 75.00 set loc marker015 marker 80.00 set loc marker016 marker 85.00 set loc marker017 marker 90.00 set loc marker018 marker 100.00 set loc marker019 marker 105.00 set loc marker020 marker 110.00 read pedigree fam300.ped set ple -1 set imp -1 set che off run set loc age qua set loc env qua set loc mal aff mal=male env=0 if (mal and rand<0.43) then env=1 if (isnon) then age=30+30*rand if (isfou) then age=60+10*rand # set loc qfac qua # factor qtl qfac # if (qfac==1) then trait=74 + sqrt(age*(0.2+mal*0.5)+48)*rnorm + env*5*rnorm # if (qfac==2) then trait=80 + sqrt(age*(0.2+mal*0.5)+48)*rnorm + env*5*rnorm # if (qfac==4) then trait=86 + sqrt(age*(0.2+mal*0.5)+48)*rnorm + env*5*rnorm dBP = 80 + env*6 + dBP + sqrt(age*(0.2+mal*0.5)+48)*rnorm + env*5*rnorm write sim1.ped write loc sib sim1.in sim1.ped $ sp2R -g sim1.in > sim1.rta order marker008 marker009 marker010 env dBP write mendel dbp.pem # write map mendel dbp.map # write loc mendel dbp.lom $ echo dBP > dbp.var $ cp one.in Control.in $ mendel und order marker008 marker009 marker010 dBP env write merlin.ped # write map merlin merlin.map # write loc merlin merlin.dat und set loc dBP2 qua if (env==1) then dBP2=dBP if (env==1) then dBP=x und order marker008 marker009 marker010 env dBP dBP2 write mendel dbp_two.pem $ echo dBP > dbp_two.var $ echo dBP2 >> dbp_two.var $ cp two.in Control.in $ mendel $ merlin --ibd $ merlin --vc $ cat dbp_one.sum