Download A Programming Language for Building Likelihood Models Version 2.1

Programming tutorial
', sig SD = ', sig_sd,
', true = ', sitsd)
PRINTLN('Absolute bias = ', sitsd - sig_mean,
', % bias = ', 100*sig_mean/sitsd)
PRINTLN('t test: (param<>0) t = ', sig_mean/sig_sd)
PRINTLN('t test: (param=', sitsd, ') t = ',
(sig_mean - sitsd)/sig_sd)
END {then}
END {model}
END
An even more complicated simulation program
This program simulates repeated datasets, each containing observations of a
bilateral morphological trait. The simulation includes the ability to add, for
example, a directional size bias. “Noise” of development is superimposed on the
underlying trait, and different variances in the noise can be specified for each
side.
159