Displaying 1 result from an estimated 1 matches for "parms1".
Did you mean:
parms
2012 Jan 19
0
Global sensitivity indices using sensitivity package: sobol, sobol2002
...<- xnews1[, 5]*I -(xnews1[, 3] + xnews1[,7])*R
dN<-(xnews1[, 2]-xnews1[, 4]*N)*N-xnews1[, 8]*I-(xnews1[, 2]-xnews1[,
4]*N)*(1-xnews1[, 6])*I
der <- c(dS, dI,dR,dN)
list(der) # the output must be returned
}) # end of 'with'
} # end of function definition
yout<-matrix(0,100,100)
parms1<-xnews1
dt<- seq(1,50,1)
inits1 <- c(S=xnews1[, 10], I=xnews1[, 11], R=xnews1[, 12],N=xnews1[, 13])
for(j in 1:100){
simulation2 <- as.data.frame(lsoda(inits=inits1, times=dt, funct=SIRdob1,
parms=parms1[j, ]))
attach(simulation2)
yout[j]<-as.numeric(simulation2[50 , 3])
}
yout
}...