function [beta_vec,gamma_vec,t_offset_vec]=my_bootstrap()

% n_bootstrap=100;
n_bootstrap=50;

temp=load('data1.mat');
I_data=temp.data1.I;
t_data=temp.data1.t;

[beta,gamma,t_offset]=get_params();
[~,I_model,~,t]=sir_solve(beta,gamma,t_offset,1);

nn=length(t)-length(t_data);
corr=0.1;
residual= log(I_model(nn+1:end)+corr)-log(I_data+corr);

beta_vec=0*(1:n_bootstrap);
gamma_vec=0*(1:n_bootstrap);
t_offset_vec=0*(1:n_bootstrap);

for i=1:n_bootstrap
    resampled_residual = randsample(residual,length(residual),'true');
    corr=0.1;
    I_data_resampled=max(exp(log(I_model(nn+1:end)+corr)-resampled_residual)-corr,corr);
    [beta,gamma,t_offset]=get_params1(t_data,I_data_resampled);

    beta_vec(i)=beta;
    gamma_vec(i)=gamma;
    t_offset_vec(i)=t_offset;
end


end