(m/s)
|
|
- こうじ しのしま
- 5 years ago
- Views:
Transcription
1 ( ) r-taka@maritime.kobe-u.ac.jp IBIS2009
2 (m/s)
3 Daily returns of the S&P 500 index Gilli & Këllezi (2006). Computational Econometrics 27(1).
4 Daily returns of the S&P 500 index
5
6 (GEV) (GP) GEV GP
7
8 block maxima) (threshold exceedances or Peaks Over Threshold)
9
10 F X F X 1, X 2,..., X n Z n = max { X 1, X 2,..., X n } = max 1 i n X i n Z n F x F Z n x F = sup{x : F (x) < 1}, n.
11 Z n F a n > 0, b n R (n = 1, 2,...) G(x) Z P ( Zn b n Z n b n a n L Z, n. a n ) x P (Z x) = G(x) G (extreme value distribution) F G (domain of attraction) F D(G) (a n, b n )
12 (1) G(x) (2) F D(G) W n = min 1 i n X i W n = min { } { } X 1, X 2,..., X n = max X1, X 2,..., X n W G W (x) = 1 G( x)
13 ( ) ( ) Zn b n P x = P (Z n a n x + b n ) = P max X i a n x + b n 1 i n a n = P (X i a n x + b n, i = 1, 2,..., n) = n i=1 P (X i a n x + b n ) = = F n (a n x + b n ). n i=1 F (a n x + b n )
14 F f(x) Exp(1) F (x) = 1 e x, f(x) = e x, x 0. F n (x + log n) = { { } n 1 e (x+log n)} n = 1 + e x n a n = 1 b n = log n e e x = exp( exp( x)), n. ( < x < ).
15 Pareto ( α > 0) F (x) = 1 1 x α, f(x) = αx α 1, x 1. F n (n 1/α x) = { 1 1 (n 1/α x) α } n = { 1 + x α n } n a n = n 1/α b n = 0 e x α = exp( x α ), n. (x > 0).
16 ( α > 0) F (x) = 1 (1 x) α, f(x) = α(1 x) α 1, 0 x 1. F n (n 1/α x + 1) = { { } n 1 ( n 1/α x) α} n = 1 + ( x)α n e ( x)α = exp( ( x) α ), n. (x 0). a n = n 1/α b n = 1 = x F
17 [Fréchet (1927) Fisher and Tippett (1928) Gnedenko (1943); Trinity Theorem ] G(x) Gumbel Λ(x) = exp( exp( x)), x R, Fréchet Φ α (x) = ( )Weibull Ψ α (x) = { 0, x 0, exp( x α ), x > 0, { exp( ( x) α ), x 0, 1, x > 0, α > 0, α > 0.
18 [Gnedenko (1943) de Haan (1970)] F D(Φ α ) x F = & lim x 1 F (tx) 1 F (x) = t α, t > 0. F D(Ψ α ) x F < & lim x xf 1 F (x F (x F x)t) 1 F (x) = t α, t > 0. 1 F (x + ts(x)) F D(Λ) s( ) > 0 s.t. ( ) lim x xf 1 F (x) ( ) = xf x = s(x) = (1 F (y))dy <, x < x F xf x = e t. (1 F (y))dy/(1 F (x)) satisfies ( ).
19 von Mises-Jenkinson ( ) G ξ (x) = exp { (1 + ξx) 1/ξ}, 1 + ξx > 0. (generalized extreme value, GEV) ξ R ξ = 0 G 0 (x) = lim ξ 0 G ξ (x) = exp{ exp( x)} = Λ(x) Φ α (x) = G 1/α (α(x 1)), Ψ α (x) = G 1/α (α(x + 1)).
20 (max stable) n A n > 0 B n R G n ξ (A nx + B n ) = G ξ (x) 3
21 F n (a n x + b n ) G ξ (x). n a n x + b n = z P (Z n z) = F n (z) G ξ ( z bn a n ). G ξ n
22 F n (a n x + b n ) = [ 1 { 1 F (a n x + b n ) }] n = [ 1 n { 1 F (a n x + b n ) }/ n ] n exp [ lim n n { 1 F (a n x + b n ) }] lim n n{1 F (a nx + b n )} = log G ξ (x) = (1 + ξx) 1/ξ. x = 0 lim n n{1 F (b n)} = 1. 2 lim n n{1 F (a n x + b n )} n{1 F (b n )} = (1 + ξx) 1/ξ.
23 P (X > a n x + b n X > b n ) = 1 F (a nx + b n ) 1 F (b n ) (1 + ξx) 1/ξ. P (X b n a n x X > b n ) 1 (1 + ξx) 1/ξ. H ξ (x) = 1 + log G ξ (x) = 1 (1 + ξx) 1/ξ (Generalized Pareto, GP) ( y b n P (X b n y X > b n ) H ξ a n ). H ξ b n
24 F F u F u (y) = P (X u y X > u), 0 y x F u [Pickands (1975)] F D(G ξ ) lim u x F F u (a(u)y) = H ξ (y), y 0, F u (a(u)y) < 1. a( ) a n n
25 G ξ H ξ ξ ξ < 0 ξ = 0 ξ > 0 Pareto t Cauchy D(G ξ ), / D(G ξ ). F (x) = 1 1/ log x, x e = F / D(G ξ ).
26 GEV (generalized extreme value) GEV(µ, σ, ξ) < µ < σ > 0 < ξ < { [ ( )] } z µ 1/ξ ( ) z µ G(z) = exp 1 + ξ = G ξ, 1 + ξ(z µ)/σ > 0. σ σ G ξ GEV(0, 1, ξ) G ξ (z) = exp [ (1 + ξz) 1/ξ], 1 + ξz > 0, µ σ ξ
27 G(z) ξ < 0 Weibull z < µ σ/ξ ξ = 0 Gumbel < z < G 0 ((z µ)/σ) = lim ξ 0 G ξ ((z µ)/σ) = exp{ exp[ (z µ)/σ]} ξ > 0 Fréchet z > µ σ/ξ G ξ (z) g ξ (z) = { } (1 + ξ z) 1/ξ 1 exp (1 + ξ z) 1/ξ, 1 + ξz > 0, ξ 0, exp { z exp( z) }, < z <, ξ = 0,
28 GEV( 2.5, 1, 0.4) 0 GEV(0, 1, 0) GEV(2.5, 1, 0.4) 0
29 GEV) {z 1, z 2,..., z n } GEV(µ, σ, ξ) { [ ( )] } z µ 1/ξ ( ) z µ G(z) = exp 1 + ξ = G ξ, 1 + ξ(z µ)/σ > 0, σ σ
30 GEV(µ, σ, ξ) ξ 0 n [ ( )] zi µ l(µ, σ, ξ) = n log σ (1 + 1/ξ) log 1 + ξ σ i=1 n [ ( )] 1/ξ zi µ 1 + ξ, σ i=1 1 + ξ(z i µ)/σ > 0, i = 1,..., n. ξ = 0 l(µ, σ) = n log σ n i=1 ( zi µ σ ) n i=1 exp { ( zi µ σ )}. ( µ, σ, ξ) ( µ, σ)
31 GEV(µ, σ, ξ) (Prescott and Walden, 1980) n σ 2 ξ 2 ξ 2 p ξ { 2 Γ(2 + ξ) p } σξ(p/ξ q) 1 2 Γ(2 + ξ) + p σ σ 2 [ Γ(2 + ξ) 1 π2 6 + ξ ( 1 γ + 1 ξ + q p ξ 1 + γ ] ) 2 2q ξ + p ξ 2 (µ, σ, ξ) Γ( ) Ψ(r) = d log Γ(r)/dr p = (1 + ξ) 2 Γ(1 + 2ξ), γ = Euler q = Γ(2 + ξ){ψ(1 + ξ) + (1 + ξ)/ξ},
32 {GEV(µ, σ, ξ), < µ <, σ > 0, < ξ < } ξ > 0.5 Smith, 1985 ξ < ξ < 0.5 (µ, σ, ξ)
33 ( ) GEV(µ, σ, ξ) 1 1/T R T ( ) RT µ G(R T ) = G ξ = 1 1 σ T µ + σ [{ log(1 1/T ) } ξ ]/ 1 ξ, ξ 0, R T = µ + σ [ log { log(1 1/T ) }], ξ = 0. R T (return period) T (return level) T = 200 R n T
34 R T R T = µ + σ { y ξ T 1}/ ξ, ξ 0, µ + σ { log y T }, ξ = 0, ( µ, σ, ξ) y T = log(1 1/T ) V ( R T ) RT T V R T V ( µ, σ, ξ) [ RT T RT = µ, R T σ, R ] T ξ [ = 1, (y ξ T 1)/ξ, σy ξ T ( log y T )/ξ σ(y ξ ] T 1)/ξ2 ( µ, σ, ξ)
35 ξ < 0 R = µ σ/ξ R = µ σ/ ξ R T = [ 1, 1/ξ, σ/ξ 2 ]
36 ξ R T ξ ξ = ξ 0 l(µ, σ, ξ 0 ) µ σ ξ 95% { ξ : 2 { l( µ, σ, ξ) } max µ, σ l(µ, σ, ξ)} χ 2 1 (0.05) = { = ξ : max µ, σ l(µ, σ, ξ) l( µ, σ, ξ) } 1.921
37 R T µ = R T σ [ y ξ T 1] /ξ (µ, σ, ξ) (R T, σ, ξ) n [ )] n [ )] 1/ξ l(r T, σ, ξ) = n log σ (1+1/ξ) i=1 log y ξ T + ξ ( zi R T σ i=1 y ξ T + ξ ( zi R T σ R T 95% { R T : 2 { l( R T, σ, ξ) max l(r T, σ, ξ) } } χ 2 1 (0.05) σ, ξ { = R T : max l(r T, σ, ξ) l( R T, σ, ξ) } χ 2 1 (0.05)/2 σ, ξ
38 (m/s) (1982) (2002) GEV (2004)
39 (m/s) (m/s)
40 GEV(µ, σ, ξ) ( µ, σ, ξ) = (15.349, 2.550, 0.111) V = ( µ, σ, ξ) PP plot QQ plot
41 Probability Plot Quantile Plot Empirical Empirical Model Model
42 Return Level Plot Density Plot Return Level f(z) Return Period z
43 ξ = > 0 Fréchet ξ 95% ξ ± = [ 0.037, 0.259]. 95% [ 0.021, 0.274]
44 Profile Log-likelihood Shape Parameter ξ ξ = > 0 95% [ 0.021, 0.274]
45 200 R % R 200 R 200 ± 1.96 R T 200 V T 200 = [25.75, 41.73]. 95% [28.29, 46.90]
46 Profile Log-likelihood Return Level 200 R 200 = % [28.29, 46.90]
47 GP (generalized Pareto) GP(σ, ξ) σ > 0 < ξ < H(y) = 1 ( 1 + ξ y ) 1/ξ ( ) y = Hξ, 1 + ξy/σ > 0. σ σ H ξ GP(1, ξ) H ξ (y) = 1 (1 + ξy) 1/ξ, 1 + ξy > 0, σ ξ
48 H(y) ξ < 0 0 < y < σ/ξ ξ = 0 0 < y < H 0 (y/σ) = lim ξ 0 H ξ (y/σ) = 1 e y/σ ξ > 0 Pareto 0 < y < H ξ (y) (1 + ξ y) 1/ξ 1, 1 + ξy > 0, ξ 0, h ξ (y) = exp( y), 0 < y <, ξ = 0,
49 GP(1, ξ) ξ = 0.4, 0, 0.4
50 (GP) {y 1, y 2,..., y n } GP(σ, ξ) H(y) = 1 ( 1 + ξ y ) 1/ξ ( ) y = Hξ, 1 + ξy/σ > 0, σ σ
51 GP(σ, ξ) ξ 0 n l(σ, ξ) = n log σ (1 + 1/ξ) log(1 + ξ y i /σ), i=1 1 + ξ y i /σ > 0, i = 1, 2,..., n. ξ = 0 l(σ) = n log σ 1 σ n i=1 y i. ( σ, ξ) σ
52 GP(σ, ξ) Fisher n [ (1 + ξ)/σ 2 1/σ ] (1 + ξ)(1 + 2ξ) 1/σ 2 ξ > 1/2 V n = 1 [ 2σ 2 (1 + ξ) σ(1 + ξ) n σ(1 + ξ) (1 + ξ) 2 Smith, 1985 ]
53 Y GP(σ, ξ) ξ < 1 E(Y ) = y+ y+ (1 H(y))dy = 0 0 y + = sup{y : H(y) < 1} Y v Y > v GP(σ + ξv, ξ) (v > 0) P (Y v > y Y > v) = = 1 H(y + v) ( 1 H(v) y 1 + ξ σ + ξv ( 1 + ξ y ) 1/ξ σ dy = σ 1 ξ. = ( 1 + ξ(y + v)/σ ) 1/ξ ) 1/ξ. ( 1 + ξv/σ ) 1/ξ
54 e(v) = E(Y v Y > v) Y (mean excess) ẽ(v) Y (median excess) P (Y v ẽ(v) Y > v) = 1/2 Y v Y > v GP(σ + ξv, ξ) e(v) = σ + ξv 1 ξ ẽ(v) = = σ 1 ξ + ξ 1 ξ v, σ(2 ξ 1)/ξ + (2 ξ 1)v, ξ 0, σ log 2, ξ = 0.
55 u u
56 (mean excess plot) u x [1], x [2],..., x [Nu ] x max u < x max u, 1 N u N u i=1 (x [i] u) : u < x max u u Y GP(σ, ξ) = E(Y v Y > v) = σ 1 ξ + ξ 1 ξ v.
57 ξ σ u < x max {x [i] u} N u i=1 GP(σ u, ξ) σ u ξ u (u, σ ) (u, ξ) σ = σ u ξu u σ ξ u Y v Y > v GP(σ + ξv, ξ) (v > 0), σ v = σ + ξv.
58 F F (VaR p ) = 1 p (Value at Risk) F (x) F (x) = (1 F (u))f u (y) + F (u), y = x u. u F u GP H ξ ( ) x u F (x) (1 F (u))h ξ + F (u), x u. σ
59 = ζ u = 1 F (u) VaR p = u + σ ( ζu ξ p ) ξ 1, u ( σ, ξ) ζ u N u /n n N u u VaR p = u + σ ( ) ξ ζ u 1 ξ p
60 VaR p ( ζ u, σ, ξ) [ ζ V = u (1 ζ u )/n 0 T ] 0 V Nu 0 T = (0, 0) VaR p V ( VaR p ) VaR T p V VaR p, VaR T p = [ VaRp ζ u, VaR p σ, VaR ] p ξ ( ζ u, σ, ξ)
61 GEV ξ 95% {ξ : max σ l(σ, ξ) l( σ, ξ) 1.921} VaR p σ = (VaR p u)ξ (ζ u /p) ξ 1 l(var p, ξ) VaR p max ξ l(var p, ξ) ζ u
62 M. Gilli & E. Këllezi An Application of Extreme Value Theory for Measuring Financial Risk. S&P 500
63 Daily returns of the S&P 500 index
64 Mean Excess u
65 Mean Excess u
66 Modified Scale Threshold Shape Threshold σ ξ
67 u = σ = (0.0359), ξ = (0.0474)
68 Probability Plot Quantile Plot Empirical Empirical Model Model S&P 500
69 Profile Log-likelihood Shape Parameter ξ 95% [0.052, 0.238]. ξ = [0.044, 0.230]. )
70 Profile Log-likelihood Return Level VaR % [2.411, 2.610].
71 Stephenson, A. and Gilleland, E. (2006). Software for the analysis of extreme events: The current state and future directions. Extremes 8, Extremes 1998 (2004) Vol. 52, No. 1 ( )
72 Beirlant, J., Goegebeur, Y., Segers, J. and Teugels, J. (2004). Statistics of Extremes, Theory and Applications. Wiley. Coles, S. G. (2001). An Introduction to Statistical Modeling of Extreme Values. Springer. Embrechts, P., Klüppelberg, C. and Mikosch, T. (2001). Modelling Extremal Events for Insurance and Finance, 3rd ed. Springer. Falk, M., Hüsler, J. and Reiss, R.-D. (2004). Laws of Small Numbers: Extremes and Rare Events, 2nd ed. Birkhäuser.
73 de Haan, L. and Ferreira, A. (2006). Extreme Value Theory: An Introduction. Springer. Kotz, S. and Nadarajah, S. (2000). Extreme Value Distributions: Theory and Applications. Imperial College Press. Leadbetter, M. R., Lindgren, G. and Rootzén, H. (1983). Extremes and Related Properties of Random Sequences and Processes. Springer-Verlag, New York. McNeil, J., Frey, R. and Embrechts, P. (2005). Quantitative Risk Management. Princeton University Press.
74 Reiss, R.-D. and Thomas, M. (2007). Statistical Analysis of Extreme Values, with Applications to Insurance, Finance, Hydrology and Other Fields, 3rd ed. Birkhäuser. Resnick, S. I. (1987). Extreme Values, Regular Variation and Point Processes. Springer-Verlag, Berlin. Resnick, S. I. (2007). Hervy-Tail Phenomena. Springer.
75 (2004) (2008) II. (2008). 21 II (2008) (5)
76 (2004). r Fisher, R. A. and Tippett, L. H. C. (1928). Limiting forms of the frequency distribution of the largest or smallest member of a sample. Proc. Cambridge Philos. Soc. 24, Fréchet, M. (1927). Sur la loi de probabilité de l écart maximum. Ann. Soc. Math. Polon. 6, Gnedenko, B. (1943). Sur la distibution limite du terme maximum d une serie aleatoire. Ann, Math. 44, Translated and
77 reprinted in: Breakthroughs in Statistics, Vol.I, 1992, eds. S. Kotz and N. L. Johnson, Springer-verlag, pp Gumbel, E. J. (1958). Statistics of Extremes. Columbia University Press. de Haan, L. (1970). On Regular Variation and Its Application to the Weak Convergence of Sample Extremes. Math. Centre Tracts, Vol.32, Amsterdam. de Haan, L. (1976). Sample extremes: an elementary introduction. Statist. Neerlandica 30,
78 Jenkinson, A. F. (1955). The frequency distribution of the annual maximum (or minimum) values of meteorological elements. Quart. J. R. Met. Soc. 81, Mises, R. von (1936). La distribution de la plus grande de n valeurs. Rev. Math. Union Interbalcanique 1, Reproduced in Selected Papers of Richard von Mises, Amer. Math. Soc. 2 (1964), pp Pickands, J. (1975). Statistical inference using extreme order statistics. Ann. Statist. 3,
79 Prescott, P. and Walden, A. T. (1980). Maximum likelihood estimation of the parameters of the generalized extreme value distribution. Biometrika 67, Smith, R. L. (1985). Maximum likelihood estimation in a class of nonregular cases. Biometrika 72, Smith, R. L. (1986). Extreme value theory based on the r largest annual events. Journal of Hydrology 86, Smith, R. L. (1987). Estimating tails of probability distributions. Ann. Statist. 15,
80 Smith, R. L. (1989). Extreme value analysis of environmental time series: an example based on ozone data (with discussion). Statistical Science 4, Smith, R. L. (1990). Extreme value theory. In Handbook of Applicable Mathematics VII, Ed. W. Ledermann, Chichester, Wiley. pp Weissman, I. (1978). Estimation of parameters and large quantiles based on the k largest observations. J. Amer. Statist. Assoc. 73,
1 Tokyo Daily Rainfall (mm) Days (mm)
( ) r-taka@maritime.kobe-u.ac.jp 1 Tokyo Daily Rainfall (mm) 0 100 200 300 0 10000 20000 30000 40000 50000 Days (mm) 1876 1 1 2013 12 31 Tokyo, 1876 Daily Rainfall (mm) 0 50 100 150 0 100 200 300 Tokyo,
More informationEinmahl and Magnus (2007) Records in athretics through extreme-value theory. 04:35
( ) r-taka@maritime.kobe-u.ac.jp Einmahl and Magnus (2007) Records in athretics through extreme-value theory. 04:35 0 2 4 6 8 10 12 14 102 104 106 108 Age Castillo, Hadi, Balakrishnan and Sarabia (2004)
More informationohpmain.dvi
fujisawa@ism.ac.jp 1 Contents 1. 2. 3. 4. γ- 2 1. 3 10 5.6, 5.7, 5.4, 5.5, 5.8, 5.5, 5.3, 5.6, 5.4, 5.2. 5.5 5.6 +5.7 +5.4 +5.5 +5.8 +5.5 +5.3 +5.6 +5.4 +5.2 =5.5. 10 outlier 5 5.6, 5.7, 5.4, 5.5, 5.8,
More informationuntitled
18 1 2,000,000 2,000,000 2007 2 2 2008 3 31 (1) 6 JCOSSAR 2007pp.57-642007.6. LCC (1) (2) 2 10mm 1020 14 12 10 8 6 4 40,50,60 2 0 1998 27.5 1995 1960 40 1) 2) 3) LCC LCC LCC 1 1) Vol.42No.5pp.29-322004.5.
More information03.Œk’ì
HRS KG NG-HRS NG-KG AIC Fama 1965 Mandelbrot Blattberg Gonedes t t Kariya, et. al. Nagahara ARCH EngleGARCH Bollerslev EGARCH Nelson GARCH Heynen, et. al. r n r n =σ n w n logσ n =α +βlogσ n 1 + v n w
More information医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.
医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009192 このサンプルページの内容は, 第 2 版 1 刷発行時のものです. i 2 t 1. 2. 3 2 3. 6 4. 7 5. n 2 ν 6. 2 7. 2003 ii 2 2013 10 iii 1987
More information2 1,2, , 2 ( ) (1) (2) (3) (4) Cameron and Trivedi(1998) , (1987) (1982) Agresti(2003)
3 1 1 1 2 1 2 1,2,3 1 0 50 3000, 2 ( ) 1 3 1 0 4 3 (1) (2) (3) (4) 1 1 1 2 3 Cameron and Trivedi(1998) 4 1974, (1987) (1982) Agresti(2003) 3 (1)-(4) AAA, AA+,A (1) (2) (3) (4) (5) (1)-(5) 1 2 5 3 5 (DI)
More informationKullback-Leibler
Kullback-Leibler 206 6 6 http://www.math.tohoku.ac.jp/~kuroki/latex/206066kullbackleibler.pdf 0 2 Kullback-Leibler 3. q i.......................... 3.2........... 3.3 Kullback-Leibler.............. 4.4
More information2 Recovery Theorem Spears [2013]Audrino et al. [2015]Backwell [2015] Spears [2013] Ross [2015] Audrino et al. [2015] Recovery Theorem Tikhonov (Tikhon
Recovery Theorem Forward Looking Recovery Theorem Ross [2015] forward looking Audrino et al. [2015] Tikhonov Tikhonov 1. Tikhonov 2. Tikhonov 3. 3 1 forward looking *1 Recovery Theorem Ross [2015] forward
More informationuntitled
2 : n =1, 2,, 10000 0.5125 0.51 0.5075 0.505 0.5025 0.5 0.4975 0.495 0 2000 4000 6000 8000 10000 2 weak law of large numbers 1. X 1,X 2,,X n 2. µ = E(X i ),i=1, 2,,n 3. σi 2 = V (X i ) σ 2,i=1, 2,,n ɛ>0
More informationseminar0220a.dvi
1 Hi-Stat 2 16 2 20 16:30-18:00 2 2 217 1 COE 4 COE RA E-MAIL: ged0104@srv.cc.hit-u.ac.jp 2004 2 25 S-PLUS S-PLUS S-PLUS S-code 2 [8] [8] [8] 1 2 ARFIMA(p, d, q) FI(d) φ(l)(1 L) d x t = θ(l)ε t ({ε t }
More informationバリュー・アット・リスクのリスク指標としての妥当性について ― 理論的サーベイによる期待ショートフォールとの比較分析―
aaaab aabab VaR VaRArtzner et al. VaR VaR VaR Artzner et al.var VaR VaR VaR ρ XY ρ (X+Y ) ρ(x) + ρ(y ) XY ρ VaRArtzner et al.1999basak and Shapiro1999Danielsson2000Rootzén and Klüppelberg VaR VaR VaRVaR
More informationL. S. Abstract. Date: last revised on 9 Feb translated to Japanese by Kazumoto Iguchi. Original papers: Received May 13, L. Onsager and S.
L. S. Abstract. Date: last revised on 9 Feb 01. translated to Japanese by Kazumoto Iguchi. Original papers: Received May 13, 1953. L. Onsager and S. Machlup, Fluctuations and Irreversibel Processes, Physical
More informationxia2.dvi
Journal of Differential Equations 96 (992), 70-84 Melnikov method and transversal homoclinic points in the restricted three-body problem Zhihong Xia Department of Mathematics, Harvard University Cambridge,
More informationCVaR
CVaR 20 4 24 3 24 1 31 ,.,.,. Markowitz,., (Value-at-Risk, VaR) (Conditional Value-at-Risk, CVaR). VaR, CVaR VaR. CVaR, CVaR. CVaR,,.,.,,,.,,. 1 5 2 VaR CVaR 6 2.1................................................
More information1
1 Borel1956 Groupes linéaire algébriques, Ann. of Math. 64 (1956), 20 82. Chevalley1956/58 Sur la classification des groupes de Lie algébriques, Sém. Chevalley 1956/58, E.N.S., Paris. Tits1959 Sur la classification
More informationwaseda2010a-jukaiki1-main.dvi
November, 2 Contents 6 2 8 3 3 3 32 32 33 5 34 34 6 35 35 7 4 R 2 7 4 4 9 42 42 2 43 44 2 5 : 2 5 5 23 52 52 23 53 53 23 54 24 6 24 6 6 26 62 62 26 63 t 27 7 27 7 7 28 72 72 28 73 36) 29 8 29 8 29 82 3
More informationkato-kuriki-2012-jjas-41-1.pdf
Vol. 41, No. 1 (2012), 1 14 2 / JST CREST T 2 T 2 2 K K K K 2,,,,,. 1. t i y i 2 2 y i = f (t i ; c) + ε i, f (t; c) = c h t h = c ψ(t), i = 1,...,N (1) h=0 c = (c 0, c 1, c 2 ), ψ(t) = (1, t, t 2 ) 3
More information(2) Fisher α (α) α Fisher α ( α) 0 Levi Civita (1) ( 1) e m (e) (m) ([1], [2], [13]) Poincaré e m Poincaré e m Kähler-like 2 Kähler-like
() 10 9 30 1 Fisher α (α) α Fisher α ( α) 0 Levi Civita (1) ( 1) e m (e) (m) ([1], [], [13]) Poincaré e m Poincaré e m Kähler-like Kähler-like Kähler M g M X, Y, Z (.1) Xg(Y, Z) = g( X Y, Z) + g(y, XZ)
More information<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>
確率的手法による構造安全性の解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/55271 このサンプルページの内容は, 初版 1 刷発行当時のものです. i 25 7 ii Benjamin &Cornell Ang & Tang Schuëller 1973 1974 Ang Mathematica
More informationカルマンフィルターによるベータ推定( )
β TOPIX 1 22 β β smoothness priors (the Capital Asset Pricing Model, CAPM) CAPM 1 β β β β smoothness priors :,,. E-mail: koiti@ism.ac.jp., 104 1 TOPIX β Z i = β i Z m + α i (1) Z i Z m α i α i β i (the
More informationφ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1)
φ 4 Minimal subtraction scheme 2-loop ε 28 University of Tokyo Atsuo Kuniba version 2/Apr/28 Formulas Γ n + ɛ = n n! ɛ + ψn + + Oɛ n =,, 2, ψn + = + 2 + + γ, 2 n ψ = γ =.5772... Euler const, log + ax x
More information1 Jensen et al.[6] GRT S&P500 GRT RT GRT Kiriu and Hibiki[8] Jensen et al.[6] GRT 3 GRT Generalized Recovery Theorem (Jensen et al.[6])
Generalized Recovery Theorem Ross[11] Recovery Theorem(RT) RT forward looking Kiriu and Hibiki[8] Generalized Recovery Theorem(GRT) Jensen et al.[6] GRT RT Kiriu and Hibiki[8] 1 backward looking forward
More informationk + (1/2) S k+(1/2) (Γ 0 (N)) N p Hecke T k+(1/2) (p 2 ) S k+1/2 (Γ 0 (N)) M > 0 2k, M S 2k (Γ 0 (M)) Hecke T 2k (p) (p M) 1.1 ( ). k 2 M N M N f S k+
1 SL 2 (R) γ(z) = az + b cz + d ( ) a b z h, γ = SL c d 2 (R) h 4 N Γ 0 (N) {( ) } a b Γ 0 (N) = SL c d 2 (Z) c 0 mod N θ(z) θ(z) = q n2 q = e 2π 1z, z h n Z Γ 0 (4) j(γ, z) ( ) a b θ(γ(z)) = j(γ, z)θ(z)
More informationA B P (A B) = P (A)P (B) (3) A B A B P (B A) A B A B P (A B) = P (B A)P (A) (4) P (B A) = P (A B) P (A) (5) P (A B) P (B A) P (A B) A B P
1 1.1 (population) (sample) (event) (trial) Ω () 1 1 Ω 1.2 P 1. A A P (A) 0 1 0 P (A) 1 (1) 2. P 1 P 0 1 6 1 1 6 0 3. A B P (A B) = P (A) + P (B) (2) A B A B A 1 B 2 A B 1 2 1 2 1 1 2 2 3 1.3 A B P (A
More informationtokei01.dvi
2. :,,,. :.... Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 3. (probability),, 1. : : n, α A, A a/n. :, p, p Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN
More informationZ[i] Z[i] π 4,1 (x) π 4,3 (x) 1 x (x ) 2 log x π m,a (x) 1 x ϕ(m) log x 1.1 ( ). π(x) x (a, m) = 1 π m,a (x) x modm a 1 π m,a (x) 1 ϕ(m) π(x)
3 3 22 Z[i] Z[i] π 4, (x) π 4,3 (x) x (x ) 2 log x π m,a (x) x ϕ(m) log x. ( ). π(x) x (a, m) = π m,a (x) x modm a π m,a (x) ϕ(m) π(x) ϕ(m) x log x ϕ(m) m f(x) g(x) (x α) lim f(x)/g(x) = x α mod m (a,
More informationPower Transformation and Its Modifications Toshimitsu HAMASAKI, Tatsuya ISOMURA, Megu OHTAKI and Masashi GOTO Key words : identity transformation, pow
Power Transformation and Its Modifications Toshimitsu HAMASAKI, Tatsuya ISOMURA, Megu OHTAKI and Masashi GOTO Key words : identity transformation, power-normal distribution, structured data, unstructured
More informationIsogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206,
H28. (TMU) 206 8 29 / 34 2 3 4 5 6 Isogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206, http://link.springer.com/article/0.007/s409-06-0008-x
More information微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.
微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. ttp://www.morikita.co.jp/books/mid/00571 このサンプルページの内容は, 初版 1 刷発行時のものです. i ii 014 10 iii [note] 1 3 iv 4 5 3 6 4 x 0 sin x x 1 5 6 z = f(x, y) 1 y = f(x)
More informationI A A441 : April 15, 2013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida )
I013 00-1 : April 15, 013 Version : 1.1 I Kawahira, Tomoki TA (Shigehiro, Yoshida) http://www.math.nagoya-u.ac.jp/~kawahira/courses/13s-tenbou.html pdf * 4 15 4 5 13 e πi = 1 5 0 5 7 3 4 6 3 6 10 6 17
More information( ) (, ) arxiv: hgm OpenXM search. d n A = (a ij ). A i a i Z d, Z d. i a ij > 0. β N 0 A = N 0 a N 0 a n Z A (β; p) = Au=β,u N n 0 A
( ) (, ) arxiv: 1510.02269 hgm OpenXM search. d n A = (a ij ). A i a i Z d, Z d. i a ij > 0. β N 0 A = N 0 a 1 + + N 0 a n Z A (β; p) = Au=β,u N n 0 A-. u! = n i=1 u i!, p u = n i=1 pu i i. Z = Z A Au
More information201711grade1ouyou.pdf
2017 11 26 1 2 52 3 12 13 22 23 32 33 42 3 5 3 4 90 5 6 A 1 2 Web Web 3 4 1 2... 5 6 7 7 44 8 9 1 2 3 1 p p >2 2 A 1 2 0.6 0.4 0.52... (a) 0.6 0.4...... B 1 2 0.8-0.2 0.52..... (b) 0.6 0.52.... 1 A B 2
More information80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = i=1 i=1 n λ x i e λ i=1 x i! = λ n i=1 x i e nλ n i=1 x
80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = n λ x i e λ x i! = λ n x i e nλ n x i! n n log l(λ) = log(λ) x i nλ log( x i!) log l(λ) λ = 1 λ n x i n =
More information1 1 1.1...................................... 1 1.2................................... 5 1.3................................... 7 1.4............................. 9 1.5....................................
More informationI, II 1, A = A 4 : 6 = max{ A, } A A 10 10%
1 2006.4.17. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 1. 1. 2. 3. 4. 5. 2. ɛ-δ 1. ɛ-n
More informationSAMA- SUKU-RU Contents p-adic families of Eisenstein series (modular form) Hecke Eisenstein Eisenstein p T
SAMA- SUKU-RU Contents 1. 1 2. 7.1. p-adic families of Eisenstein series 3 2.1. modular form Hecke 3 2.2. Eisenstein 5 2.3. Eisenstein p 7 3. 7.2. The projection to the ordinary part 9 3.1. The ordinary
More information1 Fourier Fourier Fourier Fourier Fourier Fourier Fourier Fourier Fourier analog digital Fourier Fourier Fourier Fourier Fourier Fourier Green Fourier
Fourier Fourier Fourier etc * 1 Fourier Fourier Fourier (DFT Fourier (FFT Heat Equation, Fourier Series, Fourier Transform, Discrete Fourier Transform, etc Yoshifumi TAKEDA 1 Abstract Suppose that u is
More information2.1 H f 3, SL(2, Z) Γ k (1) f H (2) γ Γ f k γ = f (3) f Γ \H cusp γ SL(2, Z) f k γ Fourier f k γ = a γ (n)e 2πinz/N n=0 (3) γ SL(2, Z) a γ (0) = 0 f c
GL 2 1 Lie SL(2, R) GL(2, A) Gelbart [Ge] 1 3 [Ge] Jacquet-Langlands [JL] Bump [Bu] Borel([Bo]) ([Ko]) ([Mo]) [Mo] 2 2.1 H = {z C Im(z) > 0} Γ SL(2, Z) Γ N N Γ (N) = {γ SL(2, Z) γ = 1 2 mod N} g SL(2,
More informationDesign of highly accurate formulas for numerical integration in weighted Hardy spaces with the aid of potential theory 1 Ken ichiro Tanaka 1 Ω R m F I = F (t) dt (1.1) Ω m m 1 m = 1 1 Newton-Cotes Gauss
More informationJorgenson F, L : L: Inada lim F =, lim F L = k L lim F =, lim F L = 2 L F >, F L > 3 F <, F LL < 4 λ >, λf, L = F λ, λl 5 Y = Const a L a < α < CES? C
27 nabe@ier.hit-u.ac.jp 27 4 3 Jorgenson Tobin q : Hayashi s Theorem Jordan Saddle Path. GDP % GDP 2. 3. 4.. Tobin q 2 2. Jorgenson F, L : L: Inada lim F =, lim F L = k L lim F =, lim F L = 2 L F >, F
More informationn ξ n,i, i = 1,, n S n ξ n,i n 0 R 1,.. σ 1 σ i .10.14.15 0 1 0 1 1 3.14 3.18 3.19 3.14 3.14,. ii 1 1 1.1..................................... 1 1............................... 3 1.3.........................
More information1 M = (M, g) m Riemann N = (N, h) n Riemann M N C f : M N f df : T M T N M T M f N T N M f 1 T N T M f 1 T N C X, Y Γ(T M) M C T M f 1 T N M Levi-Civi
1 Surveys in Geometry 1980 2 6, 7 Harmonic Map Plateau Eells-Sampson [5] Siu [19, 20] Kähler 6 Reports on Global Analysis [15] Sacks- Uhlenbeck [18] Siu-Yau [21] Frankel Siu Yau Frankel [13] 1 Surveys
More informationf(x) = f(x ) + α(x)(x x ) α(x) x = x. x = f (y), x = f (y ) y = f f (y) = f f (y ) + α(f (y))(f (y) f (y )) f (y) = f (y ) + α(f (y)) (y y ) ( (2) ) f
22 A 3,4 No.3 () (2) (3) (4), (5) (6) (7) (8) () n x = (x,, x n ), = (,, n ), x = ( (x i i ) 2 ) /2 f(x) R n f(x) = f() + i α i (x ) i + o( x ) α,, α n g(x) = o( x )) lim x g(x) x = y = f() + i α i(x )
More informationI L01( Wed) : Time-stamp: Wed 07:38 JST hig e, ( ) L01 I(2017) 1 / 19
I L01(2017-09-20 Wed) : Time-stamp: 2017-09-20 Wed 07:38 JST hig e, http://hig3.net ( ) L01 I(2017) 1 / 19 ? 1? 2? ( ) L01 I(2017) 2 / 19 ?,,.,., 1..,. 1,2,.,.,. ( ) L01 I(2017) 3 / 19 ? I. M (3 ) II,
More informationJMP V4 による生存時間分析
V4 1 SAS 2000.11.18 4 ( ) (Survival Time) 1 (Event) Start of Study Start of Observation Died Died Died Lost End Time Censor Died Died Censor Died Time Start of Study End Start of Observation Censor
More informationZ: Q: R: C: sin 6 5 ζ a, b
Z: Q: R: C: 3 3 7 4 sin 6 5 ζ 9 6 6............................... 6............................... 6.3......................... 4 7 6 8 8 9 3 33 a, b a bc c b a a b 5 3 5 3 5 5 3 a a a a p > p p p, 3,
More informationx T = (x 1,, x M ) x T x M K C 1,, C K 22 x w y 1: 2 2
Takio Kurita Neurosceince Research Institute, National Institute of Advanced Indastrial Science and Technology takio-kurita@aistgojp (Support Vector Machine, SVM) 1 (Support Vector Machine, SVM) ( ) 2
More informationx () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x
[ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),
More information第8章 位相最適化問題
8 February 25, 2009 1 (topology optimizaiton problem) ( ) H 1 2 2.1 ( ) V S χ S : V {0, 1} S (characteristic function) (indicator function) 1 (x S ) χ S (x) = 0 (x S ) 2.1 ( ) Lipschitz D R d χ Ω (Ω D)
More informationTrapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x
University of Hyogo 8 8 1 d x(t) =f(t, x(t)), dt (1) x(t 0 ) =x 0 () t n = t 0 + n t x x n n x n x 0 x i i = 0,..., n 1 x n x(t) 1 1.1 1 1 1 0 θ 1 θ x n x n 1 t = θf(t n 1, x n 1 ) + (1 θ)f(t n, x n )
More informationuntitled
3 3. (stochastic differential equations) { dx(t) =f(t, X)dt + G(t, X)dW (t), t [,T], (3.) X( )=X X(t) : [,T] R d (d ) f(t, X) : [,T] R d R d (drift term) G(t, X) : [,T] R d R d m (diffusion term) W (t)
More information, CH n. CH n, CP n,,,., CH n,,. RH n ( Cartan )., CH n., RH n CH n,,., RH n, CH n., RH n ( ), CH n ( 1.1 (v), (vi) )., RH n,, CH n,., CH n,. 1.2, CH n
( ), Jürgen Berndt,.,. 1, CH n.,,. 1.1 ([6]). CH n (n 2), : (i) CH k (k = 0,..., n 1) tube. (ii) RH n tube. (iii). (iv) ruled minimal, equidistant. (v) normally homogeneous submanifold F k tube. (vi) normally
More informationst.dvi
9 3 5................................... 5............................. 5....................................... 5.................................. 7.........................................................................
More informationばらつき抑制のための確率最適制御
( ) http://wwwhayanuemnagoya-uacjp/ fujimoto/ 2011 3 9 11 ( ) 2011/03/09-11 1 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 2 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 3 / 46 (1/2) r + Controller - u Plant y
More information5 Armitage x 1,, x n y i = 10x i + 3 y i = log x i {x i } {y i } 1.2 n i i x ij i j y ij, z ij i j 2 1 y = a x + b ( cm) x ij (i j )
5 Armitage. x,, x n y i = 0x i + 3 y i = log x i x i y i.2 n i i x ij i j y ij, z ij i j 2 y = a x + b 2 2. ( cm) x ij (i j ) (i) x, x 2 σ 2 x,, σ 2 x,2 σ x,, σ x,2 t t x * (ii) (i) m y ij = x ij /00 y
More information01.Œk’ì/“²fi¡*
AIC AIC y n r n = logy n = logy n logy n ARCHEngle r n = σ n w n logσ n 2 = α + β w n 2 () r n = σ n w n logσ n 2 = α + β logσ n 2 + v n (2) w n r n logr n 2 = logσ n 2 + logw n 2 logσ n 2 = α +β logσ
More informationver.1 / c /(13)
1 -- 11 1 c 2010 1/(13) 1 -- 11 -- 1 1--1 1--1--1 2009 3 t R x R n 1 ẋ = f(t, x) f = ( f 1,, f n ) f x(t) = ϕ(x 0, t) x(0) = x 0 n f f t 1--1--2 2009 3 q = (q 1,..., q m ), p = (p 1,..., p m ) x = (q,
More information() [REQ] 0m 0 m/s () [REQ] (3) [POS] 4.3(3) ()() () ) m/s 4. ) 4. AMEDAS
() [REQ] 4. 4. () [REQ] 0m 0 m/s () [REQ] (3) [POS] 4.3(3) ()() () 0 0 4. 5050 0 ) 00 4 30354045m/s 4. ) 4. AMEDAS ) 4. 0 3 ) 4. 0 4. 4 4.3(3) () [REQ] () [REQ] (3) [POS] () ()() 4.3 P = ρ d AnC DG ()
More information( 28 ) ( ) ( ) 0 This note is c 2016, 2017 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercial purp
( 28) ( ) ( 28 9 22 ) 0 This ote is c 2016, 2017 by Setsuo Taiguchi. It may be used for persoal or classroom purposes, but ot for commercial purposes. i (http://www.stat.go.jp/teacher/c2epi1.htm ) = statistics
More informationuntitled
c 645 2 1. GM 1959 Lindsey [1] 1960 Howard [2] Howard 1 25 (Markov Decision Process) 3 3 2 3 +1=25 9 Bellman [3] 1 Bellman 1 k 980 8576 27 1 015 0055 84 4 1977 D Esopo and Lefkowitz [4] 1 (SI) Cover and
More information日本統計学会誌, 第45巻, 第2号, 329頁-352頁
45, 2, 2016 3 329 352 Market Risk Aggregation Using Copula and Its Application to Financial Practice Toshinao Yoshiba 2009 We investigate how a copula between risk factors takes portfolio diversification
More information2 G(k) e ikx = (ik) n x n n! n=0 (k ) ( ) X n = ( i) n n k n G(k) k=0 F (k) ln G(k) = ln e ikx n κ n F (k) = F (k) (ik) n n= n! κ n κ n = ( i) n n k n
. X {x, x 2, x 3,... x n } X X {, 2, 3, 4, 5, 6} X x i P i. 0 P i 2. n P i = 3. P (i ω) = i ω P i P 3 {x, x 2, x 3,... x n } ω P i = 6 X f(x) f(x) X n n f(x i )P i n x n i P i X n 2 G(k) e ikx = (ik) n
More information.. F x) = x ft)dt ), fx) : PDF : probbility density function) F x) : CDF : cumultive distribution function F x) x.2 ) T = µ p), T : ) p : x p p = F x
203 7......................................2................................................3.....................................4 L.................................... 2.5.................................
More information(2 X Poisso P (λ ϕ X (t = E[e itx ] = k= itk λk e k! e λ = (e it λ k e λ = e eitλ e λ = e λ(eit 1. k! k= 6.7 X N(, 1 ϕ X (t = e 1 2 t2 : Cauchy ϕ X (t
6 6.1 6.1 (1 Z ( X = e Z, Y = Im Z ( Z = X + iy, i = 1 (2 Z E[ e Z ] < E[ Im Z ] < Z E[Z] = E[e Z] + ie[im Z] 6.2 Z E[Z] E[ Z ] : E[ Z ] < e Z Z, Im Z Z E[Z] α = E[Z], Z = Z Z 1 {Z } E[Z] = α = α [ α ]
More information2 (March 13, 2010) N Λ a = i,j=1 x i ( d (a) i,j x j ), Λ h = N i,j=1 x i ( d (h) i,j x j ) B a B h B a = N i,j=1 ν i d (a) i,j, B h = x j N i,j=1 ν i
1. A. M. Turing [18] 60 Turing A. Gierer H. Meinhardt [1] : (GM) ) a t = D a a xx µa + ρ (c a2 h + ρ 0 (0 < x < l, t > 0) h t = D h h xx νh + c ρ a 2 (0 < x < l, t > 0) a x = h x = 0 (x = 0, l) a = a(x,
More informationuntitled
ECOH/YG No. 6 556 754 On the Role o Spectral Width and Shape Parameters in Control o Wave Height Distribution Yoshimi GODA and Masanobu KUDAKA -4-6-4 ( -4-6-4 Distributions o individual wave heights approximately
More information1 Nelson-Siegel Nelson and Siegel(1987) 3 Nelson-Siegel 3 Nelson-Siegel 2 3 Nelson-Siegel 2 Nelson-Siegel Litterman and Scheinkman(199
Nelson-Siegel Nelson-Siegel 1992 2007 15 1 Nelson and Siegel(1987) 2 FF VAR 1996 FF B) 1 Nelson-Siegel 15 90 1 Nelson and Siegel(1987) 3 Nelson-Siegel 3 Nelson-Siegel 2 3 Nelson-Siegel 2 Nelson-Siegel
More informationyasi10.dvi
2002 50 2 259 278 c 2002 1 2 2002 2 14 2002 6 17 73 PML 1. 1997 1998 Swiss Re 2001 Canabarro et al. 1998 2001 1 : 651 0073 1 5 1 IHD 3 2 110 0015 3 3 3 260 50 2 2002, 2. 1 1 2 10 1 1. 261 1. 3. 3.1 2 1
More information確率論と統計学の資料
5 June 015 ii........................ 1 1 1.1...................... 1 1........................... 3 1.3... 4 6.1........................... 6................... 7 ii ii.3.................. 8.4..........................
More informationkubostat2017e p.1 I 2017 (e) GLM logistic regression : : :02 1 N y count data or
kubostat207e p. I 207 (e) GLM kubo@ees.hokudai.ac.jp https://goo.gl/z9ycjy 207 4 207 6:02 N y 2 binomial distribution logit link function 3 4! offset kubostat207e (https://goo.gl/z9ycjy) 207 (e) 207 4
More information: ( )
: 7 3...........................................................................................3............................................. 3.4............................ 4.5.............................
More information0. I II I II (1) linear type: GL( ), Sp( ), O( ), (2) loop type: loop current Kac-Moody affine, hyperbolic (3) diffeo t
e-mail: koyama@math.keio.ac.jp 0. I II I II (1) linear type: GL( ), Sp( ), O( ), (2) loop type: loop current Kac-Moody affine, hyperbolic (3) diffeo type: diffeo universal Teichmuller modular I. I-. Weyl
More informationD = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j
6 6.. [, b] [, d] ij P ij ξ ij, η ij f Sf,, {P ij } Sf,, {P ij } k m i j m fξ ij, η ij i i j j i j i m i j k i i j j m i i j j k i i j j kb d {P ij } lim Sf,, {P ij} kb d f, k [, b] [, d] f, d kb d 6..
More informationSiegel Hecke 1 Siege Hecke L L Fourier Dirichlet Hecke Euler L Euler Fourier Hecke [Fr] Andrianov [An2] Hecke Satake L van der Geer ([vg]) L [Na1] [Yo
Siegel Hecke 1 Siege Hecke L L Fourier Dirichlet Hecke Euler L Euler Fourier Hecke [Fr] Andrianov [An2] Hecke Satake L van der Geer ([vg]) L [Na1] [Yo] 2 Hecke ( ) 0 1n J n =, Γ = Γ n = Sp(n, Z) = {γ GL(2n,
More informationr 1 m A r/m i) t ii) m i) t B(t; m) ( B(t; m) = A 1 + r ) mt m ii) B(t; m) ( B(t; m) = A 1 + r ) mt m { ( = A 1 + r ) m } rt r m n = m r m n B
1 1.1 1 r 1 m A r/m i) t ii) m i) t Bt; m) Bt; m) = A 1 + r ) mt m ii) Bt; m) Bt; m) = A 1 + r ) mt m { = A 1 + r ) m } rt r m n = m r m n Bt; m) Aert e lim 1 + 1 n 1.1) n!1 n) e a 1, a 2, a 3,... {a n
More informationII (Percolation) ( 3-4 ) 1. [ ],,,,,,,. 2. [ ],.. 3. [ ],. 4. [ ] [ ] G. Grimmett Percolation Springer-Verlag New-York [ ] 3
II (Percolation) 12 9 27 ( 3-4 ) 1 [ ] 2 [ ] 3 [ ] 4 [ ] 1992 5 [ ] G Grimmett Percolation Springer-Verlag New-York 1989 6 [ ] 3 1 3 p H 2 3 2 FKG BK Russo 2 p H = p T (=: p c ) 3 2 Kesten p c =1/2 ( )
More informationTitle 最適年金の理論 Author(s) 藤井, 隆雄 ; 林, 史明 ; 入谷, 純 ; 小黒, 一正 Citation Issue Date Type Technical Report Text Version publisher URL
Title 最適年金の理論 Author(s) 藤井, 隆雄 ; 林, 史明 ; 入谷, 純 ; 小黒, 一正 Citation Issue 2012-06 Date Type Technical Report Text Version publisher URL http://hdl.handle.net/10086/23085 Right Hitotsubashi University Repository
More informationsakigake1.dvi
(Zin ARAI) arai@cris.hokudai.ac.jp http://www.cris.hokudai.ac.jp/arai/ 1 dynamical systems ( mechanics ) dynamical systems 3 G X Ψ:G X X, (g, x) Ψ(g, x) =:Ψ g (x) Ψ id (x) =x, Ψ gh (x) =Ψ h (Ψ g (x)) (
More information* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *
* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *1 2004 1 1 ( ) ( ) 1.1 140 MeV 1.2 ( ) ( ) 1.3 2.6 10 8 s 7.6 10 17 s? Λ 2.5 10 10 s 6 10 24 s 1.4 ( m
More informationLebesgue Fubini L p Banach, Hilbert Höld
II (Analysis II) Lebesgue (Applications of Lebesgue Integral Theory) 1 (Seiji HIABA) 1 ( ),,, ( ) 1 1 1.1 1 Lebesgue........................ 1 1.2 2 Fubini...................... 2 2 L p 5 2.1 Banach, Hilbert..............................
More information(pdf) (cdf) Matlab χ ( ) F t
(, ) (univariate) (bivariate) (multi-variate) Matlab Octave Matlab Matlab/Octave --...............3. (pdf) (cdf)...3.4....4.5....4.6....7.7. Matlab...8.7.....9.7.. χ ( )...0.7.3.....7.4. F....7.5. t-...3.8....4.8.....4.8.....5.8.3....6.8.4....8.8.5....8.8.6....8.9....9.9.....9.9.....0.9.3....0.9.4.....9.5.....0....3
More information第5章 偏微分方程式の境界値問題
October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ
More informationUntitled
II 14 14-7-8 8/4 II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ 6/ ] Navier Stokes 3 [ ] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I 1 balance law t (ρv i )+ j
More informationφ s i = m j=1 f x j ξ j s i (1)? φ i = φ s i f j = f x j x ji = ξ j s i (1) φ 1 φ 2. φ n = m j=1 f jx j1 m j=1 f jx j2. m
2009 10 6 23 7.5 7.5.1 7.2.5 φ s i m j1 x j ξ j s i (1)? φ i φ s i f j x j x ji ξ j s i (1) φ 1 φ 2. φ n m j1 f jx j1 m j1 f jx j2. m j1 f jx jn x 11 x 21 x m1 x 12 x 22 x m2...... m j1 x j1f j m j1 x
More information_TZ_4797-haus-local
1.1.................................... 3.3.................................. 4.4......................... 8.5... 10.6.................... 1.7... 14 3 16 3.1 ()........................... 16 3. 7... 17
More information(a) (b) (c) Canny (d) 1 ( x α, y α ) 3 (x α, y α ) (a) A 2 + B 2 + C 2 + D 2 + E 2 + F 2 = 1 (3) u ξ α u (A, B, C, D, E, F ) (4) ξ α (x 2 α, 2x α y α,
[II] Optimization Computation for 3-D Understanding of Images [II]: Ellipse Fitting 1. (1) 2. (2) (edge detection) (edge) (zero-crossing) Canny (Canny operator) (3) 1(a) [I] [II] [III] [IV ] E-mail sugaya@iim.ics.tut.ac.jp
More informationp *2 DSGEDynamic Stochastic General Equilibrium New Keynesian *2 2
2013 1 nabe@ier.hit-u.ac.jp 2013 4 11 Jorgenson Tobin q : Hayashi s Theorem : Jordan : 1 investment 1 2 3 4 5 6 7 8 *1 *1 93SNA 1 p.180 1936 100 1970 *2 DSGEDynamic Stochastic General Equilibrium New Keynesian
More information(2004 ) 2 (A) (B) (C) 3 (1987) (1988) Shimono and Tachibanaki(1985) (2008) , % 2 (1999) (2005) 3 (2005) (2006) (2008)
,, 23 4 30 (i) (ii) (i) (ii) Negishi (1960) 2010 (2010) ( ) ( ) (2010) E-mail:fujii@econ.kobe-u.ac.jp E-mail:082e527e@stu.kobe-u.ac.jp E-mail:iritani@econ.kobe-u.ac.jp 1 1 16 (2004 ) 2 (A) (B) (C) 3 (1987)
More information: , 2.0, 3.0, 2.0, (%) ( 2.
2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................
More informationQCD 1 QCD GeV 2014 QCD 2015 QCD SU(3) QCD A µ g µν QCD 1
QCD 1 QCD GeV 2014 QCD 2015 QCD SU(3) QCD A µ g µν QCD 1 (vierbein) QCD QCD 1 1: QCD QCD Γ ρ µν A µ R σ µνρ F µν g µν A µ Lagrangian gr TrFµν F µν No. Yes. Yes. No. No! Yes! [1] Nash & Sen [2] Riemann
More information平成 29 年度 ( 第 39 回 ) 数学入門公開講座テキスト ( 京都大学数理解析研究所, 平成 29 ~8 年月 73 月日開催 31 日 Riemann Riemann ( ). π(x) := #{p : p x} x log x (x ) Hadamard de
Riemann Riemann 07 7 3 8 4 ). π) : #{p : p } log ) Hadamard de la Vallée Poussin 896 )., f) g) ) lim f) g).. π) Chebychev. 4 3 Riemann. 6 4 Chebychev Riemann. 9 5 Riemann Res). A :. 5 B : Poisson Riemann-Lebesgue
More informationACLI-EBC-CLHIA Interim Proposal _J_ June Final.PDF
/ ACLI-EBC-CLHIA 2 5800 2004 3 3 1 2004 3 55 7 ACLI EBC / CLHIA 3 3 20047 CTE 3 ACLI-EBC-CLHIA 3 CTE CTE CTE60 CTE 1) CTE(60) CTE(80) CTE(90) 2) 100 3) 8) 3 Mercer Oliver Wyman Actuarial Consulting ACLI-EBC-CLHIA
More informationPart () () Γ Part ,
Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35
More informationK 2 X = 4 MWG(f), X P 2 F, υ 0 : X P 2 2,, {f λ : X λ P 1 } λ Λ NS(X λ ), (υ 0 ) λ : X λ P 2 ( 1) X 6, f λ K X + F, f ( 1), n, n 1 (cf [10]) X, f : X
2 E 8 1, E 8, [6], II II, E 8, 2, E 8,,, 2 [14],, X/C, f : X P 1 2 3, f, (O), f X NS(X), (O) T ( 1), NS(X), T [15] : MWG(f) NS(X)/T, MWL(f) 0 (T ) NS(X), MWL(f) MWL(f) 0, : {f λ : X λ P 1 } λ Λ NS(X λ
More information,,,17,,, ( ),, E Q [S T F t ] < S t, t [, T ],,,,,,,,
14 5 1 ,,,17,,,194 1 4 ( ),, E Q [S T F t ] < S t, t [, T ],,,,,,,, 1 4 1.1........................................ 4 5.1........................................ 5.........................................
More informationa) Extraction of Similarities and Differences in Human Behavior Using Singular Value Decomposition Kenichi MISHIMA, Sayaka KANATA, Hiroaki NAKANISHI a
a) Extraction of Similarities and Differences in Human Behavior Using Singular Value Decomposition Kenichi MISHIMA, Sayaka KANATA, Hiroaki NAKANISHI a), Tetsuo SAWARAGI, and Yukio HORIGUCHI 1. Johansson
More informationkubostat2015e p.2 how to specify Poisson regression model, a GLM GLM how to specify model, a GLM GLM logistic probability distribution Poisson distrib
kubostat2015e p.1 I 2015 (e) GLM kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2015 07 22 2015 07 21 16:26 kubostat2015e (http://goo.gl/76c4i) 2015 (e) 2015 07 22 1 / 42 1 N k 2 binomial distribution logit
More informationk2 ( :35 ) ( k2) (GLM) web web 1 :
2012 11 01 k2 (2012-10-26 16:35 ) 1 6 2 (2012 11 01 k2) (GLM) kubo@ees.hokudai.ac.jp web http://goo.gl/wijx2 web http://goo.gl/ufq2 1 : 2 2 4 3 7 4 9 5 : 11 5.1................... 13 6 14 6.1......................
More informationhttp://www.ike-dyn.ritsumei.ac.jp/ hyoo/wave.html 1 1, 5 3 1.1 1..................................... 3 1.2 5.1................................... 4 1.3.......................... 5 1.4 5.2, 5.3....................
More information