Project

General

Profile

« Previous | Next » 

Revision 7fd742b4

Added by Benoit Parmentier about 11 years ago

running gam_CAI with gam and DISTOC for daily deviation surface

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: 10/04/2013
13 
#DATE: 10/10/2013
14 14 

  		




  15
  15
  
    
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363, TASK$568--   


  		




  16
  16
  
    

  		




  ......




  67
  67
  
    
source(file.path(script_path,"sampling_script_functions_08252013.R"))


  		




  68
  68
  
    
source(file.path(script_path,"GAM_fusion_function_multisampling_10042013.R")) #Includes Fusion and CAI methods


  		




  69
  69
  
    
source(file.path(script_path,"interpolation_method_day_function_multisampling_07052013.R")) #Include GAM_day


  		




  70
  
  
    
source(file.path(script_path,"GAM_fusion_function_multisampling_validation_metrics_10042013.R"))


  		




  
  70
  
    
source(file.path(script_path,"GAM_fusion_function_multisampling_validation_metrics_10102013.R"))


  		




  71
  71
  
    

  		




  72
  72
  
    
#stages_to_run<-c(1,2,3,4,5) #May decide on antoher strategy later on...


  		




  73
  73
  
    
#stages_to_run<-c(0,2,3,4,5) #May decide on antoher strategy later on...


  		




  ......




  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_cai_lst_comb3_10042013"                #User defined output prefix


  		




  84
  
  
    
out_suffix<-"_OR_10042013"                                       #Regional suffix


  		




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


  		




  
  84
  
    
#out_prefix<-"_365d_kriging_daily_mults1_lst_comb3_10102013"                #User defined output prefix


  		




  
  85
  
    

  		




  
  86
  
    
out_suffix<-"_OR_10102013"                                       #Regional suffix


  		




  85
  87
  
    
out_suffix_modis <-"_05302013"                       #pattern to find tiles produced previously     


  		




  86
  88
  
    

  		




  87
  89
  
    
#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later


  		




  ......




  299
  301
  
    
                "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)")#,


  		




  300
  302
  
    
                #"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)")


  		




  301
  303
  
    

  		




  302
  
  
    
list_models2<-c("y_var ~ s(lat,lon) + s(elev_s)",


  		




  303
  
  
    
                "y_var ~ s(lat,lon) + s(LST)")


  		




  
  304
  
    
list_models2<-c("y_var ~ s(lat,lon) + s(DISTOC)")


  		




  304
  305
  
    

  		




  305
  306
  
    
interp_method2 <- "gam" #other options are "gwr" and "kriging"


  		




  306
  307
  
    
#list_models<-c("y_var ~ lat*lon + elev_s")


  		












Also available in:  Unified diff