Project

General

Profile

« Previous | Next » 

Revision d6fa6141

Added by Benoit Parmentier over 11 years ago

output resuults analyses, modifications to reflect general name of objects from earlier scripts

View differences:

climate/research/oregon/interpolation/results_interpolation_date_output_analyses.R
4 4 
#Part 1: Script produces plots for every selected date
5 5 
#Part 2: Examine
6 6 
#AUTHOR: Benoit Parmentier
7 
#DATE: 04/02/2013
7 
#DATE: 04/16/2013
8 8 

  		




  9
  9
  
    
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#???--   


  		




  10
  10
  
    

  		




  ......




  31
  31
  
    
#Select relevant dates and load R objects created during the interpolation step


  		




  32
  32
  
    

  		




  33
  33
  
    
##Paths to inputs and output


  		




  34
  
  
    
script_path<-"/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/"


  		




  
  34
  
    
#script_path<-"/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/"


  		




  35
  35
  
    
in_path <- "/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/input_data/"


  		




  36
  36
  
    
out_path<- "/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/output_data/"


  		




  37
  37
  
    
infile_covar<-"covariates__venezuela_region__VE_01292013.tif" #this is an output from covariate script


  		




  38
  38
  
    
date_selected<-c("20000101") ##This is for year 2000!!!


  		




  39
  
  
    
raster_prediction_obj<-load_obj("raster_prediction_obj_dailyTmin_365d_GAM_fus5_all_lstd_03292013.RData")


  		




  
  39
  
    
#raster_prediction_obj<-load_obj("raster_prediction_obj_dailyTmin_365d_GAM_fus5_all_lstd_03292013.RData")


  		




  40
  40
  
    
#out_prefix<-"_365d_GAM_fus5_all_lstd_03132013"


  		




  41
  41
  
    
#out_prefix<-"_365d_GAM_fus5_all_lstd_03142013"                #User defined output prefix


  		




  42
  42
  
    
out_prefix<-"_365d_GAM_fus5_all_lstd_03292013"                #User defined output prefix


  		




  ......




  52
  52
  
    
             "nobs_09","nobs_10","nobs_11","nobs_12")


  		




  53
  53
  
    
covar_names<-c(rnames,lc_names,lst_names)


  		




  54
  54
  
    

  		




  55
  
  
    
list_param_results_analyses<-list(in_path,out_path,script_path,raster_prediction_obj,infile_covar,covar_names,date_selected,var,out_prefix)


  		




  56
  
  
    
names(list_param_results_analyses)<-c("in_path","out_path","script_path","raster_prediction_obj",


  		




  
  55
  
    
list_param_results_analyses<-list(in_path,out_path,script_path,raster_prediction_obj,interpolation_method,infile_covar,covar_names,date_selected,var,out_prefix)


  		




  
  56
  
    
names(list_param_results_analyses)<-c("in_path","out_path","script_path","raster_prediction_obj", "interpolation_method",


  		




  57
  57
  
    
                     "infile_covar","covar_names","date_selected","var","out_prefix")


  		




  58
  58
  
    

  		




  59
  59
  
    
setwd(in_path)


  		




  ......




  83
  83
  
    
  


  		




  84
  84
  
    
  date_selected<-list_param$date_selected


  		




  85
  85
  
    
  var<-list_param$var


  		




  
  86
  
    
  interpolation


  		




  86
  87
  
    
  #gam_fus_mod<-load_obj("gam_fus_mod_365d_GAM_fus5_all_lstd_02202013.RData")


  		




  87
  88
  
    
  #validation_obj<-load_obj("gam_fus_validation_mod_365d_GAM_fus5_all_lstd_02202013.RData")


  		




  88
  89
  
    
  #clim_obj<-load_obj("gamclim_fus_mod_365d_GAM_fus5_all_lstd_02202013.RData")


  		












Also available in:  Unified diff