\section{はじめに}
感染症の流行は，毎年，百万人を越える患者を出しており，重要な国家的課題となっている\cite{国立感染症研究所2006}．特に，インフルエンザは事前に適切なワクチンを準備することにより，重篤な状態を避けることが可能なため，感染状態の把握は各国における重要なミッションとなっている\cite{Ferguson2005}．この把握は\textbf{インフルエンザ・サーベイランス}と呼ばれ，膨大なコストをかけて調査・集計が行われてきた．本邦においてもインフルエンザが流行したことによって総死亡がどの程度増加したかを示す推定値（{\bf超過死亡概念による死者数}）は毎年1万人を超えており\cite{大日2003}，国立感染症研究所を中心にインフルエンザ・サーベイランスが実施され，その結果はウェブでも閲覧することができる\footnote{https://hasseidoko.mhlw.go.jp/Hasseidoko/Levelmap/flu/index.html}．しかし，これらの従来型の集計方式は，集計に時間がかかり，また，過疎部における収集が困難だという問題が指摘されてきた\footnote{http://sankei.jp.msn.com/life/news/110112/bdy11011222430044-n1.htm}．このような背景のもと，近年，ウェブを用いた感染症サーベイランスに注目が集まっている．これらは現行の調査法と比べて，次のような利点がある．\begin{enumerate}\item{\bf大規模}：例えば，日本語単語「インフルエンザ」を含んだTwitter上での発言は平均1,000発言/日を超えている（2008年11月）．このデータのボリュームは，これまでの調査手法，例えば，本邦における医療機関の定点観測の集計を圧倒する大規模な情報収集を可能とする．\item{\bf即時性}：ユーザの情報を直接収集するため，これまでにない早い速度での情報収集が可能である．早期発見が重視される感染症の流行予測においては即時性が極めて重要な性質である．\end{enumerate}以上のように，ウェブを用いた手法は，感染症サーベイランスと相性が高い．ウェブを用いた手法は，ウェブのどのようなサービスを材料にするかで，様々なバリエーションがあるが，本研究では近年急速に広まりつつあるソーシャルメディアのひとつであるTwitterに注目する．しかしながら，実際にTwitterからインフルエンザに関する情報を収集するのは容易ではない．例えば，単語「インフルエンザ」を含む発言を収集すると，以下のような発言を抽出してしまう：\begin{enumerate}\itemカンボジアで鳥インフルエンザのヒト感染例、6歳女児が死亡（インフルエンザに関するニュース）\itemインフルエンザ怖いので予防注射してきました（インフルエンザ予防に関する発言）\itemやっと．．．インフルエンザが治った！（インフルエンザ完治後の発言）\end{enumerate}上記の例のように，単純な単語の集計では，実際に発言者がインフルエンザにかかっている本人（本稿では，{\bf当事者}と呼ぶ）かどうかが区別されない．本研究では，これを文書分類の一種とみなして，SupportVectorMachine(SVM)\cite{Vapnik1999}を用いた分類器を用いて解決する．さらに，この当事者を区別できたとしても，はたして一般の人々のつぶやきが正確にインフルエンザの流行を反映しているのかという情報の正確性の問題が残る．例えば，インフルエンザにかかった人間が，常にその病態をソーシャルメディアでつぶやくとは限らない．また，つぶやくとしても時間のずれがあるかもしれない．このように，不正確なセンサーとしてソーシャルメディアは機能していると考えられる．この不正確性は医師の診断をベースに集計する従来型のサーベイランスとの大きな違いである．実際に実験結果では，人々は流行前に過敏に反応し，流行後は反応が鈍る傾向があることが確認された．すなわち，ウェブ情報をリソースとした場合，現実の流行よりも前倒しに流行を検出してしまう恐れがある．本研究では，この時間のずれを吸収するために，感染症モデル\cite{Kermack1927}を適応し補正を行う．本論文のポイントは次の2点である：\begin{enumerate}\itemソーシャルメディアの情報はノイズを含んでいる．よって，文章分類手法にてこれを解決する．\itemソーシャルメディアのインフルエンザ報告は不正確である．これにより生じる時間的なずれを補正するためのモデルを提案する．\end{enumerate}本稿の構成は，以下のとおりである．2節では，関連研究を紹介する．3節では，構築したコーパスについて紹介する．4節では，提案する手法／モデルについて説明する．5節では，実験について報告する．6節に結論を述べる．
\section{関連研究}
これまで，インフルエンザの流行予測は，政府主導のトップダウンな集計が中心であったが（2.1節），現在では，ICT技術を用いた大規模な調査方法が検討されている（2.2節）．近年では，検索クエリやソーシャルメディアなどウェブの情報を利用した研究が行われている（2.3節）．\subsection{現在行われている調査方法}インフルエンザの流行に対しては，事前の十分なワクチン準備が必須の対応となるため，世界各国で共通の課題として対策が行われてきた．このため，多くの国で，自国の感染症調査のための機関が組織されている．例えば，アメリカではCentersforDiseaseControlandPrevention(CDC)\footnote{http://www.cdc.gov/}，EUはEUInfluenzaSurveillanceScheme(EISS)\footnote{http://ecdc.europa.eu/en/activities/surveillance/EISN/Pages/index.aspx}を運営している．本邦では，国立感染症研究所(InfectionDiseaseSurveillanceCenter;IDSC)の感染症情報センター\footnote{http://idsc.nih.go.jp/}がこの役目を担っている．これらの機関は主としてウイルスの遺伝子解析，国民の抗体保有状況や協力体制にある医療機関からの報告に頼って集計を行っている．例えば，IDSCは全国約5,000の診療所から情報を集め，集計結果を発表している．ただし，この集計には時間がかかるため，1週間前の流行状況が発表される．この遅延のため，発表時にはすでに流行入りしている可能性があり，より迅速な情報収集への期待が高まっている．\subsection{新しい調査方法}より早くインフルエンザの流行を捉えるため様々な手法が提案されている．Espinoetal.\cite{Espino2003}は電話トリアージ・アドバイス（電話を通じて医療アドバイスを行う公共サービス）に注目し，一日の電話コールの回数がインフルエンザの流行と相関していることを明らかにした．後の追試でもその有効性は確かめられている\cite{Yih2009}．Magruder\citeyear{Magruder2003}は，インフルエンザ市販薬の販売量(over-the-counterdrugsales)に注目した．ただし，本邦ではインフルエンザ薬の購入には処方が必要となっているように，この手法は万国で有効な手段ではない．このため，近年，薬事法など国ごとの制度に依存しない手法としてウェブが注目されている．\subsubsection{ログ・ベースの手法}近年，もっとも注目されている手法はGoogleのウェブ検索クエリを用いた手法\cite{Ginsberg2009}である．彼らはインフルエンザ流行と相関のある検索クエリ（相関係数の高い上位50語）を調査し，それらをモニタリングすることで，インフルエンザの予測を行っている．彼らの予測は，アメリカのCDC報告との相関係数0.97（最小値0.92；最大値0.99）という高い精度を報告している．同様のアイデアにもとづいた研究は他にも報告されている．例えば，Polgreenetal.はYahoo!のクエリを用いて同精度の予測を行った\cite{Polgreen2009}．Hulthetal.はスイス国内のウェブ検索会社のクエリを用いた\cite{Hulth2009}．Johnsonetal.は健康情報に関するウェブサイト「HealthLink」の閲覧のアクセスログを用いた\cite{Johnson2004}．これらは，それぞれ異なった情報源を用いているものの，患者の行動を直接調査するという観点からは同様のアプローチであるとみなせる．このアプローチは，多くの情報を即時的に収集することができるが，これを実現可能なのは，サービス・プロバイダのみという問題が残る．例えば，ウェブ検索クエリをリアルタイムに利用できる機関はGoogle，Yahoo!やMicrosoftといったいくつかのサービス・プロバイダに限定されてしまう．\subsubsection{ソーシャルメディアベースの手法}前述したように，ログ・ベースの手法はサービス・プロバイダに依存してしまうため，より利用が容易であるソーシャルメディアを情報源とした手法も2010年から開始されている．その中でも特にTwitterは，1.2億ユーザが参加しており，550万ツイートが毎日やりとりされている（2011年3月現在）．このデータ量と即時性はクエリログに匹敵するため，感染症の把握について多くの研究が行われてきた\cite{Lampos2010,Culotta2010,Paul2011,Aramaki2011,Tanida2011}．これらの研究は，表\ref{t_relwork}のように，単語選択，発言の分類などの観点から分類できる．\begin{table}[b]\caption{ソーシャルメディア・ベースの疾患サーベイランス研究の分類}\label{t_relwork}\input{04table01.txt}\end{table}まず，この単語選択の問題に対しては，経験則で「風邪」や「インフルエンザ」などのキーワードとなる単語を選択し，その頻度を集計することが考えられ，これまでの多くの先行研究はこの手法をとっている\cite{Culotta2010,Aramaki2011}．これに対し，最近では統計的に単語の選択を行った手法が模索されている．例えば，L1正規化を用いて単語の次元を圧縮する方法\cite{Lampos2010}，疾患をある種のトピックとみたてトピックモデルを用いる方法\cite{Paul2011}や，素性選択の手法を適応する研究\cite{Tanida2011}が試みられた．本研究では，これらの最新の単語選択手法を用いず，経験的に語を決めているが，単語選択と組み合わせることで，さらに精度を高めることが可能であり，今後の課題としたい．次に，いかに罹患した発言を集計するかという問題（分類問題）がある．この問題に対しては，高橋ら\citeyear{Takahashi2011}がBoostingを用いた文章分類，Aramakiet~al.\citeyear{Aramaki2011}がSVMによる分類手法を提案している．本研究では後者を用いた．以上の2つの単語選択問題と分類問題が，これまでの研究で扱われてきた主な問題であった．これに加え，本研究では，1節にて指摘したように，そもそもウェブから得られる情報は不正確である可能性に注目し，これを補正するために感染症モデルを用いる．
\section{インフルエンザ・コーパス}
実際に発言者がインフルエンザにかかっている当事者かどうかを分類するためには，コーパスが必要となる．本節ではこのコーパスの構築方法とその統計について述べる．まず，材料として，2008年11月から2010年7月にかけてTwitterAPIを用いて30億発言を収集した．次に，そこから「インフルエンザ」という語を含む発言を抽出した．この結果40万発言が得られ，データを次の2つの領域に分割した．\begin{enumerate}\item{\bfトレーニング・データ}：2008年11月から無作為に5,000発言．それぞれの発言について人手で，発言者が実際にインフルエンザにかかっているかどうか（当事者かどうか）を判定した．判定の結果，発言者が実際にインフルエンザにかかっている場合をポジティブな発言とみなし（本稿では{\bf陽性発言}と呼ぶ），逆にインフルエンザにかかっていない場合をネガティブな発言（{\bf陰性発言}）とみなす．\item{\bfテスト・データ}：残りのデータ．後述する実験に用いる．\end{enumerate}陽性発言と陰性発言の区別は，以下の除外基準に照らし1つでも該当するものがあれば陰性発言とみなした．\begin{enumerate}\item{\bf非当事者}：発言者（または，発言者と同一都道府県近郊の人間）のインフルエンザについてのみ扱い，それ以外の発言は除外する（陰性発言とみなす）．居住地が正確に分からない場合は陰性発言とみなす．例えば，「インフルエンザが実家で流行っている」では，「実家」の所在が不正確であるので，陰性発言とみなす．\item{\bf不適切な時制}：現在または近い過去のインフルエンザのみ扱い，それ以外の発言は除外する．ここでいう「近い過去」とは24時間以内とする．例えば，「{\bf昨年}はひどいインフルエンザで参加できなかった」は，陰性発言とみなす．\item{\bf否定／不適切なモダリティ}：「インフルエンザでなかった」等の否定の表現（ネゲーション）は除外する．また，疑問文や「かもしれない」と不確定な発言は基準を満たさないものとする．ここでいうモダリティとは狭義のモデリティの他に法（仮定法），表現類型（疑問文，命令文），価値判断（必要，許可）をも含む．例えば，「インフルにでもかかってしまえ」（命令文）や，「インフルエンザだとしても，熱がないのが不思議」（仮定法）は陰性発言とみなす．\item{\bf比喩的表現}：実際の疾患としての「インフルエンザ」を指さない表現は除外する．例えば，「インフルエンザ的な感染力じゃないですか！RT@■■■：朝からPCの調子が悪く困っています」におけるインフルエンザは，コンピューターウィルスをインフルエンザと比喩しており陰性発言とみなす．\end{enumerate}詳細な基準についてはアノテーション・ガイドラインを公開し，実際の用例とともに参照できるようにしている\footnote{http://www.mednlp.jp/resources/}．また，構築したコーパスについても同サイトにて配布している．なお，陽性発言／陰性発言のアノテーター間（3人）の平均一致率は84\%であった．トレーニング・データについてアノテートした結果，およそ半数の発言が陰性発言となった．すなわち，発言分類を行わない場合，半分近い発言はノイズとなっていることになる．ここで，陰性発言の一部（546発言）について，原因の割合を調査した（表\ref{t_内訳}）．陰性発言となる主要な要因は広告であることが分かる．\begin{table}[t]\caption{コーパスの分類結果}\label{t_内訳}\input{04table02.txt}\end{table}
\section{提案手法}
提案手法は，陽性発言と陰性発言を分類する手法（4.1節）と，集計の時間的なずれを補正する手法（4.2節）の2つからなる．\subsection{文章分類器を用いた発言の分類}前節にて述べたように，Twitterの発言の約半数は陰性発言となってしまう．そこで，前章にて述べたコーパスを用いて，発言を分類することでこれを解決する．このタスクは，スパムメール・フィルタリングや評価表現分析といった文の分類タスクと類似しており，本研究では，これらの研究で一般的な手法である機械学習による手法を用いる．\begin{description}\item{\bf【学習アルゴリズム】}学習アルゴリズムは先行研究(Aramakietal.2011)で使われていたSVMを用いた．先行研究では，モダリティごとにアノーテションを行い，これを複数のSVMを用いて学習する手法が提案されていたが，陰性発言を区別せず，単一のSVMで学習した場合と比べ，必ずしも精度を向上させず，また，コーパス作成が高コストになった．したがって，本研究では，陰性発言の種類を区別せず，一様に扱った．なお，SVMのカーネルは2次の多項カーネルとした．他のパラメータは用いたパッケージのデフォルト値\footnote{http://chasen.org/~taku/software/TinySVM/}を用いた．\item{\bf【素性の種類】}素性は注目している単語「インフルエンザ」または「インフル」とその周辺の形態素とした．この際，形態素解析器\footnote{http://nlp.ist.i.kyoto-u.ac.jp/index.php?JUMAN}を用い，原形，読み，品詞を用いた（表\ref{t_素性}）．\item{\bf【素性の範囲】}周辺文脈の範囲（ウィンドウ・サイズ）は予備実験によって調査した．予備実験は前節のトレーニング・データを10分割交差検定することによって行った．パフォーマンスはウィンドウ・サイズに依存する．表\ref{f_win_size.eps}にウィンドウ・サイズと精度（F値）の関係を示す．最もよいパフォーマンスは左右両方のコンテキストを6形態素までみる設定（表中の$\text{BOTH}=6$）で得られたため，以降の実験ではこの値を用いた．\end{description}\begin{table}[t]\caption{素性}\label{t_素性}\input{04table03.txt}\end{table}\begin{table}[t]\caption{ウィンドウ・サイズと精度（F値）}\label{f_win_size.eps}\input{04table04.txt}\end{table}素性の範囲を決める予備実験が示すように約74\%程度の精度で陰性発言を分類できる．この精度が，実際のインフルエンザ推定にとって有用かどうかは5節の実験にて検証する．\subsection{感染症モデルを用いた補正}前節では発言から陰性発言を取り除く手法について述べた．しかし，たとえ全てのノイズが取り除かれても，個々の発言が必ずしも正確にインフルエンザを報告しているわけではない．例えば，インフルエンザにかかってしまっても，あまりに症状がひどいとTwitterを使って発言している余裕がないことも考えられる．このような過誤が起こっているのならば，正しく発言を収集したとしても，実際の流行と集計結果との間にずれが生じてしまう．そこで，本研究では，この時間のずれを吸収するために，感染症の流行モデルである{\bfSIRモデル}\cite{Kermack1927}をソーシャルメディア用に修正したモデルを用いる．SIRモデル\cite{Kermack1927}は，感染症流行の古典的な計算モデルである．このモデルでは人間は下記の3つの状態のうちのいずれかをとっていると仮定し，各状態間の遷移確率を定義することで感染シミュレーションを行う．\begin{enumerate}\item{\bf状態S}(Susceptible):感染症に対して免疫を持たない者，\item{\bf状態I}(Infectious):発症者，\item{\bf状態R}(Recovered):感染症から回復し免疫を獲得した者．\end{enumerate}状態Sから状態Iの遷移は，健康な人に感染者が接触し感染が拡大する過程，状態Iから状態Rは感染者が回復する過程とみなせる．なお，感染症から回復して状態Rにいたった者は免疫を獲得したものとし再び感染しないものとする．また，各状態の総和，$\mathrm{S}(t)+\mathrm{I}(t)+\mathrm{R}(t)$,は人口全体であり定数とみなす．ここで，異なる状態にある個人同士の接触には，S-I,S-R,I-Rの3通りがあるが，感染が広がるのは，発症者と免疫を持たない者の接触であるS-Iだけである．すなわち，SとIが多いと感染は速く広がることになり，逆に少ないと感染は遅くなる．これを考慮して，状態Sと状態Iの量に比例した割合でSが感染するとみなす．\begin{equation}\frac{d\mathrm{S}}{dt}=-\beta\mathrm{IS}\end{equation}また，状態Iの者は単位時間あたり一定確率で回復し，状態Rに遷移すると考える．\begin{gather}\frac{d\mathrm{I}}{dt}=\beta\mathrm{IS}-\gamma\mathrm{I}\\\frac{d\mathrm{R}}{dt}=\gamma\mathrm{I}\end{gather}ここで$\beta$は，状態Sから状態Iへの遷移確率であり，感染率とみなせる．同様に，$\gamma$は状態Iから状態Rへの遷移確率であり，回復率とみなせる．それぞれの遷移確率は1日ごとの確率とする．つまり，感染者が100人いる状態で$\gamma$が0.1とすると翌日には，10人が回復し，90人が依然としてインフルエンザ状態のままであることを意味する．以上3つの式を用いることで，感染者が増えてからピークを迎えやがて収束するインフルエンザの流行がシミュレーション可能となる．前述のオリジナルのSIRモデルをソーシャルメディアに適応させるため，本研究では次の3点の修正を行った．\begin{enumerate}\itemインフルエンザに感染した患者（状態S--状態Iへの遷移）がウェブを通じて観測されるものとする．したがって，式(1)と式(2)で用いられる$\Delta$S(=$\beta$IS)が実数として，直接集計されるものとする．\item状態I--状態R間の遷移確率（変数$\gamma$）はインフルエンザの先行研究\cite{Anderson1979}において報告されていた値$(\gamma=0.38)$を用いる．\item各シーズンの開始時は，感染症はまったく流行していないものと仮定する$(\mathrm{I}(0)=0)$．\end{enumerate}
\section{実験}
\subsection{実験期間}テストデータのうち，インフルエンザ流行が起こりうる冬季（2009年と2010年の10月1日から3月31日まで）を実験対象とした．正解データとして，感染症情報センター(IDSC)から報告されているインフルエンザの定点当たりの患者数を用いた．これは週1回の報告であるので，各年24点（6ヶ月×4週）のデータポイントが含まれる\footnote{提案手法は日単位での集計が可能であるが，正解データの間隔に合わせた．}．\subsection{比較手法}提案手法の精度と比較するため以下の6つの手法を用意した．本研究の提案である(1)発言の分類と(2)SIRモデルによる補正のうち，後者に関しては，他のウェブベースの手法にも適応可能であり，Twitterだけでなく，Googleインフルトレンドを用いての検証を行った．\begin{description}\item{\bfTWITTER}:Twitter上での単語「インフルエンザ」の頻度を用いた手法．単なる単語頻度をそのまま出力する（ベースライン）．\item{\bfTWITTER+SVM}:上記手法TWITTERにSVMによる発言の分類を適応した手法．陽性と分類された発言のみが集計される．\item{\bfTWITTER+SIR}:ベースラインの手法TWITTERにSIRモデルを適応した手法．\item{\bfTWITTER+SVM+SIR}:ベースラインの手法TWITTERにSVMのよる発言の分類とSIRモデルの両方を適応した手法（提案手法）．\item{\bfGOOGLE}:Google検索サービスのクエリログを用いた手法\cite{Ginsberg2009}．Googleインフルトレンドのサイトにおいて，過去の推定値は公開されており\footnote{http://www.google.org/flutrends/jp/}，そのデータをそのまま比較手法として用いた．\item{\bfGOOGLE+SIR}:上記手法GOOGLEにSIRモデルを適応した手法．\end{description}SIRモデルの挙動はパラメータ$\gamma$に依存する．パラメータ$\gamma$はその季節に流行するインフルエンザによって異なるが，本実験では，先行研究\cite{Anderson1979}において報告されている値$(\gamma=0.38)$と，もっともパフォーマンスの高い値であった0.20を参考値として用いた．\subsection{評価}評価方法は正解データとシステムの出力値の相関係数（ピアソンの相関係数）を用いた．この評価法は，クエリログを用いた先行研究\cite{Ginsberg2009}やTwitterを用いた研究\cite{Aramaki2011}と共通である．なお，上記の各手法の出力する値のスケール（最小値や最大値）はそれぞれ異なるが，相関係数で評価することで，スケールの差異はキャンセルされる．\subsection{結果と考察}各種法の結果を表\ref{t_result}に示す．文章分類手法を用いたTWITTER+SVM（相関係数0.902）は，TWITTER（相関係数0.850）よりも精度を向上させている（統計的には有意ではない）．\begin{table}[b]\caption{各手法の推定精度}\label{t_result}\input{04table05.txt}\end{table}次に，SIRモデルは，すべての手法（TWITTER+SIRやGOOGLE+SIR）において精度を向上させている．特に，GOOGLE+SIRは$\gamma=0.2$と$\gamma=0.38$の両方で有意な精度向上となっている．他の場合も，いずれも精度を向上させており，SIRモデルは応用範囲が広い手法であると言える．さらに，文章分類手法とSIRモデルの両方を組み合わせて用いた場合(TWITTER+SVM+SIR)には，それぞれを個別に用いたよりも高い精度が得られた．特に，$\gamma=0.2$の時に相関係数0.850から0.919へ有意な$(p=0.10)$精度向上が見られた．この結果により，ソーシャルメディア上の情報をそのまま用いるのではなく，文章分類や疾患モデルと組み合わせて用いることで，さらにサーベイランス精度を向上できることが示された．本研究はインフルエンザを対象としたが，インフルエンザのみならず，熱中症，食中毒など，全貌の把握は困難であるが，各個人については把握が可能な疾患は数多く存在する．このような疾患においても，提案手法のパラメータを調整することで，対応できる可能性があることを強調したい．\subsection{誤り分析}発言の分類精度は，前節の実験で示されたようにF値74\%にすぎず（表\ref{f_win_size.eps}），4回に1回程度は誤っていることになる．どのように誤っているか表\ref{t_confusion}に示す．表に示されるように，陽性発言や否定や広告など定型的な表現はうまく判定できているが，不適切な時制や非当事者などは一部しか除外できていない．この理由は，時制の判定においては過去の程度（24時間以内かどうか）という判定が高度であること，また，当事者の判定においては，省略された主語の推定が困難なことなど，種々の原因が考えられる．このような課題はあるものの発言の分類はインフルエンザ流行推定という最終的な精度を向上できており，今後，より深い言語処理を導入することへの動機となりうる．\begin{table}[b]\caption{コンフュージョン・マトリクス}\label{t_confusion}\input{04table06.txt}\end{table}SIRモデルにおいても，改良の余地がある．推定値を可視化した結果を図\ref{f_result.eps}に示す．図に示されるように，TWITTERやGOOGLEの集計結果では，流行前に多めに推定し，流行後は少なめに推定するという傾向がある．これに対して，SIRモデルは，流行前の発言を後ろにずらして集計し，精度を高めることに成功している．しかし，SIRモデルにおいて，もっとも，高い推定値であったのは$\gamma=0.20$であり，先験的なパラメータ$(\gamma=0.38)$から外れていた．このベストパラメータ$(\gamma=0.20)$を用いた場合，相関係数0.919が得られるはずであり，これがモデルの上限値となる．今後，いかに，ベストパラメータを事前に設定するかが問題となる．なお，このパラメータ$\gamma$は，実験期間におけるインフルエンザ流行型の特質を示していると考えられ，シーズンを通して不変である可能性がある．この仮定が成り立つならば，シーズン開始時に正解データと誤差を最小化するようにし，最適なパラメータを得られる可能性があり，今後の課題としたい．\subsection{提案手法の有効性}本稿の冒頭にて，ウェブを用いた感染症サーベイランスの利点として，大規模性と即時性の2つの利点を挙げた．本節では，提案システムをウェブ・サービス\footnote{http://www.mednlp.jp/influ/}として運用した結果をもとに，これら2つの利点について実証的に議論する．\begin{figure}[b]\begin{center}\includegraphics{19-5ia6f1.eps}\end{center}\caption{各手法により推定された流行の可視化}\label{f_result.eps}\centerline{\footnotesize*Y軸は相対頻度を示す（最大値で正規化している）．X軸は時間軸を示す．}\end{figure}\begin{figure}[b]\begin{center}\includegraphics{19-5ia6f2.eps}\end{center}\caption{東日本災害時の福島県のサーベイランス結果（提案手法のサービスのスクリーン・ショット）}\label{f_result2}\footnotesize{*Y軸は既存手法については患者数（左目盛り），提案手法については陽性発言数を示す（右目盛り）．X軸は時間軸を示す．}\end{figure}\begin{enumerate}\item{既存手法との比較：}既存手法では全国約5,000の診療所から情報を集め，集計結果を発表している．ただし，その分布は均一でなく，過疎地については，少数の診療所からの情報に頼っている．このため，地震などの災害時には，流行把握が不可能な地域が生じてしまう事がある．例えば，東日本大震災時には，3月の中旬と下旬の2度にわたって，福島県からの報告が途絶えた（図2）．一方，提案手法は，GPS発信情報が付加された発言をもとに，福島県のインフルエンザ・サーベイランスが可能であった．さらに，既存の手法では，1週間前の流行しか把握できない．一方，提案手法では，つねに現状の流行をタイム・ラグなしに把握が可能である．\item{クエリ・ベースの手法（Googleインフルトレンド）との比較：}クエリ・ベースの手法であるGoogleインフルトレンドは，提案手法と同じくウェブを活用し，提案手法と同等の即時性と大規模性を持つ．ただし，根幹となるリソースの検索クエリは非公開であり，サービスを第三者が構築／拡張できないという問題がある．具体的には，Googleインフルトレンドの集計は，平成24年7月現在，国単位の集計にとどまっており，都道府県単位の集計を行っていない．この拡張を行う事はサービス・プロバイダ(Google)以外は不可能である．一方，提案手法のリソースであるTwitterは，公開されたAPIを通じて柔軟にサービスの開発が可能である．\end{enumerate}以上のように，提案手法は，既存手法と比較し，大規模性と即時性の両方において上回っている．また，クエリ・ベースの手法と比較し，同程度の大規模性と即時性を持つものの，サービスの拡張性という点では，ソーシャル・メディアが優位である．ただし，実用的には，両者のうち優位な方のみを利用するのではなく，いずれかが利用困難な場合でも，残った手法が稼働できる可能性があり，両者を併用することで頑健なサーベイランスを実現していくのが有効だと考えられる．
\section{おわりに}
本研究ではソーシャルメディアを材料に，インフルエンザの流行の推定に注目した．これまでの多くのシステムは，単純に単語の頻度情報をもとに患者の状態を調査するというものであった．しかし，Twitterの情報はノイズを含んでおり，実際に疾患にかかっていない場合の発言を収集してしまう恐れがある．また，そもそも，ウェブの情報が必ずしもインフルエンザの流行を反映しているとは限らない．これらの問題に対応するため，発言者が実際にインフルエンザにかかっているかどうかの分類を行った．また，現実とウェブ情報の時間の差を吸収するための感染症モデルを適応した．実験の結果，前者はTwitterベースのシステムの精度を向上させた．また，後者はTwitterベースのシステムのみならず，Googleインフルトレンドの精度も向上させた．また，これらの両手法は独立した現象を扱っており，両方を組み合わせて用いた場合，Twitterにおいて感染症情報センター報告の患者数と相関係数0.910という高い推定精度を示すことができた．本研究により，ソーシャルメディア上の単語頻度を単純に用いるのではなく，文章分類や疾患モデルを組み合わせることで，さらに精度の向上が可能であることが示された．\acknowledgment本研究は，JST戦略的創造研究推進事業さきがけ「情報環境と人」領域「自然言語処理による診断支援技術の開発」，科研費補助金若手研究A「表記ゆれ及びそれに類する現象の包括的言語処理に関する研究」，および科研費補助金挑戦的萌芽研究「ダミー診療録の構築および自動構造化に関する研究」による．本論文を書くにあたって有益な議論をいただいた東京大学知の構造化センター宮部真衣氏，東京大学医学部附属病院篠原恵美子氏に謹んで感謝の意を表する．\bibliographystyle{jnlpbbl_1.5}\begin{thebibliography}{}\bibitem[\protect\BCAY{Anderson\BBA\May}{Anderson\BBA\May}{1979}]{Anderson1979}Anderson,R.~M.\BBACOMMA\\BBA\May,R.~M.\BBOP1979\BBCP.\newblock\BBOQPopulationBiologyofInfectious-Diseases.\BBCQ\\newblock{\BemNature},{\Bbf280}(361).\bibitem[\protect\BCAY{Aramaki,Maskawa,\BBA\Morita}{Aramakiet~al.}{2011}]{Aramaki2011}Aramaki,E.,Maskawa,S.,\BBA\Morita,M.\BBOP2011\BBCP.\newblock\BBOQTwitterCatchesTheFlu:DetectingInfluenzaEpidemicsusingTwitter.\BBCQ\\newblockIn{\BemProceedingsoftheConferenceonEmpiricalMethodsinNaturalLanguageProcessing(EMNLP2011)},\mbox{\BPGS\1568--1576}.\bibitem[\protect\BCAY{Carbonell\BBA\Goldstein}{Carbonell\BBA\Goldstein}{1998}]{Carbonell1998}Carbonell,J.\BBACOMMA\\BBA\Goldstein,J.\BBOP1998\BBCP.\newblock\BBOQTheuseofMMR,diversity-basedrerankingforreorderingdocumentsandproducingsummaries.\BBCQ\\newblockIn{\BemProceedingsoftheSIGIR}.\bibitem[\protect\BCAY{Culotta}{Culotta}{2010}]{Culotta2010}Culotta,A.\BBOP2010\BBCP.\newblock\BBOQDetectinginfluenzaoutbreaksbyanalysingTwittermessages.\BBCQ\\newblock{\BemCoRR},{\Bbfabs/1007.4748}.\bibitem[\protect\BCAY{Espino,Hogan,\BBA\Wagner}{Espinoet~al.}{2003}]{Espino2003}Espino,J.,Hogan,W.,\BBA\Wagner,M.\BBOP2003\BBCP.\newblock\BBOQTelephonetriage:Atimelydatasourceforsurveillanceofinfluenza-likediseases.\BBCQ\\newblockIn{\BemProceedingsoftheAMIAAnnualSymposium},\mbox{\BPGS\215--219}.\bibitem[\protect\BCAY{Ferguson,Cummings,Cauchemez,Fraser,Riley,Meeyai,Iamsirithaworn,\BBA\Burke}{Fergusonet~al.}{2005}]{Ferguson2005}Ferguson,N.~M.,Cummings,D.A.~T.,Cauchemez,S.,Fraser,C.,Riley,S.,Meeyai,A.,Iamsirithaworn,S.,\BBA\Burke,D.~S.\BBOP2005\BBCP.\newblock\BBOQStrategiesforcontaininganemerginginfluenzapandemicinSoutheastAsia.\BBCQ\\newblock{\BemNature},{\Bbf437}(7056),\mbox{\BPGS\209--214}.\bibitem[\protect\BCAY{Ginsberg,Mohebbi,Patel,Brammer,Smolinski,\BBA\Brilliant}{Ginsberget~al.}{2009}]{Ginsberg2009}Ginsberg,J.,Mohebbi,M.~H.,Patel,R.~S.,Brammer,L.,Smolinski,M.~S.,\BBA\Brilliant,L.\BBOP2009\BBCP.\newblock\BBOQDetectinginfluenzaepidemicsusingsearchenginequerydata.\BBCQ\\newblock{\BemNature},{\Bbf457}(7232),\mbox{\BPGS\1012--1014}.\bibitem[\protect\BCAY{Hulth,Rydevik,\BBA\Linde}{Hulthet~al.}{2009}]{Hulth2009}Hulth,A.,Rydevik,G.,\BBA\Linde,A.\BBOP2009\BBCP.\newblock\BBOQWebQueriesasaSourceforSyndromicSurveillance.\BBCQ\\newblock{\BemPLoSONE},{\Bbf4}(2),\mbox{\BPG\e4378}.\bibitem[\protect\BCAY{Iwakura\BBA\Okamoto}{Iwakura\BBA\Okamoto}{2008}]{Iwakura2008}Iwakura,T.\BBACOMMA\\BBA\Okamoto,S.\BBOP2008\BBCP.\newblock\BBOQAfastboosting-basedlearnerforfeature-richtaggingandchunking.\BBCQ\\newblockIn{\BemProceedingsoftheTwelfthConferenceonComputationalNaturalLanguageLearning(CoNLL)},\mbox{\BPGS\17--24}.\bibitem[\protect\BCAY{Johnson,Wagner,Hogan,Chapman,Olszewski,Dowling,\BBA\Barnas}{Johnsonet~al.}{2009}]{Johnson2004}Johnson,H.,Wagner,M.,Hogan,W.,Chapman,W.,Olszewski,R.,Dowling,J.,\BBA\Barnas,G.\BBOP2009\BBCP.\newblock\BBOQAnalysisofWebaccesslogsforsurveillanceofinfluenza.\BBCQ\\newblock{\BemStudiesinHealthTechnologyandInformatics},{\Bbf107},\mbox{\BPGS\1202--1206}.\bibitem[\protect\BCAY{Kermack\BBA\McKendrick}{Kermack\BBA\McKendrick}{1927}]{Kermack1927}Kermack,W.~O.\BBACOMMA\\BBA\McKendrick,A.~G.\BBOP1927\BBCP.\newblock\BBOQAContributiontotheMathematicalTheoryofEpidemics.\BBCQ\\newblockIn{\BemTheRoyalSocietyofLondon},\mbox{\BPGS\700--721}.\bibitem[\protect\BCAY{国立感染症研究所}{国立感染症研究所}{2006}]{国立感染症研究所2006}国立感染症研究所\BBOP2006\BBCP.\newblock\Jem{インフルエンザ・パンデミックに関するQ＆A（2006.12改訂版）\inhibitglue}.\newblock国立感染症研究所感染症情報センター.\bibitem[\protect\BCAY{Kwak\BBA\Choi}{Kwak\BBA\Choi}{2002}]{Kwak2002}Kwak,N.\BBACOMMA\\BBA\Choi,C.-H.\BBOP2002\BBCP.\newblock\BBOQInputfeatureselectionforclassificationproblems.\BBCQ\\newblock{\BemNeuralNetworks},{\Bbf13}(2).\bibitem[\protect\BCAY{Lampos\BBA\Cristianini}{Lampos\BBA\Cristianini}{2010}]{Lampos2010}Lampos,V.\BBACOMMA\\BBA\Cristianini,N.\BBOP2010\BBCP.\newblock\BBOQTrackingtheflupandemicbymonitoringthesocialweb.\BBCQ\\newblockIn{\BemCognitiveInformationProcessing(CIP),20102ndInternationalWorkshopon},\mbox{\BPGS\411--416}.\bibitem[\protect\BCAY{Magruder}{Magruder}{2003}]{Magruder2003}Magruder,S.\BBOP2003\BBCP.\newblock\BBOQEvaluationofover-the-counterpharmaceuticalsalesasapossibleearlywarningindicatorofhumandisease.\BBCQ\\newblockIn{\BemJohnsHopkinsUniversityAPLTechnicalDigest(24)}.\bibitem[\protect\BCAY{大日\JBA重松\JBA谷口\JBA岡部}{大日\Jetal}{2003}]{大日2003}大日康史\JBA重松美加\JBA谷口清州\JBA岡部信彦\BBOP2003\BBCP.\newblockインフルエンザ超過死亡「感染研モデル」2002／03シーズン報告.\\newblock{\BemInfectiousAgentsSurveillanceReport},{\Bbf23}(11).\bibitem[\protect\BCAY{Paul\BBA\Dredze}{Paul\BBA\Dredze}{2011}]{Paul2011}Paul,M.~J.\BBACOMMA\\BBA\Dredze,M.\BBOP2011\BBCP.\newblock\BBOQYouAreWhatYouTweet:AnalysingTwitterforPublicHealth.\BBCQ\\newblockIn{\BemProcessingoftheFifthInternationalAAAIConferenceonWeblogsandSocialMedia(ICWSM)}.\bibitem[\protect\BCAY{Peng,Long,\BBA\Ding}{Penget~al.}{2005}]{Peng2005}Peng,H.,Long,F.,\BBA\Ding,C.\BBOP2005\BBCP.\newblock\BBOQFeatureselectionbasedonmutualinformation:Criteriaofmax-dependency,max-relevance,andmin-redundancy.\BBCQ\\newblock{\BemPatternAnalysisandMachineIntelligence},{\Bbf27}(18).\bibitem[\protect\BCAY{Polgreen,Chen,Pennock,Nelson,\BBA\Weinstein}{Polgreenet~al.}{2009}]{Polgreen2009}Polgreen,P.~M.,Chen,Y.,Pennock,D.~M.,Nelson,F.~D.,\BBA\Weinstein,R.~A.\BBOP2009\BBCP.\newblock\BBOQUsingInternetSearchesforInfluenzaSurveillance.\BBCQ\\newblock{\BemClinicalInfectiousDiseases},{\Bbf47}(11),\mbox{\BPGS\1443--1448}.\bibitem[\protect\BCAY{高橋\JBA野田}{高橋\JBA野田}{2011}]{Takahashi2011}高橋哲朗\JBA野田雄也\BBOP2011\BBCP.\newblock実世界のセンサーとしてのTwitterの可能性.\\newblock\Jem{電子情報通信学会情報・システムソサイエティ言語理解とコミュニケーション研究会}.\bibitem[\protect\BCAY{谷田\JBA荒牧\JBA佐藤\JBA吉田\JBA中川}{谷田\Jetal}{2011}]{Tanida2011}谷田和章\JBA荒牧英治\JBA佐藤一誠\JBA吉田稔\JBA中川裕志\BBOP2011\BBCP.\newblockTwitterによる風邪流行の推測.\\newblock\Jem{マイニングツールの統合と活用＆情報編纂研究会}.\bibitem[\protect\BCAY{Vapnik}{Vapnik}{1999}]{Vapnik1999}Vapnik,V.\BBOP1999\BBCP.\newblock{\BemTheNatureofStatisticalLearningTheory}.\newblockSpringer-Verlag.\bibitem[\protect\BCAY{Yih,Teates,Abrams,Kleinman,Kulldorff,Pinner,Harmon,Wang,\BBA\Platt}{Yihet~al.}{2009}]{Yih2009}Yih,W.~K.,Teates,K.~S.,Abrams,A.,Kleinman,K.,Kulldorff,M.,Pinner,R.,Harmon,R.,Wang,S.,\BBA\Platt,R.\BBOP2009\BBCP.\newblock\BBOQTelephoneTriageServiceDataforDetectionofInfluenza-LikeIllness.\BBCQ\\newblock{\BemPLoSONE},{\Bbf4}(4),\mbox{\BPG\e5260}.\end{thebibliography}\begin{biography}\bioauthor{荒牧英治}{2000年京都大学総合人間学部卒業．2002年京都大学大学院情報学研究科修士課程修了．2005年東京大学大学院情報理工系研究科博士課程修了（情報理工学博士）．以降，東京大学医学部附属病院企画情報運営部特任助教を経て，現在，東京大学知の構造化センター特任講師，科学技術振興機構さきがけ研究員（兼任）．自然言語処理の研究に従事．}\bioauthor{増川佐知子}{1989年お茶の水女子大学理学部物理学科卒業．1991年お茶の水女子大学大学院理学研究科物理学専攻修士課程修了．同年花王株式会社入社数理科学研究所配属．以降，国立福山病院附属看護学校，福山平成大学，福山市立女子短期大学非常勤講師を経て，2010年より，東京大学知の構造化センター学術支援専門職員．医療情報および自然言語処理の研究に従事．}\bioauthor{森田瑞樹}{2003年東京工業大学生命理工学部卒業．2005年東京工業大学大学院生命理工学研究科修士課程修了．2008年東京大学大学院農学生命科学研究科博士課程修了．同年東京大学大学院農学生命科学研究科特任助教．2009年医薬基盤研究所特任研究員．2012年東京大学知の構造化センター特任研究員．現在に至る．生命情報科学，医療分野における自然言語処理の研究に従事．}\end{biography}\biodate\end{document}