Displaying 1 result from an estimated 1 matches for "simulation2".
Did you mean:
simulation
2012 Jan 19
0
Global sensitivity indices using sensitivity package: sobol, sobol2002
...news1[, 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
}
####Input data 1
n<-10
lambda<-rtri(n,min=4,max=7,ml=5.8)
a<-rtri(n,min=0.51,max=0.87,ml=0.64)
b<-rtri(n,min=0.0001,max=0.01,m...