FEBRUARY. 2015

No.

127

The Institute of
Statistical Mathematics

02

研究室訪問
船渡川 伊久子／荻原 哲平

06

シンポジウム報告
第 9 回 Biostatisticsネットワークの開催
ワークショップ「ゲノム多様性データの統計解析」の開催報告「R-Training: Introduction to R statistics」開催報告

07

研究教育活動08

統数研トピックス「日本人の国民性第 13 次全国調査」記者発表「サイエンスアゴラ2014」に出展／リスーピア賞受賞／ SC14に出展
大学共同利用機関シンポジウムに参加・出展
科学技術コミュニケーションシンポジウムを開催
英国インペリアルカレッジ・ロンドンのGuoデータサイエンス研究所長が統計数理研究所を訪問

11

総合研究大学院大学複合科学研究科統計科学専攻関係

12

お知らせ
公開講座

14

共同利用
平成 26 年度共同利用公募追加採択課題14

外部資金・研究員等の受入れ
受託研究の受入れ／外来研究員の受入れ／寄附金の受入れ

15

人事

15

所外誌掲載論文等

16

刊行物
Research Memorandum（2014.11〜2015.1）
統計数理研究所調査研究リポート／統計数理
Annals of the Institute of Statistical Mathematics

18

コラム

情報・システム研究機構

2014 年 11月―12月の公開講座実施状況
統計数理セミナー実施報告（2014 年 11月〜2015 年 1月）
2014 年公開講演会「ビッグデータを活かせ：スパコン利用の最前線」

大学共同利用機関法人

CONTENTS

統計数理研究所ニュースNews研 究 室 訪 問

船渡川 伊久子
データ科学研究系
計量科学グループ准教授

02

The Institute of Statistical Mathematics

既存データを新しい解釈で読み解くtResearch Fron

Line
「データを見ながら考えていると、『あ！』と気づく瞬間があります。数字がそれまでとは違う意味を持って見え始めるのです」。
手持ちのデータから推論を立て、それを裏付けるために足りないデータを探し出して立証する。このプロセスには、謎解きのような面白さがあると船渡川は言う。

経済学の自己回帰モデルを個体差のある医学に応用
専門は経時データ解析。経時データとは、複数の対象者それぞれから複数の時点で繰り返し反応を測定したデータのことだ。
モデルを使った数理的な解析方法と、既存のデータからわかりやすい概念を導き出す方法の両方を研究している。
東京大学の学生時代は疫学・生物統計学を専攻した。1990
年代初頭、まだ日本では新しい分野であり、医学系のデータを扱うデータ解析手法の開発も黎明期だった。船渡川は未知の世界に惹かれ、同大学の大橋靖雄教授が日本で初めて開設した生物統計の教室に飛び込んだ。卒業論文の作成にあたり紹介された経時測定データ解析が、今につながる研究テーマになった。
統計数理研究所に入ったのは2014 年 4月。大学時代から取り組む統計モデルの研究については、前任所の帝京大学医学部を経て現在、特に自己回帰線形混合効果モデルに注力している。このモデルを使い、副甲状腺機能亢進症の薬物療法における経時データ解析に携わった。投薬と生体反応の推移を見て、適切な投与量を探るのが目的だ。それまで多く使われていた
線形混合効果モデルに、経済学で用いられる自己回帰の概念を取り入れた点が斬新だった。1 時点前の反応数値を考慮する
この方法により、過去の投与履歴が反映され、時間の経緯とともに薬効が均衡する過程を合理的に、また、より的確に当てはめられるようになった。

カラクリを解くカギは「出生コホート」
医学部では、公衆衛生に関する研究も開始した。まず、肥満の問題に取り組み、肥満度の指標であるBMI（ボディマスイン

図 1．線形混合効果モデルと自己回帰線形混合効果モデル

Research Front Line
他分野と連携して、社会に役立つ研究を続けたいと思っていますデックス）のデータを見直した。公表されているデータは、年齢

出生年を考慮しながら肺がん死亡率と対比した。

別の横断調査によるものが一般的だ。日本人の成人男性のグラ

喫煙は肺がんのみならず多くの疾患のリスク要因だが、喫煙

フは逆 U 字カーブを描き、一旦上昇したのち年を取るにつれ痩を開始してから死亡までの期間が長いために、両者の関連につせていた。また、日本人女性では20 歳代は年々痩せていくのにいては不明な点が少なくない。例えば、1980 年代前半まで、日対し、60歳代では逆に年々太っていく。船渡川は「出生コホート」

本人の男性は喫煙率が高いにもかかわらず肺がん死亡率が低に着目して、これらの現象の原因を解き明かした。

かったことから、「日本人は喫煙の影響を受けにくい国民である」

コホートとは共通した因子を持つ集団を意味する。出生年を共という見方があった。船渡川は日英米の出生コホートから、日本通項とする出生コホートでBMIを調べた結果、世代間で比較すの喫煙開始年齢が他国より遅かったことを突き止めた。また女ると男性では新しい世代ほどBMI が高いが、女性では1930 年

性では、1920 年以前に生まれた喫煙率の高い集団と、高い肺代生まれの集団が最も高く、それより前あるいは後に生まれた世

がん死亡率を示した集団が一致している。

代ほど低いことがわかった。世代ごとの年齢による推移は、いず

このように、既存のデータから新しい解釈を見つけ出す手腕

れの世代も6 歳頃から60 歳代までBMIは上昇する。ただし、日が、船渡川の真骨頂だ。ただ、喫煙と疾病の関連のように長期

本人女性では、10 代後半頃から減少する時期がある。

スパンで追う必要のある研究テーマでは、昔のデータが見つか

ところが、こうした出生コホートの経過は、年齢別グラフには表らないことも多い。この研究では、イギリスの喫煙開始年齢のれない。もともとBMI が低い世代が差し掛かった年齢ではBMIデータを探すのに苦労した。
「推論であやふやな部分を確認すが下降し、もともと高い世代が差し掛かった年齢では上昇しているためには、どうしても必要なデータでした。目星をつけたイギリ

く。男性の年齢別グラフでは60 歳代に他世代よりもBMI が低い

スの古本でやっと見つけ、やっぱり、と確信しました」。ミッシング世代が分布しているため、見かけ上は「いったん太って再び痩リンクがつながった瞬間だった。

せる」ように表されていた。実際は、同じ年代に生まれた男性たちが加齢により痩せていたわけではない。

これまで、医学や健康学の分野で統計を扱ってきた船渡川。
「統数研に着任し、人口学や経済学、生態学など様々な分野でさらに、日本人女性は、最近のコホートほど10 代の数値が下

統計が活用されているのを見て新鮮な気持ちです。他分野と連がり始める年齢が若く、BMI が低くなっていることもわかった。

携して、社会に役立つ研究を続けたいと思っています。よりわか

1970 年代生まれでは、子どもの頃は40 年代生まれより太っていりやすく説明するために、データの可視化などにも取り組んでみるのに、16 歳以降になると逆転し、痩せている。船渡川は肥満

たい」。そう話す眼差しは、まっすぐに先へ向かっている。

