A Statistical Variance Components Framework for Mapping Imprinted Quantitative Trait Locus in Experimental Crosses

Table 4

The power, MLEs of the QTL position, and effect parameters estimated based on 100 simulation replicates for data simulated with Mendelian and imprinting models based on the 20×20 sampling design.

Simulation	Analysis	Position	𝜇1	𝜇2	𝜇3	ℎ^2𝑎			ℎ^2𝑔	𝜎^2𝑒
model	model	48 cm	10	8	6	ℎ^2𝑚	ℎ^2𝑓	ℎ^2𝑚𝑓	0.083	2	Power^1	Power^2
M	M	48.08	10.123	7.890	6.023	0.544			0.093	1.982	1.00	1.00
		(2.981)	(0.887)	(0.730)	(1.046)	(0.121)			(0.137)	(0.250)
	I	48.260	9.922	8.081	6.074	0.263	0.272	0.239	0.049	1.983	1.00	1.00
		(3.015)	(1.010)	(0.720)	(0.949)	(0.097)	(0.110)	(0.124)	(0.080)	(0.249)
I1	M	48.940	10.079	7.913	6.023	0.323			0.111	2.022	0.86	0.96
		(11.660)	(0.588)	(0.511)	(0.712)	(0.287)			(0.134)	(0.205)
	I	47.86	10.077	7.900	6.018	0.008	0.473	0.089	0.091	2.015	0.95	0.98
		(7.524)	(0.579)	(0.575)	(0.665)	(0.022)	(0.151)	(0.099)	(0.053)	(0.86)
I2	M	48.640	10.070	7.917	6.028	0.331			0.103	2.023	0.91	0.98
		(8.389)	(0.613)	(0.517)	(0.706)	(0.278)			(0.129)	(0.195)
	I	49.200	10.059	7.926	6.008	0.079	0.396	0.089	0.089	2.023	0.93	0.98
		(9.818)	(0.617)	(0.518)	(0.673)	(0.075)	(0.156)	(0.097)	(0.050)	(0.184)

M and I refer to Mendelian and imprinting models, respectively. Simulated parameters for model M:(ℎ²_𝑎=0.583); I₁∶(ℎ²_𝑚,ℎ²_𝑓,ℎ²_𝑚𝑓)=(0, 0.5, 0.083); I₂∶(ℎ²_𝑚,ℎ²_𝑓,ℎ²_𝑚𝑓)=( 0.083, 0.417, 0.083). Power¹ and Power² correspond to the power calculated using the empirical cutoff and the theoretical threshold, respectively. The numbers given in the parenthesis with normal and italic fonts correspond to the RMSEs and standard errors of the parameter estimates, respectively. See Table 2 for other explanations.