統計数学セミナー

過去の記録 ~12/08次回の予定今後の予定 12/09~

担当者 吉田朋広、増田弘毅、荻原哲平、小池祐太
セミナーURL http://www.sigmath.es.osaka-u.ac.jp/~kamatani/statseminar/
目的 確率統計学およびその関連領域に関する研究発表, 研究紹介を行う.

過去の記録

2007年10月31日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
深澤 正彰 氏 (東京大学大学院数理科学研究科)
最尤推定量の漸近展開とその応用:とくに拡散過程の場合について
[ 講演概要 ]
最尤推定量とそのスチューデント化統計量の漸近展開公式について、 スキューネス修正の観点から考察し、AR過程や、あるクラスの拡散過程モデルへの応用について述べる。 一般の対称拡散過程モデルにおける最尤推定量のバイアス推定量、 スキューネス推定量も提案する。
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/07.html

2007年10月10日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
清 智也 氏 (東大情報理工)
勾配モデルの摂動解析と許容領域の評価
[ 講演概要 ]
多変量標準正規分布を凸関数の勾配写像によって 引き戻すと, 一つの確率分布が得られる. さらにパラメトリックな勾配写像を考えれば, 統計モデルが得られる. この統計モデルを勾配モデルと呼ぶことにする. 本講演は二つの内容からなる. 第一に, 恒等写像に摂動を加えた勾配写像を考え, 対応する密度関数, キュムラント母関数, ダイバージェンスなどの摂動展開を求める. 第二に, より具体的な勾配モデルに対して, パラメータが定義域に属すための十分条件を示す. このような考察の必要性は, 定義域が無限個の 制約式で与えられることによる.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/06.html

2007年07月25日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
小林 景 氏 (統計数理研究所, 学振特別研究員)
大規模ランダム行列のスペクトル理論とデータ解析への応用(Review)
[ 講演概要 ]
The empirical spectral distribution of random matrices have been studied since Wigner's pioneering work on the semicircular law in the 1950's. The result says that the empirical spectral distribution of a symmetric matrix with i.i.d. random elements converges to the semicircular law as the size of the matrix becomes large. Though this result is beautiful in theory, its application has been limited to a few problems in nuclear physics and coding theory. The next breakthrough was the Marcenko-Pastur (M-P) law for the asymptotic spectral distribution of sample covariance matrices. The M-P law has found more applications, in particular high dimensional statistical data analysis, than the semicircular law.

In this talk I will first review these two significant results. Each of them has three completely different proofs. Then I will explain several other theoretical results that have mostly been studied this decade. Finally, I will present some of the applications of these results. This review is partly based on lectures on random matrices given by P. Bickel, N. El-Karoui and A. Guionnet, and also some seminars at UC Berkeley.