に関するこうした研究成果をBMJ（イギリス医師会雑誌）や疫（広報室）

学の国際ジャーナルなどに発表した。

根気よくデータを探して ミッシングリンク をつなぐ
肥満の次に選んだテーマが喫煙だ。日本、イギリス、アメリカ
における喫煙開始年齢と喫煙率の長期の加齢変化を推定し、

図 2．Body Mass Index（BMI）の年齢による変化。1940 年代生まれと1970年代生まれの男女。国民健康・栄養調査より作成。

図 3．喫煙開始割合、喫煙率、肺癌死亡率の年齢変化。1925 年生まれ日本
人男女と米国男、非喫煙者。Funatogawaら(WHO Bulletin, 2014) 改変。

The Institute of Statistical Mathematics

03

研 究 室 訪 問

荻原 哲平
数理・推論研究系
統計基礎数理グループ助教

04

The Institute of Statistical Mathematics

高頻度観測金融データからリスク管理を研究

tResearch FronLine
1984 年（昭和 59 年）生まれの若手研究者である。千葉、神奈川で育ち、千葉の渋谷教育学園幕張高校から東大数学科へ進み、修士課程を終了した。数学は早くから得意で、中学で高校の、高校で大学教養課程の教科書をマスターした。幼少期に父親から数のナゾナゾの本を与えられたりしたことが刺激になったという。
2009 年、三菱 UFJトラスト投資工学研究所へ入った。世界経済が大混乱したリーマンショックの翌年だった。投資する株を数式のモデルを立てて選ぶ仕事に携わった。人間のカンだけではなく数式でも投資し、リスク分散をはかる。発展途上国の株式市場を対象に利益をあげることに貢献したという。
ここでの最大の体験は2011 年 3月11日の東日本大震災で株式市場の激変ぶりを目の当たりにしたことだった。出向していた銀行で、ファンドマネージャーの隣席にいた。多くの株価が激しく変動し、それまで優良株の代表だった東電株はその瞬間からガラッと動きを変えた。株価は秒単位で動いた。この時の体験がその後の研究に影響することになる。

秒単位の高頻度データの中で複数株式の連動性を見る
勤めながら東大の博士課程に通い、2012 年 7月、大阪大学
金融・保険教育研究センターの特任助教になった。経済、理、
基礎工学部など専門の枠を超え、金融の関連を研究している研究者が横断的に集まっているところで、実務家との協力関係も盛んだった。高頻度観測金融データの研究に取り組み始めた。
いまの株売買はコンピューターを利用するため秒単位で成立する。大型株だと1日数千回の売買があり、市場全銘柄のその日の取引はすべて記録される。１日の終値しか記録されなかった時代とは比較にならないデータ量である。この高頻度データを読み取って有意な情報をつかみ、リスク管理に役立てるのが研究の目的である。
統計数理研究所には2014 年 7月に移り、研究はより本格化し

図 1．高頻度観測金融データには、株価の観測時刻が一致しない「非同期観測」
や観測ノイズの混入の問題がある。

Research Front Lineパラダイムを変えるような局面のきっかけとなる数学をつくればいいかなと思っていますている。秒単位の高頻度データの中で主に複数株式の連動性

実はこれらの研究はハードルが高い。厖大な取引データを収に関する研究を行っている。例えば、株の値動きは、同じ産業

容するハードディスクや証券市場から購入するデータそのものが内では同じように動くことが多い。その際、より大企業の方が、そ高額だ。高性能コンピューターも不可欠。大規模な研究機関でれ以下の企業よりも早く動く傾向がある。先行、遅行の差が出ないと取り組めないこともあり、国内の研究者はあまりいない。欧る。まったく産業の違う企業の株が連動して動くこともある。これ

米諸国と比べると人数、研究環境の整備ともに遅れている。

らの連動性の特徴を秒単位・分単位のデータからとらえ、日々

現在よりもっと正確なリスク管理手法をそう遠くない先にのリスク管理の参考とする。

ノイズ除去など新しい統計手法の開発にも取り組む

しかし、難しかった高頻度データの解析も研究段階ではクリア
されつつある。荻原の研究が高いレベルに達した時は、連動性高頻度観測金融データでは、株価の変動の大きさや株取引の少ない複数の株を持つなどの方法で安定的なリスク管理が十の活発さ、つまり株を売りたい時にすぐに買いの注文が入るかど分に期待できる。
「リーマンショックのような急激な変化の時に最うかという流動性もより正確に見ることができる。年金などの機関初の一発を避けるのは無理ですが、市場が不安定になった時

投資家にとってはリスク管理上、大きな関心事である。

に現在よりもっと正確にリスク管理できる手法はできるかなと思っ

わずか秒単位の変化を統計的にとらえることで、リスク管理にています。そう遠くない先です」

実際に対処できるのだろうか。荻原は「数秒の差でも役立ちま2015 年 3月には、この研究が進んでいるフランスの大学へ 1

す。統計的に有効な情報を取り出すには数百、数千程度のデー

か月ほど行くことになっている。フランスでは研究機関と実務家のタは必要なので、１日の中で有効な情報を取り出すためには秒提携が進んでいる。ここで研究者仲間のネットワークを広げ、共単位の変化のデータを集める必要があります。市場には数千銘

同研究に進めることも視野に入れている。

柄もあるので、連動性の高さ、大きさは人間の判断で対応でき

当面はこの研究に力を注いでいくつもりだ。
「高頻度データのず、数値化して管理しなければならない。銘柄間の連動性をも

リスク評価法は、いまは理論的なところと実証に少し手を出したとにリスク管理を行うので、連動性が一番、重要な情報です」とところなので、もう少し金融機関が使える実用的なレベルにする語る。

ことを近い目標としていこうというところです」。まだ30 歳、前途は研究所での専門は数理統計学と確率解析である。高頻度観洋々である。
「遠い先の目標は」と尋ねると「社会の中で科学が測データの特性に応じた新しい統計手法の開発にも取り組んで

１つ、レベルが上がり、パラダイムを変えるような局面のきっかけといる。その１つはデータの中に含まれているノイズを薄めること。

なる数学をつくればいいかなと思っています」と頼もしい答えがデータを自動的に平均化する計算方式をつくった。もう１つは、返ってきた。具体的には機械学習や人工知能の研究を金融や取引時期の違う複数株の連動性を見るために非同期観測の方

関連の分野へ応用することを考えているという。これからが大い法を研究すること。阪大の時に取引時期が違っても連動の大きに楽しみな研究者である。

（広報室）

さを効率的に推定する方法を見つけた。昔の１日ごとのデータといまの高頻度データを比較することができるようになった。

図 2．最新の統計手法により、株価の高頻度データから非同期観測や観測ノイズを除去し、株価相関が計算される。

図 3．株式注文情報から、取引量に応じた価格変動率を示した図。取引コストを抑える上で価格変動率の予測は重要。

The Institute of Statistical Mathematics

05

シンポジウム報告
Report

