Project

General

Profile

« Previous | Next » 

Revision 34598372

Added by Benoit Parmentier about 11 years ago

running gam fusion with 20-30% monthly holdout, multitimescale paper, 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: 09/09/2013                                                                                 
13
#DATE: 09/10/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_gam_fus_lst_comb3_09092013"                #User defined output prefix
84
out_suffix<-"_OR_09092013"                                       #Regional suffix
83
out_prefix<-"_365d_gam_fus_lst_comb3_09102013"                #User defined output prefix
84
out_suffix<-"_OR_09102013"                                       #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
......
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 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...
257
prop_minmax_month <-c(0.2,0.3)  #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")
......
318 318
                                "list_models","lst_avg","out_path","script_path","use_clim_image","join_daily",
319 319
                                "interpolation_method")
320 320

  
321
debug(raster_prediction_fun)
321
#debug(raster_prediction_fun)
322 322
raster_prediction_obj <-raster_prediction_fun(list_param_raster_prediction)
323 323

  
324 324
############## STAGE 5: OUTPUT ANALYSES ##################

Also available in: Unified diff