function PlottMWelch(Datum, NormZP, RepGr1, RepGr2)

load(sprintf('%i/%s/%i_norm_%s.mat', Datum, addZeros(NormZP, 5), Datum, addZeros(NormZP, 5)));
iNorm = find(Zeit(:,4)==NormZP);

iRepGr = [RepGr1 RepGr2];
[iRepGr  Names] = benennung(LComp, iRepGr);

L = length(ZPListe); % Anzahl Zeitpunkte

for t = 1:L                         
    eval(sprintf('n=n%s;', addZeros(Zeit(t,4), 5)));
    p(t) = WelchTest(M(t,RepGr1), S(t,RepGr1), N(t,RepGr1), M(t,RepGr2), S(t,RepGr2), N(t,RepGr2));
end

figure
hold on

for g = 1:2
    y(:,g) = M(:,iRepGr(g));
    fehler(:,g) = S(:, iRepGr(g));    
    RGB(g,:) = LFarbe(iRepGr(g),:);
    FehlerRGB(g,:) = 0.2+RGB(g,:);
    FehlerRGB(find(FehlerRGB>1))=1;
    
    xmax = Zeit(length(Zeit),1);
    
       
    text(xmax+2, y(length(Zeit),g), LComp(iRepGr(g)), 'Color', RGB(g,:), 'FontSize', 10, 'FontWeight', 'bold')
    text(xmax+8, y(length(Zeit),g), LConc(iRepGr(g)), 'Color', RGB(g,:), 'FontSize', 10, 'FontWeight', 'bold')   
    
    %nur jeder vierte Wert plotten

    for i=1:floor(length(Zeit)/4)
        k=i*4;
        FehlerZeit(i)=Zeit(k,1);
        FehlerY(i,g)=y(k,g);
        Fehler(i,g)=fehler(k,g);
     
    end
    errorbar(FehlerZeit, FehlerY(:,g), Fehler(:,g), 'Color', FehlerRGB(g,:))     
end

xlim=xmax+20;
ylim=ceil(max(max(y+fehler)));

% for i=1:ylim
%     ylimlog(i)=0.1/10^i;
%     ylimi(i)=i;
% end
ylimi=[0 0.5 1 1.5 2];
ylimlog=[0.1 0.01 0.001 0.0001 0.00001];
axis([0 xlim 0 ylim])

title(sprintf('Date: %i / Normalized to %i h %i min',Datum, Zeit(iNorm,2), Zeit(iNorm,3)))


plot(Zeit(:,1), y(:,1),'LineWidth', 2, 'Color', RGB(1,:));
plot(Zeit(:,1), y(:,2),'LineWidth', 2, 'Color', RGB(2,:));

[AX,H1,H2]=plotyy(Zeit(:,1), y(:,2), Zeit(:,1), (p));

%fuer Zoya: macht Linie bei p=0.05
yLinie=(-log10(0.05)-1)/2;
plot([0 xlim],[yLinie yLinie], 'r', 'LineWidth', 2)

xlabel('Time [h]');
set(AX, 'XLimMode', 'manual', 'XLim', [0 xlim], 'YColor', 'k')

ylabel(AX(1), 'Cell Index')
set(AX(1), 'YLimMode', 'manual', 'YLim', [0 ylim], 'YTick', ylimi, 'YGrid','on', 'XGrid', 'on')
set (H1, 'Color', RGB(2,:))

ylabel(AX(2), 'p-value (Welch)')
set(AX(2), 'YDir', 'rev', 'Yscale', 'log')
set(AX(2), 'YLimMode', 'manual', 'YLim', [min(ylimlog) max(ylimlog)])
set(AX(2), 'YMinorTick', 'on', 'YTick', [10^-5 10^-4 10^-3 10^-2 0.05 0.1])%fliplr(ylimlog))
set(get(AX(2),'YLabel'),'Color','k'); 


hold off
% 
% speichern = input('Do you want to save the figure? (Yes = 1) ');
% if (speichern == 1)
%     Dateiname = sprintf('Bilder/PlotMWelch/jpg/%i-%s_%s.jpg', Datum, addZeros(NormZP, 5), Names)
%     saveas(gcf, Dateiname, 'jpg');
%     Dateiname = sprintf('Bilder/PlotMWelch/fig/%i-%s_%s.fig', Datum, addZeros(NormZP, 5), Names)
%     saveas(gcf, Dateiname, 'fig');
%     
%     close
% elseif (speichern == 0)
%     close
% end