addpath(genpath('/Users/anhjaeil/work/CNV/bayes_rnaseq/isopure/'))
    

    dname1='t_input.txt';

    
    [alldata,~] = readtext(dname1,'[,\s]','','','numeric');

    tumordata=alldata(:, 5:16);
    normaldata=alldata(:, 1:4);

    addpath(genpath('/Users/anhjaeil/Work/cnv/ISOpureMatlab/ISOpureS1/'))
    [ISOpureS1model, loglikelihood]=learnmodel(tumordata, normaldata);

    
    addpath(genpath('/Users/anhjaeil/Work/cnv/ISOpureMatlab/ISOpureS2/'))
    [ISOpureS2model, loglikelihood]=learnmodel(tumordata, normaldata, ISOpureS1model);

    iter=0
    s0='Ritchie_t'
    s2='_alphapurities.txt'
    s3='_cc_cancerprofiles.txt'
    aa=[s0  s2]  
    cc	=[s0  s3]
    
    A=ISOpureS2model.alphapurities
    C=ISOpureS2model.cc_cancerprofiles

    addpath(genpath('/Users/anhjaeil/Work/cnv/bayes_rnaseq/isopure/'))
    save(aa,  'A', '-ASCII')
    save(cc,  'C', '-ASCII')



	addpath(genpath('/Users/anhjaeil/work/CNV/bayes_rnaseq/isopure/'))
    

    dname1='mrna_input.txt';

    
    [alldata,~] = readtext(dname1,'[,\s]','','','numeric');

    tumordata=alldata(:, 5:16);
    normaldata=alldata(:, 1:4);

    addpath(genpath('/Users/anhjaeil/Work/cnv/ISOpureMatlab/ISOpureS1/'))
    [ISOpureS1model, loglikelihood]=learnmodel(tumordata, normaldata);

    
    addpath(genpath('/Users/anhjaeil/Work/cnv/ISOpureMatlab/ISOpureS2/'))
    [ISOpureS2model, loglikelihood]=learnmodel(tumordata, normaldata, ISOpureS1model);

    iter=0
    s0='manRitchie'
    s2='_alphapurities.txt'
    s3='_cc_cancerprofiles.txt'
    aa=[s0  s2]  
    cc	=[s0  s3]
    
    A=ISOpureS2model.alphapurities
    C=ISOpureS2model.cc_cancerprofiles

    addpath(genpath('/Users/anhjaeil/Work/cnv/bayes_rnaseq/isopure/'))
    save(aa,  'A', '-ASCII')
    save(cc,  'C', '-ASCII')




    
    
    
    


    tumordata=alldata(1:100, 51:100);
    normaldata=alldata(1:100, 1:50);