This functions creates a stochastic Occurence Multi-Site
Model for the variable x
(PrecipitationOccurenceMultiSiteModel S3 object)
through a calibration from observed data.
data frame (each column is a site) of variable
utilized for the auto-regression of its occurrence, e.g.
daily precipitaton
exogen
exogenous predictors
station
character string vectors containing the
codes of the station used for model calibration
valmin
minimum admitted value for daily
precipitation amount
multisite_type
string indicating the utilized
approach for spatial multi-site dependence description.
Default is "wilks".
tolerance_wilks
see tolerance used by
omega_inv through CCGamma
origin
character string (yyyy-dd-mm) indicating
the date of the first row of "x".
p
auto-regression order
...
further arguments
Value
The function returns a
PrecipitationOccurenceModel-class S3 object
containing the following elements:
... PrecipitationOccurenceModel S3 class
objects for each analyzed site. The name is the site (or
station) code
ccgamaCCGammaObjectListPerEachMonth object,
i.e. matices of Gaussian Inter-Site Correlation returned by
CCGamma;
type string indicating the utilized approach for
spatial multi-site dependence description, only
"wilks" type is implemented;
station character string vectors containing the
codes of the station used in
PrecipitationMultiSiteOccurenceModel.
See Also
PrecipitationOccurenceModel,CCGamma
Examples
library(RGENERATEPREC)
data(trentino)
year_min <- 1961
year_max <- 1990
origin <- paste(year_min,1,1,sep="-")
period <- PRECIPITATION$year>=year_min & PRECIPITATION$year<=year_max
period_temp <- TEMPERATURE_MAX$year>=year_min & TEMPERATURE_MAX$year<=year_max
prec_mes <- PRECIPITATION[period,]
Tx_mes <- TEMPERATURE_MAX[period_temp,]
Tn_mes <- TEMPERATURE_MIN[period_temp,]
accepted <- array(TRUE,length(names(prec_mes)))
names(accepted) <- names(prec_mes)
for (it in names(prec_mes)) {
acc <- TRUE
acc <- (length(which(!is.na(Tx_mes[,it])))==length(Tx_mes[,it]))
acc <- (length(which(!is.na(Tn_mes[,it])))==length(Tn_mes[,it])) & acc
accepted[it] <- (length(which(!is.na(prec_mes[,it])))==length(prec_mes[,it])) & acc
}
valmin <- 1.0
prec_mes <- prec_mes[,accepted]
Tx_mes <- Tx_mes[,accepted]
Tn_mes <- Tn_mes[,accepted]
prec_occurence_mes <- prec_mes>=valmin
station <- names(prec_mes)[!(names(prec_mes) %in% c("day","month","year"))]
station <- station[1:2] # to save example elapsed time!!
exogen <- Tx_mes-Tn_mes
months <- factor(prec_mes$month)
#' ### Not Run!!
# The following lines are commented to save example elapsed time!!
#model_multisite <- PrecipitationOccurenceMultiSiteModel(x=prec_mes,exogen=exogen,
#origin=origin,multisite_type="wilks")
### Not Run!!
# The following lines are commented to save example elapsed time!!
#model_multisite_logit <- PrecipitationOccurenceMultiSiteModel(x=prec_mes,exogen=exogen,
#origin=origin,multisite_type="logit")
###
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(RGENERATEPREC)
Loading required package: copula
Loading required package: RGENERATE
Loading required package: RMAWGEN
Loading required package: chron
Loading required package: date
Loading required package: vars
Loading required package: MASS
Loading required package: strucchange
Loading required package: zoo
Attaching package: 'zoo'
The following objects are masked from 'package:base':
as.Date, as.Date.numeric
Loading required package: sandwich
Loading required package: urca
Loading required package: lmtest
Attaching package: 'vars'
The following object is masked from 'package:copula':
A
Loading required package: blockmatrix
Loading required package: Matrix
Loading required package: stringr
Attaching package: 'stringr'
The following object is masked from 'package:strucchange':
boundary
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/RGENERATEPREC/PrecipitationOccurenceMultiSiteModel.Rd_%03d_medium.png", width=480, height=480)
> ### Name: PrecipitationOccurenceMultiSiteModel
> ### Title: Precipitation Occurence Multi-Site Model
> ### Aliases: PrecipitationOccurenceMultiSiteModel
>
> ### ** Examples
>
> library(RGENERATEPREC)
>
> data(trentino)
>
> year_min <- 1961
> year_max <- 1990
> origin <- paste(year_min,1,1,sep="-")
>
> period <- PRECIPITATION$year>=year_min & PRECIPITATION$year<=year_max
> period_temp <- TEMPERATURE_MAX$year>=year_min & TEMPERATURE_MAX$year<=year_max
>
> prec_mes <- PRECIPITATION[period,]
> Tx_mes <- TEMPERATURE_MAX[period_temp,]
> Tn_mes <- TEMPERATURE_MIN[period_temp,]
> accepted <- array(TRUE,length(names(prec_mes)))
> names(accepted) <- names(prec_mes)
> for (it in names(prec_mes)) {
+ acc <- TRUE
+ acc <- (length(which(!is.na(Tx_mes[,it])))==length(Tx_mes[,it]))
+ acc <- (length(which(!is.na(Tn_mes[,it])))==length(Tn_mes[,it])) & acc
+ accepted[it] <- (length(which(!is.na(prec_mes[,it])))==length(prec_mes[,it])) & acc
+
+ }
>
> valmin <- 1.0
> prec_mes <- prec_mes[,accepted]
>
>
>
> Tx_mes <- Tx_mes[,accepted]
> Tn_mes <- Tn_mes[,accepted]
> prec_occurence_mes <- prec_mes>=valmin
>
> station <- names(prec_mes)[!(names(prec_mes) %in% c("day","month","year"))]
> station <- station[1:2] # to save example elapsed time!!
> exogen <- Tx_mes-Tn_mes
> months <- factor(prec_mes$month)
> #' ### Not Run!!
> # The following lines are commented to save example elapsed time!!
> #model_multisite <- PrecipitationOccurenceMultiSiteModel(x=prec_mes,exogen=exogen,
> #origin=origin,multisite_type="wilks")
>
> ### Not Run!!
> # The following lines are commented to save example elapsed time!!
> #model_multisite_logit <- PrecipitationOccurenceMultiSiteModel(x=prec_mes,exogen=exogen,
> #origin=origin,multisite_type="logit")
> ###
>
>
>
>
>
> dev.off()
null device
1
>