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