Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/Predictability.R
\name{Predictability}
\alias{Predictability}
\title{Computing scores of predictability using two dinamical proxies
based on dynamical systems theory.}
\usage{
Predictability(dim, theta)
}
\arguments{
\item{dim}{data to create the attractor. Must be a matrix with the timesteps in nrow and the grids in ncol
for instance: dat(time,grids)=(1:24418,1:1060), where time=24418 timesteps and grids=1060 gridpoints}
\item{theta}{list of arbitrary length of secondary grids. Each secondary grid is to be provided as a list of length 2 with longitudes and latitudes}
}
\value{
A list of length 2:
\itemize{
\item\code{pred.dim} {a list of two lists 'qdim' and 'pos.d'. The 'qdim' list
contains values of local dimension 'dim' divided by terciles:
d1: lower tercile and more predictability,
d2: middle tercile,
d3: higher tercile and less predictability
The 'pos.d' list contains the position of each tercile in parameter 'dim'}
\item\code{pred.theta} {a list of two lists 'qtheta' and 'pos.t'.
The 'qtheta' list contains values of the inverse of persistence 'theta'
divided by terciles:
th1: lower tercile and more predictability,
th2: middle tercile,
th3: higher tercile and less predictability
The 'pos.t' list contains the position of each tercile in parameter 'theta'}
}
dyn_scores values from 0 to 1. A dyn_score of 1 indicates higher
predictability.
}
\description{
This function divides in terciles the two dynamical proxies
computed with CST_ProxiesAttractor or ProxiesAttractor. These terciles will
be used to measure the predictability of the system in dyn_scores. When the
local dimension 'dim' is small and the inverse of persistence 'theta' is
small the predictability is higher, and viceversa.
}
\examples{
# Creating an example of matrix dat(time,grids):
m <- matrix(rnorm(2000) * 10, nrow = 50, ncol = 40)
# imposing a threshold
quanti <- 0.90
# computing dyn_scores from parameters dim and theta of the attractor
attractor <- ProxiesAttractor(dat = m, quanti = 0.60, iplot = FALSE)
predyn <- Predictability(dim = attractor$dim, theta = attractor$theta)
}
\references{
Faranda, D., Alvarez-Castro, M.C., Messori, G., Rodriguez, D.,
and Yiou, P. (2019). The hammam effect or how a warm ocean enhances large
scale atmospheric predictability.Nature Communications, 10(1), 1316.
DOI = https://doi.org/10.1038/s41467-019-09305-8 "
Faranda, D., Gabriele Messori and Pascal Yiou. (2017).
Dynamical proxies of North Atlantic predictability and extremes.
Scientific Reports, 7-41278, 2017.
}
\author{
Carmen Alvarez-Castro, \email{carmen.alvarez-castro@cmcc.it}
Maria M. Chaves-Montero, \email{mdm.chaves-montero@cmcc.it}
Veronica Torralba, \email{veronica.torralba@cmcc.it}
Davide Faranda, \email{davide.faranda@lsce.ipsl.fr}
}