(# This talk is almost the same as the talk I gave at ISM on June 1.)
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/05.html

2007年07月18日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
増田 弘毅 氏 (九州大学大学院数理学研究院)
Easy full-joint estimators of stable processes
[ 講演概要 ]
安定レヴィ過程からの高頻度観測に基づいた同時推定に関しては,尤度比確率場の漸近挙動は安易に二次項の観測情報量まで見ただけでは解明されず,最尤推定量の典型的な“良い漸近挙動”が保証されないことが分かっている.本発表では,標本メディアンおよび標本メディアンのプラグイン型統計量の漸近挙動に基づき,モデルに入る全パラメータの同時推定を可能とする計算容易な推定量の構成法を紹介し,正則な漸近共分散行列を有する漸近正規性を導出する.推定量の有限標本での挙動を数値実験で検証する.時間があれば,非対称安定レヴィ過程の場合に関して分かっている事柄についても触れる.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/04.html

2007年07月06日(金)

15:00-16:10   数理科学研究科棟(駒場) 122号室
Arturo KOHATSU-HIGA 氏 (大阪大学大学院基礎工学研究科)
Estimating multidimensional densities through the Malliavin-Thalmaier formula
[ 講演概要 ]
TBA
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/03.html

2007年06月27日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
小方 浩明 氏 (早稲田大学, 国際教養学部)
Empirical likelihood method for time series analysis
[ 講演概要 ]
For a class of vector-valued non-Gaussian stationary processes with unkown parameters, we develop the empirical likelihood approach which was proposed in the i.i.d. setting. In the time series analysis it is known that Whittle likelihood is one of fundamental tools to get a good estimator of unknown parameters and that the score functions are asymptotically normal. Motivated by the Whittle likelihood, we take its score as an estimating function and obtain the asymptotic distribution of our test statistic. Since the fitted spectral model may be different from true spectral structure, the results enable us to construct confidence rigions for various important time series parameters without knowing true spectral structure. We also consider the approach to a minimum contrast estimation and Cressie-Read power-divergence statistic. Numerical studies are introduced and illuminate some interesting features of the empirical approach.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/02.html

2007年06月06日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
小池 健一 氏 (筑波大学大学院数理物質科学研究科)
非正則な位置尺度母数分布族における位置母数の逐次点推定について
[ 講演概要 ]
有界な台をもつ非正則な位置尺度母数分布族に対して,その位置母数の逐次点推定を考える.ここでは,平均二乗誤差に費用も加えてリスクを考える.レンジに基づく停止則を提案し,これが漸近有効であることを示す.また,良く知られているRobbinsの逐次推定方式との比較を行い,密度関数の台の端点で密度関数が急激に変化する場合には,提案する逐次推定方式が標本数やリスクの意味で優れていることを示す.この結果は,逐次区間推定に関するKoike (2007)のものと同様であることが分かる
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/01.html

2007年05月23日(水)

16:20-17:30   数理科学研究科棟(駒場) 122号室
沖本 竜義 氏 (横浜国立大学経済学部・大学院国際社会科学研究科)
New Evidence of Asymmetric Dependence Structures in International Equity Markets
[ 講演概要 ]
A number of recent studies found two asymmetries in dependence structures in international equity markets; specifically, dependence tends to be high in (1) highly volatile markets and (2) bear markets. In this paper, a further investigation on asymmetric dependence structures in international equity markets is performed under the use of the Markov switching model and copula theory. Combining these two theories enables us to model dependence structures with sufficient flexibility. Using this flexible framework we indeed found that there are two distinct regimes in the US-UK market. We also showed that, for the US-UK market, the bear regime is better described by an asymmetric copula with lower tail dependence with clear rejection of the Markov switching multivariate Normal model. In addition, we showed ignorance of this further asymmetry in bear markets is very costly for risk management. Lastly, we conducted similar analysis for other G7 countries, where we found other c ases where the use of a Markov switching multivariate Normal model would be inappropriate.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2007/00.html

2007年01月31日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
西山 慶彦 氏 (京都大学経済研究所)
A Sequential Unit Root Test
[ 講演概要 ]
It is well known that conventional unit root tests such as Dickey=Fuller and its variants do not have good power properties when sample size is not large. Lai and Siegmund (1983, AS) proved that OLS estimator of the AR(1) coefficient is asymptotically normally distributed in a sequential framework even if the time series has a unit root unlike the OLS estimator under a standard sampling scheme. We pursue this direction to propose a unit root test under a sequential sampling. The proposed test uses not only the OLS estimator of the AR(1) coefficient, which is asymptotically normal, but also the stopping time to construct the critical region, anticipating a better power property. We obtain analytic expressions of the joint distribution of the two statistics as well as its marginals under the null. We also consider the distribution of the statistics under local alternatives. The properties of the stopping time, to the best of our knowledge, have not been studied in the unit root literature. We calculate its expectation and variance.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/18.html

2007年01月17日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
玉置 健一郎 氏 (早稲田大学)
Second order optimality for estimators in time series regression models
[ 講演概要 ]
We consider the second order asymptotic properties of an efficient frequency domain regression coefficient estimator $\\hat{\\beta}$ proposed by Hannan (1963). This estimator is a semiparametric estimator based on nonparametric spectral estimators. We derive the second order Edgeworth expansion of the distribution of $\\hat{\\beta}$. Then it is shown that the second order asymptotic properties are independent of the bandwidth choice for residual spectral estimator, which implies that $\\hat{\\beta}$ has the same rate of convergence as in regular parametric estimation. This is a sharp contrast with the general semiparametric estimation theory. We also examine the second order Gaussian efficiency of $\\hat{\\beta}$. Numerical studies are given to confirm the theoretical results.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/17.html

2006年12月06日(水)

15:00-16:10   数理科学研究科棟(駒場) 128号室
Stefano IACUS 氏 (Department of Economics Business and Statistics, University of Milan, Italy)
Inference problems for the standard and geometric telegraph process
[ 講演概要 ]
The telegraph process {X(t), t>0}, has been introduced (see Goldstein, 1951) as an alternative model to the Brownian motion B(t). This process describes a motion of a particle on the real line which alternates its velocity, at Poissonian times, from +v to -v. The density of the distribution of the position of the particle at time t solves the hyperbolic differential equation called telegraph equation and hence the name of the process. Contrary to B(t) the process X(t) has finite variation and continuous and differentiable paths. At the same time it is mathematically challenging to handle.

In this talk we will discuss inference problems for the estimation of the intensity of the Poisson process, either homogeneous and non homogeneous, from continuous and discrete time observations of X(t). We further discuss estimation problems for the geometric telegraph process S(t) = S(0) * exp{m - 0.5 * s^2) * t + s X(t)} where m is a known constant and s>0 and the intensity of the underlying Poisson process are two parameter of interest to be estimated. The geometric telegraph process has been recently introduced in Mathematical Finance to describe the dynamics of assets as an alternative to the usual geometric Brownian motion.

