library(XML) # Load the XML packages
library(sensitivity)
source("Yield_APSIM.R")
yield_sa <- fast99(model = Yield_APSIM, factors=c("grains_per_gram_stem",
"potential_grain_filling_rate",
"potential_grain_growth_rate",
"max_grain_size",
"tt_emergence",
"tt_end_of_juvenile",
"tt_floral_initiation",
"tt_flowering",
"tt_start_grain_fill",
"tt_end_grain_fill",
"tt_maturity",
"tt_end_crop",
"tt_harvest_ripe",
"startgf_to_mat",
"est_days_emerg_to_init",
"vern_sens",
"photop_sens",
"tt_end_grain_to_maturity",
"tt_maturity_to_ripe"),
n = 1000,q ="qunif", q.arg =
list("grains_per_gram_stem"=list(min=12.5,max=37.5),
"potential_grain_filling_rate"=list(min=0.001,max=0.005),
"potential_grain_growth_rate"=list(min=0.0005,max=0.0015),
"max_grain_size"=list(min=0.02,max=0.06),
"tt_emergence"=list(min=1,max=1000),
"tt_end_of_juvenile"=list(min=200,max=600),
"tt_floral_initiation"=list(min=277.5,max=832.5),
"tt_flowering"=list(min=60,max=180),
"tt_start_grain_fill"=list(min=400,max=1000),
"tt_end_grain_fill"=list(min=25,max=40),
"tt_maturity"=list(min=1,max=1.5),
"tt_end_crop"=list(min=500,max=1500),
"tt_harvest_ripe"=list(min=500,max=1500),
"startgf_to_mat"=list(min=290,max=870),
"est_days_emerg_to_init"=list(min=41.5,max=124.5),
"vern_sens"=list(min=1,max=5),
"photop_sens"=list(min=1,max=5),
"tt_end_grain_to_maturity"=list(min=30,max=70),
"tt_maturity_to_ripe"=list(min=1,max=1.5))
)
print(yield_sa)
plot(yield_sa)
the results as follow:
Model runs: 1
Estimations of the indices:
first order total order
grains_per_gram_stem NA NA
potential_grain_filling_rate NA NA
potential_grain_growth_rate NA NA
max_grain_size NA NA
tt_emergence NA NA
tt_end_of_juvenile NA NA
tt_floral_initiation NA NA
tt_flowering NA NA
tt_start_grain_fill NA NA
tt_end_grain_fill NA NA
tt_maturity NA NA
tt_end_crop NA NA
tt_harvest_ripe NA NA
startgf_to_mat NA NA
est_days_emerg_to_init NA NA
vern_sens NA NA
photop_sens NA NA
tt_end_grain_to_maturity NA NA
tt_maturity_to_ripe NA NA
why model runs only 1?
--
ºÎÁÁ
Öйú¿ÆѧԺ´óѧ µØÀí¿ÆѧÓë×ÊÔ´Ñо¿Ëù ½µØˮѻ·ÓëµØ±í¹ý³ÌÖصãʵÑéÊÒ
±±¾©ÊÐ ³¯ÑôÇø ´óÍÍ·¼×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.
