Revision 0e76845c
Added by Benoit Parmentier over 11 years ago
| 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: 07/11/2013
|
|
| 13 |
#DATE: 07/15/2013
|
|
| 14 | 14 |
|
| 15 | 15 |
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363, TASK$568-- |
| 16 | 16 |
|
| ... | ... | |
| 56 | 56 |
#source(file.path(script_path,"download_and_produce_MODIS_LST_climatology_06112013.R")) |
| 57 | 57 |
source(file.path(script_path,"covariates_production_temperatures_07052013.R")) |
| 58 | 58 |
source(file.path(script_path,"Database_stations_covariates_processing_function_06112013.R")) |
| 59 |
source(file.path(script_path,"GAM_fusion_analysis_raster_prediction_multisampling_07052013.R"))
|
|
| 59 |
source(file.path(script_path,"GAM_fusion_analysis_raster_prediction_multisampling_07142013.R"))
|
|
| 60 | 60 |
source(file.path(script_path,"results_interpolation_date_output_analyses_06112013.R")) |
| 61 | 61 |
#source(file.path(script_path,"results_covariates_database_stations_output_analyses_04012013.R")) #to be completed |
| 62 | 62 |
|
| ... | ... | |
| 71 | 71 |
stages_to_run<-c(0,2,3,4,5) #May decide on antoher strategy later on... |
| 72 | 72 |
|
| 73 | 73 |
var<-"TMAX" # variable being interpolated |
| 74 |
out_prefix<-"_365d_kriging_day_lst_comb3_07112013" #User defined output prefix
|
|
| 75 |
out_suffix<-"_OR_07112013"
|
|
| 74 |
out_prefix<-"_365d_gam_day_lst_comb4_07152013" #User defined output prefix
|
|
| 75 |
out_suffix<-"_OR_07152013"
|
|
| 76 | 76 |
out_suffix_modis <-"_05302013" #use tiles produce previously |
| 77 | 77 |
|
| 78 | 78 |
#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later
|
| 79 | 79 |
#interpolation_method<-c("gam_CAI") #other otpions to be added later
|
| 80 | 80 |
#interpolation_method<-c("gam_fusion") #other otpions to be added later
|
| 81 |
#interpolation_method<-c("gam_daily") #other otpions to be added later
|
|
| 82 |
interpolation_method<-c("kriging_daily") #other otpions to be added later
|
|
| 81 |
interpolation_method<-c("gam_daily") #other otpions to be added later
|
|
| 82 |
#interpolation_method<-c("kriging_daily") #other otpions to be added later
|
|
| 83 | 83 |
#interpolation_method<-c("gwr_daily") #other otpions to be added later
|
| 84 | 84 |
|
| 85 | 85 |
#out_path <- paste("/home/parmentier/Data/IPLANT_project/Venezuela_interpolation/Venezuela_01142013/output_data",
|
| ... | ... | |
| 238 | 238 |
|
| 239 | 239 |
#Models to run...this can be change for each run |
| 240 | 240 |
#LC1: Evergreen/deciduous needleleaf trees |
| 241 |
|
|
| 242 |
#Combination 4: for paper baseline=s(lat,lon) |
|
| 243 |
list_models<-c("y_var ~ s(lat,lon)",
|
|
| 244 |
"y_var ~ s(lat,lon) + s(elev_s)", |
|
| 245 |
"y_var ~ s(lat,lon) + s(N_w)", |
|
| 246 |
"y_var ~ s(lat,lon) + s(E_w)", |
|
| 247 |
"y_var ~ s(lat,lon) + s(LST)", |
|
| 248 |
"y_var ~ s(lat,lon) + s(DISTOC)", |
|
| 249 |
"y_var ~ s(lat,lon) + s(LC1)", |
|
| 250 |
"y_var ~ s(lat,lon) + s(CANHGHT)") |
|
| 251 |
# "y_var ~ s(lat,lon) + s(LST) + ti(LST,LC1)", |
|
| 252 |
# "y_var ~ s(lat,lon) + s(LST) + ti(LST,CANHGHT)") |
|
| 253 |
|
|
| 241 | 254 |
#Combination 3: for paper baseline=s(lat,lon)+s(elev) |
| 242 | 255 |
#list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
|
| 243 | 256 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w)", |
| ... | ... | |
| 250 | 263 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)") |
| 251 | 264 |
|
| 252 | 265 |
# list_models<-c("y_var ~ s(elev_s)",
|
| 253 |
# "y_var ~ s(LST)", |
|
| 254 |
# "y_var ~ s(elev_s,LST)", |
|
| 255 |
# "y_var ~ s(lat) + s(lon)+ s(elev_s)", |
|
| 256 |
# "y_var ~ s(lat,lon,elev_s)", |
|
| 257 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST)", |
|
| 258 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(LC6)", |
|
| 259 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + ti(LC6,LST)", |
|
| 260 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(DISTOC)") |
|
| 261 |
|
|
| 262 |
#testing new combinations and tensor for paper
|
|
| 263 |
list_models<-c("y_var ~ lat*lon + elev_s",
|
|
| 264 |
"y_var ~ lat*lon + elev_s + N_w", |
|
| 265 |
"y_var ~ lat*lon + elev_s + E_w",
|
|
| 266 |
"y_var ~ lat*lon + elev_s + LST", |
|
| 267 |
"y_var ~ lat*lon + elev_s + DISTOC",
|
|
| 268 |
"y_var ~ lat*lon + elev_s + LC1",
|
|
| 269 |
"y_var ~ lat*lon + elev_s + CANHGHT",
|
|
| 270 |
"y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
|
|
| 271 |
"y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)") |
|
| 272 |
|
|
| 273 |
|
|
| 274 |
#list_models<-c("y_var~ lat + lon + elev_",
|
|
| 275 |
# "y_var~ I(lat*lon) + elev_s",
|
|
| 276 |
# "y_var~ lat + lon + elev_s + N + E + DISTOC",
|
|
| 277 |
# "y_var~ I(lat*lon) + elev_s + I(N*E) + DISTOC + LST",
|
|
| 278 |
# "y_var~ lat + lon + ELEV_SRTM + N_w + E_w + DISTOC + LST",
|
|
| 279 |
# "y_var~ lat + lon + ELEV_SRTM + N_w + E_w + DISTOC + LST + LC2",
|
|
| 280 |
# "y_var~ lat + lon + ELEV_SRTM + Northness_w + Eastness_w + DISTOC + LST + LC6",
|
|
| 281 |
# "y_var~ lat + lon + ELEV_SRTM + Northness_w + Eastness_w + DISTOC + LST + I(LST*LC2)", |
|
| 266 |
# "y_var ~ s(LST)",
|
|
| 267 |
# "y_var ~ s(elev_s,LST)",
|
|
| 268 |
# "y_var ~ s(lat) + s(lon)+ s(elev_s)",
|
|
| 269 |
# "y_var ~ s(lat,lon,elev_s)",
|
|
| 270 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST)",
|
|
| 271 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(LC6)",
|
|
| 272 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + ti(LC6,LST)",
|
|
| 273 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(DISTOC)")
|
|
| 274 |
|
|
| 275 |
#testing combination 3, part 1 for GWR day
|
|
| 276 |
#list_models<-c("y_var ~ lat*lon + elev_s",
|
|
| 277 |
# "y_var ~ lat*lon + elev_s + N_w",
|
|
| 278 |
# "y_var ~ lat*lon + elev_s + E_w")
|
|
| 279 |
|
|
| 280 |
#testing combination 3, part 2 for GWR day
|
|
| 281 |
#list_models<-c("y_var ~ lat*lon + elev_s + LST",
|
|
| 282 |
# "y_var ~ lat*lon + elev_s + DISTOC",
|
|
| 283 |
# "y_var ~ lat*lon + elev_s + LC1")
|
|
| 284 |
|
|
| 285 |
#list_models<-c("y_var ~ lat*lon + elev_s",
|
|
| 286 |
# "y_var ~ lat*lon + elev_s + N_w", |
|
| 287 |
# "y_var ~ lat*lon + elev_s + E_w",
|
|
| 288 |
# "y_var ~ lat*lon + elev_s + LST",
|
|
| 289 |
# "y_var ~ lat*lon + elev_s + DISTOC",
|
|
| 290 |
# "y_var ~ lat*lon + elev_s + LC1",
|
|
| 291 |
# "y_var ~ lat*lon + elev_s + CANHGHT",
|
|
| 292 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
|
|
| 293 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)")
|
|
| 294 |
|
|
| 282 | 295 |
|
| 283 | 296 |
#Default name of LST avg to be matched |
| 284 | 297 |
lst_avg<-c("mm_01","mm_02","mm_03","mm_04","mm_05","mm_06","mm_07","mm_08","mm_09","mm_10","mm_11","mm_12")
|
Also available in: Unified diff
gam day, predictions for paper testing covariates combination 4