For discrete time observations we consider the "high frequency" approach, which means that data are collected at n+1 equidistant time points Ti=i * Dn, i=0,1,..., n, n*Dn = T, T fixed and such that Dn shrinks to 0 as n increases.

The process X(t) in non Markovian, non stationary and not ergodic thus we use approximation arguments to derive estimators. Given the complexity of the equations involved only estimators on the standard telegraph process can be studied analytically. We will also present a Monte Carlo study on the performance of the estimators for small sample size, i.e. Dn not shrinking to 0.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/16.html

2006年11月22日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
鎌谷 研吾 氏 (東京大学大学院数理科学研究科)
A Note on Haplotype Estimation
[ 講演概要 ]
Haplotype information is important for many analyses but it is not always possible to obtain. This work is motivated to seek haplotype information from diploid population data. We present a new approach to know the haplotype information using classical methods. We do not intend to say that our method is better than the well-known EM based approache for practical purposes, but our way is attractive in some sense.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/15.html

2006年11月17日(金)

15:00-16:10   数理科学研究科棟(駒場) 118号室
清水 泰隆 氏 (大阪大学大学院基礎工学研究科)
Functional estimation of L'evy measure for jump-type processes
[ 講演概要 ]
Recently, stochastic processes with Poissonian jumps are frequently used in finance and insurance. In their applications, it often becomes important to estimate some functionals of integral types with respect to L'evy measures. In this talk, we propose a nonparametric estimator of their functionals based on both continuous and discrete observations. If time permits, we shall also mention the application to the mathematical insurance, in particular, the estimates of ruin probabilities for genelarized risk processes.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/13.html

2006年11月08日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
蛭川 潤一 氏 (早稲田大学)
LAN Theorem for Non-Gaussian Locally Stationary Processes and Their Discriminant and Cluster Analyses
[ 講演概要 ]
This talk is concerned with asymptotic inference for non-Gaussian locally stationary processes. Lucien LeCam established the most important and sophisticated foundation of the general statistical asymptotic theory. He introduced the concept of local asymptotic normality (LAN) for the likelihood ratio of general statistical models. Once LAN is proved, the asymptotic optimality of estimators and tests is described in terms of the LAN property. The techniques of statistical inference for stationary time series have been well established. However, stationary time series model is not plausible to describe the real world. One of the difficult problem when we deal with nonstationary processes is how to set up an adequate model. Otherwise, the observation in the future will bring no information for the present structure. Recently, Dahalhaus has proposed an important class of nonstationary processes, called locally stationary processes. Locally stationary processes have the time varying densities whose spectral structures smoothly change in time. In this talk, we first show the LAN results for locally stationary processes under the assumption of the non-Gaussianity. Then, we apply the LAN theorem to estimation and testing theory, non-Gaussian robustness and adaptive estimation. Our LAN theorem elucidates various non-Gaussian asymptotics. Next, we develop asymptotic theory for discriminant and cluster analyses of non-Gaussian locally stationary processes. We discuss about non-Gaussian robustness of our classification statistic. Furthermore, we execute the clustering of stock returns in Tokyo Stock Exchanges. Consequently, we observe that the clustering results well extract features of relationships among companies.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/13.html

