Revision 7145d96b
Added by Benoit Parmentier over 11 years ago
| climate/research/oregon/interpolation/master_script_temp.R | ||
|---|---|---|
| 78 | 78 |
met_stations_outfiles_obj_file<-"/data/project/layers/commons/data_workflow/output_data_365d_gam_fus_lst_test_run_07172013/met_stations_outfiles_obj_gam_fusion__365d_gam_fus_lst_test_run_07172013.RData" |
| 79 | 79 |
|
| 80 | 80 |
var<-"TMAX" # variable being interpolated |
| 81 |
out_prefix<-"_365d_kriging_cai_lst_comb3_07312013" #User defined output prefix
|
|
| 82 |
out_suffix<-"_OR_07312013" #Regional suffix
|
|
| 81 |
out_prefix<-"_365d_gam_cai_lst_comb3_08012013" #User defined output prefix
|
|
| 82 |
out_suffix<-"_OR_08012013" #Regional suffix
|
|
| 83 | 83 |
out_suffix_modis <-"_05302013" #pattern to find tiles produced previously |
| 84 | 84 |
|
| 85 | 85 |
#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later
|
| 86 |
#interpolation_method<-c("gam_CAI") #other otpions to be added later
|
|
| 86 |
interpolation_method<-c("gam_CAI") #other otpions to be added later
|
|
| 87 | 87 |
#interpolation_method<-c("gam_fusion") #other otpions to be added later
|
| 88 | 88 |
#interpolation_method<-c("kriging_fusion") #other otpions to be added later
|
| 89 | 89 |
#interpolation_method<-c("gwr_fusion") #other otpions to be added later
|
| 90 | 90 |
#interpolation_method<-c("gwr_CAI") #other otpions to be added later
|
| 91 |
interpolation_method<-c("kriging_CAI")
|
|
| 91 |
#interpolation_method<-c("kriging_CAI")
|
|
| 92 | 92 |
|
| 93 | 93 |
#interpolation_method<-c("gam_daily") #other otpions to be added later
|
| 94 | 94 |
#interpolation_method<-c("kriging_daily") #other otpions to be added later
|
| ... | ... | |
| 253 | 253 |
#LC1: Evergreen/deciduous needleleaf trees |
| 254 | 254 |
|
| 255 | 255 |
#Combination 3: for paper baseline=s(lat,lon)+s(elev) |
| 256 |
# list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
|
|
| 257 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w)",
|
|
| 258 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(E_w)",
|
|
| 259 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST)",
|
|
| 260 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)",
|
|
| 261 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LC1)",
|
|
| 262 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)",
|
|
| 263 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)",
|
|
| 264 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)")
|
|
| 265 |
|
|
| 266 |
list_models<-c("y_var ~ lat*lon + elev_s",
|
|
| 267 |
"y_var ~ lat*lon + elev_s + N_w", |
|
| 268 |
"y_var ~ lat*lon + elev_s + E_w", |
|
| 269 |
"y_var ~ lat*lon + elev_s + LST", |
|
| 270 |
"y_var ~ lat*lon + elev_s + DISTOC", |
|
| 271 |
"y_var ~ lat*lon + elev_s + LC1", |
|
| 272 |
"y_var ~ lat*lon + elev_s + CANHGHT", |
|
| 273 |
"y_var ~ lat*lon + elev_s + LST + I(LST*LC1)", |
|
| 274 |
"y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)") |
|
| 256 |
list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
|
|
| 257 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w)", |
|
| 258 |
"y_var ~ s(lat,lon) + s(elev_s) + s(E_w)", |
|
| 259 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST)", |
|
| 260 |
"y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)", |
|
| 261 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LC1)", |
|
| 262 |
"y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)", |
|
| 263 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)", |
|
| 264 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)") |
|
| 265 |
|
|
| 266 |
# list_models<-c("y_var ~ lat*lon + elev_s",
|
|
| 267 |
# "y_var ~ lat*lon + elev_s + N_w",
|
|
| 268 |
# "y_var ~ lat*lon + elev_s + E_w",
|
|
| 269 |
# "y_var ~ lat*lon + elev_s + LST",
|
|
| 270 |
# "y_var ~ lat*lon + elev_s + DISTOC",
|
|
| 271 |
# "y_var ~ lat*lon + elev_s + LC1",
|
|
| 272 |
# "y_var ~ lat*lon + elev_s + CANHGHT",
|
|
| 273 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
|
|
| 274 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)")
|
|
| 275 | 275 |
|
| 276 | 276 |
#Default name of LST avg to be matched |
| 277 | 277 |
lst_avg<-c("mm_01","mm_02","mm_03","mm_04","mm_05","mm_06","mm_07","mm_08","mm_09","mm_10","mm_11","mm_12")
|
Also available in: Unified diff
running GAM_CAI with comb3 for comparison, OR temp