Displaying 1 result from an estimated 1 matches for "xnews1".
Did you mean:
xnew1
2012 Jan 19
0
Global sensitivity indices using sensitivity package: sobol, sobol2002
...)
u[s2] <- 1-sqrt((1-mode)*(1-u[s2]))
min+(max-min)*u
}
###Loading required packages
#install.packages("odesolve")
library(odesolve)
#install.packages("sensitivity",dependencies=T)
library(sensitivity)
###function that generates responses
simfunzcom<-function(xnews1,So=4250, I=250, R=0, No=4500){
SIRdob1<- function(t, x, parms){
with(as.list(c(parms,x)),{
dS <- (xnews1[, 2]-xnews1[, 4]*N)*(S+R+I*xnews1[, 6]*(1-xnews1[,
9]))-xnews1[, 3]*S + xnews1[,7]*R-(xnews1[, 1]*I*S)/N
dI <- (xnews1[, 1]*I*S)/N + I*(xnews1[, 2]-xnews1[, 3]*N)*xnews1[,
9]*xnews1[, 6...