4
社区成员




原始代码如下,问题出现在最后一条,报错信息为“number of items to replace is not a multiple of replacement length”,求教如何根据beta.rseq的不同取值生成lambda.rseq的序列。感激不尽!!!
R = seq(1,5,by=0.1);
durInf = 6;
dt = 0.1;
tmax = 365;
numSteps = tmax/dt;
S.r = rep(NA,numSteps);
I.r = rep(NA,numSteps);
R.r = rep(NA,numSteps);
lambda.rseq = rep(NA,numSteps);
incidence.rseq = rep(NA,numSteps);
cumulativeIncidence.rseq = rep(NA,numSteps);
time.r = rep(NA,numSteps);
beta.rseq = c(R/durInf/N);
I.r[1] = 100;
R.r[1] = 0;
S.r[1] = N-I.r[1]-R.r[1];
lambda.rseq[stepIndex] = c(beta.rseq*I.r[stepIndex]);
没用过R语言,帮顶!