第9回Biostatisticsネットワークの開催
平成 26 年 8月3日、京都大学大学院医学研究科におい間を通して、参加大学の学生同士が協同して企画・準備をて、統 計 思 考 院の公 募 型 人 材 育 成 事 業の一 環として、進めることにより、交流を深められるという点からも大きな収「Biostatisticsネットワーク」が開催されました。Biostatistics

穫が得られています。本年度は、生存時間解析や診断法

ネットワークは、国内の医学統計学を専攻する大学院生の研の開発における統計手法、観察研究における傾向スコアを

究交流のための集会であり、今回で9 度目の開催を迎えま用いた統計解析、ゲノムワイド研究におけるスパース推定をす。本年度は、京都大学、総合研究大学院大学、東京理用いた遺伝子選択方法など、大学ごとに特色のある、個性科大学、久留米大学、大分大学の5 大学が参加し、45 名的な研究報告が行われました。

の参加者を得て、盛況のもと、開催されました。
本研究集会は、参加大学の学生の手によって、企画・運本研究会が、有意義な交流の場となり、参加者の皆さんの将来に資する機会となれば幸いです。

（野間久史）

営が行われるということも大きな特徴であり、1 年間の準備期

Report

ワークショップ「ゲノム多様性データの統計解析」の開催報告
ワークショップ「ゲノム多様性データの統計解析」は、統計

イズ統計を組み合わせてパラメータを推定する方法について

思考院公募型人材育成事業の補助により2014 年 12月4日の実践的な解説が提供されました。ワークショップ後にアン（木）と12月5日（金）の二日にわたり、九州大学西新プラザ

ケートを行なったところ、おおむね好意的な反応をいただきま

中会議室で開催されました。北は北海道から南は沖縄まで、した。

（九州大学理学部・手島康介／間野修平）

45 名の学生および研究者のみなさまにご参加いただきました。初日は解析に必要となる集団遺伝および統計的推測の基礎に関する話が、総合研究院大学院大学、統計数理研
究所、九州大学の演者らによって、数理的かつ直感的に提供されました。二日目はデータが得られた後の処理および解析に関する話題が中心でした。統計数理研究所、国立遺伝学研究所、理化学研究所の演者らにより、ゲノム多様性
データの処理に関する基礎的話題および実際の研究現場における話題や、ゲノム多様性データから集団遺伝の理論とベReport「R-Training: Introduction to R statistics」開催報告
平成 26 年 12月15‑17日の三日間に渡り、フィリピン・ロス

R statistics」を開催しました。本ワークショップでは、フィリピ

バニョスにおいて、統計数理研究所リスク解析戦略研究センにおける土地利用や森林資源管理に携わる研究者、技ンターとInstitute of Biological Sciences（フィリピン大学ロ

術者、および学生を対象として、統計ソフトRによるデータハ

スバニョス校）の共催により、「R-Training: Introduction to

ンドリングおよび基礎統計解析手法の習得を目的としました。
一日目は、Rの導入作業、Rによる外部ファイルの入出力
およびデータ編集について、二日目は、グラフ作成、要約統
計量の計算、および単回帰・重回帰分析について、三日目
は、非線形回帰モデル、および AICを使ったモデル選択に関する講義と演習を行いました。最終日の閉会式では、修了証書の授与が行われました。
ほとんどの参加者が R 未経験という状況でしたが、参加者たちは非常に熱心で、ワークショップ中にも知識と技術の06

The Institute of Statistical Mathematics

目覚しい向上が見て取れました。次回開催に対する要望もワークショップには総勢３７名が参加しました。ワークショッ

あり、参加者・開催者共に実りの多いワークショップとなりま

プの運 営は、Dr. Damasa B. Magcale-Macandogらが担した。

当しました。ここに、感謝の意を表します。

（吉本 敦）

研究教育活動Report2014年11月―12月の公開講座実施状況
平 成 26 年 度 後 期 公 開 講 座は11月11日（火）の「欠 測れました。成城大学の塚原英敦教授、日本銀行金融研究

データの統計解析：理論と応用」で始まりました。当研究所所の吉羽要直氏、そして当研究所の加藤昇吾准教授が、の野間久史助教が、初級〜中級程度の確率・統計に関す複数の確率変数間の依存構造を理解するためのツールとしる知識を前提として、調査・実験における欠測データの扱いて近年着目されているコピュラの、理論と応用について、概方と具体的な統計解析の方法について、独りで解説をしまし

説をしました。会場は満席となりました。

（情報資源室）た。静岡、兵庫、岡山など、首都圏の他からの参加が多数ありました。
次の講座は12月9日（火）と10日（水）の「バイオイメージ
データ解析」でした。1日目は、当研究所の徳永旭将特任助
教と吉田亮准教授が、午前と午後に分けて講義をしました。
2日目の講義も１日目同様 2コマで、午前は金沢大学の広瀬修助教が、午後は講師全員 3 名が入れ替わり、担当しました。受講生の中には、宮城、石川、兵庫などの遠方からの参加者が含まれていました。
続いて12月19日（金）に「コピュラの理論と応用」が行わ

公開講座「バイオイメージデータ解析」資料より

Report

統計数理セミナー実施報告（2014年11月〜2015年1月）
毎週水曜 16 時から所内研究教育職員及び外部の方が 1 人 40 分ずつ、1日に2 人の講演を行っています。11月から2015
年 1月のセミナーは以下の通り行われました。
日 程
2014年11月12日

氏 名
栗木

哲

タイトル
時空間スキャン統計量の多重性調整p値の正確計算

11月12日芝井 清久

アジア・太平洋価値観国際比較調査による南北ベトナム人の比較分析

11月19日園 信太郎

ベイズ統計学の基礎をなす確率概念

11月19日庄

Inverting earthquake fault geometry through seismicity triggering

11月26日福水 健次

無限次元指数分布族を用いたスコアマッチングによるノンパラメ
トリック推定11月26日小林

高次元データにおける近傍構造の統計解析

12月3日

吉田 亮

建倉
景

Repulsive parallel MCMC algorithm for discovering diverse motifs from large sequence sets

12月3日

斎藤 正也

連結SIRモデルのパラメータ推定

12月10日吉野 諒三アジア・太平洋価値観国際比較̶普遍的価値観に基づく信頼醸成へ

12月10日鈴木 香寿恵ISM-確率台風モデル開発の紹介と流跡線解析モデルによる大気移流のデモンストレーション

12月17日松井 知子

ガウス過程状態空間モデルによる音楽感情認識

12月17日柏木 宣久

環境データの統計解析

2015年 1月14日

金藤 浩司

Web of Scienceのデータを活用した研究評価と研究戦略についての一考察

1月14日

椿

計数値データと品質管理

1月28日

藤澤 洋徳

適応的な主成分スコアに基づくスパース回帰

1月28日

伊藤

ヒルベルト空間における局所リプシッツ最適化

広計
聡

セミナーの開催予定はホームページにてご案内しています。http://www.ism.ac.jp/（メディア開発室）
The Institute of Statistical Mathematics

07

