Project

General

Profile

« Previous | Next » 

Revision 12b8c55d

Added by Benoit Parmentier about 11 years ago

running gam cai, 0 to 70% monthly holdout comb5 for multi-time scale paper, max temp OR

View differences:

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