Re: [R] wrong when use the function fast99 in package "Sensitivity"

2013-07-09 Thread David Winsemius
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",

[R] wrong when use the function fast99 in package "Sensitivity"

2013-07-09 Thread 何亮
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