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