8.5.1 How do I export a numerator relationship matrix so that I can read it into Octave?

Numerator relationship matrices may be exported to text files for use with, e.g., Octave using the NewAMatrix.save method:

example = pyp_newclasses.loadPedigree(optionsfile='denny.ini')
amatrix = pyp_newclasses.NewAMatrix(example.kw)
amatrix.form_a_matrix(example.pedigree)
amatrix.tofile('Ainv.txt')

When matrices are written to text files array elements are separated by sepchar (Table 3.1).

Matrices may also be written to a binary format. The default value of the nrm_format pedigree option is text. To write files in binary format you must either specify the value of the nrm_format option as binary before you load your pedigree file or use the nrm_format keyword when you call NewAMatrix.save:

amatrix.tofile('Ainv.bin',nrm_format='binary')

Once you've saved the NRM to a file, say 'Ainv.txt', in text format you can easily read it into Octave:

octave:1> myfile = fopen ("Ainv.txt", "r");
octave:2> ainv = fscanf(myfile,'%f',[18,18])

This has been verified with Octave 2.1.57 under RedHat Enterprise Linux on small matrices.

See About this document... for information on suggesting changes.