Report2014年公開講演会「ビッグデータを活かせ：スパコン利用の最前線」
統計数理研究所では毎年、教育・文化週間（11月1日〜7日）関連行事の一つとして公開講演会を開催しています。この講演会では、一般の方を対象に、統計数理に関連したテーマをわかりやすくお話ししています。今回のテーマ
はスーパーコンピュータで、11月4日の午後、統数研 2 階の大会議室で開催されました。このテーマは、本年度から統数研で3 台の新しいスーパーコンピュータ（愛称 A、I、C、
写真右）が稼働し始めたことを記念したものです。以下のようなプログラムで開催されました。
０．開会挨拶 樋口知之所長
１．スーパーコンピュータと統計数理研究所
中野純司（統計数理研究所 モデリング研究系教授／統計科学技術センター長）
２．Extreme Big Data: 次世代ビッグデータとスーパーコトによると「スーパーコンピュータがいろいろな分野に活用されていることがよく理解できた」「大変分かりやすく、理解がンピューティングの必然的な統合深まった。また来たい。」「今後の研究方向・内容の説明を

松岡 聡（東京工業大学 学術国際情報センター 教授）

いただけ、大変有益だった。」など、多くの好意的意見をい３．シミュレーションで手術をアシスト
大島まり（東京大学 大学院情報学環・生産技術研究所 教授）

ただきました。
また、統数研のスーパーコンピュータ見学ツアーも同時に行われ、こちらもご好評をいただきました。

（中野純司）

司会 川崎能典（統計数理研究所 統計科学技術センター）
最初の講演は、スーパーコンピュータの基礎的な話、統計数理研究所におけるこれまでのスパコン導入の歴史、そして最新の3 台のスパコンについての概略を説明しました。
2 番目の講演では東工大のスパコンTSUBAMEの話題を中心にスパコンの過去・現在・未来が語られました。
3 番目の講演ではスパコンの応用分野のひとつである医療分野のコンピュータシミュレーション、特に複雑な形状の血管内の血流シミュレーションが説明されました。
136 名の方が熱心に聴講されましたが、講演後のアンケー

統数研トピックス
Report「日本人の国民性第13次全国調査」記者発表
2014 年 10月30日（木）午後 2 時から、文部科学省記者会室において、日本人の国民性第 13 次全国調査の結果を発表いたしました。当日夕刻から、ネットニュースで流れ始め、翌日の新聞各紙朝刊の紙面を飾りました。新聞コラムも含めると2 週間ほどで68 件以上の参照がありました。
日本人の国民性調査は、統計数理研究所によって1953
年（昭和 28 年）に第１次調査が行われ、以降 5 年ごとに継続して実施、最初の調査から60 年たった一昨年（2013 年）
秋に、第 13 次調査を実施しました。

08

The Institute of Statistical Mathematics

記者発表に臨む吉野教授、中村教授、前田准教授（左から）調査結果の主なポイントは以下の通りです。日本人の長穏無事に」が拡大、再び楽観に転じ始めた将来の見通し。

所として「礼儀正しい」「親切」が過去最高、もう一度生ま

記者発表の内容と集計結果は、下記の統計数理研究所

れかわるとしたら「日本」に、生活水準〜日本を再評価も自
身は「変らない」が最多、「努力しても報われない」が増加、蔓延する「いらいら」、若年層で「わずらわしさを避けて、平のホームページから見ることができます。
http://www.ism.ac.jp/kokuminsei/index.html（中村 隆）

Report「サイエンスアゴラ2014」に出展
平成 26 年 11月9日（日）に江東区のお台場地区で行なわ講演会は砂田利一氏（明治大学）に総合司会をお願いれた「サイエンスアゴラ2014」に、当研究所が受託している

し、3 人の研究者赤木和人氏（東北大）、若野友一郎氏（明文部科学省委託事業「数学協働プログラム」から、講演会

治大学）、大内正己氏（東京大学）にそれぞれ聞き手の寺「科学における発見、数学における発見」と展示「『創って本紫織氏（フリーランス）、横山広美氏（東京大学）、江田

動かす』生物研究̶数理科学とロボット工学からのアプロー

慧子氏（信州大学）とのかけあいでご研究を紹介して頂き、チ̶」の二つの企画を出展致しました。

最後は司会者を含めて7 人のパネルディスカッション形式で科学談義を行って頂きました。聞き手にも若手研究者を含め、「数理を中心に据えた科学の異文化交流」を目指しての企画でしたが、聴衆の方々にも楽しんでいただけたようです。
展示の方は、動物の運動のしくみを数理的に理解するためのロボットの紹介を意図したもので、東北大学石黒章夫研究室からの御協力も頂き、ヘビ型、4 足歩行型、クモヒトデ型、ヒラムシ型の4 種が実際に動く様子を来場者に見て頂きました。数学の魅力・有用性を感じて頂けたと思います。
（丸山直昌）

Reportリスーピア賞受賞
平成 26 年 11月9日(日)に江東区のお台場地区で行なわれた「サイエンスアゴラ2014」に当研究所の受託事業「数学協働プログラム」から出展した「『創って動かす』生物研
究̶数理科学とロボット工学からのアプローチ̶」が「リスーピア賞」を受賞致しました。12月24日に行なわれた表彰式には風間俊哉特任助教が出席し、JST 理事長より賞状を受けました。

（丸山直昌）

ReportSC14に出展
11月17日から20日までの4日間、アメリカ合衆国ニューオ
リンズで開催されたSC14に参加し、ブース展示を行いました。SCはスーパーコンピュータ（スパコン）に関して、ハード
ウェアから計算手法、応用までを幅広く扱う国際会議です。
スパコン関連企業や大規模計算機を所有する研究機関などによる展示も大規模に行われており、統数研も毎年ブース
出展しています。今回の展示では、今年度新しく導入した
スパコンの紹介をはじめ、スパコンを活用した研究の成果や国民性調査の紹介を行いました。特にデータ同化スーパー
The Institute of Statistical Mathematics

09コンピュータシステムには、64テラバイトの仮想共有メモリとTOP500の301 位に、データ同化スーパーコンピュータシス

いう稀少な構成ということもあって、多くの人が関心を寄せてテムが Graph500の40 位に格付けされました。今後、このよいました。この会議は、TOP500など、スパコンの速さの世

うな充実した計算機資源を活用し、挑戦的な研究課題に取界ランキングが発表される会議としても知られています。今り組んでいきたいと考えています。

（中野慎也）

回、新しい統計科学スーパーコンピュータシステム「I」がReport

大学共同利用機関シンポジウムに参加・出展
2014 年 11月22日（土）に東京国際フォーラムに於いて、

藤森俊彦教授（自然科学研究機構基礎生物学研究所）、大学共同利用機関大学共同利用機関シンポジウム2014

阿部健一教授（人間文化研究機構総合地球環境学研究「研究者に会い行こう！日本の学術研究を支える大学共同利所）、千田俊哉教授（高エネルギー加速器研究機構物質構用機関の研究者博覧会」が開催されました。プログラムは、造科学研究所）、新井紀子教授（情報・システム研究機構
国立情報学研究所）によるトークセッションがメイン会場で行
われ、また、隣のブースセッション会場では、各研究所の展
示及び研究所の紹介が行われました。本研究所のブースには、研究者、一般社会人、高校生、小・中学生 約 370
名の方が訪れました。すっきりしたデザインのポスター顔写
真はなかなか好評であったのと、DVDビデオによる「研究所の紹介」の中身を一新したこと、さらに今年は、「日本人の国民性第 13 次全国調査」に関する最新の研究成果を発表することが出来て、一区切りがついた1日でもありました。
（広報室）

