SAS Output

Obs	animal	sire	dam	afa	afb	sfa	sfb	dfa	dfb	mgsfa	mgsfb	mgdfa	mgdfb	sex	bw	ww
1	1	0	0	1.00	0.00	1	0	1.0	0.0	1	0	1	0	1	33	289
2	2	0	0	0.00	1.00	0	1	0.0	1.0	0	1	0	1	2	29	245
3	3	0	2	0.50	0.50	1	0	0.0	1.0	0	1	0	1	2	32	256
4	4	1	0	0.50	0.50	1	0	0.0	1.0	0	1	0	1	2	30	261
5	5	1	2	0.50	0.50	1	0	0.0	1.0	0	1	0	1	1	38	292
6	6	1	3	0.75	0.25	1	0	0.5	0.5	1	0	0	1	1	35	286

MABM_19M_DMMAM_ALL MODELS_Example_November 05,2014_a November 5, 2014

ANIMAL BREEDING NOTES

CHAPTER 19M DMMAM ALL MODELS

DIRECT MATERNAL MULTIBREED ANIMAL MODELS WITH:

1) UNEQUAL RESIDUAL VARIANCE, UNEQUAL ADDITIVE AND NONADDITIVE GENETIC COVARIANCE MATRICES

2) EQUAL RESIDUAL VARIANCE, UNEQUAL ADDITIVE AND NONADDITIVE GENETIC COVARIANCE MATRICES

3) EQUAL RESIDUAL VARIANCE, EQUAL ADDITIVE GENETIC COVARIANCE MATRICES, UNEQUAL NONADDITIVE GENETIC COVARIANCE MATRICES

4) EQUAL RESIDUAL VARIANCE, EQUAL ADDITIVE GENETIC COVARIANCE MATRICES, NO RANDOM NONADDITIVE GENETIC EFFECTS

Mauricio A. Elzo, University of Florida, maelzo@ufl.edu

Read input dataset (SAS file)

OBS=     1  animal=   1.0000  sire=        0  dam=        0  afa=   1.0000  afb=        0  sfa=   1.0000  sfb=        0  dfa=   1.0000  dfb=        0 
mgsfa=   1.0000  mgsfb=        0  mgdfa=   1.0000  mgdfb=        0  sex=   1.0000  bw=  33.0000  ww= 289.0000                                         
OBS=     2  animal=   2.0000  sire=        0  dam=        0  afa=        0  afb=   1.0000  sfa=        0  sfb=   1.0000  dfa=        0  dfb=   1.0000 
mgsfa=        0  mgsfb=   1.0000  mgdfa=        0  mgdfb=   1.0000  sex=   2.0000  bw=  29.0000  ww= 245.0000                                         
OBS=     3  animal=   3.0000  sire=        0  dam=   2.0000  afa=   0.5000  afb=   0.5000  sfa=   1.0000  sfb=        0  dfa=        0  dfb=   1.0000 
mgsfa=        0  mgsfb=   1.0000  mgdfa=        0  mgdfb=   1.0000  sex=   2.0000  bw=  32.0000  ww= 256.0000                                         
OBS=     4  animal=   4.0000  sire=   1.0000  dam=        0  afa=   0.5000  afb=   0.5000  sfa=   1.0000  sfb=        0  dfa=        0  dfb=   1.0000 
mgsfa=        0  mgsfb=   1.0000  mgdfa=        0  mgdfb=   1.0000  sex=   2.0000  bw=  30.0000  ww= 261.0000                                         
OBS=     5  animal=   5.0000  sire=   1.0000  dam=   2.0000  afa=   0.5000  afb=   0.5000  sfa=   1.0000  sfb=        0  dfa=        0  dfb=   1.0000 
mgsfa=        0  mgsfb=   1.0000  mgdfa=        0  mgdfb=   1.0000  sex=   1.0000  bw=  38.0000  ww= 292.0000                                         
OBS=     6  animal=   6.0000  sire=   1.0000  dam=   3.0000  afa=   0.7500  afb=   0.2500  sfa=   1.0000  sfb=        0  dfa=   0.5000  dfb=   0.5000 
mgsfa=   1.0000  mgsfb=        0  mgdfa=        0  mgdfb=   1.0000  sex=   1.0000  bw=  35.0000  ww= 286.0000

datmat = matrix of input data

datmat
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16
ROW1	1	0	0	1	0	1	0	1	0	1	0	1	0	1	33	289
ROW2	2	0	0	0	1	0	1	0	1	0	1	0	1	2	29	245
ROW3	3	0	2	0.5	0.5	1	0	0	1	0	1	0	1	2	32	256
ROW4	4	1	0	0.5	0.5	1	0	0	1	0	1	0	1	2	30	261
ROW5	5	1	2	0.5	0.5	1	0	0	1	0	1	0	1	1	38	292
ROW6	6	1	3	0.75	0.25	1	0	0.5	0.5	1	0	0	1	1	35	286

Enter Parameters for Current Run

Enter restronsol = 1 to impose restrictions on solutions to solve the MME, else = 0 if not

restronsol
0

No restrictions imposed on solutions to solve MME

Enter nrec = Number of records

nrec
6

Make nanim = nrec for this example

nanim
6

Enter nf = Number of fixed effects in the MME

nf
9

Enter nga = Number of random additive genetic effects in the MME

nga
12

ngn = Number of random nonadditive genetic effects in the MME

ngn
12

Compute neq = nf+nga+ngn = total number of MME

neq
33

Define pedigf = pedigree file with breed composition of animals, sires, and dams

pedigf
1	0	0	1	0	1	0	1	0	1	0	1	0
2	0	0	0	1	0	1	0	1	0	1	0	1
3	0	2	0.5	0.5	1	0	0	1	0	1	0	1
4	1	0	0.5	0.5	1	0	0	1	0	1	0	1
5	1	2	0.5	0.5	1	0	0	1	0	1	0	1
6	1	3	0.75	0.25	1	0	0.5	0.5	1	0	0	1

Construct xf = matrix of fixed and random effects

Construct fixed effects in matrix xf

xf
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	1	1	1	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW2	1	0	0	1	1	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW3	1	0.5	0	0.5	1	1	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW4	1	0.5	0	0.5	1	1	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW5	1	0.5	0	0.5	1	1	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW6	1	0.75	0.5	0.25	0.5	0.5	1	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Construct random additive genetic effects in matrix xf

xf
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	1	1	1	0	0	0	0	1	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW2	1	0	0	1	1	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW3	1	0.5	0	0.5	1	1	0	0	1	0	0	0	1	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW4	1	0.5	0	0.5	1	1	0	0	1	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW5	1	0.5	0	0.5	1	1	0	1	0	0	0	0	1	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW6	1	0.75	0.5	0.25	0.5	0.5	1	1	0	0	0	0	0	0	1	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0

Construct random nonadditive genetic effects in matrix xf

xf
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	1	1	1	0	0	0	0	1	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW2	1	0	0	1	1	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW3	1	0.5	0	0.5	1	1	0	0	1	0	0	0	1	1	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0
ROW4	1	0.5	0	0.5	1	1	0	0	1	0	0	0	0	0	0	1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0
ROW5	1	0.5	0	0.5	1	1	0	1	0	0	0	0	1	0	0	0	0	1	0	0	0	1	0	1	0	0	0	0	0	0	0	0	0
ROW6	1	0.75	0.5	0.25	0.5	0.5	1	1	0	0	0	0	0	0	1	0	0	0	0	1	0	0.5	0	0	0	0.5	1	0	0	0	0	0	0

Make x = xf, i.e., use computed xf

x
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	1	1	1	0	0	0	0	1	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW2	1	0	0	1	1	0	0	0	1	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW3	1	0.5	0	0.5	1	1	0	0	1	0	0	0	1	1	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0
ROW4	1	0.5	0	0.5	1	1	0	0	1	0	0	0	0	0	0	1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0
ROW5	1	0.5	0	0.5	1	1	0	1	0	0	0	0	1	0	0	0	0	1	0	0	0	1	0	1	0	0	0	0	0	0	0	0	0
ROW6	1	0.75	0.5	0.25	0.5	0.5	1	1	0	0	0	0	0	0	1	0	0	0	0	1	0	0.5	0	0	0	0.5	1	0	0	0	0	0	0

Enter intrabreed and interbreed environmental variances

veaa	vebb	veab
49	16	25

Compute vef = vector of multibreed environmental variances

pedigf
1	0	0	1	0	1	0	1	0	1	0	1	0
2	0	0	0	1	0	1	0	1	0	1	0	1
3	0	2	0.5	0.5	1	0	0	1	0	1	0	1
4	1	0	0.5	0.5	1	0	0	1	0	1	0	1
5	1	2	0.5	0.5	1	0	0	1	0	1	0	1
6	1	3	0.75	0.25	1	0	0.5	0.5	1	0	0	1

vef
49
16
32.5
32.5
32.5
47

Make ve = vef, i.e., use computed ve

r = matrix of residual covariances

r
49	0	0	0	0	0
0	16	0	0	0	0
0	0	32.5	0	0	0
0	0	0	32.5	0	0
0	0	0	0	32.5	0
0	0	0	0	0	47

invr = inverse of matrix of residual covariances

invr
0.0204082	0	0	0	0	0
0	0.0625	0	0	0	0
0	0	0.0307692	0	0	0
0	0	0	0.0307692	0	0
0	0	0	0	0.0307692	0
0	0	0	0	0	0.0212766

Read yf = vector of records

yf
289
245
256
261
292
286

Make y = yf, i.e., use read yf

y
289
245
256
261
292
286

Compute xtinvr = x transpose times r

xtinvr
0.0204082	0.0625	0.0307692	0.0307692	0.0307692	0.0212766
0.0204082	0	0.0153846	0.0153846	0.0153846	0.0159574
0.0204082	0	0	0	0	0.0106383
0	0.0625	0.0153846	0.0153846	0.0153846	0.0053191
0	0.0625	0.0307692	0.0307692	0.0307692	0.0106383
0	0	0.0307692	0.0307692	0.0307692	0.0106383
0	0	0	0	0	0.0212766
0.0204082	0	0	0	0.0307692	0.0212766
0	0.0625	0.0307692	0.0307692	0	0
0.0204082	0	0	0	0	0
0	0	0	0	0	0
0	0.0625	0	0	0	0
0	0	0.0307692	0	0.0307692	0
0	0	0.0307692	0	0	0
0	0	0	0	0	0.0212766
0	0	0	0.0307692	0	0
0	0	0	0	0	0
0	0	0	0	0.0307692	0
0	0	0	0	0	0
0	0	0	0	0	0.0212766
0	0	0	0	0	0
0	0	0	0.0307692	0.0307692	0.0106383
0	0	0	0	0	0
0	0	0.0307692	0	0.0307692	0
0	0	0	0	0	0
0	0	0	0	0	0.0106383
0	0	0	0	0	0.0212766
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0
0	0	0	0	0	0

Compute xtinvr = x transpose times r times x

