function [data_por]=sim_doby_do_por(t,X,n,p_0,censoring,n_c)

tic
k=1;
k_c=1;
%pocet=0;
data_por=zeros(1,200);
data_por_c=zeros(1,200);
while k<=n
    nn=25+round(rand(1,1)*10);
    j=round(nn/2+rand(1,1)*5);
    clear p pp
    [Y,H,mi_sigma]=tran_sim_data(t,X,nn);
    m=size(Y,2);
    l=1;
    pp(1)=0;
    for i=j:size(Y,1)
        p(i-j+1)=0;
        [H,p(i-j+1)]=pravd_prekrocenimodely([30 2 60 50 20 400 3800],Y(1:i-1,1),Y(1:i-1,1),Y(1:i-1,1),Y(1:i-1,1),Y(1:i-1,1),Y(1:i-1,1),Y(1:i-1,1),Y(1:i-1,2),Y(1:i-1,3),Y(1:i-1,4),Y(1:i-1,5),Y(1:i-1,6),Y(1:i-1,7),Y(1:i-1,8),Y(i,1),1,3);
        if p(i-j+1)>p_0;
            pp(l)=i;
            l=l+1;
        end;
    end;
    if pp(1)||0;
        zapis=1;
        if k>1;
            for i=1:k-1
                if data_por(i,1) == Y(pp(1),1);
                    zapis=0;
                end;
            end;
        end;
    
        if zapis==1;
            data_por(k,1:m)=Y(pp(1),:);
            k=k+1;
        end;
        
        if censoring==1;
            c=1+round(rand(1,1)*5);
            if k_c<=n_c;
                zapis_c=1;
                if k_c>1;
                    for i=1:k_c-1
                        if data_por_c(i,1) == Y(pp(1)-c,1);
                            zapis_c=0;
                        end;
                    end;
                end;
                
                if zapis_c==1;
                    data_por_c(k_c,1:m)=Y(pp(1)-c,:);
                    k_c=k_c+1;
                end;
                
            end;
        
        end;
    %pocet=pocet+1;
    end;
end;
data_por=data_por(:,1:m);
data_por=[data_por(:,1) zeros(n,1) data_por(:,2:end)];

for i=1:n
     data_por(i,1)
    end    
end    
data_por_c(:,1)=T;

if censoring==1;
    data_por_c=data_por_c(:,1:m);
    data_por_c=[data_por_c(:,1) ones(n_c,1) data_por_c(:,2:end)];
    data_por=[data_por;data_por_c];
    %data_por_2=[zeros(n,1);ones(n_c,1)];
    %data_por=[data_por(:,1) data_por_2 data_por(:,2:end)];
end;
t=toc
data_por