

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: xhl860728@163.com
Tel: +86(10) 6488 9513
Mob:+86 0 13651318417






	[[alternative HTML version deleted]]

