Optional: Rscript plot_postgibbsf90.R can use postind file as input

[bonifazi]

postind file has the index of the G and R effects in order that gibbs_samples follows. It can be used to make it easier to plot G and R effects directly in one file. Now the user has to provide the G and R subset of columns (with the first 3 cols) separately. Thus, providing postind file would make life easier.

[bonifazi]

Also need to account for user-defined se_covar_se parameters.
Example of posting file in such case:

run4 $ cat postind
G matrix for effect = 6
1 2 3 4 5 6
0 7 8 9 10 11
0 0 12 13 14 15
0 0 0 16 17 18
0 0 0 0 19 20
0 0 0 0 0 21
R matrix
22 0 0 23 0 0
0 24 0 0 25 0
0 0 26 0 0 27
0 0 0 28 0 0
0 0 0 0 29 0
0 0 0 0 0 30

Functions of (co)variances
31 rg12 - g_6_6_1_2/(g_6_6_1_1g_6_6_2_2)^0.5
32 rg13 - g_6_6_1_3/(g_6_6_1_1g_6_6_3_3)^0.5
33 rg14 - g_6_6_1_4/(g_6_6_1_1g_6_6_4_4)^0.5
34 rg15 - g_6_6_1_5/(g_6_6_1_1g_6_6_5_5)^0.5
35 rg16 - g_6_6_1_6/(g_6_6_1_1g_6_6_6_6)^0.5
36 rg23 - g_6_6_2_3/(g_6_6_2_2g_6_6_3_3)^0.5
37 rg24 - g_6_6_2_4/(g_6_6_2_2g_6_6_4_4)^0.5
38 rg25 - g_6_6_2_5/(g_6_6_2_2g_6_6_5_5)^0.5
39 rg26 - g_6_6_2_6/(g_6_6_2_2g_6_6_6_6)^0.5
40 rg34 - g_6_6_3_4/(g_6_6_3_3g_6_6_4_4)^0.5
41 rg35 - g_6_6_3_5/(g_6_6_3_3g_6_6_5_5)^0.5
42 rg36 - g_6_6_3_6/(g_6_6_3_3g_6_6_6_6)^0.5
43 rg45 - g_6_6_4_5/(g_6_6_4_4g_6_6_5_5)^0.5
44 rg46 - g_6_6_4_6/(g_6_6_4_4g_6_6_6_6)^0.5
45 rg56 - g_6_6_5_6/(g_6_6_5_5g_6_6_6_6)^0.5
46 h2t1 - g_6_6_1_1/(g_6_6_1_1+r_1_1)
47 h2t2 - g_6_6_2_2/(g_6_6_2_2+r_2_2)
48 h2t3 - g_6_6_3_3/(g_6_6_3_3+r_3_3)
49 h2t4 - g_6_6_4_4/(g_6_6_4_4+r_4_4)
50 h2t5 - g_6_6_5_5/(g_6_6_5_5+r_5_5)
51 h2t6 - g_6_6_6_6/(g_6_6_6_6+r_6_6)
52 ret1t4 - r_1_4/(r_1_1r_4_4)^0.5
53 ret2t5 - r_2_5/(r_2_2r_5_5)^0.5
54 ret2t6 - r_3_6/(r_3_3r_6_6)^0.5