an approach to the extreme value distribution of non ... · * falk et al., laws of small numbers:...
TRANSCRIPT
![Page 1: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/1.jpg)
25~26 Feb. 2005 Cherry Bud Workshop 2005 (Keio University)
An approach to the extreme value distribution of non-stationary process
Hang CHOI & Jun KANDA
Institute of Environmental StudiesGraduate School of Frontier Sciences
The University of Tokyo
![Page 2: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/2.jpg)
01 Theoretical Frameworks of EVA
1) Stationary Random Process (Sequence)- Distribution Ergodicity- (In)dependent and Identically Distributedrandom variables (I.I.D. assumption)
→ strict stationarity* Conventional approach
2) Non-stationary Random Process (Sequence)- (In)dependent but non-identically Distributedrandom variables (non-I.I.D. assumption)→ weak stationarity, non-stationarity
3) Ultimate (Asymptotic) and penultimate forms
![Page 3: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/3.jpg)
1) I.I.D. random variable Approach
Ultimate form : Fisher-Tippett theorem (1928)
lim lim ( ) ( )nn nn nn n
n
Z aP x F a b x G xb→+∞ →+∞
⎛ ⎞−≤ = + = ∈⎜ ⎟
⎝ ⎠F
{ }1max ,...,n nZ X X=1, , ( ) ,nX X F x∈K
=F { Gumbel, Fréchet, Weibull }
GEVD, GPD, POT-GPD + MLM, PWM, MOM etc.
R.A. Fisher & L.H.C. Tippett (1928), Limiting forms of the frequency distribution of the largest or smallest member of a sample, Proc. Cambridge Philosophical Society, Vol 24, p180~190
![Page 4: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/4.jpg)
reduced variates y=-log(-log G(x))
xFréchet
Weibull
Gumbel
![Page 5: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/5.jpg)
2) non-I.I.D. random variable Approach
{ }1 1 1( ), , ( ), max , ,n n n nX F x X F x Z X X∈ ∈ =K K
1
lim lim ( ) ( )n
n nj j jn n
jn
Z cP x F d c x Q xd→+∞ →+∞
=
⎛ ⎞−≤ = + = ∈⎜ ⎟
⎝ ⎠∏ S
⊂F S=S { Gumbel, Fréchet, Weibull,…..}
* Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994
The class of EVD for non-i.i.d. case is much larger.
![Page 6: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/6.jpg)
3) Ultimate / penultimate form andfinite epoch T in engineering practice
In engineering practice, the epoch of interest, T is finite.e.g. annual maximum value, monthly maximum value and
maximum/minimum pressure coefficients in 10min etc.
As such, the number of independent random variables min the epoch T becomes a finite integer, i.e. m<∞, and consequently, the theoretical framework for ultimate formis no longer available regardless i.i.d. or non-i.i.d.case.
Following discussions are restricted on the penultimate d.f.for the extremes of non-stationary random process in a finite epoch T.
![Page 7: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/7.jpg)
02 Statement of the Problem
( ), : [0, )x t t +∈ = ∞。Let be the continuous observation record ofa non-stationary continuous stochastic process X(t) and assumethat X(t) is a mean square differentiable process and hold thefollowing condition.
[ ]( , ) : ( ) ( ) ( , ) log 0 asr t E X t X t r t Tτ τ τ τ τ= + ⇒ → →
As such, how to estimate the extreme value distribution of X(t) in an epoch T<∞ from the continuous observation record x(t)?
![Page 8: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/8.jpg)
03 Assumptions and Formulation
According to the conventional approach in wind engineering, let assume the non-stationary process X(t) can be partitioned with a finite epoch T, in which the partitioned process Xi(t),(i-1)T≦t < iT, can be assumed as an independent stationary random process, and define the d.f. FZ(x) as follows.
( ) ( ; : sup ( ), [0, ))ZF x P Z x Z X t t T= ≤ = ∈ < ∞
Then, by the Glivenko-Cantelli theorem and the block maximaapproach
( )( , ]1 1
1 1( ) lim : sup ( ) lim ( )n n
Z x i i in ni iF x I Z X t P Z x
n n−∞→∞ →∞= =
= = = ≤∑ ∑where ( ) 1 if else 0cI x x c= ∈
![Page 9: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/9.jpg)
04 i.i.d. random sequence (EQRS) approach of ( )iP Z x≤
By partitioning the interval [(i-1)T,iT) into finite subpartitions[( 1) , ), 1 [ / ] ij h jh j T h m− ≤ ≤ = in the manner of that
( )( )* *
1
( ) ; : sup ( ), [( 1) , ( )i
i
i
mm
i j j i Zj
P Z x P Z x Z X t t j h jh F x<∞
=
≤ = ≤ = ∈ − =∏
the required d.f. FZ(x) can be defined as follows.
1
1( ) : lim ( )i
i
nm
Z Zn iF x F x
n→∞=
= ∑
If it is possible to assume that all mi≈m, then
1
1( ) : lim ( )i
nm
Z Zn iF x F x
n→∞=
= ∑
![Page 10: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/10.jpg)
05 lower bound of FZ(x)
From the inequality (geometric mean)<(arithmetic mean), a lowerbound of FZ(x) can be defined as follows:
/
11
1( ) : lim ( ) lim ( ) ( )i i
n nm n m
Z Z Z Zn n ii
F x F x F x F xn→∞ →∞
==
= < =∑∏%
06 alternative definition of FZ(x)
1
1ˆ ( ) : ( ) ( )i
mnm
Z Z Zi
F x F x F xn =
⎛ ⎞= =⎜ ⎟⎝ ⎠∑
This definition can be found easily in engineering applications andmay be reasonable for the case of
1 nZ ZF F; L ; .
![Page 11: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/11.jpg)
07 Inequality of quantile functions Let define quantile functions of each definition as follows.
1 1 1ˆ ˆ( ) : ( ), ( ) : ( ), ( ) : ( )Z Z ZQ F Q F Q Fα α α α α α− − −= = =% %
Then, by the inequality for the means and the comparison of the distribution of order statistics, i.e. Zn:n and Z1:n,
ˆ( ) ( ), ( ) ( ) (0,1)Q Q Q Q for allα α α α α≤ < ∈% %
ˆ ( ) ( ) ( ) for largeˆ( ) ( ) ( ) for small
Q Q Q
Q Q Q
α α α α
α α α α
⎧ < →⎪⎨
< <⎪⎩
%
%
Therefore, the alternative definition results in smaller variance ofextremes.
![Page 12: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/12.jpg)
Complement for the inequalities of quantile functions
/: : : :
1 1 1
~ , ~ ( ) ( )nn n n
m m n mn n i n n i i n n n n
i i i
Z F Z F F Z Zα α= = =
⎛ ⎞= ⇒ =⎜ ⎟
⎝ ⎠∏ ∏ ∏% %
/: : : :
1 11 1
1 1~ , ~ , ( ) ( )nmn mnn nn n
m m nn n i n n i i i n n n n
i ii i
Z F Z F F F Z Zn n
α α= == =
⎛ ⎞⎛ ⎞ ⎛ ⎞< ⇒ >⎜ ⎟⎜ ⎟ ⎜ ⎟⎝ ⎠ ⎝ ⎠⎝ ⎠∑ ∑∏ ∏% %
( ) ( )2 21: 1: 1
1 11
1: 1:
1~ 1 1 ( ), ~ 1 1 ( )
( ) ( )
n n nnm m m m m mn i i n im
i ii
n n
Z F F o F Z F F o Fn
Z Zα α
−= ==
− − − − − −
⇒ <
∑ ∑∏ ; ;
( )
/1: 1:
1 1
/1: 1:
1 1
~ 1 (1 ), ~ 1 1 ,
(1 ) 1 1 ( ) ( )
nn nm m n
n i n ii i
nn nnm m m ni i n n
i i
Z F Z F
F F F Z Zα α
= =
= =
⎛ ⎞− − − −⎜ ⎟
⎝ ⎠
⎛ ⎞− < − < − ⇒ <⎜ ⎟
⎝ ⎠
∏ ∏
∏ ∏
%
%
①
②
③
④
![Page 13: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/13.jpg)
08 Numerical example ( ) (0,1)F x N=( m=10,000, n=5,000, iteration=100 )
( )Q α%
( )Q α
ˆ ( )Q α
~ (0, ), ~ (1,0.05)i i iX N Nσ σ
~ (0, ), ~ (1,0.1)i i iX N Nσ σ ~ (0, ), ~ (1,0.2)i i iX N Nσ σ
~ (0,1)iX N
![Page 14: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/14.jpg)
09 practical application:Annual maximum wind speed in Japan
1) Observation records and Historical annual maximum wind speeds at 155 sites in JapanObservation Records: JMA records (CSV format)
1961~1990 : 10 min average wind speed per 3 hours1991~2002 : 10 min average wind speed per hour
Historical annual maximum wind speeds record:1929~1999 : A historical annual max. wind speed data set
compiled by Ishihara et al.(2002)2000~2002 : extracted from JMA records (CSV format)
T. Ishihara et al. (2002), A database of annual maximum wind speed and corrections for anemometers in Japan,Wind Engineers, JAWE, No.92, p5~54 (in Japanese)
![Page 15: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/15.jpg)
09 practical application:Annual maximum wind speed in Japan
2) The effect of different observation recording format on the basic statistics
Base on the recently opened continuous 1 min average wind speed records (1997. 3~2002. 2), calculating every 10 min average wind speeds, 10 min average per hour and 3 hours, and comparing the basic statistics for each recording format, then the effect of different recording format becomes to be clear.
![Page 16: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/16.jpg)
µ σ
1 1γ β= 2 2 3γ β= −
![Page 17: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/17.jpg)
3) Examples of non-stationarity : the basic statistics (1961~2002)
Coefficient of Variations (C.O.V)Site
Abashiri 5.1% 5.0% 11.7% 14.4%
Katsuura 5.81% 8.87% 24.06% 26.71%
Kobe 5.71% 8.25% 12.55% 16.31%
Kumamoto 7.62% 5.93% 14.19% 32.36%
Makurazaki 4.60% 6.07% 25.97% 47.98%
Morioka 6.99% 4.43% 12.09% 11.55%
Oita 3.42% 8.41% 16.49% 27.41%
Shionomisaki 4.34% 4.07% 19.36% 30.10%
Tokyo 5.90% 8.06% 18.90% 17.88%
Minimum 3.42% 4.07% 11.72% 11.55%
Maximum 7.62% 8.87% 25.97% 47.98%
µ σ 1γ 2γ
![Page 18: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/18.jpg)
Morioka
Kobe
Shionomisaki
Makurazaki
Naha
![Page 19: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/19.jpg)
Morioka (1961~2002)
µ σ
1 1γ β= 2 2 3γ β= −
![Page 20: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/20.jpg)
Kobe (1961~2002)
![Page 21: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/21.jpg)
Shionomisaki (1961~2002)
![Page 22: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/22.jpg)
Makurazaki (1961~2002)
![Page 23: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/23.jpg)
Naha (1964~2002)
![Page 24: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/24.jpg)
09 practical application:Annual maximum wind speed in Japan
4) Approximation of the annual wind speed distributionBased on the probability integral transformation,
( )
( )2 3
1
( ) ( ) ( )
( )
( ) ( )
i i
i
X X
X
z F x F g z
x g z a bz cz dz
F x g x
α α α
α α α α α
α α
α
−
Φ = = =
= = + + +
⇒ = Φ
The coefficients a,b,c and d can be estimated from the givenbasic statistics of annual wind speed, i.e. mean, standard deviation, skewness and kurtosis (Choi & Kanda 2003).
H. Choi and J. Kanda (2003), Translation Method: a historical review and its application to simulation on non-Gaussianstationary processes, Wind and Struct., 6(5), p357~386
![Page 25: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/25.jpg)
5) Estimation by Monte Carlo Simulation (MCS)
5.1) Simulation methods
① Based on the Spectral representation theorem for stationary stochastic process→Using a given spectral density function, discrete
stationary stochastic process is simulated.→time consuming method
② Based on equivalent i.i.d. random sequence (EQRS)→A stochastic process, which can be approximated
by Poisson process, is modeled as an i.i.d. random sequence having same quantile function.→time effective method
![Page 26: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/26.jpg)
5.2) Required information for MCS based on EQRS
① m : the number of Independent rv→ approximated by mean zero crossing rate(Normal process)From Rice theorem and Poisson approximation, normal quantile function is given as follows:
{ }02 log Tz T yα µ= +
0 : mean zero crossinin whcih g rate, log( log )Tyµ α= − −
From ( )m zαΦ ,
( ){ }2 log / 4 Tz m yα π +;
By comparison,04m Tπµ;
* Choi & Kanda (2004), A new method of the extreme value distributions based on the translation method, Summaries of Tech. Papers of Ann. Meet. of AIJ, Vol. B1, p23~24 (in Japanese)
*
![Page 27: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/27.jpg)
{ }2 log(1/ 4 )z yπ= +
![Page 28: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/28.jpg)
{ }2 log( / 4 ) 0.57722z mµ π +;
![Page 29: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/29.jpg)
(non-Normal process)
To Rice formula for the expected number of crossings, i.e.
0 0( ) ( | ) ( ) ( )x x f x X x dx f x x f x dxν
∞ ∞= ⋅ = = ⋅ ⋅∫ ∫& & & & & &
( )( ){ }2
1
0
1 2 1 2
0
exp / ( ) / 21( ) ( )( ) 2 ( )
{ ( )} { ( )}exp exp2 2 2
z
z
z
x x g zx g x dx
g z g z
g x g x
σν φ
π σ
σ νπ
∞−
− −
′⋅ − ⋅= ⋅
′ ′⋅ ⋅
⎛ ⎞ ⎛ ⎞= ⋅ − = ⋅ −⎜ ⎟ ⎜ ⎟
⎝ ⎠ ⎝ ⎠
∫&
&
&
& &&
applying translation function g(z)
From Poisson approximation
{ }( )0( ) 2 log( ) Tx g z g T yα α ν= = +
With the same manner04m Tπν;
* The distribution of dx/dt is assumed as normal distribution and the assumption is reasonable.e.g. H. Choi (1988), Characteristics of natural wind for wind load estimation, Master Thesis, Univ. of Tokyo (in Japanese)
*
![Page 30: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/30.jpg)
Practical example (T=1 year)
• estimated from long term observation records in Tokyo (1985~1987, Choi 1988)
0Tν
2384311525863124
187937963Height(m)
23072739254526652883
5848464545Height(m)
0Tν
0Tν
04 4 (2300 ~ 3000) 30,000m Tπν π= = ⋅ →
![Page 31: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/31.jpg)
5.2) Required information for MCS based on EQRS
② parent distribution function for each year→Generalized bootstrap method+
Translation method (Probability Integral Transform)
For the year i,
{ }
2 3
1 2
1 1, , 1 1, ,
( ( )) ( ), ( ), , , ( , , , )
max( , , ) (max( , , ) )
i i i i i i i
i i i ii
m i n i m i n
F x g z z g z a b z c z d za b c d
x x g z z
µ σ γ γ
= =
= = Φ = + + +
= Ψ
=K KK K
![Page 32: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/32.jpg)
Example : Tokyo (1961~2002)
µ σ
1γ 2γ
Estimated from historical records Monte Carlo Simulation
![Page 33: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/33.jpg)
Correlation between the basic statistics (Tokyo)
1.0000.9460.433-0.3150.9461.0000.541-0.3720.4330.5411.0000.242-0.315-0.3720.2421.000
µ σ 1γ 2γµσ
1γ2γ
Such correlation characteristics between the basic statistics are regeneratedby Cholesky decomposition of correlation matrix.
![Page 34: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/34.jpg)
Regeneration of correlation characteristics
correlation coefficients of simulated ones and given values
1.0000.958(0.946)
0.476(0.433)
-0.343(-0.315)
0.958(0.946)
1.0000.565(0.541)
-0.389(-0.372)
0.476(0.433)
0.565(0.541)
1.0000.242(0.242)
-0.343(-0.315)
-0.389(-0.372)
0.242(0.242)
1.000
(・): given correlation coefficient
µ σ 1γ 2γ
µ
σ
1γ
2γ
![Page 35: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/35.jpg)
No. of Simulation:n=100 year x 1000 times
○:annual mean and standard dev. from historical records
![Page 36: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/36.jpg)
No. of Simulation:n=100 years x 1000 times
○:skewness and kurtosis from historical records
![Page 37: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/37.jpg)
Aomori
Sendai
TokyoKobe
Shionomisaki
Makurazaki
Naha
![Page 38: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/38.jpg)
6) Comparison of the quantile functions from MCS and historical records in normalized form
, ,( ) /n m n m nZ a b−
Aomori95% confidence intervals
Mean of 1000 samples
○:before 1961 ●:after 1962
(Z-am,n)/bm,n
![Page 39: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/39.jpg)
Sendai
![Page 40: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/40.jpg)
Tokyo
![Page 41: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/41.jpg)
Kobe
![Page 42: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/42.jpg)
Shionomisaki
![Page 43: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/43.jpg)
Makurazaki
![Page 44: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/44.jpg)
Naha
![Page 45: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/45.jpg)
7) Comparison of the quantile functions from MCS and historical records in full scale
Aomori
![Page 46: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/46.jpg)
Sendai
![Page 47: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/47.jpg)
Tokyo
![Page 48: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/48.jpg)
Kobe
![Page 49: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/49.jpg)
Shionomisaki
![Page 50: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/50.jpg)
Makurazaki
![Page 51: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/51.jpg)
Naha
![Page 52: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/52.jpg)
8) Comparison of the attraction coefficients from MCS and the historical records
○:n>50 (136 sites), ■:n≦50 (19 sites)
,m na ,m nb
, ,lim , limm n m m n mn na a b b
→∞ →∞= =
from
his
toric
al d
ata
from
his
toric
al d
ata
from Monte Carlo Simulation from Monte Carlo Simulation
![Page 53: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/53.jpg)
9) Defect of the alternative definition ˆ ( ) ( )mZF x F x=
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Aomori
![Page 54: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/54.jpg)
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Sendai
![Page 55: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/55.jpg)
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Tokyo
![Page 56: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/56.jpg)
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Kobe
![Page 57: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/57.jpg)
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Shionomisaki
![Page 58: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/58.jpg)
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Makurazaki
![Page 59: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/59.jpg)
Mean of 1000 MCS
ˆ ( ) ( )mZF x F x=
Naha
![Page 60: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/60.jpg)
10) Comparison of the attraction coefficients from the alternative definition and the
historical records
,m na ,m nbfrom
his
toric
al re
cord
s
from
his
toric
al re
cord
s
from alternative definition from alternative definition
![Page 61: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/61.jpg)
10 The cult of isolated statistics and The law of large number
How many extreme values should be used to estimate an extreme value distribution?
( case study for max/min pressure coefficients)*
* Choi & Kanda (2004), Stability of extreme quantile function estimation from relatively short records having different parent distributions, Proc. 18th Natl. Symp. Wind Engr., p455~460 (in Japanese)
![Page 62: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/62.jpg)
B/5, B/5, B/10
H=400mm
D=100mm
B=100mmB/2
H/10H/10
H/10H/10
H/20
H/20
B=100mm
D=100mmD/2
D/10, D/5, D/5, D/5, D/5, D/10
1 2 3 4 5 6 7 8
wind
910
wind
Windward face Side and Leeward faces
Tap #11 Tap #15 Tap #17
![Page 63: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/63.jpg)
Which one is the best estimation?
2
( )( )1/ 2
sp
H
p t pC tUρ−
= ( ) : instant total pressure, : static pressure: air density, : wind velocity at model height
s
H
p t pUρ
* 50 blocks contain about 480,000 discrete data
![Page 64: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/64.jpg)
Which one is the best estimation?
![Page 65: An approach to the extreme value distribution of non ... · * Falk et al., Laws of Small Numbers: Extremes and Rare Events, Birkhäuser, 1994 The class of EVD for non-i.i.d. case](https://reader036.vdocument.in/reader036/viewer/2022071006/5fc3a962421a664ebc2d8dde/html5/thumbnails/65.jpg)
10 The cult of isolated statistics and The law of large number
We never be free from the law of large number.
The statistician and the scientists/technologist need to understand that models are necessarily simplifications of the system being modelled; that they are , an absolute sense, wrong; that they are certainly provisional, but nonetheless are useful and necessary for successful quantitative thinking.
from J.A. Nelder (1986), Statistics, Science and Technology – The address of the president, delivered to theRoyal Statistical Society on Wednesday, April 16th, 1986, J. R. Statist. Soc. A 149(2), p109~121