Report

科学技術コミュニケーションシンポジウムを開催
独立行政法人科学技術振興機構（JST）の平成 26 年

県屋代高等学校）による「中等教育におけるデータ分析人

度科学技術コミュニケーション推進事業の支援を受け、去る材育成に向けた統計カリキュラム」、宍倉 剛氏、橋本武彦

11月29日（土）及び 12月13日（土）に秋葉原 UDXに於い氏（データサイエンティスト協会）による「産業界におけるて、統計数理研究所が主催する科学技術コミュニケーショ

データ活用の実態と、データサイエンティストの定義の紹介」

ンシンポジウム「ビッグデータは社会に何をもたらすのか」をと題した講演がそれぞれ行われました。参加者は、講師・開催しました。11月29日（土）のシンポジウムテーマは「デースタッフを含め両日とも74 名で延べ 148 名有り、パネル討論

タ秘匿化の現状」で、当日展開されたプログラムは、佐藤一

では、参加者から事前に受け付けた質問に対して講演者

郎氏（国立情報学研究所）による「パーソナルデータ法制が回答する形式で進められ、フロアからの発言を含めて、度変更」、佐久間 淳氏（筑波大学）による「ビッグデータ

両日とも熱気のある活発な議論が展開されました。

解析におけるプライバシー保護の技術」、生貝直人氏（東京大学）による「ビッグデータと共同規制」と題した講演がそれぞれ行われました。一方、12月13日（土）の方は、テーマは「データ分析人材の育成」で、田村義保副所長（統計数理研究所）による「高等教育機関としての統計数理研究所」、長尾篤志氏（文部科学省）による「学習指導要領と数学教育・統計教育」、渡辺美智子氏（慶応義塾大学）
による「データサイエンス力の高い人材へのニーズと大学教育の役割」、遠藤寛和氏（熊谷女子高等学校）による「高等教育機関としての統計数理研究所」、横澤克彦氏（長野10The Institute of Statistical Mathematics（広報室）

Report

英国インペリアルカレッジ・ロンドンのGuoデータサイエンス研究所長が統計数理研究所を訪問
2014 年 12月16日（火）、英国インペリアルカレッジ・ロンと実践

新研究分野の創生データに基づく意思決定ドン（ICL=Imperial College London）のデータサイエンができる人材育成 等とも非常にマッチした内容であり、聴ス 研 究 所（DSI＝Data Science Institute）所 長 の Yike

講した所内教員や大学院生とも活発な質疑応答がなされGuo 教授が統計数理研究所を訪問されました。

ました。

所長室での会談後、会議室 3における研究活動の紹介

その後の中野教授と川崎准教授の案内による所内施設に移り、統計数理研究所からは、樋口所長、丸山副所長、見学も含め、Guo 教授は、日本国内で唯一統計数理を専松井教授、中野（純）教授、持橋准教授、有吉特任研究門として研究教育を行っている本研究所の活動や設備に員が参加。研究所の概要、データ同化研究開発センター、

高い関心を示されました。今回の Guo 教授の訪問により、サービス科学研究センター、統計的機械学習研究センター

DSIと本研究所の交流の活性化が今後期待されます。

の各センターの紹介、スーパーコンピュータシステムの紹介、（NOE 推進室）

そして、持 橋 准 教 授 による Nonparametric Bayesian
Methods in Audio and Language Processing の研究発表が行われました。
また、Guo 教授による45 分間の講演では、DSI の紹介
や研究活動、今後の動静および課題等について述べられ、DSI のヴィジョンとして「よりよい世界のため、大きな科
学イノベーションの土壌：データを作ること」とまとめられました。本研究所が目標に掲げる データ中心科学の確立

総合研究大学院大学複合科学研究科統計科学専攻関係
Report大学院説明会について
平成 26 年 11月7日（金）会議室 1において、平成 26 年
度第 2 回大学院説明会を開催し、20 名の参加者がありました。内容は、今年度の入試ガイダンス、5 年の課程と後期 3 年の課程のカリキュラムの説明、在学生による研究
テーマの紹介、修了後の進路の紹介などについての説明を行い、説明会終了後には、希望者に対して教員との面談を行いました。

The Institute of Statistical Mathematics

11

Report

大学院生の受賞
総研大生の佃康司さんが、論文「信号因子が複数存在

日本品質管理学会第 43 年研究奨励賞を受賞しました。

するシステムに対する動特性ロバストパラメータ設計」により、（企画グループ・研究支援担当）

お 知 ら せI n fo rma t i o n●公開講座
一般社会人・学生を対象に、下記の公開講座を開催します。

心に解説し、時変係数 AR/VARモデルで代表される分散（共分散）非定常なモデルにも触れる。取り扱う予定の内容は以下の通り。
［1］定常時系列モデル［2］時系列の変換とモデリング

A 統計学概論［3］状態空間モデル［4］トレンドの推定と季節調整［5］日時：5月12日（火）〜5月15日（金）10時〜16時（20時間）

制約付き最小二乗法とベイズモデリング［6］循環変動を講師：山下智志・小林景・野間久史・荻原哲平

含む時系列の季節調整［7］説明変数・干渉変数の導（統計数理研究所）
申込受付：3月30日（月）
10時〜4月6日（月）10時

入［8］時変係数 AR（VAR）モデル［9］局所定常 AR

受講料:20,000円

モデル

定員：100名（応募者多数の場合は抽選）
講義レベル：初級学部卒業程度の統計学・統計解析の知識（例えば線形回帰モデルに関する推測統計的な議論になじみが統計学の入門編として基礎的な講義です。

あるレベル）を前提に講義を行う。

これからデータ分析・モデリングを行う初心者に必要

参考書：●北川源四郎「時系列解析入門」岩波書店●姜興起「ベイズ統計データ解析」共立出版

な統計学の知識を、網羅的に紹介します。教科書的な数学による統計学の解説だけでなく、簡易なデータ分析
例を適宜示すことにより、直感的な理解を深めることを目的としています。
なお、高校数学Ⅲ程度の微積分と初歩的な線形代数の知識があることを前提としています。

C サンプリング法入門
日時：7月2日（木）〜7月3日（金）10時〜16時（10時間）
講師：中村隆（統計数理研究所）
申込受付：5月25日（月）
10時〜6月1日（月）
10時

・統計のための基礎数学

受講料:10,000円

・記述統計、確率分布

定員：60名（応募者多数の場合は抽選）

・統計的推定・統計的検定講義レベル：初級

サンプリングの入門的内容を、簡単な実習を行いなが・回帰分析と統計モデル

らやさしく紹介します。これから調査に携わろうとしてい・重回帰分析とモデル選択

る人、調査を始めて間もない人向けの講座。2日間の午・多変量解析の概要と用例