xtinvrx
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.1964925	0.0825195	0.0310465	0.113973	0.165446	0.102946	0.0212766	0.072454	0.1240385	0.0204082	0	0.0625	0.0615385	0.0307692	0.0212766	0.0307692	0	0.0307692	0	0.0212766	0	0.0721768	0	0.0615385	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW2	0.0825195	0.0554532	0.0283869	0.0270663	0.0541326	0.0541326	0.0159574	0.0517502	0.0307692	0.0204082	0	0	0.0307692	0.0153846	0.0159574	0.0153846	0	0.0153846	0	0.0159574	0	0.038748	0	0.0307692	0	0.0079787	0.0159574	0	0	0	0	0	0
ROW3	0.0310465	0.0283869	0.0257273	0.0026596	0.0053191	0.0053191	0.0106383	0.0310465	0	0.0204082	0	0	0	0	0.0106383	0	0	0	0	0.0106383	0	0.0053191	0	0	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW4	0.113973	0.0270663	0.0026596	0.0869067	0.1113134	0.0488134	0.0053191	0.0207038	0.0932692	0	0	0.0625	0.0307692	0.0153846	0.0053191	0.0153846	0	0.0153846	0	0.0053191	0	0.0334288	0	0.0307692	0	0.0026596	0.0053191	0	0	0	0	0	0
ROW5	0.165446	0.0541326	0.0053191	0.1113134	0.1601268	0.0976268	0.0106383	0.0414075	0.1240385	0	0	0.0625	0.0615385	0.0307692	0.0106383	0.0307692	0	0.0307692	0	0.0106383	0	0.0668576	0	0.0615385	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW6	0.102946	0.0541326	0.0053191	0.0488134	0.0976268	0.0976268	0.0106383	0.0414075	0.0615385	0	0	0	0.0615385	0.0307692	0.0106383	0.0307692	0	0.0307692	0	0.0106383	0	0.0668576	0	0.0615385	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW7	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	0	0	0	0	0	0.0212766	0	0	0	0	0.0212766	0	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW8	0.072454	0.0517502	0.0310465	0.0207038	0.0414075	0.0414075	0.0212766	0.072454	0	0.0204082	0	0	0.0307692	0	0.0212766	0	0	0.0307692	0	0.0212766	0	0.0414075	0	0.0307692	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW9	0.1240385	0.0307692	0	0.0932692	0.1240385	0.0615385	0	0	0.1240385	0	0	0.0625	0.0307692	0.0307692	0	0.0307692	0	0	0	0	0	0.0307692	0	0.0307692	0	0	0	0	0	0	0	0	0
ROW10	0.0204082	0.0204082	0.0204082	0	0	0	0	0.0204082	0	0.0204082	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW11	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW12	0.0625	0	0	0.0625	0.0625	0	0	0	0.0625	0	0	0.0625	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW13	0.0615385	0.0307692	0	0.0307692	0.0615385	0.0615385	0	0.0307692	0.0307692	0	0	0	0.0615385	0.0307692	0	0	0	0.0307692	0	0	0	0.0307692	0	0.0615385	0	0	0	0	0	0	0	0	0
ROW14	0.0307692	0.0153846	0	0.0153846	0.0307692	0.0307692	0	0	0.0307692	0	0	0	0.0307692	0.0307692	0	0	0	0	0	0	0	0	0	0.0307692	0	0	0	0	0	0	0	0	0
ROW15	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	0	0	0	0	0	0.0212766	0	0	0	0	0.0212766	0	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW16	0.0307692	0.0153846	0	0.0153846	0.0307692	0.0307692	0	0	0.0307692	0	0	0	0	0	0	0.0307692	0	0	0	0	0	0.0307692	0	0	0	0	0	0	0	0	0	0	0
ROW17	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW18	0.0307692	0.0153846	0	0.0153846	0.0307692	0.0307692	0	0.0307692	0	0	0	0	0.0307692	0	0	0	0	0.0307692	0	0	0	0.0307692	0	0.0307692	0	0	0	0	0	0	0	0	0
ROW19	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW20	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	0	0	0	0	0	0.0212766	0	0	0	0	0.0212766	0	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW21	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW22	0.0721768	0.038748	0.0053191	0.0334288	0.0668576	0.0668576	0.0106383	0.0414075	0.0307692	0	0	0	0.0307692	0	0.0106383	0.0307692	0	0.0307692	0	0.0106383	0	0.0668576	0	0.0307692	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW23	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW24	0.0615385	0.0307692	0	0.0307692	0.0615385	0.0615385	0	0.0307692	0.0307692	0	0	0	0.0615385	0.0307692	0	0	0	0.0307692	0	0	0	0.0307692	0	0.0615385	0	0	0	0	0	0	0	0	0
ROW25	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW26	0.0106383	0.0079787	0.0053191	0.0026596	0.0053191	0.0053191	0.0106383	0.0106383	0	0	0	0	0	0	0.0106383	0	0	0	0	0.0106383	0	0.0053191	0	0	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW27	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	0	0	0	0	0	0.0212766	0	0	0	0	0.0212766	0	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW28	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW29	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW30	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW31	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW32	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW33	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Enter intrabreed and interbreed additive genetic covariance matrices

vaaa		vabb		vaab
36	8	9	2	4	1
8	24	2	4	1	3

Compute vaf = multibreed additive genetic covariance matrices for individual animals

vaf
36	8	0	0	0	0	0	0	0	0	0	0
8	24	0	0	0	0	0	0	0	0	0	0
0	0	9	2	0	0	0	0	0	0	0	0
0	0	2	4	0	0	0	0	0	0	0	0
0	0	0	0	22.5	5	0	0	0	0	0	0
0	0	0	0	5	14	0	0	0	0	0	0
0	0	0	0	0	0	22.5	5	0	0	0	0
0	0	0	0	0	0	5	14	0	0	0	0
0	0	0	0	0	0	0	0	22.5	5	0	0
0	0	0	0	0	0	0	0	5	14	0	0
0	0	0	0	0	0	0	0	0	0	30.25	6.75
0	0	0	0	0	0	0	0	0	0	6.75	19.75

Compute daf = block-diagonal matrix of residual additive genetic covariance matrices

