S=S+delay_fun(i-j)*Hospitalization[i]/(Din[i]+0.00001)# I added a small number after Din to avoid dividing by 0 (happens when there is a long sequence of 0 in the hospitalization data)
}
lambda1[j]=lambda0[j]/qj*S
}
lambda0=lambda1
Ein=rep(0,L)# expected number of hospitalization to occur on day j
for(iin1:L){
s=0
for(jin1:i){
s=s+delay_fun(i-j)*lambda0[j]
}
Ein[i]=s
}
w=which(Ein>0)
KHI=1/N0*sum((Ein[w]-Din[w])^2/Ein[w])
}
# lambda1[1]=0
lambda1[L]=0
}
lambda1=lambda1[-seq(1,N1)]
return(lambda1)
}
# Compute the number of infection for different age class and region ----