Revision b337838a
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/05/2013
|
|
| 13 |
#DATE: 07/06/2013
|
|
| 14 | 14 |
|
| 15 | 15 |
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363, TASK$568-- |
| 16 | 16 |
|
| ... | ... | |
| 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_gam_day_lst_comb2_07052013" #User defined output prefix
|
|
| 75 |
out_suffix<-"_OR_07052013"
|
|
| 74 |
out_prefix<-"_365d_gam_day_lst_comb1_07062013" #User defined output prefix
|
|
| 75 |
out_suffix<-"_OR_07062013"
|
|
| 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
|
| ... | ... | |
| 239 | 239 |
#Models to run...this can be change for each run |
| 240 | 240 |
|
| 241 | 241 |
#Combination 1 |
| 242 |
# list_models<-c("y_var ~ s(elev_s)",
|
|
| 243 |
# "y_var ~ s(LST)",
|
|
| 244 |
# "y_var ~ s(elev_s,LST)",
|
|
| 245 |
# "y_var ~ s(lat) + s(lon)+ s(elev_s)",
|
|
| 246 |
# "y_var ~ s(lat,lon,elev_s)",
|
|
| 247 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST)",
|
|
| 248 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(LC6)",
|
|
| 249 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + ti(LC6,LST)",
|
|
| 250 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(DISTOC)")
|
|
| 242 |
list_models<-c("y_var ~ s(elev_s)",
|
|
| 243 |
"y_var ~ s(LST)", |
|
| 244 |
"y_var ~ s(elev_s,LST)", |
|
| 245 |
"y_var ~ s(lat) + s(lon)+ s(elev_s)", |
|
| 246 |
"y_var ~ s(lat,lon,elev_s)", |
|
| 247 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST)", |
|
| 248 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(LC6)", |
|
| 249 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + ti(LC6,LST)", |
|
| 250 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w,E_w) + s(LST) + s(DISTOC)") |
|
| 251 | 251 |
|
| 252 | 252 |
#testing new combinations and tensor for paper |
| 253 |
list_models<-c("y_var ~ s(x,y)",
|
|
| 254 |
"y_var ~ s(x)+ s(y) + s(LST) + ti(x,y) + ti(x,y,LST)", |
|
| 255 |
"y_var ~ s(elev_s) + s(LST) + ti(elev_s,LST)", |
|
| 256 |
"y_var ~ s(lat) + s(lon) + s(LST)", |
|
| 257 |
"y_var ~ s(lat) + s(lon) + s(elev_s)+ ti(lat,lon) + ti(lat,lon,elev_s)", |
|
| 258 |
"y_var ~ s(lat,lon) + s(LST)", |
|
| 259 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST)", |
|
| 260 |
"y_var ~ s(x) + s(y) + s(elev_s) + ti(x,y) + ti(x,y,elev_s)") |
|
| 253 |
# list_models<-c("y_var ~ s(x,y)",
|
|
| 254 |
# "y_var ~ s(x)+ s(y) + s(LST) + ti(x,y) + ti(x,y,LST)",
|
|
| 255 |
# "y_var ~ s(elev_s) + s(LST) + ti(elev_s,LST)",
|
|
| 256 |
# "y_var ~ s(lat) + s(lon) + s(LST)",
|
|
| 257 |
# "y_var ~ s(lat) + s(lon) + s(elev_s)+ ti(lat,lon) + ti(lat,lon,elev_s)",
|
|
| 258 |
# "y_var ~ s(lat,lon) + s(LST)",
|
|
| 259 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST)",
|
|
| 260 |
# "y_var ~ s(x) + s(y) + s(elev_s) + ti(x,y) + ti(x,y,elev_s)")
|
|
| 261 | 261 |
|
| 262 | 262 |
|
| 263 | 263 |
#list_models<-c("y_var~ lat + lon + elev_",
|
Also available in: Unified diff
gam day testing covariates combination for paper OR