前と午後で、サンプリングの考え方と非復元・復元単純

・時系列モデルの初歩

ランダムサンプリング、層別サンプリングと2 段サンプリング
について順に学び、統計数理研究所が 1953 年（昭和B 非定常時系列解析
日時：6月11日（木）〜6月12日（金）10時〜16時（10時間）
講師：川﨑能典（統計数理研究所）、
姜興起（帯広畜産大学）
申込受付：4月20日（月）
10時〜4月27日（月）10時
受講料:10,000円

査」における実際例を知ることにより、サンプリング法の基礎固めをします。実習のために電卓を用意してきてください。

定員：70名（応募者多数の場合は抽選）

D ポアソン分布・ポアソン回帰・ポアソン過程

講義レベル：中級

日時：7月28日（火）10時〜16時（5時間）

本講座では、平均が非定常な（すなわち時間とともに変動する）時系列に対する統計的モデリングの方法を中1228 年）から5 年ごとに実施している「日本人の国民性調

The Institute of Statistical Mathematics

講師：島谷健一郎（統計数理研究所）
申込受付：6月15日（月）
10時〜6月22日（月）10時
受講料:5,000円

定員：100名（応募者多数の場合は抽選）
講義レベル：初級学部卒業程度の情報科学（特にアルゴリズム論）、及び平均、分散、標準偏差等、統計学の基礎的な概念を知っていることを前提とする。

0、1、2、… という整数しかとらないカウントデータには、
今日、いわゆる線形回帰でなくポアソン回帰という統計手

参考書：●Benjamin Fung, etc.「Introduction to Privacy-Preserving
Data Publishing」CRC Press.

お 知 ら せInformati on

法が広く使われるようになってきている。そこではポアソン分布という確率分布が基本となるが、確率・統計の授F 多変量解析法【社会調査士資格E科目対応（予定）】

業や教科書でポアソン分布を学んだとき、何らかの違和

日時：9月15日（火）〜9月18日（金）
10時〜16時（20時間）感を抱いた人は少なくないだろう。2 項分布のような直観

講師：馬場康維・清水信夫（統計数理研究所）、

的な意味付けがなく、ややこしそうな数式が天下り的に与えられ、「稀なイベントの起こる回数がこの分布に従うことが知られている」といった説明が添えられている。そ大森拓哉（多摩大学）
申込受付：7月27日（月）
10時〜8月3日（月）10時
受講料:20,000円
定員：100名（応募者多数の場合は抽選）

れっきり忘れていた確率分布が、カウントデータの解析で講義レベル：初級

突然、必要となる。そして、当然のことながら、ポアソン

多くの現象は一つの変数で観測されるものではなく多分布に関する理解なしに、どのようなカウントデータにはポ数の変数の観測によって把握される。すなわち多次元

アソン回帰が適切で、結果をどう解釈するかの判断もでのデータによって現象が表現される。得られた多次元のきず、誤用・濫用も数多くみられる。

データを用いて、数量の推測をする、判別をする、ある本講座では、（ポアソン過程として定式化される）ラン

いは尺度を作る等の手法の総称が多変量解析法であダムに発生するイベントからポアソン分布が導かれるとこる。重回帰分析、判別分析、数量化など、多変量解析

ろから始め、カウントデータ解析の基本であるポアソン回の古典的・標準的な手法の解説をする。平均、分散、帰の適用と結果の解釈に最低限必要な基本事項を解説

標準偏差等、統計学の基礎的な概念を知っていることをする。

前提とする。手法の数学的な説明よりは、用い方に重点

E ビッグデータのプライバシー保護技術
日時：9月3日（木）10時〜16時（5時間）
講師：南和宏
（統計数理研究所）
申込受付：7月13日（月）
10時〜7月21日（火）
10時をおいた解説を行うが、理解をたやすくする点から、大学初級程度の線形代数と微分積分の知識があることが望ましい。
資料を配布する。

受講料:5,000円
定員：70名（応募者多数の場合は抽選）
講義レベル：中級

ビッグデータの安全な流通と利活用には個人情報の漏

詳細は、以下のweb サイトをご覧ください。
http://www.ism.ac.jp/lectures/kouza.html（情報資源室）洩を防ぐ匿名化技術が不可欠である。本講座では、k匿名性、l- 多様性、差分プライバシー等の代表的な安全性指標とその実現手法を概説する。特に、ライフログのような個人の行動履歴データに対しては通常の匿名化技術の適用が困難であり、多次元データの匿名化に関する新しい提案手法も紹介する。
・プライバシー保護データパブリッシングとは？
・k- 匿名化・l- 多様性・t- 近似性・匿名化アルゴリズム
・差分プライバシー
・多次元データの匿名化

The Institute of Statistical Mathematics

13

共同利用
Report

平成26年度共同利用公募追加採択課題【共同利用登録】3 件
研究課題名
異質的主体マクロモデルの構造推定

分野
e7

砂川

研究代表者（所属）武貴（東京大学・講師）
穣（名古屋大学・研究員）

e3

臨床・全ゲノム・オミックスのビッグデータの解析に基づく疾患の原因探索・亜病態分類とリスク予測

西野

e7

Relationship Between Term Structure of Local Currency
Sovereign Bond Yield and Term Structure of Sovereign
CDS Spread鶴田 大（一橋大学・大学院生 博士課程）【一般研究 1】1 件
研究課題名
融合研究プロジェクトにおけるNGSデータ解析

分野
h3

研究代表者（所属）堀内 陽子（情報･システム研究機構・特任研究員）（企画グループ・研究支援担当）

外部資金・研究員等の受入れ
Report

受託研究の受入れ
研究題目

研究期間

研究経費（円）受入担当研究教育職員

国立大学法人東京大学大学院
医学系研究科研究科長
宮園 浩平

委託者の名称

感染症対策における政策判断のための数理モデル研究基盤の構築と発展

H26.4.1〜
H27.3.31

1,000,000

データ同化研究開発センター
斎藤 正也 特任助教

独立行政法人
国立精神•神経医療研究センター総長 樋口 輝彦

自殺対策のための自殺死亡の地域統計の更新版作成

H26.12.1〜
H27.3.31

990,000

リスク解析戦略研究センター
椿 広計 教授
竹林 由武 特任助教

独立行政法人科学技術振興機構
分任研究契約担当者
執行役 小原 英雄

多層オミックスデータから医学的に有用な知見を抽出するためのデータ解析の方法論の深化と発展・秘匿性確保の研究

H26.10.1〜
H27.3.31

936,000

数理・推論研究系
間野 修平 准教授

独立行政法人科学技術振興機構
分任研究契約担当者
執行役 小原 英雄

疫学と遺伝学の統合モデルを利用したデータ同化による流行予測

H26.10.1〜
H27.3.31

8,450,000

樋口 知之 所長

独立行政法人科学技術振興機構
分任研究契約担当者
執行役 小原 英雄

統計学に基づく情報処理に関する研究H26.10.1〜
H27.3.31

650,000

数理・推論研究系
池田 思朗 准教授