2006年11月01日(水)

15:00-16:10   数理科学研究科棟(駒場) 128号室
Ilia NEGRI 氏 (Department of Management and Information Technology, University of Bergamo, Italy)
Some problems related to the estimation of the invariant measure of an ergodic diffusion.
[ 講演概要 ]
We consider a one dimensional ergodic diffusion process solution of a stochastic differential equation. We suppose that the diffusion coefficient is known whereas the drift coefficient $S$ is unknown. Our interest is the invariant measure of the process denoted as $\\mu $. We denote by $f_S$ and $F_S$ the invariant density and the invariant distribution function of $\\mu$ respectively. We present the problems of finding efficient estimators when we observe the trajectory of the diffusion in continuos time over $[0,T]$ and we study asymptotic properties of the estimators when $T$ goes to infinity.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/12.html

2006年10月18日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
松田 安昌 氏 (東北大学大学院経済学研究科)
Fourier analysis of irregularly spaced data on R^d
[ 講演概要 ]
1970年代にはじまった時系列解析の研究は、1変量時系列からはじまり、多変量時系列、空間系列、時空間系列へと対象が拡張されてきている。近年の空間系列、時空間系列の文献を調べてみれば、時系列と同じく等間隔に観測されたデータ(regularly spaced data) を前提としている場合が多い。しかし時系列とは異なり、空間系列ではランダムに散らばった地点で観測されたデータ (irregularly spaced data)を仮定する方が、応用範囲がひろい。そこで本発表では、irregularly spaced dataによる空間相関の推定問題を扱う。具体的にはデータを周波数領域に変換してスペクトルを推定する方法を提案し、推定量が一致性、漸近正規性を持つための条件を示す。本方法による東京地価データ分析例も紹介する。本研究は、矢島美寛教授(東京大学大学院経済学研究科)との共同研究です。
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/11.html

2006年10月11日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
石川保志 氏 (愛媛大学)
Optimal stopping problem associated with jump-diffusion processes
[ 講演概要 ]
We study an optimal stopping problem of some performance function with respect to a jump-diffusion process.

2006年09月25日(月)

13:00-14:10   数理科学研究科棟(駒場) 128号室
西山 陽一 氏 (統計数理研究所)
Nonparametric testing time-homogeneity for L'evy processes
[ 講演概要 ]
First, a review about uniform central limit theorems for martingales is given. The main part of the talk is concerned with a change point problem for L'evy processes. The null hypothesis is that the L'evy process is time-homogeneous, and the alternative is that the L'evy measure changes at a certain time point of the observation period. We propose an empirical process type statistics, and derive its asymptotic behaviour under the null and the alternative hypotheses. The limiting distribution under the null hypothesis is a functional of the standard Brownian motion.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/09.html

2006年08月22日(火)

15:30-16:40   数理科学研究科棟(駒場) 128号室
Jeannette H.C. WOERNER 氏 (University of Gottingen)
A unifying approach to inference in semimartingale and long-memory models
[ 講演概要 ]
Over the recent years classical stochastic volatility models based on Brownian motion have been generalized in different ways, either replacing the Brownian motion by a pure jump Levy process, which leads to a pure jump model, or by a fractional Brownian motion, which makes it possible to model both long memory or turbulent behaviour. We consider robust and easily computable estimators for the inte- grated volatility, providing insight in the level of volatility, as needed for risk measurement and pricing of variance and volatility swaps. We discuss consistency and distributional results for the power and multi- power variation estimates based on high frequency data. Furthermore, we consider robustness against additive components and compare the results for the different classes of semimartingale and fractional Brow- nian motion models.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/07.html

2006年08月22日(火)

16:50-18:00   数理科学研究科棟(駒場) 128号室
Delphine DAVID 氏 (Departement de Mathematiques, Universite de La Rochelle)
A computation of Theta in a jump diffusion model by integration by parts
[ 講演概要 ]
Using Malliavin weights in a jump-diffusion model we obtain an expression for Theta (the sensitivity of an option price with respect to the time remaining until exercise), with application to non-smooth payoff functions. Optimal weights are computed by minimization of variance and numerical simulations are presented for digital and European options. Some results are also presented for Asian options. Our representation formula for Theta differs in general from the one obtained from the Black-Scholes PDE in terms of Delta and Gamma.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/08.html

