Matvec Users’ Guide
Matvec Users' Guide
Matvec Users' Guide
- No tags were found...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
52 CHAPTER 9. MACRO PACKAGES<br />
9.1.1 demo() function<br />
Here is the macro script of demo() function:<br />
function demo(program_name)<br />
{<br />
#<br />
# run a demo program in $(MATVECHOME)/examples/interface/txt<br />
#<br />
# Tianlin Wang at UIUC, Sat Jan 28 16:13:05 CST 1995<br />
#<br />
}<br />
local(retval,filename);<br />
if (program_name.class() == "String") {<br />
filename = getenv("MATVECHOME") + "/examples/interface/txt/";<br />
filename += program_name;<br />
retval = input(filename);<br />
}<br />
else {<br />
error("demo(program_name): program_name must be a string");<br />
retval = 0;<br />
}<br />
return retval;<br />
There are a number of demo examples coming with the <strong>Matvec</strong> standard distribution. Usually, these<br />
examples are residing in MATVECHOME/examples/interface/txt directory, where MATVECHOME can be<br />
obtained using <strong>Matvec</strong> function getenv(”MATVECHOME”). These demo examples are listed and briefly<br />
described below:<br />
hello: greetings to every body.<br />
gianola: matrix version of BLUP.<br />
try: variance component estimation based on a five-animal pedigree<br />
VG: single trait variance component estimation.<br />
susan: single trait BLUP.<br />
For examples<br />
> package("demo");<br />
> demo("try")<br />
.matvec.799959823.0<br />
.matvec.799959823.1<br />
iteration sigma_1 ..... sigma_e log_likelihood<br />
0: 1 2 -6.10457<br />
1: 1.74667 2.87176 -6.48199<br />
2: 1.65664 2.91128 -4.73543<br />
3: 2.23916 2.68895 -5.76399<br />
4: 2.22348 2.6937 -5.76309<br />
5: 2.21558 2.69613 -5.76264<br />
6: 2.21156 2.69737 -5.76241<br />
7: 2.20951 2.69801 -5.76229<br />
8: 2.20846 2.69833 -5.76223<br />
9: 2.20792 2.6985 -5.7622