独立行政法人科学技術振興機構
分任研究契約担当者
執行役 小原 英雄

データ空間の幾何学的特徴を活用す
る解析手法と統計理論

H26.10.1〜
H27.3.31

6,058,000

数理・推論研究系
小林 景 助教（企画グループ・研究支援担当）

Report14

外来研究員の受入れ
職 名

研究題目

研究期間

Ashad Md
Alam

氏 名

ハジャ モハメド ダニッシュ科学技術大学
助教

カーネル法による次元削減のロバスト化の研究

H26.9.30 〜
H26.10.10

福水

郭北京大学 地球宇宙科学研究科
博士後期課程断層形状と震源機構の影響を含める高精度な地震活動の統計モデルの開発

H26.11.1 〜
H26.12.29

庄

Aaditya
Kumar
Ramdasカーネギーメロン大学
博士課程 5 年在学中

カーネル平均埋め込みの高次元における挙動の解析

H26.11.24〜
H26.11.28

福水

健次 教授

大谷

晋一

ジョンズホプキンス大学
応用物理学研究所 主任研究員

磁気圏MHDモデルと衛星及び地上観測との同化による電離層電気伝導度の解析

H26.11.17〜
H26.11.21

樋口

知之 所長

Mieczyslaw
Pokorski

ポーランド科学アカデミー医学研究センター 呼吸器研究部門 教授

中枢性呼吸調節における鉄分子の役割：摘出脳幹脊髄標本に動的膜電位イメージング法を適用した解析

H26.12.4 〜
H27.3.31

三分一

一村

The Institute of Statistical Mathematics

受入担当研究教育職員
健次 教授
建倉 准教授史和 准教授

職 名

研究題目

研究期間

岡村

氏 名
毅

地方行政独立法人東京都健康長寿
研究センター 研究員

時空間統計解析によるホームレスの疫学研
究H26.12.4 〜
H27.3.31

受入担当研究教育職員
椿

広計 教授

周防

節雄

兵庫県立大学 名誉教授

自殺地域統計をビジュアライズするソフトウェ
アの開発

H26.12.4 〜
H27.3.31

椿

広計 教授

河口

朋子

株式会社筑波総合研究所
主任研究員

自殺の地域統計作成のための政府統計を用いたデータベースの構築とその利活用

H26.12.4 〜
H27.3.31

椿広計 教授（企画グループ・研究支援担当）

Report

寄附金の受入れ
受入決定年月日
H26.11.5

寄附者

寄附金額（円）グラクソ・スミスクライン
株式会社

1,000,000

担当教員
野間

久史

寄附目的
統計数理研究所の研究助成のため（企画グループ・財務担当）

人事

平成26年11月1日所内異動兼務（研究教育職員）
異動内容

氏

名

兼務先

兼

荻原

哲平

リスク解析戦略研究センター

務

本 務
統計思考院助教

平成26年11月30日転出者（研究教育職員）
異動内容

氏

名

辞

河村

敏彦

職

新職名等

旧職名等

島根大学医学部附属病院医療情報部准教授

データ科学研究系データ設計グループ助教

外国人研究員（客員）氏 名
ネグリ イリア
Negri Ilia

現 職

所 属

職名

研究課題

期 間

受入教員

ベルガモ大学情報工学と数理 数理・推論研究系
客員 離散的観測に基づく H26.11.10〜 西山
方法学科准教授
統計基礎数理グル 准教授 拡散過程の統計的推 H26.12.5
イタリア共和国測ープ

陽一 准教授（企画グループ・人事担当）

Report

外国人客員紹介●Ilia Negri

客員准教授

I would like to thank Prof. Nishiyama for inviting me once again at the ISM. It is a great pleasure
for me to visit ISM and to have the opportunity to work in this institute. The ISM has always
offered to me remarkable working conditions and it is always a very stimulating and rewarding
experience to be visiting professor here. During my stay I will continue my scientificcollaboration with Prof. Nishiyama on inference for stochastic processes. In particular this year
our research will focus on moment convergence for Z-estimators.

所外誌掲載論文等
本研究所の教員、研究員、総研大（統計科学専攻）大学院生によって発表された論文等を前号に引き続き紹介します。
Cheng, M. -Y., Honda, T., Li, J. and Peng, H., Nonparametric independence screening and structure identification for
ultra-high dimensional longitudinal data, The Annals of Statistics, 42(5), 1819-1849, doi:10.1214/14-AOS1236, 2014.10
Fukaya, K., Okuda, T., Nakaoka, M. and Noda, T., Effects of spatial structure of population size on the population dynamics
of barnacles across their elevational range, Journal of Animal Ecology, 83, 1334-1343, doi:10.1111/1365-2656.12234,
2014.11Imoto, T., A generalized Conway-Maxwell-Poisson distribution which includes the negative binomial distribution,
Applied Mathematics and Computation, 247, 824, doi:10.1016/j.amc.2014.09.052, 2014.09
The Institute of Statistical Mathematics

15川畑 拓矢, 上野 玄太, 中野 慎也, 小守 信正, 増田 周平, 茂木 耕作, 三好 建正, 藤井 陽介, 多田 英夫, 吉田 隆, 第７回気象庁数値モデル研究会・第４回データ同化ワークショップの報告, 天気, 61(8), 668-670, 2014.08
Nagao, H., What is required for data assimilation that is applicable to big data in the solid Earth science?, The Proceedings of 17th International Conference on Information Fusion, 1-6, 2014.07
Brouste, A., Fukasawa, M., Hino, H., Iacus, S., Kamatani, K., Koike, Y., Masuda, H., Nomura, R., Ogihara, T., Shimizu,
Y., Uchida, M. and Yoshida, N., The yuima project: A computational framework for simulation and inference ofstochastic differential equations, Journal of Statistical Software, 57(4), 1-51, 2014.04
Ogihara, T. and Yoshida, N., Quasi-likelihood analysis for nonsynchronously observed diffusion processes, Stochastic
Processes and their Applications, 124(9), 2954-3008, doi:10.1016/j.spa.2014.03.014, 2014.09
芝井 清久, 欧州の核不拡散と東アジアの核拡散の因果関係―西ドイツをめぐる核不拡散交渉とその影響―, 国際政治, 180,2015.03
Uchida, M. and Yoshida, N., Adaptive Bayes type estimators of ergodic diffusion processes from discrete observations,
Statistical Inference for Stochastic Processes, 17(2), 181-219, doi:10.1007/s11203-014-9095-4, 2014.07（情報資源室）刊行物
ReportResearch Memorandum（2014.11〜2015.1）
No.1189: Ikeda, S., A min-max problem for prediction
No.1190: Takahashi, J., Yamashita, S., Imputing missing values using the k-NN method for extremely large-scale financial（メディア開発室）

Report

統計数理研究所調査研究リポート
No.115: 朴 堯星、
土屋 隆裕、
多摩地域住民意識調査―八王子市郵送調査(2014)― (2014.9)（メディア開発室）

Report