2006年07月19日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
深澤 正彰 氏 (東京大学大学院数理科学研究科)
Edgeworth Expansion for Likelihood Analysis on Ergodic Diffusions with applications to Bootstrap
[ 講演概要 ]
We shall consider the maximal lilelihood estimator for the drift coefficient of a given one-dimensional diffusion. An Edgeworth expansion formula will be presented and verify a second-order correct confidence interval we shall newly propose. We are also going to mention the likelihood ratio statistic, which enjoys second-order correctness. There are Bootstrap methods closely related to the subject and introduced recently by the author. Some generalized results on those methods will be also introduced in this talk.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/06.html

2006年06月21日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
石川 保志 氏 (愛媛大学理学部)
Malliavin calculus applied to mathematical finance and a new formulation of the intgration-by-parts
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/05.html

2006年05月31日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
津田 美幸 氏 (統計数理研究所)
Bhattacharyya inequality for quantum state estimation II
[ 講演概要 ]
前回導出した三種類(S型, R型, L型)の量子Bhattacharyya不等式を量子ガウス状態の複素振幅θの多項式g(θ)の推定問題に応用する. 量子ガウス状態は, レーザ光の量子状態の典型的なモデルであり, 量子光学や量子情報で重要な研究対象である. 未知の複素振幅θを推定する方法としては, θが実軸にある場合はホモダイン測定, 一般の複素数の場合はヘテロダイン測定が知られており, それぞれS型とR型の量子Cramer-Rao不等式の下限を達成するUMVUEである. さらにここでは, θが実数の場合に (1), (2)を示し, θが複素数の場合に (3), (4)を示す.

(1) g(θ)=θ^2に対するUMVUEは存在してS型Bhattacharyya下限を達成する. その推定量は, 物理系にスクイジングと呼ばれる操作を施した後の個数測定によって与えられる.
(2) g(θ)=θ^3に対するUMVUEは, 生成消滅作用素の多項式の形では存在しない.
(3) g(θ)が正則, 或いは反正則, ならば, ヘテロダイン測定によってUMVUEが与えられ, それぞれR型, L型のBhattacharyya下限を達成する.
(4) g(θ)が実数値ならば, ある測定によりUMVUEが与えられ, R型, L型両方の下限を達成する.

量子ガウス状態は古典の正規分布に似ている. しかし, 古典では, 平均の多項式は Hermite多項式により常にUMVUEを構成できるが, 量子では上記のように事情が異なる.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/04.html

2006年05月24日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
津田 美幸 氏 (統計数理研究所)
Bhattacharyya inequality for quantum state estimation I
[ 講演概要 ]
量子状態推定のBhattacharyya不等式の導出とその応用例を前後二回に分けて紹介する. 今回は, 一般的な形で問題設定を行い, 量子Cramer-Rao不等式と量子 Bhattacharyya不等式について述べる.
量子状態推定は量子力学系の未知の状態に関する統計的推定問題である. 古典的な統計モデルの推定問題との違いは, データを観測するための測定を, 量子力学的に可狽ネ範囲で, 選択する点にある. 実数または複素数でパラメトライズされたモデルに対しては, 不偏推定量の分散の最小化が基本的な問題である. ただしここでは, 複素パラメータz=x+iyの分散とは, (x,y)の二次元の共分散行列のトレースをさす. この問題に対しては Yuen and Lax (1973) 等により, パラメータの一階微分に基づいた量子Cramer-Rao不等式が導出されており, 量子ガウス状態の複素振幅θのUMVUEが知られている. 二階以上の微分に基づくBhattacharyya型の不等式は, Brody and Hughston (1998) により, ある特殊なモデルにおいて導入され, 漸近論へ応用された. ここではより一般的なモデルに適用可能な形で量子Bhattacharyya不等式を三種類定式化する.
[ 参考URL ]
https://www.ms.u-tokyo.ac.jp/~kengok/statseminar/2006/03.html

2006年05月17日(水)

16:20-17:30   数理科学研究科棟(駒場) 128号室
Arnak DALALYAN 氏 (Universite Paris 6, France)
Second-order efficiency in the semiparametric problem of estimating the shift of a signal

< 前へ 12345678910 次へ >