plotFlux.m

function C=plotFlux(sample,varargin)
%%----------------------------------------
%% plot the raw data
%%----------------------------------------
    %
    opt = 1;
    if (numel(varargin)>0)
        opt = varargin{1};
    end
    if (ishandle(opt))
        close(opt);
    end
    h1=figure(opt);
    set(h1,'Position',[200,300,1800,900]);
    %
    C = getPmts('Y:\ctrl13\csv');
    C.samp = sample;
    C.D = readtable([C.csv,'\',C.samp,'.csv']);
    C.batch = unique(C.D.batch);    
    %
    annotation('textbox',[0,0.95,1,0.05],...
       'String',C.samp,'EdgeColor','none',...
       'HorizontalAlignment','center','FontSize',12);
    %
    P.axis = [min(C.batch)-1 max(C.batch)+1 -5 5];
    for m=2:7
        subplot(2,3,m-1);
            hold on;
            plotOneMarker(C,m,P);
    end
    shg
    %
end


function plotOneMarker(C,m,P)
    %
    D=C.D;
    marks   = {'o','^','s','<','>','d'};
    color   = colormap(lines);
    batch   = unique(D.batch);
    %
    title(C.markers.marker{m});
    %
    for t=1:numel(C.tissues.ntissue)
        cores  = C.cores(C.cores.ntissue==t,:);
        for i=1:numel(cores.core)
            c = cores.core(i);
            ix = D.nmarker==m & D.core==c;
            d  = D(ix,:);
            d  = sortrows(d,'batch');
            %
            if (numel(d)==0)
                continue
            end
            %
            clr = color(t,:);
            p(t) = plot(d.batch,log(d.f0),'-',...
                'MarkerSize',5, 'MarkerFaceColor',clr,...
                'MarkerEdgeColor',clr*0.4,...
                'Marker',marks{t},'Color',clr);
        end
    end
    %
    axis(P.axis);
    box on
    grid on
    xlabel('batch');
    ylabel('log(flux)');
    xticks((min(batch)-1:max(batch)+1));
    yticks(-5:5);
    legend(p, C.tissues.tissue);
    %
end