Recall: (Ga)-1 = (I - 1/2 P') (Block-diagonal Da)-1 (I - 1/2 P) for [dai]-1 blocks

daf
36	8	0	0	0	0	0	0	0	0	0	0
8	24	0	0	0	0	0	0	0	0	0	0
0	0	9	2	0	0	0	0	0	0	0	0
0	0	2	4	0	0	0	0	0	0	0	0
0	0	0	0	20.25	4.5	0	0	0	0	0	0
0	0	0	0	4.5	13	0	0	0	0	0	0
0	0	0	0	0	0	13.5	3	0	0	0	0
0	0	0	0	0	0	3	8	0	0	0	0
0	0	0	0	0	0	0	0	11.25	2.5	0	0
0	0	0	0	0	0	0	0	2.5	7	0	0
0	0	0	0	0	0	0	0	0	0	15.625	3.5
0	0	0	0	0	0	0	0	0	0	3.5	10.25

Make da = daf, i.e., use computed da

Compute dainv = inverse of da

dainv = inverse of block-diagonal matrix of residual additive genetic covariance matrices

dainv
0.03	-0.01	0	0	0	0	0	0	0	0	0	0
-0.01	0.045	0	0	0	0	0	0	0	0	0	0
0	0	0.125	-0.0625	0	0	0	0	0	0	0	0
0	0	-0.0625	0.28125	0	0	0	0	0	0	0	0
0	0	0	0	0.0534979	-0.018519	0	0	0	0	0	0
0	0	0	0	-0.018519	0.0833333	0	0	0	0	0	0
0	0	0	0	0	0	0.0808081	-0.030303	0	0	0	0
0	0	0	0	0	0	-0.030303	0.1363636	0	0	0	0
0	0	0	0	0	0	0	0	0.0965517	-0.034483	0	0
0	0	0	0	0	0	0	0	-0.034483	0.1551724	0	0
0	0	0	0	0	0	0	0	0	0	0.0693007	-0.023664
0	0	0	0	0	0	0	0	0	0	-0.023664	0.1056412

Compute gainv = inverse of the matrix of multibreed additive genetic covariances

Using algorithm to compute gainv directly; Elzo (1990a),JAS 68:1215-1228

gainv
0.0916651	-0.032112	0.0241379	-0.008621	0.0173252	-0.005916	-0.040404	0.0151515	-0.048276	0.0172414	-0.03465	0.0118318
-0.032112	0.1442943	-0.008621	0.0387931	-0.005916	0.0264103	0.0151515	-0.068182	0.0172414	-0.077586	0.0118318	-0.052821
0.0241379	-0.008621	0.1625124	-0.07575	-0.026749	0.0092593	0	0	-0.048276	0.0172414	0	0
-0.008621	0.0387931	-0.07575	0.3408764	0.0092593	-0.041667	0	0	0.0172414	-0.077586	0	0
0.0173252	-0.005916	-0.026749	0.0092593	0.0708231	-0.024434	0	0	0	0	-0.03465	0.0118318
-0.005916	0.0264103	0.0092593	-0.041667	-0.024434	0.1097436	0	0	0	0	0.0118318	-0.052821
-0.040404	0.0151515	0	0	0	0	0.0808081	-0.030303	0	0	0	0
0.0151515	-0.068182	0	0	0	0	-0.030303	0.1363636	0	0	0	0
-0.048276	0.0172414	-0.048276	0.0172414	0	0	0	0	0.0965517	-0.034483	0	0
0.0172414	-0.077586	0.0172414	-0.077586	0	0	0	0	-0.034483	0.1551724	0	0
-0.03465	0.0118318	0	0	-0.03465	0.0118318	0	0	0	0	0.0693007	-0.023664
0.0118318	-0.052821	0	0	0.0118318	-0.052821	0	0	0	0	-0.023664	0.1056412

gainv
0.092	-0.032	0.024	-0.009	0.017	-0.006	-0.040	0.015	-0.048	0.017	-0.035	0.012
-0.032	0.144	-0.009	0.039	-0.006	0.026	0.015	-0.068	0.017	-0.078	0.012	-0.053
0.024	-0.009	0.163	-0.076	-0.027	0.009	0.000	0.000	-0.048	0.017	0.000	0.000
-0.009	0.039	-0.076	0.341	0.009	-0.042	0.000	0.000	0.017	-0.078	0.000	0.000
0.017	-0.006	-0.027	0.009	0.071	-0.024	0.000	0.000	0.000	0.000	-0.035	0.012
-0.006	0.026	0.009	-0.042	-0.024	0.110	0.000	0.000	0.000	0.000	0.012	-0.053
-0.040	0.015	0.000	0.000	0.000	0.000	0.081	-0.030	0.000	0.000	0.000	0.000
0.015	-0.068	0.000	0.000	0.000	0.000	-0.030	0.136	0.000	0.000	0.000	0.000
-0.048	0.017	-0.048	0.017	0.000	0.000	0.000	0.000	0.097	-0.034	0.000	0.000
0.017	-0.078	0.017	-0.078	0.000	0.000	0.000	0.000	-0.034	0.155	0.000	0.000
-0.035	0.012	0.000	0.000	-0.035	0.012	0.000	0.000	0.000	0.000	0.069	-0.024
0.012	-0.053	0.000	0.000	0.012	-0.053	0.000	0.000	0.000	0.000	-0.024	0.106

Enter intrabreed and interbreed nonadditive genetic variances

vnab
16	6
6	12

Compute vnf = multibreed nonadditive genetic covariance matrices for individual animals

vnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	16	6	0	0	0	0	0	0
0	0	0	0	6	12	0	0	0	0	0	0
0	0	0	0	0	0	16	6	0	0	0	0
0	0	0	0	0	0	6	12	0	0	0	0
0	0	0	0	0	0	0	0	16	6	0	0
0	0	0	0	0	0	0	0	6	12	0	0
0	0	0	0	0	0	0	0	0	0	16	6
0	0	0	0	0	0	0	0	0	0	6	12

Compute dnf = block-diagonal matrix of residual nonadditive genetic covariance matrices

Recall: (Gn)-1 = (I - 1/2 P') (Block-diagonal Dn)-1 (I - 1/2 P) for [dni]-1 blocks

dnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	12	4.5	0	0	0	0	0	0
0	0	0	0	4.5	9	0	0	0	0	0	0
0	0	0	0	0	0	12	4.5	0	0	0	0
0	0	0	0	0	0	4.5	9	0	0	0	0
0	0	0	0	0	0	0	0	8	3	0	0
0	0	0	0	0	0	0	0	3	6	0	0
0	0	0	0	0	0	0	0	0	0	8	3
0	0	0	0	0	0	0	0	0	0	3	6

Make dn = dnf, i.e., use computed dn

Compute dainv = inverse of dn

dninv = inverse of block-diagonal matrix of residual nonadditive genetic covariance matrices

dninv
0.0769231	-0.038462	0	0	0	0	0	0	0	0	0	0
-0.038462	0.1025641	0	0	0	0	0	0	0	0	0	0
0	0	0.0769231	-0.038462	0	0	0	0	0	0	0	0
0	0	-0.038462	0.1025641	0	0	0	0	0	0	0	0
0	0	0	0	0.1025641	-0.051282	0	0	0	0	0	0
0	0	0	0	-0.051282	0.1367521	0	0	0	0	0	0
0	0	0	0	0	0	0.1025641	-0.051282	0	0	0	0
0	0	0	0	0	0	-0.051282	0.1367521	0	0	0	0
0	0	0	0	0	0	0	0	0.1538462	-0.076923	0	0
0	0	0	0	0	0	0	0	-0.076923	0.2051282	0	0
0	0	0	0	0	0	0	0	0	0	0.1538462	-0.076923
0	0	0	0	0	0	0	0	0	0	-0.076923	0.2051282

Compute gninv = inverse of the matrix of regression nonadditive genetic covariances

Using algorithm to compute gninv directly; Elzo (1990b),JAS 68:4079-4099

gninv
0.1794872	-0.089744	0.0384615	-0.019231	0.0384615	-0.019231	-0.051282	0.025641	-0.076923	0.0384615	-0.076923	0.0384615
-0.089744	0.2393162	-0.019231	0.0512821	-0.019231	0.0512821	0.025641	-0.068376	0.0384615	-0.102564	0.0384615	-0.102564
0.0384615	-0.019231	0.1410256	-0.070513	-0.051282	0.025641	0	0	-0.076923	0.0384615	0	0
-0.019231	0.0512821	-0.070513	0.1880342	0.025641	-0.068376	0	0	0.0384615	-0.102564	0	0
0.0384615	-0.019231	-0.051282	0.025641	0.1410256	-0.070513	0	0	0	0	-0.076923	0.0384615
-0.019231	0.0512821	0.025641	-0.068376	-0.070513	0.1880342	0	0	0	0	0.0384615	-0.102564
-0.051282	0.025641	0	0	0	0	0.1025641	-0.051282	0	0	0	0
0.025641	-0.068376	0	0	0	0	-0.051282	0.1367521	0	0	0	0
-0.076923	0.0384615	-0.076923	0.0384615	0	0	0	0	0.1538462	-0.076923	0	0
0.0384615	-0.102564	0.0384615	-0.102564	0	0	0	0	-0.076923	0.2051282	0	0
-0.076923	0.0384615	0	0	-0.076923	0.0384615	0	0	0	0	0.1538462	-0.076923
0.0384615	-0.102564	0	0	0.0384615	-0.102564	0	0	0	0	-0.076923	0.2051282

gninv
0.179	-0.090	0.038	-0.019	0.038	-0.019	-0.051	0.026	-0.077	0.038	-0.077	0.038
-0.090	0.239	-0.019	0.051	-0.019	0.051	0.026	-0.068	0.038	-0.103	0.038	-0.103
0.038	-0.019	0.141	-0.071	-0.051	0.026	0.000	0.000	-0.077	0.038	0.000	0.000
-0.019	0.051	-0.071	0.188	0.026	-0.068	0.000	0.000	0.038	-0.103	0.000	0.000
0.038	-0.019	-0.051	0.026	0.141	-0.071	0.000	0.000	0.000	0.000	-0.077	0.038
-0.019	0.051	0.026	-0.068	-0.071	0.188	0.000	0.000	0.000	0.000	0.038	-0.103
-0.051	0.026	0.000	0.000	0.000	0.000	0.103	-0.051	0.000	0.000	0.000	0.000
0.026	-0.068	0.000	0.000	0.000	0.000	-0.051	0.137	0.000	0.000	0.000	0.000
-0.077	0.038	-0.077	0.038	0.000	0.000	0.000	0.000	0.154	-0.077	0.000	0.000
0.038	-0.103	0.038	-0.103	0.000	0.000	0.000	0.000	-0.077	0.205	0.000	0.000
-0.077	0.038	0.000	0.000	-0.077	0.038	0.000	0.000	0.000	0.000	0.154	-0.077
0.038	-0.103	0.000	0.000	0.038	-0.103	0.000	0.000	0.000	0.000	-0.077	0.205

Compute lhs = left hand side of the MME

Add gainv to lhs

Add gninv to lhs

lhs
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.1964925	0.0825195	0.0310465	0.113973	0.165446	0.102946	0.0212766	0.072454	0.1240385	0.0204082	0	0.0625	0.0615385	0.0307692	0.0212766	0.0307692	0	0.0307692	0	0.0212766	0	0.0721768	0	0.0615385	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW2	0.0825195	0.0554532	0.0283869	0.0270663	0.0541326	0.0541326	0.0159574	0.0517502	0.0307692	0.0204082	0	0	0.0307692	0.0153846	0.0159574	0.0153846	0	0.0153846	0	0.0159574	0	0.038748	0	0.0307692	0	0.0079787	0.0159574	0	0	0	0	0	0
ROW3	0.0310465	0.0283869	0.0257273	0.0026596	0.0053191	0.0053191	0.0106383	0.0310465	0	0.0204082	0	0	0	0	0.0106383	0	0	0	0	0.0106383	0	0.0053191	0	0	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW4	0.113973	0.0270663	0.0026596	0.0869067	0.1113134	0.0488134	0.0053191	0.0207038	0.0932692	0	0	0.0625	0.0307692	0.0153846	0.0053191	0.0153846	0	0.0153846	0	0.0053191	0	0.0334288	0	0.0307692	0	0.0026596	0.0053191	0	0	0	0	0	0
ROW5	0.165446	0.0541326	0.0053191	0.1113134	0.1601268	0.0976268	0.0106383	0.0414075	0.1240385	0	0	0.0625	0.0615385	0.0307692	0.0106383	0.0307692	0	0.0307692	0	0.0106383	0	0.0668576	0	0.0615385	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW6	0.102946	0.0541326	0.0053191	0.0488134	0.0976268	0.0976268	0.0106383	0.0414075	0.0615385	0	0	0	0.0615385	0.0307692	0.0106383	0.0307692	0	0.0307692	0	0.0106383	0	0.0668576	0	0.0615385	0	0.0053191	0.0106383	0	0	0	0	0	0
ROW7	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	0	0	0	0	0	0.0212766	0	0	0	0	0.0212766	0	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW8	0.072454	0.0517502	0.0310465	0.0207038	0.0414075	0.0414075	0.0212766	0.072454	0	0.0204082	0	0	0.0307692	0	0.0212766	0	0	0.0307692	0	0.0212766	0	0.0414075	0	0.0307692	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW9	0.1240385	0.0307692	0	0.0932692	0.1240385	0.0615385	0	0	0.1240385	0	0	0.0625	0.0307692	0.0307692	0	0.0307692	0	0	0	0	0	0.0307692	0	0.0307692	0	0	0	0	0	0	0	0	0
ROW10	0.0204082	0.0204082	0.0204082	0	0	0	0	0.0204082	0	0.1120733	-0.032112	0.0241379	-0.008621	0.0173252	-0.005916	-0.040404	0.0151515	-0.048276	0.0172414	-0.03465	0.0118318	0	0	0	0	0	0	0	0	0	0	0	0
ROW11	0	0	0	0	0	0	0	0	0	-0.032112	0.1442943	-0.008621	0.0387931	-0.005916	0.0264103	0.0151515	-0.068182	0.0172414	-0.077586	0.0118318	-0.052821	0	0	0	0	0	0	0	0	0	0	0	0
ROW12	0.0625	0	0	0.0625	0.0625	0	0	0	0.0625	0.0241379	-0.008621	0.2250124	-0.07575	-0.026749	0.0092593	0	0	-0.048276	0.0172414	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW13	0.0615385	0.0307692	0	0.0307692	0.0615385	0.0615385	0	0.0307692	0.0307692	-0.008621	0.0387931	-0.07575	0.4024149	0.0400285	-0.041667	0	0	0.0480106	-0.077586	0	0	0.0307692	0	0.0615385	0	0	0	0	0	0	0	0	0
ROW14	0.0307692	0.0153846	0	0.0153846	0.0307692	0.0307692	0	0	0.0307692	0.0173252	-0.005916	-0.026749	0.0400285	0.1015923	-0.024434	0	0	0	0	-0.03465	0.0118318	0	0	0.0307692	0	0	0	0	0	0	0	0	0
ROW15	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	-0.005916	0.0264103	0.0092593	-0.041667	-0.024434	0.1310202	0	0	0	0	0.0331084	-0.052821	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW16	0.0307692	0.0153846	0	0.0153846	0.0307692	0.0307692	0	0	0.0307692	-0.040404	0.0151515	0	0	0	0	0.1115773	-0.030303	0	0	0	0	0.0307692	0	0	0	0	0	0	0	0	0	0	0
ROW17	0	0	0	0	0	0	0	0	0	0.0151515	-0.068182	0	0	0	0	-0.030303	0.1363636	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW18	0.0307692	0.0153846	0	0.0153846	0.0307692	0.0307692	0	0.0307692	0	-0.048276	0.0172414	-0.048276	0.0480106	0	0	0	0	0.127321	-0.034483	0	0	0.0307692	0	0.0307692	0	0	0	0	0	0	0	0	0
ROW19	0	0	0	0	0	0	0	0	0	0.0172414	-0.077586	0.0172414	-0.077586	0	0	0	0	-0.034483	0.1551724	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW20	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	-0.03465	0.0118318	0	0	-0.03465	0.0331084	0	0	0	0	0.0905773	-0.023664	0.0106383	0	0	0	0.0106383	0.0212766	0	0	0	0	0	0
ROW21	0	0	0	0	0	0	0	0	0	0.0118318	-0.052821	0	0	0.0118318	-0.052821	0	0	0	0	-0.023664	0.1056412	0	0	0	0	0	0	0	0	0	0	0	0
ROW22	0.0721768	0.038748	0.0053191	0.0334288	0.0668576	0.0668576	0.0106383	0.0414075	0.0307692	0	0	0	0.0307692	0	0.0106383	0.0307692	0	0.0307692	0	0.0106383	0	0.2463448	-0.089744	0.0692308	-0.019231	0.0437807	-0.008592	-0.051282	0.025641	-0.076923	0.0384615	-0.076923	0.0384615
ROW23	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	-0.089744	0.2393162	-0.019231	0.0512821	-0.019231	0.0512821	0.025641	-0.068376	0.0384615	-0.102564	0.0384615	-0.102564
ROW24	0.0615385	0.0307692	0	0.0307692	0.0615385	0.0615385	0	0.0307692	0.0307692	0	0	0	0.0615385	0.0307692	0	0	0	0.0307692	0	0	0	0.0692308	-0.019231	0.2025641	-0.070513	-0.051282	0.025641	0	0	-0.076923	0.0384615	0	0
ROW25	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	-0.019231	0.0512821	-0.070513	0.1880342	0.025641	-0.068376	0	0	0.0384615	-0.102564	0	0
ROW26	0.0106383	0.0079787	0.0053191	0.0026596	0.0053191	0.0053191	0.0106383	0.0106383	0	0	0	0	0	0	0.0106383	0	0	0	0	0.0106383	0	0.0437807	-0.019231	-0.051282	0.025641	0.1463448	-0.059875	0	0	0	0	-0.076923	0.0384615
ROW27	0.0212766	0.0159574	0.0106383	0.0053191	0.0106383	0.0106383	0.0212766	0.0212766	0	0	0	0	0	0	0.0212766	0	0	0	0	0.0212766	0	-0.008592	0.0512821	0.025641	-0.068376	-0.059875	0.2093108	0	0	0	0	0.0384615	-0.102564
ROW28	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	-0.051282	0.025641	0	0	0	0	0.1025641	-0.051282	0	0	0	0
ROW29	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.025641	-0.068376	0	0	0	0	-0.051282	0.1367521	0	0	0	0
ROW30	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	-0.076923	0.0384615	-0.076923	0.0384615	0	0	0	0	0.1538462	-0.076923	0	0
ROW31	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.0384615	-0.102564	0.0384615	-0.102564	0	0	0	0	-0.076923	0.2051282	0	0
ROW32	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	-0.076923	0.0384615	0	0	-0.076923	0.0384615	0	0	0	0	0.1538462	-0.076923
ROW33	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.0384615	-0.102564	0	0	0.0384615	-0.102564	0	0	0	0	-0.076923	0.2051282

lhs
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.196	0.083	0.031	0.114	0.165	0.103	0.021	0.072	0.124	0.020	0.000	0.063	0.062	0.031	0.021	0.031	0.000	0.031	0.000	0.021	0.000	0.072	0.000	0.062	0.000	0.011	0.021	0.000	0.000	0.000	0.000	0.000	0.000
ROW2	0.083	0.055	0.028	0.027	0.054	0.054	0.016	0.052	0.031	0.020	0.000	0.000	0.031	0.015	0.016	0.015	0.000	0.015	0.000	0.016	0.000	0.039	0.000	0.031	0.000	0.008	0.016	0.000	0.000	0.000	0.000	0.000	0.000
ROW3	0.031	0.028	0.026	0.003	0.005	0.005	0.011	0.031	0.000	0.020	0.000	0.000	0.000	0.000	0.011	0.000	0.000	0.000	0.000	0.011	0.000	0.005	0.000	0.000	0.000	0.005	0.011	0.000	0.000	0.000	0.000	0.000	0.000
ROW4	0.114	0.027	0.003	0.087	0.111	0.049	0.005	0.021	0.093	0.000	0.000	0.063	0.031	0.015	0.005	0.015	0.000	0.015	0.000	0.005	0.000	0.033	0.000	0.031	0.000	0.003	0.005	0.000	0.000	0.000	0.000	0.000	0.000
ROW5	0.165	0.054	0.005	0.111	0.160	0.098	0.011	0.041	0.124	0.000	0.000	0.063	0.062	0.031	0.011	0.031	0.000	0.031	0.000	0.011	0.000	0.067	0.000	0.062	0.000	0.005	0.011	0.000	0.000	0.000	0.000	0.000	0.000
ROW6	0.103	0.054	0.005	0.049	0.098	0.098	0.011	0.041	0.062	0.000	0.000	0.000	0.062	0.031	0.011	0.031	0.000	0.031	0.000	0.011	0.000	0.067	0.000	0.062	0.000	0.005	0.011	0.000	0.000	0.000	0.000	0.000	0.000
ROW7	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	0.000	0.000	0.000	0.000	0.000	0.000	0.021	0.000	0.000	0.000	0.000	0.021	0.000	0.011	0.000	0.000	0.000	0.011	0.021	0.000	0.000	0.000	0.000	0.000	0.000
ROW8	0.072	0.052	0.031	0.021	0.041	0.041	0.021	0.072	0.000	0.020	0.000	0.000	0.031	0.000	0.021	0.000	0.000	0.031	0.000	0.021	0.000	0.041	0.000	0.031	0.000	0.011	0.021	0.000	0.000	0.000	0.000	0.000	0.000
ROW9	0.124	0.031	0.000	0.093	0.124	0.062	0.000	0.000	0.124	0.000	0.000	0.063	0.031	0.031	0.000	0.031	0.000	0.000	0.000	0.000	0.000	0.031	0.000	0.031	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW10	0.020	0.020	0.020	0.000	0.000	0.000	0.000	0.020	0.000	0.112	-0.032	0.024	-0.009	0.017	-0.006	-0.040	0.015	-0.048	0.017	-0.035	0.012	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW11	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.032	0.144	-0.009	0.039	-0.006	0.026	0.015	-0.068	0.017	-0.078	0.012	-0.053	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW12	0.063	0.000	0.000	0.063	0.063	0.000	0.000	0.000	0.063	0.024	-0.009	0.225	-0.076	-0.027	0.009	0.000	0.000	-0.048	0.017	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW13	0.062	0.031	0.000	0.031	0.062	0.062	0.000	0.031	0.031	-0.009	0.039	-0.076	0.402	0.040	-0.042	0.000	0.000	0.048	-0.078	0.000	0.000	0.031	0.000	0.062	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW14	0.031	0.015	0.000	0.015	0.031	0.031	0.000	0.000	0.031	0.017	-0.006	-0.027	0.040	0.102	-0.024	0.000	0.000	0.000	0.000	-0.035	0.012	0.000	0.000	0.031	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW15	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	0.000	-0.006	0.026	0.009	-0.042	-0.024	0.131	0.000	0.000	0.000	0.000	0.033	-0.053	0.011	0.000	0.000	0.000	0.011	0.021	0.000	0.000	0.000	0.000	0.000	0.000
ROW16	0.031	0.015	0.000	0.015	0.031	0.031	0.000	0.000	0.031	-0.040	0.015	0.000	0.000	0.000	0.000	0.112	-0.030	0.000	0.000	0.000	0.000	0.031	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW17	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.015	-0.068	0.000	0.000	0.000	0.000	-0.030	0.136	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW18	0.031	0.015	0.000	0.015	0.031	0.031	0.000	0.031	0.000	-0.048	0.017	-0.048	0.048	0.000	0.000	0.000	0.000	0.127	-0.034	0.000	0.000	0.031	0.000	0.031	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW19	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.017	-0.078	0.017	-0.078	0.000	0.000	0.000	0.000	-0.034	0.155	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW20	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	0.000	-0.035	0.012	0.000	0.000	-0.035	0.033	0.000	0.000	0.000	0.000	0.091	-0.024	0.011	0.000	0.000	0.000	0.011	0.021	0.000	0.000	0.000	0.000	0.000	0.000
ROW21	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.012	-0.053	0.000	0.000	0.012	-0.053	0.000	0.000	0.000	0.000	-0.024	0.106	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW22	0.072	0.039	0.005	0.033	0.067	0.067	0.011	0.041	0.031	0.000	0.000	0.000	0.031	0.000	0.011	0.031	0.000	0.031	0.000	0.011	0.000	0.246	-0.090	0.069	-0.019	0.044	-0.009	-0.051	0.026	-0.077	0.038	-0.077	0.038
ROW23	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.090	0.239	-0.019	0.051	-0.019	0.051	0.026	-0.068	0.038	-0.103	0.038	-0.103
ROW24	0.062	0.031	0.000	0.031	0.062	0.062	0.000	0.031	0.031	0.000	0.000	0.000	0.062	0.031	0.000	0.000	0.000	0.031	0.000	0.000	0.000	0.069	-0.019	0.203	-0.071	-0.051	0.026	0.000	0.000	-0.077	0.038	0.000	0.000
ROW25	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.019	0.051	-0.071	0.188	0.026	-0.068	0.000	0.000	0.038	-0.103	0.000	0.000
ROW26	0.011	0.008	0.005	0.003	0.005	0.005	0.011	0.011	0.000	0.000	0.000	0.000	0.000	0.000	0.011	0.000	0.000	0.000	0.000	0.011	0.000	0.044	-0.019	-0.051	0.026	0.146	-0.060	0.000	0.000	0.000	0.000	-0.077	0.038
ROW27	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	0.000	0.000	0.000	0.000	0.000	0.000	0.021	0.000	0.000	0.000	0.000	0.021	0.000	-0.009	0.051	0.026	-0.068	-0.060	0.209	0.000	0.000	0.000	0.000	0.038	-0.103
ROW28	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.051	0.026	0.000	0.000	0.000	0.000	0.103	-0.051	0.000	0.000	0.000	0.000
ROW29	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.026	-0.068	0.000	0.000	0.000	0.000	-0.051	0.137	0.000	0.000	0.000	0.000
ROW30	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.077	0.038	-0.077	0.038	0.000	0.000	0.000	0.000	0.154	-0.077	0.000	0.000
ROW31	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.038	-0.103	0.038	-0.103	0.000	0.000	0.000	0.000	-0.077	0.205	0.000	0.000
ROW32	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.077	0.038	0.000	0.000	-0.077	0.038	0.000	0.000	0.000	0.000	0.154	-0.077
ROW33	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.038	-0.103	0.000	0.000	0.038	-0.103	0.000	0.000	0.000	0.000	-0.077	0.205

Compute rhs = right hand side of the MME

rhs
52.187873
22.907943
8.9405124
29.27993
43.247361
27.934861
6.0851064
20.967681
31.220192
5.8979592
0
15.3125
16.861538
7.8769231
6.0851064
8.0307692
0
8.9846154
0
6.0851064
0
20.057938
0
16.861538
0
3.0425532
6.0851064
0
0
0
0
0
0

rhs
52.19
22.91
8.94
29.28
43.25
27.93
6.09
20.97
31.22
5.90
0.00
15.31
16.86
7.88
6.09
8.03
0.00
8.98
0.00
6.09
0.00
20.06
0.00
16.86
0.00
3.04
6.09
0.00
0.00
0.00
0.00
0.00
0.00

Compute ginvlhs = generalized inverse of the left hand side of the MME

ginvlhs
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	4.364257	3.506757	5.471267	0.857500	-1.107010	-3.929020	-6.418328	3.377990	0.986267	-6.920270	-1.537838	-1.901014	-0.477703	-1.265203	-0.308784	-3.250338	-0.722297	-4.410642	-1.007770	-4.092736	-0.923311	0.248649	0.093243	-0.248649	-0.093243	-0.124324	-0.046622	0.124324	0.046622	-0.000000	-0.000000	0.062162	0.023311
ROW2	3.506757	15.356757	12.679392	-11.85000	-9.172635	5.354730	-2.015203	-6.772635	10.279392	-10.07027	-2.237838	2.436486	0.397297	-3.652703	-0.883784	-7.187838	-1.597297	-1.566892	-0.420270	-6.861486	-1.560811	0.648649	0.243243	-0.648649	-0.243243	-0.324324	-0.121622	0.324324	0.121622	-0.000000	-0.000000	0.162162	0.060811
ROW3	5.471267	12.679392	16.961448	-7.208125	-11.49018	-8.564113	-3.913788	-6.423932	11.895198	-9.088176	-2.019595	1.331841	1.024493	-1.259882	0.084291	-4.610220	-1.024493	-0.503167	0.252449	-5.174029	-0.967652	4.721622	1.770608	3.278378	1.229392	1.639189	0.614696	2.360811	0.885304	4.000000	1.500000	3.180405	1.192652
ROW4	0.857500	-11.85000	-7.208125	12.707500	8.065625	-9.283750	-4.403125	10.150625	-9.293125	3.150000	0.700000	-4.337500	-0.875000	2.387500	0.575000	3.937500	0.875000	-2.843750	-0.587500	2.768750	0.637500	-0.400000	-0.150000	0.400000	0.150000	0.200000	0.075000	-0.200000	-0.075000	0.000000	-0.000000	-0.100000	-0.037500
ROW5	-1.107010	-9.172635	-11.49018	8.065625	10.383171	4.635093	-2.504540	9.801921	-10.90893	2.167905	0.481757	-3.232855	-1.502196	-0.005321	-0.393074	1.359882	0.302196	-3.907475	-1.260220	1.081292	0.044341	-4.472973	-1.677365	-3.527027	-1.322635	-1.763514	-0.661318	-2.236486	-0.838682	-4.000000	-1.500000	-3.118243	-1.169341
ROW6	-3.929020	5.354730	-8.564113	-9.283750	4.635093	27.837686	3.797171	-0.697407	-3.231613	-1.964189	-0.436486	2.209291	-1.254392	-4.785642	-1.936149	-5.155236	-1.145608	-2.127449	-1.345439	-3.374916	-1.186318	-8.145946	-3.054730	-7.854054	-2.945270	-3.927027	-1.472635	-4.072973	-1.527365	-8.000000	-3.000000	-6.036486	-2.263682
ROW7	-6.418328	-2.015203	-3.913788	-4.403125	-2.504540	3.797171	127.40530	-18.94829	12.529962	6.795608	1.510135	0.796030	0.612331	-11.74704	-14.39274	1.744510	0.387669	9.420819	2.311233	-18.10072	-9.941301	-1.959459	-0.734797	2.959459	-3.765203	-9.020270	-13.13260	-0.979730	-0.367399	0.500000	-2.250000	-5.489865	-6.933699
ROW8	3.377990	-6.772635	-6.423932	10.150625	9.801921	-0.697407	-18.94829	22.450671	-19.07268	-7.732095	-1.718243	-2.657855	-1.552196	3.319679	0.256926	-0.215118	-0.047804	-10.81997	-2.885220	-2.206208	-0.730659	-3.672973	-1.377365	-4.327027	-1.622635	-2.163514	-0.811318	-1.836486	-0.688682	-4.000000	-1.500000	-2.918243	-1.094341
ROW9	0.986267	10.279392	11.895198	-9.293125	-10.90893	-3.231613	12.529962	-19.07268	20.058948	0.811824	0.180405	0.756841	1.074493	-4.584882	-0.565709	-3.035220	-0.674493	6.409333	1.877449	-1.886529	-0.192652	3.921622	1.470608	4.078378	1.529392	2.039189	0.764696	1.960811	0.735304	4.000000	1.500000	2.980405	1.117652
ROW10	-6.920270	-10.07027	-9.088176	3.150000	2.167905	-1.964189	6.795608	-7.732095	0.811824	33.810811	7.513514	0.790541	0.608108	2.858108	0.851351	15.263514	3.391892	17.300676	4.060811	18.334459	4.182432	-1.945946	-0.729730	1.945946	0.729730	0.972973	0.364865	-0.972973	-0.364865	0.000000	0.000000	-0.486486	-0.182432
ROW11	-1.537838	-2.237838	-2.019595	0.700000	0.481757	-0.436486	1.510135	-1.718243	0.180405	7.513514	23.891892	0.175676	0.135135	0.635135	0.189189	3.391892	11.864865	3.844595	12.013514	4.074324	12.040541	-0.432432	-0.162162	0.432432	0.162162	0.216216	0.081081	-0.216216	-0.081081	0.000000	0.000000	-0.108108	-0.040541
ROW12	-1.901014	2.436486	1.331841	-4.337500	-3.232855	2.209291	0.796030	-2.657855	0.756841	0.790541	0.175676	8.714527	1.780405	3.467905	0.692568	0.988176	0.219595	4.752534	0.978041	2.129223	0.434122	0.702703	0.263514	-0.702703	-0.263514	-0.351351	-0.131757	0.351351	0.131757	0.000000	0.000000	0.175676	0.065878
ROW13	-0.477703	0.397297	1.024493	-0.875000	-1.502196	-1.254392	0.612331	-1.552196	1.074493	0.608108	0.135135	1.780405	3.831081	0.206081	1.763514	0.760135	0.168919	1.194257	1.983108	0.407095	0.949324	0.540541	0.202703	-0.540541	-0.202703	-0.270270	-0.101351	0.270270	0.101351	0.000000	-0.000000	0.135135	0.050676
ROW14	-1.265203	-3.652703	-1.259882	2.387500	-0.005321	-4.785642	-11.74704	3.319679	-4.584882	2.858108	0.635135	3.467905	0.206081	18.768581	3.888514	3.572635	0.793919	3.163007	0.420608	10.813345	2.261824	2.540541	0.952703	-2.540541	-0.952703	-1.270270	-0.476351	1.270270	0.476351	0.000000	0.000000	0.635135	0.238176
ROW15	-0.308784	-0.883784	0.084291	0.575000	-0.393074	-1.936149	-14.39274	0.256926	-0.565709	0.851351	0.189189	0.692568	1.763514	3.888514	13.668919	1.064189	0.236486	0.771959	0.976351	2.369932	6.929054	0.756757	0.283784	-0.756757	-0.283784	-0.378378	-0.141892	0.378378	0.141892	0.000000	-0.000000	0.189189	0.070946
ROW16	-3.250338	-7.187838	-4.610220	3.937500	1.359882	-5.155236	1.744510	-0.215118	-3.035220	15.263514	3.391892	0.988176	0.760135	3.572635	1.064189	19.079392	4.239865	8.125845	2.076014	9.418074	2.228041	-2.432432	-0.912162	2.432432	0.912162	1.216216	0.456081	-1.216216	-0.456081	-0.000000	-0.000000	-0.608108	-0.228041
ROW17	-0.722297	-1.597297	-1.024493	0.875000	0.302196	-1.145608	0.387669	-0.047804	-0.674493	3.391892	11.864865	0.219595	0.168919	0.793919	0.236486	4.239865	13.831081	1.805743	6.016892	2.092905	6.050676	-0.540541	-0.202703	0.540541	0.202703	0.270270	0.101351	-0.270270	-0.101351	-0.000000	0.000000	-0.135135	-0.050676
ROW18	-4.410642	-1.566892	-0.503167	-2.843750	-3.907475	-2.127449	9.420819	-10.81997	6.409333	17.300676	3.844595	4.752534	1.194257	3.163007	0.771959	8.125845	1.805743	22.276605	5.019426	10.231841	2.308277	-0.621622	-0.233108	0.621622	0.233108	0.310811	0.116554	-0.310811	-0.116554	0.000000	-0.000000	-0.155405	-0.058277
ROW19	-1.007770	-0.420270	0.252449	-0.587500	-1.260220	-1.345439	2.311233	-2.885220	1.877449	4.060811	12.013514	0.978041	1.983108	0.420608	0.976351	2.076014	6.016892	5.019426	13.998311	2.240709	6.494932	0.054054	0.020270	-0.054054	-0.020270	-0.027027	-0.010135	0.027027	0.010135	-0.000000	0.000000	0.013514	0.005068
ROW20	-4.092736	-6.861486	-5.174029	2.768750	1.081292	-3.374916	-18.10072	-2.206208	-1.886529	18.334459	4.074324	2.129223	0.407095	10.813345	2.369932	9.418074	2.092905	10.231841	2.240709	30.198902	6.722128	0.297297	0.111486	-0.297297	-0.111486	-0.148649	-0.055743	0.148649	0.055743	0.000000	0.000000	0.074324	0.027872
ROW21	-0.923311	-1.560811	-0.967652	0.637500	0.044341	-1.186318	-9.941301	-0.730659	-0.192652	4.182432	12.040541	0.434122	0.949324	2.261824	6.929054	2.228041	6.050676	2.308277	6.494932	6.722128	19.734797	0.162162	0.060811	-0.162162	-0.060811	-0.081081	-0.030405	0.081081	0.030405	0.000000	0.000000	0.040541	0.015203
ROW22	0.248649	0.648649	4.721622	-0.400000	-4.472973	-8.145946	-1.959459	-3.672973	3.921622	-1.945946	-0.432432	0.702703	0.540541	2.540541	0.756757	-2.432432	-0.540541	-0.621622	0.054054	0.297297	0.162162	14.270270	5.351351	1.729730	0.648649	0.864865	0.324324	7.135135	2.675676	8.000000	3.000000	7.567568	2.837838
ROW23	0.093243	0.243243	1.770608	-0.150000	-1.677365	-3.054730	-0.734797	-1.377365	1.470608	-0.729730	-0.162162	0.263514	0.202703	0.952703	0.283784	-0.912162	-0.202703	-0.233108	0.020270	0.111486	0.060811	5.351351	11.756757	0.648649	0.243243	0.324324	0.121622	2.675676	5.878378	3.000000	6.000000	2.837838	5.939189
ROW24	-0.248649	-0.648649	3.278378	0.400000	-3.527027	-7.854054	2.959459	-4.327027	4.078378	1.945946	0.432432	-0.702703	-0.540541	-2.540541	-0.756757	2.432432	0.540541	0.621622	-0.054054	-0.297297	-0.162162	1.729730	0.648649	14.270270	5.351351	7.135135	2.675676	0.864865	0.324324	8.000000	3.000000	4.432432	1.662162
ROW25	-0.093243	-0.243243	1.229392	0.150000	-1.322635	-2.945270	-3.765203	-1.622635	1.529392	0.729730	0.162162	-0.263514	-0.202703	-0.952703	-0.283784	0.912162	0.202703	0.233108	-0.020270	-0.111486	-0.060811	0.648649	0.243243	5.351351	11.756757	2.675676	5.878378	0.324324	0.121622	3.000000	6.000000	1.662162	3.060811
ROW26	-0.124324	-0.324324	1.639189	0.200000	-1.763514	-3.927027	-9.020270	-2.163514	2.039189	0.972973	0.216216	-0.351351	-0.270270	-1.270270	-0.378378	1.216216	0.270270	0.310811	-0.027027	-0.148649	-0.081081	0.864865	0.324324	7.135135	2.675676	15.567568	5.837838	0.432432	0.162162	4.000000	1.500000	8.216216	3.081081
ROW27	-0.046622	-0.121622	0.614696	0.075000	-0.661318	-1.472635	-13.13260	-0.811318	0.764696	0.364865	0.081081	-0.131757	-0.101351	-0.476351	-0.141892	0.456081	0.101351	0.116554	-0.010135	-0.055743	-0.030405	0.324324	0.121622	2.675676	5.878378	5.837838	11.939189	0.162162	0.060811	1.500000	3.000000	3.081081	6.030405
ROW28	0.124324	0.324324	2.360811	-0.200000	-2.236486	-4.072973	-0.979730	-1.836486	1.960811	-0.972973	-0.216216	0.351351	0.270270	1.270270	0.378378	-1.216216	-0.270270	-0.310811	0.027027	0.148649	0.081081	7.135135	2.675676	0.864865	0.324324	0.432432	0.162162	15.567568	5.837838	4.000000	1.500000	3.783784	1.418919
ROW29	0.046622	0.121622	0.885304	-0.075000	-0.838682	-1.527365	-0.367399	-0.688682	0.735304	-0.364865	-0.081081	0.131757	0.101351	0.476351	0.141892	-0.456081	-0.101351	-0.116554	0.010135	0.055743	0.030405	2.675676	5.878378	0.324324	0.121622	0.162162	0.060811	5.837838	11.939189	1.500000	3.000000	1.418919	2.969595
ROW30	0.000000	0.000000	4.000000	-0.000000	-4.000000	-8.000000	0.500000	-4.000000	4.000000	0.000000	0.000000	-0.000000	-0.000000	-0.000000	-0.000000	0.000000	0.000000	-0.000000	0.000000	-0.000000	-0.000000	8.000000	3.000000	8.000000	3.000000	4.000000	1.500000	4.000000	1.500000	16.000000	6.000000	6.000000	2.250000
ROW31	-0.000000	-0.000000	1.500000	-0.000000	-1.500000	-3.000000	-2.250000	-1.500000	1.500000	0.000000	0.000000	-0.000000	-0.000000	-0.000000	-0.000000	0.000000	0.000000	-0.000000	0.000000	-0.000000	-0.000000	3.000000	6.000000	3.000000	6.000000	1.500000	3.000000	1.500000	3.000000	6.000000	12.000000	2.250000	4.500000
ROW32	0.062162	0.162162	3.180405	-0.100000	-3.118243	-6.036486	-5.489865	-2.918243	2.980405	-0.486486	-0.108108	0.175676	0.135135	0.635135	0.189189	-0.608108	-0.135135	-0.155405	0.013514	0.074324	0.040541	7.567568	2.837838	4.432432	1.662162	8.216216	3.081081	3.783784	1.418919	6.000000	2.250000	15.891892	5.959459
ROW33	0.023311	0.060811	1.192652	-0.037500	-1.169341	-2.263682	-6.933699	-1.094341	1.117652	-0.182432	-0.040541	0.065878	0.050676	0.238176	0.070946	-0.228041	-0.050676	-0.058277	0.005068	0.027872	0.015203	2.837838	5.939189	1.662162	3.060811	3.081081	6.030405	1.418919	2.969595	2.250000	4.500000	5.959459	11.984797

ginvlhs
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	4.364	3.507	5.471	0.857	-1.107	-3.929	-6.418	3.378	0.986	-6.920	-1.538	-1.901	-0.478	-1.265	-0.309	-3.250	-0.722	-4.411	-1.008	-4.093	-0.923	0.249	0.093	-0.249	-0.093	-0.124	-0.047	0.124	0.047	-0.000	-0.000	0.062	0.023
ROW2	3.507	15.357	12.679	-11.85	-9.173	5.355	-2.015	-6.773	10.279	-10.07	-2.238	2.436	0.397	-3.653	-0.884	-7.188	-1.597	-1.567	-0.420	-6.861	-1.561	0.649	0.243	-0.649	-0.243	-0.324	-0.122	0.324	0.122	-0.000	-0.000	0.162	0.061
ROW3	5.471	12.679	16.961	-7.208	-11.49	-8.564	-3.914	-6.424	11.895	-9.088	-2.020	1.332	1.024	-1.260	0.084	-4.610	-1.024	-0.503	0.252	-5.174	-0.968	4.722	1.771	3.278	1.229	1.639	0.615	2.361	0.885	4.000	1.500	3.180	1.193
ROW4	0.858	-11.85	-7.208	12.707	8.066	-9.284	-4.403	10.151	-9.293	3.150	0.700	-4.337	-0.875	2.388	0.575	3.937	0.875	-2.844	-0.587	2.769	0.638	-0.400	-0.150	0.400	0.150	0.200	0.075	-0.200	-0.075	0.000	-0.000	-0.100	-0.038
ROW5	-1.107	-9.173	-11.49	8.066	10.383	4.635	-2.505	9.802	-10.91	2.168	0.482	-3.233	-1.502	-0.005	-0.393	1.360	0.302	-3.907	-1.260	1.081	0.044	-4.473	-1.677	-3.527	-1.323	-1.764	-0.661	-2.236	-0.839	-4.000	-1.500	-3.118	-1.169
ROW6	-3.929	5.355	-8.564	-9.284	4.635	27.838	3.797	-0.697	-3.232	-1.964	-0.436	2.209	-1.254	-4.786	-1.936	-5.155	-1.146	-2.127	-1.345	-3.375	-1.186	-8.146	-3.055	-7.854	-2.945	-3.927	-1.473	-4.073	-1.527	-8.000	-3.000	-6.036	-2.264
ROW7	-6.418	-2.015	-3.914	-4.403	-2.505	3.797	127.41	-18.95	12.530	6.796	1.510	0.796	0.612	-11.75	-14.39	1.745	0.388	9.421	2.311	-18.10	-9.941	-1.959	-0.735	2.959	-3.765	-9.020	-13.13	-0.980	-0.367	0.500	-2.250	-5.490	-6.934
ROW8	3.378	-6.773	-6.424	10.151	9.802	-0.697	-18.95	22.451	-19.07	-7.732	-1.718	-2.658	-1.552	3.320	0.257	-0.215	-0.048	-10.82	-2.885	-2.206	-0.731	-3.673	-1.377	-4.327	-1.623	-2.164	-0.811	-1.836	-0.689	-4.000	-1.500	-2.918	-1.094
ROW9	0.986	10.279	11.895	-9.293	-10.91	-3.232	12.530	-19.07	20.059	0.812	0.180	0.757	1.074	-4.585	-0.566	-3.035	-0.674	6.409	1.877	-1.887	-0.193	3.922	1.471	4.078	1.529	2.039	0.765	1.961	0.735	4.000	1.500	2.980	1.118
ROW10	-6.920	-10.07	-9.088	3.150	2.168	-1.964	6.796	-7.732	0.812	33.811	7.514	0.791	0.608	2.858	0.851	15.264	3.392	17.301	4.061	18.334	4.182	-1.946	-0.730	1.946	0.730	0.973	0.365	-0.973	-0.365	0.000	0.000	-0.486	-0.182
ROW11	-1.538	-2.238	-2.020	0.700	0.482	-0.436	1.510	-1.718	0.180	7.514	23.892	0.176	0.135	0.635	0.189	3.392	11.865	3.845	12.014	4.074	12.041	-0.432	-0.162	0.432	0.162	0.216	0.081	-0.216	-0.081	0.000	0.000	-0.108	-0.041
ROW12	-1.901	2.436	1.332	-4.337	-3.233	2.209	0.796	-2.658	0.757	0.791	0.176	8.715	1.780	3.468	0.693	0.988	0.220	4.753	0.978	2.129	0.434	0.703	0.264	-0.703	-0.264	-0.351	-0.132	0.351	0.132	0.000	0.000	0.176	0.066
ROW13	-0.478	0.397	1.024	-0.875	-1.502	-1.254	0.612	-1.552	1.074	0.608	0.135	1.780	3.831	0.206	1.764	0.760	0.169	1.194	1.983	0.407	0.949	0.541	0.203	-0.541	-0.203	-0.270	-0.101	0.270	0.101	0.000	-0.000	0.135	0.051
ROW14	-1.265	-3.653	-1.260	2.387	-0.005	-4.786	-11.75	3.320	-4.585	2.858	0.635	3.468	0.206	18.769	3.889	3.573	0.794	3.163	0.421	10.813	2.262	2.541	0.953	-2.541	-0.953	-1.270	-0.476	1.270	0.476	0.000	0.000	0.635	0.238
ROW15	-0.309	-0.884	0.084	0.575	-0.393	-1.936	-14.39	0.257	-0.566	0.851	0.189	0.693	1.764	3.889	13.669	1.064	0.236	0.772	0.976	2.370	6.929	0.757	0.284	-0.757	-0.284	-0.378	-0.142	0.378	0.142	0.000	-0.000	0.189	0.071
ROW16	-3.250	-7.188	-4.610	3.938	1.360	-5.155	1.745	-0.215	-3.035	15.264	3.392	0.988	0.760	3.573	1.064	19.079	4.240	8.126	2.076	9.418	2.228	-2.432	-0.912	2.432	0.912	1.216	0.456	-1.216	-0.456	-0.000	-0.000	-0.608	-0.228
ROW17	-0.722	-1.597	-1.024	0.875	0.302	-1.146	0.388	-0.048	-0.674	3.392	11.865	0.220	0.169	0.794	0.236	4.240	13.831	1.806	6.017	2.093	6.051	-0.541	-0.203	0.541	0.203	0.270	0.101	-0.270	-0.101	-0.000	0.000	-0.135	-0.051
ROW18	-4.411	-1.567	-0.503	-2.844	-3.907	-2.127	9.421	-10.82	6.409	17.301	3.845	4.753	1.194	3.163	0.772	8.126	1.806	22.277	5.019	10.232	2.308	-0.622	-0.233	0.622	0.233	0.311	0.117	-0.311	-0.117	0.000	-0.000	-0.155	-0.058
ROW19	-1.008	-0.420	0.252	-0.587	-1.260	-1.345	2.311	-2.885	1.877	4.061	12.014	0.978	1.983	0.421	0.976	2.076	6.017	5.019	13.998	2.241	6.495	0.054	0.020	-0.054	-0.020	-0.027	-0.010	0.027	0.010	-0.000	0.000	0.014	0.005
ROW20	-4.093	-6.861	-5.174	2.769	1.081	-3.375	-18.10	-2.206	-1.887	18.334	4.074	2.129	0.407	10.813	2.370	9.418	2.093	10.232	2.241	30.199	6.722	0.297	0.111	-0.297	-0.111	-0.149	-0.056	0.149	0.056	0.000	0.000	0.074	0.028
ROW21	-0.923	-1.561	-0.968	0.637	0.044	-1.186	-9.941	-0.731	-0.193	4.182	12.041	0.434	0.949	2.262	6.929	2.228	6.051	2.308	6.495	6.722	19.735	0.162	0.061	-0.162	-0.061	-0.081	-0.030	0.081	0.030	0.000	0.000	0.041	0.015
ROW22	0.249	0.649	4.722	-0.400	-4.473	-8.146	-1.959	-3.673	3.922	-1.946	-0.432	0.703	0.541	2.541	0.757	-2.432	-0.541	-0.622	0.054	0.297	0.162	14.270	5.351	1.730	0.649	0.865	0.324	7.135	2.676	8.000	3.000	7.568	2.838
ROW23	0.093	0.243	1.771	-0.150	-1.677	-3.055	-0.735	-1.377	1.471	-0.730	-0.162	0.264	0.203	0.953	0.284	-0.912	-0.203	-0.233	0.020	0.111	0.061	5.351	11.757	0.649	0.243	0.324	0.122	2.676	5.878	3.000	6.000	2.838	5.939
ROW24	-0.249	-0.649	3.278	0.400	-3.527	-7.854	2.959	-4.327	4.078	1.946	0.432	-0.703	-0.541	-2.541	-0.757	2.432	0.541	0.622	-0.054	-0.297	-0.162	1.730	0.649	14.270	5.351	7.135	2.676	0.865	0.324	8.000	3.000	4.432	1.662
ROW25	-0.093	-0.243	1.229	0.150	-1.323	-2.945	-3.765	-1.623	1.529	0.730	0.162	-0.264	-0.203	-0.953	-0.284	0.912	0.203	0.233	-0.020	-0.111	-0.061	0.649	0.243	5.351	11.757	2.676	5.878	0.324	0.122	3.000	6.000	1.662	3.061
ROW26	-0.124	-0.324	1.639	0.200	-1.764	-3.927	-9.020	-2.164	2.039	0.973	0.216	-0.351	-0.270	-1.270	-0.378	1.216	0.270	0.311	-0.027	-0.149	-0.081	0.865	0.324	7.135	2.676	15.568	5.838	0.432	0.162	4.000	1.500	8.216	3.081
ROW27	-0.047	-0.122	0.615	0.075	-0.661	-1.473	-13.13	-0.811	0.765	0.365	0.081	-0.132	-0.101	-0.476	-0.142	0.456	0.101	0.117	-0.010	-0.056	-0.030	0.324	0.122	2.676	5.878	5.838	11.939	0.162	0.061	1.500	3.000	3.081	6.030
ROW28	0.124	0.324	2.361	-0.200	-2.236	-4.073	-0.980	-1.836	1.961	-0.973	-0.216	0.351	0.270	1.270	0.378	-1.216	-0.270	-0.311	0.027	0.149	0.081	7.135	2.676	0.865	0.324	0.432	0.162	15.568	5.838	4.000	1.500	3.784	1.419
ROW29	0.047	0.122	0.885	-0.075	-0.839	-1.527	-0.367	-0.689	0.735	-0.365	-0.081	0.132	0.101	0.476	0.142	-0.456	-0.101	-0.117	0.010	0.056	0.030	2.676	5.878	0.324	0.122	0.162	0.061	5.838	11.939	1.500	3.000	1.419	2.970
ROW30	0.000	0.000	4.000	-0.000	-4.000	-8.000	0.500	-4.000	4.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	8.000	3.000	8.000	3.000	4.000	1.500	4.000	1.500	16.000	6.000	6.000	2.250
ROW31	-0.000	-0.000	1.500	-0.000	-1.500	-3.000	-2.250	-1.500	1.500	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	3.000	6.000	3.000	6.000	1.500	3.000	1.500	3.000	6.000	12.000	2.250	4.500
ROW32	0.062	0.162	3.180	-0.100	-3.118	-6.036	-5.490	-2.918	2.980	-0.486	-0.108	0.176	0.135	0.635	0.189	-0.608	-0.135	-0.155	0.014	0.074	0.041	7.568	2.838	4.432	1.662	8.216	3.081	3.784	1.419	6.000	2.250	15.892	5.959
ROW33	0.023	0.061	1.193	-0.038	-1.169	-2.264	-6.934	-1.094	1.118	-0.182	-0.041	0.066	0.051	0.238	0.071	-0.228	-0.051	-0.058	0.005	0.028	0.015	2.838	5.939	1.662	3.061	3.081	6.030	1.419	2.970	2.250	4.500	5.959	11.985

Compute gl = ginvlhs*lhs = matrix of expectations of solutions

gl
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.400	0.200	0.200	0.200	0.200	-0.000	0.000	0.200	0.200	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW2	0.200	0.400	0.300	-0.200	-0.100	0.200	-0.000	0.100	0.100	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW3	0.200	0.300	0.400	-0.100	-0.200	-0.200	0.000	0.100	0.100	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW4	0.200	-0.200	-0.100	0.400	0.300	-0.200	0.000	0.100	0.100	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000
ROW5	0.200	-0.100	-0.200	0.300	0.400	0.200	-0.000	0.100	0.100	0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW6	-0.000	0.200	-0.200	-0.200	0.200	0.800	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW7	0.000	-0.000	0.000	0.000	0.000	-0.000	1.000	0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000
ROW8	0.200	0.100	0.100	0.100	0.100	0.000	-0.000	0.600	-0.400	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000
ROW9	0.200	0.100	0.100	0.100	0.100	-0.000	-0.000	-0.400	0.600	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW10	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	0.000	-0.000	1.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000
ROW11	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	1.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW12	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	1.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW13	0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	0.000	1.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000
ROW14	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	1.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW15	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	1.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000
ROW16	0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	1.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000
ROW17	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	1.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000
ROW18	0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	0.000	1.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000
ROW19	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	1.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000
ROW20	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	1.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW21	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	1.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000
ROW22	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	1.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000
ROW23	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	1.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW24	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	1.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW25	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	1.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW26	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	1.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000
ROW27	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	1.000	-0.000	-0.000	0.000	0.000	0.000	-0.000
ROW28	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	1.000	-0.000	0.000	0.000	0.000	0.000
ROW29	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	1.000	0.000	0.000	0.000	0.000
ROW30	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	1.000	0.000	-0.000	0.000
ROW31	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	1.000	0.000	0.000
ROW32	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	1.000	0.000
ROW33	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	1.000

Notice that lg = gl (i.e., lhs*ginvlhs = lhs*ginvlhs)

lg
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.400	0.200	0.200	0.200	0.200	-0.000	0.000	0.200	0.200	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW2	0.200	0.400	0.300	-0.200	-0.100	0.200	0.000	0.100	0.100	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW3	0.200	0.300	0.400	-0.100	-0.200	-0.200	0.000	0.100	0.100	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW4	0.200	-0.200	-0.100	0.400	0.300	-0.200	0.000	0.100	0.100	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW5	0.200	-0.100	-0.200	0.300	0.400	0.200	0.000	0.100	0.100	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW6	-0.000	0.200	-0.200	-0.200	0.200	0.800	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW7	-0.000	-0.000	0.000	0.000	-0.000	0.000	1.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW8	0.200	0.100	0.100	0.100	0.100	0.000	0.000	0.600	-0.400	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW9	0.200	0.100	0.100	0.100	0.100	-0.000	0.000	-0.400	0.600	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW10	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	1.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000
ROW11	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	0.000	1.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW12	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	1.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW13	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	0.000	1.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000
ROW14	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	1.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW15	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	1.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000
ROW16	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	1.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW17	0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	1.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000
ROW18	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	1.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000
ROW19	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	1.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000
ROW20	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	1.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW21	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	1.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	0.000
ROW22	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	1.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW23	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	1.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW24	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	0.000	1.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW25	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	1.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW26	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	1.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW27	0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	1.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000
ROW28	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	1.000	-0.000	0.000	0.000	0.000	0.000
ROW29	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	1.000	0.000	0.000	0.000	0.000
ROW30	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	1.000	0.000	0.000	0.000
ROW31	-0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	1.000	0.000	0.000
ROW32	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	1.000	0.000
ROW33	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	1.000

Verify that lgl = lhs (i.e., lhs*ginvlhs*lhs = lhs => generalized inverse is correct)

lgl
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.196	0.083	0.031	0.114	0.165	0.103	0.021	0.072	0.124	0.020	-0.000	0.063	0.062	0.031	0.021	0.031	-0.000	0.031	0.000	0.021	0.000	0.072	-0.000	0.062	-0.000	0.011	0.021	0.000	0.000	-0.000	0.000	-0.000	0.000
ROW2	0.083	0.055	0.028	0.027	0.054	0.054	0.016	0.052	0.031	0.020	-0.000	-0.000	0.031	0.015	0.016	0.015	0.000	0.015	0.000	0.016	0.000	0.039	-0.000	0.031	-0.000	0.008	0.016	0.000	0.000	-0.000	0.000	-0.000	0.000
ROW3	0.031	0.028	0.026	0.003	0.005	0.005	0.011	0.031	-0.000	0.020	-0.000	-0.000	0.000	0.000	0.011	-0.000	-0.000	0.000	-0.000	0.011	0.000	0.005	0.000	-0.000	-0.000	0.005	0.011	0.000	0.000	0.000	-0.000	-0.000	-0.000
ROW4	0.114	0.027	0.003	0.087	0.111	0.049	0.005	0.021	0.093	-0.000	-0.000	0.063	0.031	0.015	0.005	0.015	0.000	0.015	0.000	0.005	0.000	0.033	-0.000	0.031	-0.000	0.003	0.005	0.000	0.000	-0.000	0.000	-0.000	0.000
ROW5	0.165	0.054	0.005	0.111	0.160	0.098	0.011	0.041	0.124	0.000	-0.000	0.063	0.062	0.031	0.011	0.031	0.000	0.031	0.000	0.011	0.000	0.067	-0.000	0.062	0.000	0.005	0.011	-0.000	0.000	-0.000	0.000	-0.000	0.000
ROW6	0.103	0.054	0.005	0.049	0.098	0.098	0.011	0.041	0.062	0.000	0.000	-0.000	0.062	0.031	0.011	0.031	0.000	0.031	0.000	0.011	-0.000	0.067	-0.000	0.062	0.000	0.005	0.011	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW7	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.021	-0.000	0.000	0.000	0.000	0.021	0.000	0.011	-0.000	0.000	0.000	0.011	0.021	-0.000	0.000	0.000	-0.000	0.000	0.000
ROW8	0.072	0.052	0.031	0.021	0.041	0.041	0.021	0.072	-0.000	0.020	-0.000	-0.000	0.031	0.000	0.021	-0.000	0.000	0.031	0.000	0.021	0.000	0.041	-0.000	0.031	-0.000	0.011	0.021	0.000	0.000	-0.000	0.000	-0.000	0.000
ROW9	0.124	0.031	0.000	0.093	0.124	0.062	0.000	0.000	0.124	-0.000	0.000	0.063	0.031	0.031	0.000	0.031	-0.000	-0.000	-0.000	0.000	0.000	0.031	-0.000	0.031	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000
ROW10	0.020	0.020	0.020	-0.000	0.000	0.000	0.000	0.020	-0.000	0.112	-0.032	0.024	-0.009	0.017	-0.006	-0.040	0.015	-0.048	0.017	-0.035	0.012	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000
ROW11	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.032	0.144	-0.009	0.039	-0.006	0.026	0.015	-0.068	0.017	-0.078	0.012	-0.053	-0.000	-0.000	0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000
ROW12	0.063	-0.000	-0.000	0.063	0.063	0.000	-0.000	-0.000	0.063	0.024	-0.009	0.225	-0.076	-0.027	0.009	0.000	-0.000	-0.048	0.017	-0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	0.000
ROW13	0.062	0.031	0.000	0.031	0.062	0.062	0.000	0.031	0.031	-0.009	0.039	-0.076	0.402	0.040	-0.042	-0.000	0.000	0.048	-0.078	0.000	-0.000	0.031	-0.000	0.062	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000
ROW14	0.031	0.015	-0.000	0.015	0.031	0.031	-0.000	-0.000	0.031	0.017	-0.006	-0.027	0.040	0.102	-0.024	0.000	-0.000	-0.000	-0.000	-0.035	0.012	-0.000	-0.000	0.031	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000
ROW15	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	-0.000	-0.006	0.026	0.009	-0.042	-0.024	0.131	-0.000	0.000	0.000	0.000	0.033	-0.053	0.011	-0.000	-0.000	0.000	0.011	0.021	0.000	-0.000	0.000	-0.000	-0.000	0.000
ROW16	0.031	0.015	0.000	0.015	0.031	0.031	0.000	-0.000	0.031	-0.040	0.015	0.000	-0.000	0.000	0.000	0.112	-0.030	-0.000	0.000	0.000	0.000	0.031	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW17	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.015	-0.068	-0.000	-0.000	0.000	0.000	-0.030	0.136	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	0.000	0.000	-0.000
ROW18	0.031	0.015	0.000	0.015	0.031	0.031	0.000	0.031	-0.000	-0.048	0.017	-0.048	0.048	-0.000	-0.000	-0.000	0.000	0.127	-0.034	-0.000	0.000	0.031	-0.000	0.031	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000
ROW19	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.017	-0.078	0.017	-0.078	0.000	0.000	0.000	0.000	-0.034	0.155	-0.000	-0.000	0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW20	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	0.000	-0.035	0.012	-0.000	-0.000	-0.035	0.033	-0.000	0.000	0.000	0.000	0.091	-0.024	0.011	-0.000	0.000	0.000	0.011	0.021	-0.000	0.000	0.000	0.000	0.000	0.000
ROW21	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.012	-0.053	0.000	-0.000	0.012	-0.053	-0.000	-0.000	-0.000	-0.000	-0.024	0.106	0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000
ROW22	0.072	0.039	0.005	0.033	0.067	0.067	0.011	0.041	0.031	0.000	-0.000	-0.000	0.031	-0.000	0.011	0.031	0.000	0.031	0.000	0.011	0.000	0.246	-0.090	0.069	-0.019	0.044	-0.009	-0.051	0.026	-0.077	0.038	-0.077	0.038
ROW23	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.090	0.239	-0.019	0.051	-0.019	0.051	0.026	-0.068	0.038	-0.103	0.038	-0.103
ROW24	0.062	0.031	0.000	0.031	0.062	0.062	0.000	0.031	0.031	-0.000	0.000	0.000	0.062	0.031	0.000	0.000	-0.000	0.031	-0.000	0.000	-0.000	0.069	-0.019	0.203	-0.071	-0.051	0.026	-0.000	0.000	-0.077	0.038	-0.000	-0.000
ROW25	-0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	0.000	-0.019	0.051	-0.071	0.188	0.026	-0.068	0.000	0.000	0.038	-0.103	-0.000	0.000
ROW26	0.011	0.008	0.005	0.003	0.005	0.005	0.011	0.011	0.000	0.000	-0.000	-0.000	0.000	0.000	0.011	-0.000	0.000	0.000	0.000	0.011	0.000	0.044	-0.019	-0.051	0.026	0.146	-0.060	0.000	-0.000	0.000	-0.000	-0.077	0.038
ROW27	0.021	0.016	0.011	0.005	0.011	0.011	0.021	0.021	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.021	0.000	-0.000	0.000	-0.000	0.021	-0.000	-0.009	0.051	0.026	-0.068	-0.060	0.209	-0.000	-0.000	0.000	-0.000	0.038	-0.103
ROW28	-0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.051	0.026	0.000	0.000	0.000	0.000	0.103	-0.051	0.000	0.000	0.000	0.000
ROW29	-0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.026	-0.068	0.000	0.000	0.000	0.000	-0.051	0.137	-0.000	0.000	0.000	0.000
ROW30	0.000	0.000	0.000	0.000	0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	0.000	-0.000	0.000	-0.077	0.038	-0.077	0.038	-0.000	-0.000	0.000	-0.000	0.154	-0.077	0.000	0.000
ROW31	-0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.038	-0.103	0.038	-0.103	0.000	0.000	-0.000	-0.000	-0.077	0.205	0.000	0.000
ROW32	-0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	-0.000	0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.000	-0.000	-0.000	0.000	-0.000	-0.077	0.038	0.000	-0.000	-0.077	0.038	-0.000	-0.000	0.000	0.000	0.154	-0.077
ROW33	0.000	0.000	-0.000	0.000	0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	-0.000	0.000	-0.000	-0.000	0.000	0.000	-0.000	0.000	-0.000	0.000	0.038	-0.103	0.000	0.000	0.038	-0.103	-0.000	0.000	0.000	0.000	-0.077	0.205

Compute ranklhs = rank of the MME = trace of ginvlhs*lhs

ranklhs
29

Compute yhat = vector of solutions for the MME

yhat
106.7223
57.997297
53.624493
48.725
53.097804
8.7456081
-3.887669
70.047804
36.674493
0.6081081
0.1351351
-0.219595
-0.168919
-0.793919
-0.236486
0.7601351
0.1689189
0.1942568
-0.016892
-0.092905
-0.050676
0.5405405
0.2027027
-0.540541
-0.202703
-0.27027
-0.101351
0.2702703
0.1013514
1.386E-13
-1.99E-13
0.1351351
0.0506757

yhat
106.72
58.00
53.62
48.73
53.10
8.75
-3.89
70.05
36.67
0.61
0.14
-0.22
-0.17
-0.79
-0.24
0.76
0.17
0.19
-0.02
-0.09
-0.05
0.54
0.20
-0.54
-0.20
-0.27
-0.10
0.27
0.10
0.00
-0.00
0.14
0.05

Compute sesol = standard error of solutions

sesol
2.09
3.92
4.12
3.56
3.22
5.28
11.29
4.74
4.48
5.81
4.89
2.95
1.96
4.33
3.70
4.37
3.72
4.72
3.74
5.50
4.44
3.78
3.43
3.78
3.43
3.95
3.46
3.95
3.46
4.00
3.46
3.99
3.46

Computation of Additive, Nonadditive, and Total Genetic Predictions

Using matrix computations

Define ka = coefficient matrix of direct and maternal additive genetic predictions deviated from B

ka
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0	1	0	-1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW2	0	0	1	0	-1	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW3	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW4	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW5	0	0.5	0	-0.5	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW6	0	0	0.5	0	-0.5	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW7	0	0.5	0	-0.5	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW8	0	0	0.5	0	-0.5	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW9	0	0.5	0	-0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW10	0	0	0.5	0	-0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW11	0	0.75	0	-0.75	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0
ROW12	0	0	0.75	0	-0.75	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0

ka
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.00	1.00	0.00	-1.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW2	0.00	0.00	1.00	0.00	-1.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW3	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW4	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW5	0.00	0.50	0.00	-0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.00	0.50	0.00	-0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW7	0.00	0.50	0.00	-0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW8	0.00	0.00	0.50	0.00	-0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW9	0.00	0.50	0.00	-0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW10	0.00	0.00	0.50	0.00	-0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW11	0.00	0.75	0.00	-0.75	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW12	0.00	0.00	0.75	0.00	-0.75	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00

Compute kagl = ka*ginvlhs*lhs to check if functions in matrix ka are estimable

(kagl = ka if functions in ka are estimable)

kagl
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	-0.00	0.60	0.40	-0.60	-0.40	0.40	-0.00	0.00	-0.00	1.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00
ROW2	0.00	0.40	0.60	-0.40	-0.60	-0.40	0.00	0.00	-0.00	-0.00	1.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW3	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	1.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW4	0.00	-0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	1.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	0.00
ROW5	-0.00	0.30	0.20	-0.30	-0.20	0.20	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	1.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.20	0.30	-0.20	-0.30	-0.20	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	1.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00
ROW7	0.00	0.30	0.20	-0.30	-0.20	0.20	0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	1.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW8	0.00	0.20	0.30	-0.20	-0.30	-0.20	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	1.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW9	0.00	0.30	0.20	-0.30	-0.20	0.20	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	1.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	0.00
ROW10	0.00	0.20	0.30	-0.20	-0.30	-0.20	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	1.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW11	0.00	0.45	0.30	-0.45	-0.30	0.30	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	1.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW12	0.00	0.30	0.45	-0.30	-0.45	-0.30	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	1.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00

difkaglka
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	-0.00	-0.40	0.40	0.40	-0.40	0.40	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00
ROW2	0.00	0.40	-0.40	-0.40	0.40	-0.40	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW3	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW4	0.00	-0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	0.00
ROW5	-0.00	-0.20	0.20	0.20	-0.20	0.20	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.20	-0.20	-0.20	0.20	-0.20	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00
ROW7	0.00	-0.20	0.20	0.20	-0.20	0.20	0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW8	0.00	0.20	-0.20	-0.20	0.20	-0.20	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW9	0.00	-0.20	0.20	0.20	-0.20	0.20	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	0.00
ROW10	0.00	0.20	-0.20	-0.20	0.20	-0.20	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW11	0.00	-0.30	0.30	0.30	-0.30	0.30	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW12	0.00	0.30	-0.30	-0.30	0.30	-0.30	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00

Compute uaka = vector of multibreed additive genetic predictions

uaka
9.88
0.66
-0.22
-0.17
3.84
0.03
5.40
0.43
4.83
0.25
6.86
0.34

Compute vepuaka = matrix of variance of errors of additive genetic predictions

vepuaka
59.13	30.45	7.56	1.88	16.09	12.33	23.41	13.85	37.85	17.16	37.61	21.39
30.45	69.21	4.74	2.66	16.47	24.58	14.52	34.45	24.34	37.44	23.46	46.90
7.56	4.74	8.71	1.78	6.85	2.97	4.38	2.50	8.14	3.26	7.21	3.86
1.88	2.66	1.78	3.83	0.84	3.03	1.40	1.43	1.83	3.25	1.36	2.84
16.09	16.47	6.85	0.84	25.67	11.81	7.93	8.21	13.72	9.16	20.88	14.14
12.33	24.58	2.97	3.03	11.81	26.73	6.63	12.39	11.03	14.55	12.07	25.65
23.41	14.52	4.38	1.40	7.93	6.63	20.90	9.30	16.14	8.46	15.67	10.57
13.85	34.45	2.50	1.43	8.21	12.39	9.30	25.09	11.55	18.69	11.03	23.42
37.85	24.34	8.14	1.83	13.72	11.03	16.14	11.55	36.49	16.09	25.79	17.68
17.16	37.44	3.26	3.25	9.16	14.55	8.46	18.69	16.09	28.09	13.16	26.00
37.61	23.46	7.21	1.36	20.88	12.07	15.67	11.03	25.79	13.16	44.87	21.27
21.39	46.90	3.86	2.84	14.14	25.65	10.57	23.42	17.68	26.00	21.27	46.52

Compute sepuaka = vector of standard errors of additive genetic predictions

sepuaka
7.69
8.32
2.95
1.96
5.07
5.17
4.57
5.01
6.04
5.30
6.70
6.82

Define kn = coefficient matrix of direct and maternal nonadditive genetic predictions

Assume that males will be mated to (1/2A 1/2B) females and viceversa

kn
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0
ROW2	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0
ROW3	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0
ROW4	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0
ROW5	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0
ROW6	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0
ROW7	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0
ROW8	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0
ROW9	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0
ROW10	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0
ROW11	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5	0
ROW12	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0.5

kn
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW2	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW3	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW4	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW5	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00
ROW7	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00
ROW8	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00
ROW9	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00
ROW10	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00
ROW11	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00
ROW12	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50

Compute kngl = kn*ginvlhs*lhs to check if functions in matrix kn are estimable

(kngl = kn if functions in kn are estimable)

kngl
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.50	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW2	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.50	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.50	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW3	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.50	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW4	0.00	-0.00	0.00	0.00	0.00	-0.00	0.50	0.00	-0.00	0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	0.00	0.00	0.00	0.50	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00
ROW5	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.50	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.00	-0.00	0.00	0.00	-0.00	0.50	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	0.00	0.00	0.50	-0.00	0.00	-0.00	-0.00	-0.00	0.00
ROW7	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.50	-0.00	0.00	0.00	0.00	0.00
ROW8	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.50	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	0.50	-0.00	0.00	0.00	0.00
ROW9	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.50	0.00	0.00	0.00
ROW10	0.00	-0.00	-0.00	0.00	0.00	0.00	0.50	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.50	0.00	0.00
ROW11	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.50	0.00
ROW12	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.50	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.50

difknglkn
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW2	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW3	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW4	0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00
ROW5	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00
ROW7	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW8	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	0.00
ROW9	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW10	0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW11	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW12	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00

Compute uakn = vector of multibreed nonadditive genetic predictions

uakn
4.64
-1.84
4.10
-2.05
4.24
-1.99
4.51
-1.89
4.37
-1.94
4.44
-1.92

Compute vepuakn = matrix of variance of errors of nonadditive genetic predictions

vepuakn
6.45	1.03	3.39	-0.11	4.16	0.17	5.69	0.75	4.92	0.46	5.31	0.60
1.03	34.42	1.09	30.79	-1.99	28.41	0.61	33.05	1.06	32.61	-0.48	31.42
3.39	1.09	6.60	2.29	5.80	1.99	4.19	1.39	5.00	1.69	4.59	1.54
-0.11	30.79	2.29	32.91	-1.37	29.10	0.05	30.85	1.09	31.85	-0.74	29.94
4.16	-1.99	5.80	-1.37	8.89	-0.21	5.07	-1.65	4.98	-1.68	6.52	-1.10
0.17	28.41	1.99	29.10	-0.21	28.27	0.38	28.49	1.08	28.76	-0.02	28.34
5.69	0.61	4.19	0.05	5.07	0.38	8.81	1.78	4.94	0.33	5.38	0.49
0.75	33.05	1.39	30.85	-1.65	28.49	1.78	34.65	1.07	31.95	-0.45	30.77
4.92	1.06	5.00	1.09	4.98	1.08	4.94	1.07	6.96	1.82	4.95	1.07
0.46	32.61	1.69	31.85	-1.68	28.76	0.33	31.95	1.82	33.73	-0.61	30.68
5.31	-0.48	4.59	-0.74	6.52	-0.02	5.38	-0.45	4.95	-0.61	7.91	0.50
0.60	31.42	1.54	29.94	-1.10	28.34	0.49	30.77	1.07	30.68	0.50	31.38

Compute sepuakn = vector of standard errors of nonadditive genetic predictions

sepuakn
2.54
5.87
2.57
5.74
2.98
5.32
2.97
5.89
2.64
5.81
2.81
5.60

Define kt = coefficient matrix of total direct and maternal genetic predictions

Assume that males will be mated to (1/2A 1/2B) females and viceversa

kt
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0	1	0	-1	0	0.5	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0	0
ROW2	0	0	1	0	-1	0	0.5	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0	0
ROW3	0	0	0	0	0	0.5	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0	0
ROW4	0	0	0	0	0	0	0.5	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0	0
ROW5	0	0.5	0	-0.5	0	0.5	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0	0
ROW6	0	0	0.5	0	-0.5	0	0.5	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0	0
ROW7	0	0.5	0	-0.5	0	0.5	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0	0
ROW8	0	0	0.5	0	-0.5	0	0.5	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0	0
ROW9	0	0.5	0	-0.5	0	0.5	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0	0
ROW10	0	0	0.5	0	-0.5	0	0.5	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0	0
ROW11	0	0.75	0	-0.75	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5	0
ROW12	0	0	0.75	0	-0.75	0	0.5	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0.5

kt
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	0.00	1.00	0.00	-1.00	0.00	0.50	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW2	0.00	0.00	1.00	0.00	-1.00	0.00	0.50	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW3	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW4	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW5	0.00	0.50	0.00	-0.50	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.00	0.50	0.00	-0.50	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00
ROW7	0.00	0.50	0.00	-0.50	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00	0.00
ROW8	0.00	0.00	0.50	0.00	-0.50	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00	0.00
ROW9	0.00	0.50	0.00	-0.50	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00	0.00
ROW10	0.00	0.00	0.50	0.00	-0.50	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00	0.00
ROW11	0.00	0.75	0.00	-0.75	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50	0.00
ROW12	0.00	0.00	0.75	0.00	-0.75	0.00	0.50	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	1.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.50

Compute ktgl = kt*ginvlhs*lhs to check if functions in matrix kt are estimable

(ktgl = kt if functions in kt are estimable)

ktgl
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	-0.00	0.70	0.30	-0.70	-0.30	0.80	-0.00	-0.00	-0.00	1.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.50	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW2	0.00	0.40	0.60	-0.40	-0.60	-0.40	0.50	0.00	-0.00	-0.00	1.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	0.00	0.50	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW3	-0.00	0.10	-0.10	-0.10	0.10	0.40	-0.00	-0.00	-0.00	0.00	0.00	1.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.50	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW4	0.00	-0.00	0.00	0.00	0.00	-0.00	0.50	0.00	-0.00	0.00	-0.00	-0.00	1.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	0.00	0.00	0.00	0.50	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00
ROW5	-0.00	0.40	0.10	-0.40	-0.10	0.60	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	1.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.50	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.20	0.30	-0.20	-0.30	-0.20	0.50	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	1.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.50	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW7	-0.00	0.40	0.10	-0.40	-0.10	0.60	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	1.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	0.50	-0.00	0.00	0.00	0.00	0.00
ROW8	0.00	0.20	0.30	-0.20	-0.30	-0.20	0.50	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	1.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.50	-0.00	-0.00	-0.00	-0.00
ROW9	-0.00	0.40	0.10	-0.40	-0.10	0.60	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	1.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.50	0.00	0.00	0.00
ROW10	0.00	0.20	0.30	-0.20	-0.30	-0.20	0.50	0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	0.00	-0.00	-0.00	0.00	1.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.50	-0.00	-0.00
ROW11	-0.00	0.55	0.20	-0.55	-0.20	0.70	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	1.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.50	0.00
ROW12	0.00	0.30	0.45	-0.30	-0.45	-0.30	0.50	0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	1.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	0.50

difktglkt
	COL1	COL2	COL3	COL4	COL5	COL6	COL7	COL8	COL9	COL10	COL11	COL12	COL13	COL14	COL15	COL16	COL17	COL18	COL19	COL20	COL21	COL22	COL23	COL24	COL25	COL26	COL27	COL28	COL29	COL30	COL31	COL32	COL33
ROW1	-0.00	-0.30	0.30	0.30	-0.30	0.30	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW2	0.00	0.40	-0.40	-0.40	0.40	-0.40	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW3	-0.00	0.10	-0.10	-0.10	0.10	-0.10	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00
ROW4	0.00	-0.00	0.00	0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00
ROW5	-0.00	-0.10	0.10	0.10	-0.10	0.10	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW6	0.00	0.20	-0.20	-0.20	0.20	-0.20	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW7	-0.00	-0.10	0.10	0.10	-0.10	0.10	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW8	0.00	0.20	-0.20	-0.20	0.20	-0.20	0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW9	-0.00	-0.10	0.10	0.10	-0.10	0.10	-0.00	0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW10	0.00	0.20	-0.20	-0.20	0.20	-0.20	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00
ROW11	-0.00	-0.20	0.20	0.20	-0.20	0.20	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00	-0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	-0.00	0.00	0.00	0.00	0.00
ROW12	0.00	0.30	-0.30	-0.30	0.30	-0.30	0.00	0.00	-0.00	-0.00	0.00	-0.00	0.00	-0.00	-0.00	-0.00	-0.00	0.00	-0.00	0.00	0.00	0.00	0.00	-0.00	0.00	-0.00	-0.00	0.00	0.00	-0.00	-0.00	-0.00	-0.00

Compute uakt = vector of multibreed total genetic predictions

uakt
14.52
-1.18
3.88
-2.21
8.08
-1.97
9.90
-1.46
9.20
-1.70
11.30
-1.57

Compute vepuakt = matrix of variance of errors of total genetic predictions

vepuakt
77.37	33.47	19.20	6.17	29.61	15.58	35.34	17.26	51.66	20.57	53.49	24.53
33.47	107.02	3.16	35.26	4.78	47.18	11.61	68.97	27.87	73.27	11.31	75.35
19.20	3.16	16.82	3.44	13.32	2.35	11.89	2.45	16.88	3.05	16.26	2.76
6.17	35.26	3.44	37.15	-7.14	25.33	2.78	33.22	7.62	36.83	-8.30	28.55
29.61	4.78	13.32	-7.14	35.56	2.43	16.59	-1.31	22.59	-0.93	32.58	3.60
15.58	47.18	2.35	25.33	2.43	40.40	5.73	34.05	15.16	37.63	1.19	42.04
35.34	11.61	11.89	2.78	16.59	5.73	30.92	9.51	24.74	7.44	25.96	8.67
17.26	68.97	2.45	33.22	-1.31	34.05	9.51	60.18	16.05	52.47	0.16	49.76
51.66	27.87	16.88	7.62	22.59	15.16	24.74	16.05	48.65	21.25	37.12	21.52
20.57	73.27	3.05	36.83	-0.93	37.63	7.44	52.47	21.25	64.93	2.00	53.70
53.49	11.31	16.26	-8.30	32.58	1.19	25.96	0.16	37.12	2.00	60.66	10.50
24.53	75.35	2.76	28.55	3.60	42.04	8.67	49.76	21.52	53.70	10.50	68.69

Compute sepuakt = vector of standard errors of total genetic predictions

sepuakt
8.80
10.35
4.10
6.09
5.96
6.36
5.56
7.76
6.97
8.06
7.79
8.29

vnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	16	6	0	0	0	0	0	0
0	0	0	0	6	12	0	0	0	0	0	0
0	0	0	0	0	0	16	6	0	0	0	0
0	0	0	0	0	0	6	12	0	0	0	0
0	0	0	0	0	0	0	0	16	6	0	0
0	0	0	0	0	0	0	0	6	12	0	0
0	0	0	0	0	0	0	0	0	0	16	6
0	0	0	0	0	0	0	0	0	0	6	12

dnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	12	4.5	0	0	0	0	0	0
0	0	0	0	4.5	9	0	0	0	0	0	0
0	0	0	0	0	0	12	4.5	0	0	0	0
0	0	0	0	0	0	4.5	9	0	0	0	0
0	0	0	0	0	0	0	0	8	3	0	0
0	0	0	0	0	0	0	0	3	6	0	0
0	0	0	0	0	0	0	0	0	0	8	3
0	0	0	0	0	0	0	0	0	0	3	6

vnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	16	6	0	0	0	0	0	0
0	0	0	0	6	12	0	0	0	0	0	0
0	0	0	0	0	0	16	6	0	0	0	0
0	0	0	0	0	0	6	12	0	0	0	0
0	0	0	0	0	0	0	0	16	6	0	0
0	0	0	0	0	0	0	0	6	12	0	0
0	0	0	0	0	0	0	0	0	0	16	6
0	0	0	0	0	0	0	0	0	0	6	12

dnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	12	4.5	0	0	0	0	0	0
0	0	0	0	4.5	9	0	0	0	0	0	0
0	0	0	0	0	0	12	4.5	0	0	0	0
0	0	0	0	0	0	4.5	9	0	0	0	0
0	0	0	0	0	0	0	0	8	3	0	0
0	0	0	0	0	0	0	0	3	6	0	0
0	0	0	0	0	0	0	0	0	0	8	3
0	0	0	0	0	0	0	0	0	0	3	6

vnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	16	6	0	0	0	0	0	0
0	0	0	0	6	12	0	0	0	0	0	0
0	0	0	0	0	0	16	6	0	0	0	0
0	0	0	0	0	0	6	12	0	0	0	0
0	0	0	0	0	0	0	0	16	6	0	0
0	0	0	0	0	0	0	0	6	12	0	0
0	0	0	0	0	0	0	0	0	0	16	6
0	0	0	0	0	0	0	0	0	0	6	12

dnf
16	6	0	0	0	0	0	0	0	0	0	0
6	12	0	0	0	0	0	0	0	0	0	0
0	0	16	6	0	0	0	0	0	0	0	0
0	0	6	12	0	0	0	0	0	0	0	0
0	0	0	0	12	4.5	0	0	0	0	0	0
0	0	0	0	4.5	9	0	0	0	0	0	0
0	0	0	0	0	0	12	4.5	0	0	0	0
0	0	0	0	0	0	4.5	9	0	0	0	0
0	0	0	0	0	0	0	0	8	3	0	0
0	0	0	0	0	0	0	0	3	6	0	0
0	0	0	0	0	0	0	0	0	0	8	3
0	0	0	0	0	0	0	0	0	0	3	6