clear;
addpath 'C:\Program files\MATLAB\R2008a\work\SELAPver3\Utilities';

tmp=importdata('C:\Program files\MATLAB\R2008a\work\SELAPver3\MYC_orig_predict_r1_eran.txt');
dat=tmp.data;
dat=logit(dat);
sids=tmp.textdata(2:end,2);
pathways=tmp.textdata(1,3:end);
clear tmp;

use=1:size(dat,2);

penalty=-33;

[mu sig prob]=selapMixMap(dat(:,use),penalty);
componentProbabilities=mvnMixtureProb(dat(:,use)',mu,sig,prob);
