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.arg = 
list(alpha=list(min=0,max=1),
                           
beta=list(min=0,max=1),gamma=list(min=0,max=1),delta=list(min=0,max=1),
                           lammda=list(min=0,max=1),omiga=list(min=0,max=1), 
Rd=list(min=0,max=4)))


print(Pn_I_sensitivity)
plot(Pn_I_sensitivity)


the results as follow

Error in response(x, ...) : 
  (list) object cannot be coerced to type 'double'


does anybody know where is wrong in my code?
--
ºÎÁÁ
Öйú¿ÆѧԺ´óѧ µØÀí¿ÆѧÓë×ÊÔ´Ñо¿Ëù ½µØˮѭ»·ÓëµØ±í¹ý³ÌÖصãʵÑéÊÒ
±±¾©ÊÐ ³¯ÑôÇø ´óÍÍ·¼×11ºÅ
µç»°£º13651318417
Liang He
Institute of Geographic Sciences and Natural resource research
Chinese Academy of Science,
No.11, Datun Road, Chaoyang District,
Beijing 100101, China
Email: xhl860...@163.com
Tel: +86(10) 6488 9513
Mob:+86 0 13651318417






        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to