統計数理 第62巻 第2号
特集「生体高分子の揺らぎとダイナミクス̶シミュレーションと実験の統計解析̶」「特集 生体高分子の揺らぎとダイナミクス̶シミュレーションと実験の統計解析̶」について
伊庭 幸人、藤崎 弘士、松永 康佑 …………………………………………………………………………………163
相関から眺める生体分子運動の解析［研究詳解］
桜庭 俊 ………………………………………………………………………………………………………………171
生体分子シミュレーションの摂動解析［研究詳解］小山 洋平 ……………………………………………………………………………………………………………185
タンパク質分子の構造ダイナミクス：ウェーブレット変換による解析［研究詳解］
鎌田 真由美、戸田 幹人 ……………………………………………………………………………………………203
生体高分子系の緩和モード解析［研究詳解］高野 宏 ………………………………………………………………………………………………………………221
独立成分分析 tICAでタンパク質の複雑な運動を解きほぐす［研究詳解］
渕上 壮太郎 ……………………………………………………………………………………………………………243
１分子計測実験から分子状態を識別する統計的データ解析［研究詳解］
岡本 憲二 ……………………………………………………………………………………………………………257
カスケード型超並列シミュレーションによるタンパク質構造遷移のパスウェイ探索［研究詳解］
西原 泰孝、原田 隆平、北尾 彰朗 …………………………………………………………………………………273
ストリング法によるタンパク質構造変化解析［研究詳解］
松永 康佑 ……………………………………………………………………………………………………………285
パスサンプリングを使った分子動力学とベイズ推定［研究詳解］
藤崎 弘士 ……………………………………………………………………………………………………………301

16

The Institute of Statistical Mathematics

Individual Participant Dataに基づくメタアナリシス［研究詳解］
野間 久史 ……………………………………………………………………………………………………………313

ReportAnnals of the Institute of Statistical Mathematics
Volume 67, Number 1 (February 2015)
Benjamin Colling, Cédric Heuchenne, Rawane Samb and Ingrid Van Keilegom
Estimation of the error density in a semiparametric transformation model ……………………………………………1
Krzysztof Podgórski and Jonas Wallin
Maximizing leave-one-out likelihood for the location parameter of unbounded densities ……………………………19
Kentaro Akashi and Naoto Kunitomo
The limited information maximum likelihood approach to dynamic panel structural equation models ……………39
Nenghui Kuang and Huantian Xie
Maximum likelihood estimator for the sub-fractional Brownian motion approximated by a random walk …………75
Guang Cheng, Hao Helen Zhang and Zuofeng Shang
Sparse and efficient estimation for partial spline models with increasing dimension …………………………………93
Aerambamoorthy Thavaneswaran, Nalini Ravishanker and You Liang
Generalized duration models and optimal estimation using estimating functions ……………………………………129
Ting-Li Chen
On the convergence and consistency of the blurring mean-shift process ………………………………………………157
T. Hotz and S. Huckemann
Intrinsic means on the circle: uniqueness, locus and asymptotics ………………………………………………………177ˇˇ
A. Elijio and V. Cekanavicius
Compound Poisson approximation to weighted sums of symmetric discrete variables ………………………………195（メディア開発室）

The Institute of Statistical Mathematics

17

コ ラム

ColumnNO.

127

走るということ：クーパーって知ってますか？
上野 玄太

モデリング研究系

その日の仕事を終えて運動着に着替え、市役所前の横断歩道を渡ると、昭和記念公園の外周・8キロのジョグのスタートである。これを書いている12月は手袋がいる。
速く走るには、つま先で走ればい
いのだ。そうわかったのは保育園のころである。運動会も手つなぎ鬼もハンカチ落としもそれで十分だった。

雨上がりのルンド駅にてしかし小学校に入ってサッカーを始めると、それだけでは足りないことが分かってきた。15分ハー

心拍数が毎分

フの試合でも、ずっとつま先で走り続けるのは不可能だった

130 回を超えないのである。

というのは本当にアルヘルゴナ教会前にて

そんなサッカー部の練習メニューに、「クーパー」というもゆっくりで、厳密に守って鴨川でウォーキングの人に抜かれのがあった。
「クーパー30 分」などと使う。あらかじめ決めた

たこともある。途中で給水もするし、心拍の上限設定があ

時間をずっと走り続けるという練習である。
「クーパー30 分」

るために息が上がることもないので、きつい練習ではない。

ならば 30 分間走り続ける。30 分間校庭のトラックを周回す

LSDを行うことで、今までの走り方よりも楽な走り方があるこるのだが、これはきつい練習だった。小学生の心理として、とに気づかされた。正しい姿勢を維持し楽な呼吸を確保し1 周 100mのトラックを何周もし、ときに周回遅れが起こる状

たうえで、正しい動作を繰り返すことが大事だ。
「つま先」

況は闘争心をかき立てる。いきおい、スタート直後に得意の以来の開眼だった。

つま先で飛ばして大幅リードを狙ってしまうのだが、2 周目く

本来は「クーパー」の小学生時代にLSDの効果が得ら

らいでばて、以降ペースは上がらず抜かれていく。そう苦れたはずであった。しかし、「走る」ことと「速く走る」ことがしんでからは、30 分間の持久走なのだから、初めから飛ば同義で闘争心溢れる小学生に、ゆっくりしたペースで走る

さず先頭集団についていくのだ、早く飛ばしても後で抜かことの意義を理解させるには情報が少なすぎた。練習時間

れるぞ、飛ばすのはラスト30 秒だ、と心を改めた。

中には水を飲んではいけない時代である。

誤りである。これは「クーパー」のあるべき姿ではない。と市役所前をスタートした昭和記念公園の外周コースのお気づいたのは、大学院を終えて統数研でポスドクを始めた

気に入りの景色は、昭島口を過ぎて公園を左手に見ながらころである。港区スポーツセンターでエアロビクスに興味をの左カーブである。都道 153 号線の道なりに、夜の中に青持ったおり、エアロビクスを有酸素運動のプログラムとして梅線の明かりが見えてくる。暗い中のことなので、写真でお提唱したケネス・クーパー博士の名をふと思い出した。そ見せするのは技術的に難しい。そのかわりに示すのは、海ういえばそんな名前のきつい練習があったなあ、とぼんやり

外出張先で早朝や昼休みに走ったときに出会った風景で思った直後にあっと声が出た。まさにその博士の名前からある。
走る時間のある生活では、研究活動がうまく流れている。

来ているに違いない。
そのとき有酸素運動の練習法の代表格は、大学のス
キー 部 で 経 験した LSD（long slow distance）である。

走り終わって高松駅に向かうのは、明日への期待を予感させる至福の時間である。

LSD とは、ゆっくりしたペースで長時間走る練習法である。

統計数理研究所ニュース No.127（平成 27 年 2 月 25 日発行）
発行：大学共同利用機関法人 情報・システム研究機構
統計数理研究所 広報委員会〒190-8562 東京都立川市緑町 10-3
TEL 050-5533-8500（代） http://www.ism.ac.jp/

制作：統計科学技術センター
送付先変更は、Fax（042-526-4334）
にてお知らせ下さい。
本誌に対するご感想・ご意見はism̲news@ism.ac.jp宛にお送り下さい。