Revision 12b8c55d
Added by Benoit Parmentier about 11 years ago
climate/research/oregon/interpolation/master_script_temp.R | ||
---|---|---|
10 | 10 |
#STAGE 5: Output analyses: assessment of results for specific dates... |
11 | 11 |
# |
12 | 12 |
#AUTHOR: Benoit Parmentier |
13 |
#DATE: 11/27/2013
|
|
13 |
#DATE: 11/19/2013
|
|
14 | 14 |
|
15 | 15 |
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363, TASK$568-- |
16 | 16 |
|
... | ... | |
81 | 81 |
|
82 | 82 |
var<-"TMAX" # variable being interpolated |
83 | 83 |
#out_prefix<-"_365d_gam_cai_lst_comb3_10102013" #User defined output prefix |
84 |
out_prefix<-"_365d_kriging_cai_lst_mults_0_70_comb5_11272013" #User defined output prefix
|
|
84 |
out_prefix<-"_365d_gam_cai_lst_mults_0_70_comb5_11192013" #User defined output prefix
|
|
85 | 85 |
|
86 |
out_suffix<-"_OR_11272013" #Regional suffix
|
|
86 |
out_suffix<-"_OR_11192013" #Regional suffix
|
|
87 | 87 |
out_suffix_modis <-"_05302013" #pattern to find tiles produced previously |
88 | 88 |
|
89 | 89 |
#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later |
90 |
#interpolation_method<-c("gam_CAI") #other otpions to be added later
|
|
90 |
interpolation_method<-c("gam_CAI") #other otpions to be added later |
|
91 | 91 |
#interpolation_method<-c("gam_fusion") #other otpions to be added later |
92 | 92 |
#interpolation_method<-c("kriging_fusion") #other otpions to be added later |
93 | 93 |
#interpolation_method<-c("gwr_fusion") #other otpions to be added later |
94 | 94 |
#interpolation_method<-c("gwr_CAI") #other otpions to be added later |
95 |
interpolation_method<-c("kriging_CAI") |
|
95 |
#interpolation_method<-c("kriging_CAI")
|
|
96 | 96 |
|
97 | 97 |
#interpolation_method<-c("gam_daily") #other otpions to be added later |
98 | 98 |
#interpolation_method<-c("kriging_daily") #other otpions to be added later |
... | ... | |
255 | 255 |
#prop_minmax<-c(0.1,0.7) #if prop_min=prop_max and step=0 then predictions are done for the number of dates... |
256 | 256 |
|
257 | 257 |
seed_number_month <- 100 |
258 |
nb_sample_month <- 1 #number of time random sampling must be repeated for every hold out proportion
|
|
258 |
nb_sample_month <-1 #number of time random sampling must be repeated for every hold out proportion |
|
259 | 259 |
#step_month <-0 |
260 | 260 |
step_month <-0.1 |
261 | 261 |
constant_month <- 0 #if value 1 then use the same samples as date one for the all set of dates |
... | ... | |
274 | 274 |
#LC1: Evergreen/deciduous needleleaf trees |
275 | 275 |
|
276 | 276 |
#Combination 5: for paper multi-timescale paper |
277 |
#list_models<-c("y_var ~ s(lat,lon)",
|
|
278 |
# "y_var ~ s(lat,lon) + s(LST)",
|
|
279 |
# "y_var ~ s(lat,lon) + s(elev_s)",
|
|
280 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w)",
|
|
281 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)",
|
|
282 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST)",
|
|
283 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)")
|
|
277 |
list_models<-c("y_var ~ s(lat,lon)", |
|
278 |
"y_var ~ s(lat,lon) + s(LST)", |
|
279 |
"y_var ~ s(lat,lon) + s(elev_s)", |
|
280 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w)", |
|
281 |
"y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)", |
|
282 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST)", |
|
283 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)") |
|
284 | 284 |
|
285 | 285 |
#Combination 5: for paper multi-timescale paper |
286 |
list_models<-c("y_var ~ lat*lon", |
|
287 |
"y_var ~ lat*lon + LST", |
|
288 |
"y_var ~ lat*lon + elev_s", |
|
289 |
"y_var ~ lat*lon + elev_s + N_w*E_w", |
|
290 |
"y_var ~ lat*lon + elev_s + DISTOC", |
|
291 |
"y_var ~ lat*lon + elev_s + LST", |
|
292 |
"y_var ~ lat*lon + elev_s + LST + I(LST*LC1)") |
|
286 |
#list_models<-c("y_var ~ lat*lon",
|
|
287 |
# "y_var ~ lat*lon + LST",
|
|
288 |
# "y_var ~ lat*lon + elev_s",
|
|
289 |
# "y_var ~ lat*lon + elev_s + N_w*E_w",
|
|
290 |
# "y_var ~ lat*lon + elev_s + DISTOC",
|
|
291 |
# "y_var ~ lat*lon + elev_s + LST",
|
|
292 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*LC1)")
|
|
293 | 293 |
|
294 | 294 |
#Combination 3: for paper baseline=s(lat,lon)+s(elev) |
295 | 295 |
#list_models<-c("y_var ~ s(lat,lon) + s(elev_s)", |
Also available in: Unified diff
running gam cai, 0 to 70% monthly holdout comb5 for multi-time scale paper, max temp OR