On Jul 8, 2013, at 11:57 PM, 何亮 wrote:
>
>
> Pn_I<-function(x){
>
> I<-1000
> Rp<-x[4]*(1-x[5]*I)*I/(1+x[6]*I)
> Pn<-x[1]*(1-x[2]*I)*I/(1+x[3]*I)-Rp-x[7]
> return(Pn)
>
> }
>
>
>
>
> Pn_I_sensitivity <- fast99(model = Pn_I,
> factors=c("alpha","beta","gamma","delta","lammda","omiga",
Pn_I<-function(x){
I<-1000
Rp<-x[4]*(1-x[5]*I)*I/(1+x[6]*I)
Pn<-x[1]*(1-x[2]*I)*I/(1+x[3]*I)-Rp-x[7]
return(Pn)
}
Pn_I_sensitivity <- fast99(model = Pn_I,
factors=c("alpha","beta","gamma","delta","lammda","omiga","Rd"),
n = 3000,q =rep("qunif",7), q.ar
2 matches
Mail list logo