clear all; close all; runs = 3000; % Sample size N Para_settings;% LHS matrix rm_LHS=LHS_Call(0.036, rm, 42.5, 0 ,runs,'unif'); delta_LHS=LHS_Call(0.051, delta, 0.093, 0 ,runs,'unif'); dl_LHS=LHS_Call(0.213, dl, 16.9, 0, runs,'unif'); kappa_LHS=LHS_Call(0.0001, kappa ,1, 0 ,runs,'unif'); dm_LHS=LHS_Call(0.016, dm, 0.07, 0, runs,'unif'); bm_LHS=LHS_Call(0.2, bm, 0.75, 0 ,runs,'unif'); beta_m_LHS=LHS_Call(0.02, beta_m, 0.24, 0, runs,'unif'); beta_b_LHS= LHS_Call(0.8 , beta_b , 1 , 0 , runs , 'unif'); k_LHS=LHS_Call(0.087,k,0.125,0 ,runs,'unif'); mu_LHS=LHS_Call(0.12,mu,0.2, 0 ,runs,'unif'); gamma_LHS=LHS_Call(0,gamma,0.2, 0 ,runs,'unif'); % LHS matrix and para. labels LHSmatrix=[rm_LHS delta_LHS dl_LHS kappa_LHS dm_LHS bm_LHS beta_m_LHS beta_b_LHS k_LHS mu_LHS gamma_LHS]; for x=1:runs f=@ODE_WNV; x LHSmatrix(x,:) [t,y]=ode15s(@(t,y)f(t,y,LHSmatrix,x,runs),tspan,y0,[]); A=[t y]; if size(A,1)>=max(time_points+1) Time_lhs(:,x)=A(time_points+1,1); L_lhs(:,x)=A(time_points+1,2); Mm_lhs(:,x)=A(time_points+1,3); Sm_lhs(:,x)=A(time_points+1,4); Em_lhs(:,x)=A(time_points+1,5); Imm_lhs(:,x)=A(time_points+1,6); Sb_lhs(:,x)=A(time_points+1,7); Ib_lhs(:,x)=A(time_points+1,8); Rb_lhs(:,x)=A(time_points+1,9); Xb_lhs(:,x)=A(time_points+1,10); end R00 = R0(LHSmatrix,x,t); if length(R00)>=max(time_points+1) R0_lhs(:,x) = R00(time_points+1); end end alpha = 0.05; %PRCC for adult mosquito abundance [prcc sign sign_label]=PRCCM(LHSmatrix,(Mm_lhs+Em_lhs+Sm_lhs+Imm_lhs),1:length(time_points),PRCC_var,alpha); %PRCC for R_0 [prcc sign sign_label]=PRCCR0(LHSmatrix,R0_lhs,1:length(time_points),PRCC_var,alpha);