Project

General

Profile

« Previous | Next » 

Revision 1ca79db5

Added by Benoit Parmentier about 11 years ago

running gwr daily OR temp 2010 with combination 5 for multi timescale methods

View differences:

climate/research/oregon/interpolation/master_script_temp.R
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_gam_fss_lst_comb5_11062013"                #User defined output prefix
84
out_prefix<-"_365d_gwr_daily_lst_comb5p1_3_11062013"                #User defined output prefix
85 85

  
86 86
out_suffix<-"_OR_11062013"                                       #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 90
#interpolation_method<-c("gam_CAI") #other otpions to be added later
91
interpolation_method<-c("gam_fusion") #other otpions to be added later
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
......
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
99
#interpolation_method<-c("gwr_daily") #other otpions to be added later
99
interpolation_method<-c("gwr_daily") #other otpions to be added later
100 100

  
101 101
#out_path<-"/home/parmentier/Data/IPLANT_project/Oregon_interpolation/Oregon_03142013/output_data"
102 102
out_path <- "/data/project/layers/commons/Oregon_interpolation/output_data"
......
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",
286
list_models<-c("y_var ~ lat*lon",
287
               "y_var ~ lat*lon + LST",
288
               "y_var ~ lat*lon + elev_s")
289 289
#               "y_var ~ lat*lon + elev_s + N_w*E_w",
290 290
#               "y_var ~ lat*lon + elev_s + DISTOC",
291 291
#               "y_var ~ lat*lon + elev_s + LST",

Also available in: Unified diff