Project

General

Profile

« Previous | Next » 

Revision 2cf9c858

Added by Benoit Parmentier about 11 years ago

running gwr CAI with monthly hold out from 0.4 to 0.7 in OR, paper analyses

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/16/2013
13 
#DATE: 09/18/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_gwr_CAI_lst_comb3_09172013"                #User defined output prefix


  		




  84
  
  
    
out_suffix<-"_OR_09172013"                                       #Regional suffix


  		




  
  83
  
    
out_prefix<-"_365d_gwr_CAI_lst_comb3_09182013"                #User defined output prefix


  		




  
  84
  
    
out_suffix<-"_OR_09182013"                                       #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.2,0.3)  #if prop_min=prop_max and step=0 then predictions are done for the number of dates...


  		




  
  257
  
    
prop_minmax_month <-c(0.4,0.7)  #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")


  		












Also available in:  Unified diff