Matvec Users’ Guide

62 CHAPTER 10. LINEAR MODEL ANALYSES
• In statement
D.input("try.dat","animal\$ herd _skip y");
the dollar sign (\$) after animal is used to indicate that this data-column is string column. The _skip
is one of reserved data-column name meaning to skip or ignore the column.
• In statement
M.equation("y = intercept + herd + animal");
the intercept is a reserved model-term meaning an intercept or a grand-mean. The + between model
terms is optional. The interaction term or the nested term can be expressed in the same way as SAS.
For instance, A*B and x(A). The above statement about model can also be re-written as
M.equation("y = intercept +",
"herd + animal");
This should give you enough hint about how to break down a model with dozens of terms.
If you add a semicolon(;) after statement M.blup(), then BLUP solution would not be displayed on
screen.
Now, at the Matvec prompt, I type
> input try2
i = 1 i = 2 i = 3 i = 4 i = 5
7.42 0 1.69 0.26 -0.26
i = 6 i = 7 i = 8
0.26 -0.67 0.67
If you want a more descriptive output for BLUP, then use statement M.save("try.out") right after
M.blup(). The following could be the possible contents in "try.out":
some extra information in the model
------------------------------------------------
variance for animal =
Col 1
Row 1 1
residual variance =
Col 1
Row 1 2
MME dimension : 8
non-zeros in MME: 23
basic statistics for dependent variables
------------------------------------------------
trait-name n mean std
y 6 8.33333 2.06559
--------------------------------------------------
BLUP (BLUE, Posterior_Mean)
-----------------------------------------------