matthiasm@8: function mk_tetrad_data_file(filename, samples, sig)
matthiasm@8: % MK_TETRAD_DATA_FILE Make a file containing raw discrete data for input to TETRAD
matthiasm@8: % mk_tetrad_data_file(filename, samples, sig)
matthiasm@8: %
matthiasm@8: % samples(i,j) is the value for case i, variable j
matthiasm@8: % The resulting file can be used for the 'build' part of Tetrad.
matthiasm@8: % For details on tetrad, see hss.cmu.edu/html/departments/philosophy/TETRAD/tetrad.html
matthiasm@8: 
matthiasm@8: [nsamples N] = size(samples);
matthiasm@8: 
matthiasm@8: fid = fopen(filename, 'w');
matthiasm@8: fprintf(fid, '/Raw\n');
matthiasm@8: fprintf(fid, '%d\n', nsamples);
matthiasm@8: for i=1:N
matthiasm@8:   fprintf(fid, 'x%d ', i);
matthiasm@8: end
matthiasm@8: fprintf(fid, '\n');
matthiasm@8: for i=1:nsamples
matthiasm@8:   fprintf(fid, '%d ', samples(i,:)-1); % tetrad counts from 0
matthiasm@8:   fprintf(fid, '\n');
matthiasm@8: end
matthiasm@8: %fprintf(fid, '/Knowledge\n');
matthiasm@8: %fprintf(fid, 'Significance %4.2f\n', sig);
matthiasm@8: fclose(fid);
matthiasm@8: