14.09.2015 Views

Matvec Users’ Guide

Matvec Users' Guide

Matvec Users' Guide

SHOW MORE
SHOW LESS
  • 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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!