Project

General

Profile

« Previous | Next » 

Revision b516681e

Added by Benoit Parmentier over 11 years ago

running kriging fusion comb3 OR with revised screening for paper analyses

View differences:

climate/research/oregon/interpolation/master_script_temp.R
80 80
#met_stations_outfiles_obj_file<-"met_stations_outfiles_obj_kriging_daily__365d_kriging_daily_mults10_lst_comb3_08062013.RData"
81 81

  
82 82
var<-"TMAX" # variable being interpolated
83
out_prefix<-"_365d_gam_cai_lst_comb3_08122013"                #User defined output prefix
83
out_prefix<-"_365d_kriging_fus_lst_comb3_08122013"                #User defined output prefix
84 84
out_suffix<-"_OR_08122013"                                       #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
interpolation_method<-c("gam_CAI") #other otpions to be added later
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 92
#interpolation_method<-c("gwr_CAI") #other otpions to be added later
93 93
#interpolation_method<-c("kriging_CAI") 
......
256 256
#LC1: Evergreen/deciduous needleleaf trees
257 257

  
258 258
#Combination 3: for paper baseline=s(lat,lon)+s(elev)
259
list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
260
              "y_var ~ s(lat,lon) + s(elev_s) + s(N_w)",
261
              "y_var ~ s(lat,lon) + s(elev_s) + s(E_w)",
262
              "y_var ~ s(lat,lon) + s(elev_s) + s(LST)",
263
              "y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)",
264
              "y_var ~ s(lat,lon) + s(elev_s) + s(LC1)",
265
              "y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)",
266
              "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)",
267
              "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)")
259
# list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
260
#                "y_var ~ s(lat,lon) + s(elev_s) + s(N_w)",
261
#                "y_var ~ s(lat,lon) + s(elev_s) + s(E_w)",
262
#                "y_var ~ s(lat,lon) + s(elev_s) + s(LST)",
263
#                "y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)",
264
#                "y_var ~ s(lat,lon) + s(elev_s) + s(LC1)",
265
#                "y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)",
266
#                "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)",
267
#                "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)")
268 268

  
269 269
#list_models<-c("y_var ~ lat*lon + elev_s")
270 270

  
271 271
#list_models<-c("y_var ~ s(lat,lon) + s(elev_s)")
272 272

  
273
#list_models<-c("y_var ~ lat*lon + elev_s",
274
#              "y_var ~ lat*lon + elev_s + N_w",
275
#              "y_var ~ lat*lon + elev_s + E_w",
276
#              "y_var ~ lat*lon + elev_s + LST",
277
#              "y_var ~ lat*lon + elev_s + DISTOC",
278
#              "y_var ~ lat*lon + elev_s + LC1",
279
#             "y_var ~ lat*lon + elev_s + CANHGHT",
280
#              "y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
281
#              "y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)")
273
list_models<-c("y_var ~ lat*lon + elev_s",
274
               "y_var ~ lat*lon + elev_s + N_w",
275
               "y_var ~ lat*lon + elev_s + E_w",
276
               "y_var ~ lat*lon + elev_s + LST",
277
               "y_var ~ lat*lon + elev_s + DISTOC",
278
               "y_var ~ lat*lon + elev_s + LC1",
279
               "y_var ~ lat*lon + elev_s + CANHGHT",
280
               "y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
281
               "y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)")
282 282

  
283 283
#Default name of LST avg to be matched               
284 284
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