Model functions
Context: Model line
and() function
and(t,r)
and
a(t,r)
overlay (add) r times the design matrix for model
term t to the existing design matrix.
Specifically, if the model up to this point has p effects
and t
has a effects,
the a columns of the design matrix for t are multiplied by the scalar r
(default value 1.0) and added to the last a of the p columns already defined.
The overlaid term must agree in size with the term it overlays.
This can be used to force a correlation of 1 between two terms as in a diallel analysis
male and(female)
assuming the i th male is the same individual
as the i th female.
Note that if the overlaid term is complex, it must be predefined; e.g.
Simple Diallel
female !A
male !AS female
ht
dbh
diallel.dat !skip 1
ht dbh ~ Trait !r Tr.male -Tr.female and(Tr.female)
1 2 2
0 # Represents number of trees
Trait 0 US !GP
3*0 # use default initial values
Tr.male 2 # G structure
Trait 0 US !GP
3*0
male
Return to start