Revision 26f78228
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: 09/14/2013
|
|
13 |
#DATE: 09/16/2013
|
|
14 | 14 |
|
15 | 15 |
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363, TASK$568-- |
16 | 16 |
|
... | ... | |
80 | 80 |
#met_stations_outfiles_obj_file<-"met_stations_outfiles_obj_gam_CAI__365d_gam_CAI_lst_comb3_08252013.RData" |
81 | 81 |
|
82 | 82 |
var<-"TMAX" # variable being interpolated |
83 |
out_prefix<-"_365d_kriging_fus_lst_comb3_09142013" #User defined output prefix
|
|
83 |
out_prefix<-"_365d_gwr_CAI_lst_comb3_09162013" #User defined output prefix
|
|
84 | 84 |
out_suffix<-"_OR_09142013" #Regional suffix |
85 | 85 |
out_suffix_modis <-"_05302013" #pattern to find tiles produced previously |
86 | 86 |
|
87 | 87 |
#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later |
88 | 88 |
#interpolation_method<-c("gam_CAI") #other otpions to be added later |
89 | 89 |
#interpolation_method<-c("gam_fusion") #other otpions to be added later |
90 |
interpolation_method<-c("kriging_fusion") #other otpions to be added later |
|
90 |
#interpolation_method<-c("kriging_fusion") #other otpions to be added later
|
|
91 | 91 |
#interpolation_method<-c("gwr_fusion") #other otpions to be added later |
92 |
#interpolation_method<-c("gwr_CAI") #other otpions to be added later
|
|
92 |
interpolation_method<-c("gwr_CAI") #other otpions to be added later |
|
93 | 93 |
#interpolation_method<-c("kriging_CAI") |
94 | 94 |
|
95 | 95 |
#interpolation_method<-c("gam_daily") #other otpions to be added later |
... | ... | |
253 | 253 |
seed_number_month <- 100 |
254 | 254 |
nb_sample_month <-1 #number of time random sampling must be repeated for every hold out proportion |
255 | 255 |
step_month <-0.1 |
256 |
constant_month <-0 #if value 1 then use the same samples as date one for the all set of dates |
|
257 |
prop_minmax_month <-c(0.4,0.7) #if prop_min=prop_max and step=0 then predictions are done for the number of dates...
|
|
256 |
constant_month <- 0 #if value 1 then use the same samples as date one for the all set of dates
|
|
257 |
prop_minmax_month <-c(0,0.1) #if prop_min=prop_max and step=0 then predictions are done for the number of dates...
|
|
258 | 258 |
|
259 | 259 |
#dates_selected<-c("20100101","20100102","20100103","20100901") # Note that the dates set must have a specific format: yyymmdd |
260 | 260 |
#dates_selected<-c("20100101","20100102","20100301","20100302","20100501","20100502","20100701","20100702","20100901","20100902","20101101","20101102") |
Also available in: Unified diff
runnin gwr CAI, monthly hold out 0 to 0.1 for OR, paper analyses