请问大牛，由多变量的微分方程组构成的模型，其中多个自变量都具有不同的时滞，能够用延迟微分方程DDE在R语言环境中实现吗？
具体如下模型中的两个微分方程，代码中的E和In两个变量，E需要采用t-6的变量值，In需要采用t-5的变量值，
#dE/dt=beta_1*S*In/N+beta_2*S*E/N-k*E
#dIn/dt=k*E-z*In

可能描述的不是很准确，敬请不吝赐教。