Revision a1721e6a
Added by Benoit Parmentier over 11 years ago
| climate/research/oregon/interpolation/master_script_temp.R | ||
|---|---|---|
| 80 | 80 |
#met_stations_outfiles_obj_file<-"met_stations_outfiles_obj_kriging_daily__365d_kriging_daily_mults10_lst_comb3_08062013.RData" |
| 81 | 81 |
|
| 82 | 82 |
var<-"TMAX" # variable being interpolated |
| 83 |
out_prefix<-"_365d_gam_day_lst_comb3_08132013" #User defined output prefix
|
|
| 84 |
out_suffix<-"_OR_08132013" #Regional suffix
|
|
| 83 |
out_prefix<-"_365d_gam_day_lst_comb4_08152013" #User defined output prefix
|
|
| 84 |
out_suffix<-"_OR_08152013" #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
|
| ... | ... | |
| 256 | 256 |
#LC1: Evergreen/deciduous needleleaf trees |
| 257 | 257 |
|
| 258 | 258 |
#Combination 3: for paper baseline=s(lat,lon)+s(elev) |
| 259 |
list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
|
|
| 260 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w)", |
|
| 261 |
"y_var ~ s(lat,lon) + s(elev_s) + s(E_w)", |
|
| 262 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST)", |
|
| 263 |
"y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)", |
|
| 264 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LC1)", |
|
| 265 |
"y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)", |
|
| 266 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)", |
|
| 267 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)") |
|
| 259 |
# list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
|
|
| 260 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w)", |
|
| 261 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(E_w)", |
|
| 262 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST)", |
|
| 263 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)", |
|
| 264 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LC1)", |
|
| 265 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)", |
|
| 266 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)", |
|
| 267 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)") |
|
| 268 |
|
|
| 269 |
#Combination 4: for paper baseline=s(lat,lon) |
|
| 270 |
list_models<-c("y_var ~ s(lat,lon)",
|
|
| 271 |
"y_var ~ s(lat,lon) + s(elev_s)", |
|
| 272 |
"y_var ~ s(lat,lon) + s(N_w)", |
|
| 273 |
"y_var ~ s(lat,lon) + s(E_w)", |
|
| 274 |
"y_var ~ s(lat,lon) + s(LST)", |
|
| 275 |
"y_var ~ s(lat,lon) + s(DISTOC)", |
|
| 276 |
"y_var ~ s(lat,lon) + s(LC1)", |
|
| 277 |
"y_var ~ s(lat,lon) + s(CANHGHT)", |
|
| 278 |
"y_var ~ s(lat,lon) + s(LST) + ti(LST,LC1)", |
|
| 279 |
"y_var ~ s(lat,lon) + s(LST) + ti(LST,CANHGHT)") |
|
| 268 | 280 |
|
| 269 | 281 |
#list_models<-c("y_var ~ lat*lon + elev_s")
|
| 270 | 282 |
|
Also available in: Unified diff
running gam daily comb3 OR with revised screening for paper analyses