Project

General

Profile

« Previous | Next » 

Revision 92c36e33

Added by Benoit Parmentier over 12 years ago

  • ID 92c36e334b16f2d8e77dce99fab4aeae1d442881

made minor changes to PRISM variables

View differences:

climate/research/oregon/interpolation/Linear_reg.R
7 7

  
8 8
###Loading r library and packages
9 9
library(raster)                                                                        # loading the raster package
10
library(gtools)                                                                       l# loading ...
10
library(gtools)                                                                        # loading ...
11 11
library(sp)
12 12
library(mgcv)
13 13

  
......
53 53

  
54 54
###Regression part 1: linear models
55 55

  
56
lm_ANUSPLIN1=lm(tmax~lat+lon+ELEV_SRTM, data=ghcn1507_s)
57
lm_PRISM1=lm(tmax~lat+lon+ELEV_SRTM+ASPECT+DISTOC, data=ghcn1507_s) #Note that a variable on inversion is missing
56
lm_ANUSPLIN1<-lm(tmax~lat+lon+ELEV_SRTM, data=ghcn1507_s)
57
lm_PRISM1<-lm(tmax~lat+lon+ELEV_SRTM+ASPECT+DISTOC, data=ghcn1507_s) #Note that a variable on inversion is missing
58 58
summary(lm_ANUSPLIN1)
59 59
summary(lm_PRISM1)
60 60

  
61 61
###Regression part 2: GAM models
62
GAM_PRISM1<-gam(tmax~ s(lat) + s (lon) + s (ELEV_SRTM), data=ghcn1507_s)
63 62
GAM_ANUSPLIN1<-gam(tmax~ s(lat) + s (lon) + s (ELEV_SRTM), data=ghcn1507_s)
63
GAM_PRISM1<-gam(tmax~ s(lat) + s (lon) + s (ELEV_SRTM) + s (ASPECT)+ s(DISTOC), data=ghcn1507_s)
64 64
#use the s() for smoothing function
65 65

  
66 66
###Compare the models
......
105 105
AIC_all<-AIC(lm_ANUSPLIN1,GAM_ANUSPLIN1,lm_PRISM1,GAM_PRISM1)
106 106
results<-data.frame(model=mod_name,RMSE=RMSE_all,df=AIC_all$df,AIC=AIC_all$AIC)
107 107

  
108

  
109 108
#dump("results", file= paste(path,"/","results_reg1507.txt",sep=""))
110 109
write.csv(results, file= paste(path,"/","results_reg1507.txt",sep=""))
111 110

  

Also available in: Unified diff