JP4498666B2 - 予測装置、予測プログラム、および記録媒体 - Google Patents
予測装置、予測プログラム、および記録媒体 Download PDFInfo
- Publication number
- JP4498666B2 JP4498666B2 JP2002241139A JP2002241139A JP4498666B2 JP 4498666 B2 JP4498666 B2 JP 4498666B2 JP 2002241139 A JP2002241139 A JP 2002241139A JP 2002241139 A JP2002241139 A JP 2002241139A JP 4498666 B2 JP4498666 B2 JP 4498666B2
- Authority
- JP
- Japan
- Prior art keywords
- prediction
- distribution
- equation
- logistic
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000009826 distribution Methods 0.000 claims description 64
- 238000000034 method Methods 0.000 claims description 26
- 230000001186 cumulative effect Effects 0.000 claims description 23
- 230000009385 viral infection Effects 0.000 claims description 14
- 238000000611 regression analysis Methods 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 claims description 8
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 7
- 238000009792 diffusion process Methods 0.000 description 7
- 238000007796 conventional method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 241000700605 Viruses Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- PVCRZXZVBSCCHH-UHFFFAOYSA-N ethyl n-[4-[benzyl(2-phenylethyl)amino]-2-(4-phenoxyphenyl)-1h-imidazo[4,5-c]pyridin-6-yl]carbamate Chemical compound N=1C(NC(=O)OCC)=CC=2NC(C=3C=CC(OC=4C=CC=CC=4)=CC=3)=NC=2C=1N(CC=1C=CC=CC=1)CCC1=CC=CC=C1 PVCRZXZVBSCCHH-UHFFFAOYSA-N 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
Images
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、ロジスティック曲線モデルを用いて推定するパラメータの推定値に基づいて、注目する現象の中で発生する事象の予測を行う予測方法、予測装置、予測プログラム、および予測プログラムを記録した記録媒体に関する。
【０００２】
【従来の技術】
コンピュータウィルスの感染数、ソフトウェアの信頼度、人口や生物個体数の変化、並びに新製品、技術、およびサービスの普及等の現象の予測を行うとき、従来さまざまな手法が用いられている。その中で、ロジスティック曲線等に代表される曲線モデルを用いる統計的手法が知られている。この統計的手法では、曲線モデルのパラメータを推定することにより、注目している現象の中で発生する特定の事象の予測値、発生期間、予測曲線などを得ることができる。
【０００３】
図３は、上述した各種現象を記述するモデルとして好適であると考えられているロジスティック曲線モデルを用いたときに、実績データから得られる予測曲線の例を示した説明図である。同図において、予測曲線は初期時点で求められ、この予測曲線に基づいて予測値の飽和値等の推定が行われる。
【０００４】
ロジスティック曲線モデルを用いたパラメータ推定方法について、より詳細に説明する。ロジスティック曲線モデルは、次の微分方程式（以後、ロジスティック方程式と呼ぶ）で定義される。
【数１】
【０００５】
ロジスティック方程式(1)は、厳密解
【数２】
を有する（ここで、ｋ＞０，ｍ＞０，α＞０）。ここでＬ(ｔ）は、注目する現象（モデル）で予測対象となる事象が時間（期間）ｔまでに発生する事象の累積件数（累積事象数）、ｍおよびαは事象数の実績データからそれぞれ求められるパラメータである。
【０００６】
式(1)より
Ｌ(ｔ) → ｋ （ｔ→∞） …（3）
となるので、このｋが累積事象数の飽和値を意味する。
【０００７】
パラメータｍ、ｋ、αを求めるために、式(1)を次のように書き直す。
【数３】
【０００８】
ここで改めて
【数４】
とおくと、式(4)は
Ｙ＝Ａ＋ＢＸ …（5）
と表される。ここで現実には、微分値
【数５】
を求めることはできないので、δをデータ集計期間（事象数を累積して集計する所定の期間）として、
【数６】
とおくことにより、次式(9)に示す回帰式を使用する。
Ｙn ＝Ａ＋ＢＬn …（9）
なお、Ｙn として式(8)の代わりに
【数７】
を利用することもある。
【０００９】
【数８】
【００１０】
式(9)を回帰式として回帰分析を行うことにより、Ａ，Ｂの推定値Ａ＾，Ｂ＾を求め、これらの推定値により、パラメータｍ，ｋ，αの推定値ｍ＾，ｋ＾，α＾が、それぞれ次のように得られる。
【数９】
【００１１】
しかしながら、上述した試験工程初期のデータによるパラメータ推定結果は精度が悪いといわれており、少なくとも適用するロジスティック曲線モデルの変曲点を超えた時点のデータでパラメータ推定を行うことが必要である。三觜武著「ソフトウェアの品質評価法」（日科技連、1981）によれば、
【数１０】
となれば、その時点までのデータでパラメータ推定を行うとしている。ここで、潜在する事象数の予測値をｋ＾、実績累積事象数をｙ−としている。このｋ＾は経験的または統計的に予測する。また、ｗは経験的にｗ＝０．６〜０．８の値をとるものとしている。
【００１２】
以上説明したパラメータ推定方法（以後「従来技術１」と呼ぶ）の問題点、すなわち、初期のデータによるパラメータ推定結果の精度が悪い、という問題点を解決するために、ロジスティック方程式(1)を通常とは異なる差分方程式に書き直してパラメータ推定を行う方法も提案されている（詳細については、特開２０００−１２２８６０公報を参照。以後、この技術を「従来技術２」と呼ぶ）。
【００１３】
以下、この従来技術２をロジスティック曲線モデルに適用した場合について説明する。ロジスティック方程式(1)を、δを差分間隔として次式のように差分化する。
【数１１】
【００１４】
式(16)の厳密解は、差分間隔０の連続極限（δ→０）でロジスティック方程式(1)の厳密解（式(2)参照）に一致する。さらに、
｜１−δα｜＜１
という条件の下で、
Ｌn →ｋ （ｔn → ∞） …（17）
となり、式(1)が有している性質（式(3)参照）を保存していることがわかる。
【００１５】
パラメータｍ，ｋ，αを求めるために、ｔn ＝ｎδ，δ＝１とおいて式(15)を次のように書き替える。
Ｙn ＝Ａ＋ＢＬn+1 …（18）
ここで、Ｙn ，Ａ，Ｂは、それぞれ
【数１２】
である。
【００１６】
式(18)を回帰式として回帰分析を行うことによって得られるＡ，Ｂの推定値をＡ＾，Ｂ＾とすると、パラメータｍ，ｋ，αの推定値ｍ＾，ｋ＾，α＾はそれぞれ次のように求められる。
【数１３】
【００１７】
この従来技術２の方法によれば、従来技術１よりも早期の段階で正確なパラメータ推定が可能となり、試験工程初期のデータを用いて正確な予測を行うことができる。
【００１８】
【発明が解決しようとする課題】
上述した従来技術１では、初期でのパラメータ推定精度が低く、ある程度の精度を要求するには、少なくとも適用している曲線モデルの変曲点を超えた時点までのデータが必要となる。どの時点までのデータを利用してパラメータ推定を行うかの判定基準の一つとして式(14)があるが、この式におけるｋ＾は経験に頼って予測せざるを得ないため、ｋ＾の精度そのものが低いという問題点があった。
【００１９】
また、従来技術１では、実績データと曲線モデルによる分析が一致しなかった場合、その原因が曲線モデルの選択が不適切であることによるものなのか、あるいはパラメータの推定精度の低さによるものなのかがはっきりしないという問題点があった。
【００２０】
さらに、従来技術１では、使用している方程式が決定論に従うものであるため、確率論的な議論によって分布等の統計的性質を予測することは不可能であった。すなわち、注目する現象で発生する特定の事象の予測値の分布や、予測値がある値を下回る確率（分布）や上回る確率（補分布）を求めることができなかった。
【００２１】
一方、上記従来技術２では、パラメータの推定値は精度よく求められるものの、従来技術１と同様、使用している方程式が決定論に従うものであるため、確率論的な議論によって分布等の統計的性質を予測することは不可能であった。
【００２２】
本発明は、上記に鑑みてなされたものであり、その目的は、ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定するとともに、この推定したパラメータに基づいて、前記現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測方法、予測装置、予測プログラム、および予測プログラムを記録した記録媒体を提供することにある。
【００２３】
【課題を解決するための手段】
本発明は、従来技術２におけるロジスティック差分方程式が厳密解を有するという性質を保持しながら、パラメータの一つを確率変数として扱うことによって、注目する現象を記述するモデルの中で予測対象となる事象の累積事象数の予測値について、その統計的な分布を推定すること主要な特徴とする。
【００２４】
そこで、上記目的を達成するために、本願の装置発明は、
ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定し、この推定したパラメータに基づいて前記現象の中で発生する事象の予測を行う予測装置であって、前記現象の所定の期間ごとの発生件数を記憶する事象内容記憶手段と、予測する時点と予測目標値からなる予測条件を記憶する予測条件記憶手段と、を有し、この２つの記憶手段を活用する下記に示す５つの処理手段を有する。
【００２５】
１つ目は、前記事象内容記憶手段から前記発生件数を読み出し、注目している事象について所定の時間間隔で累積件数を集計し、実績値を得る累積件数集計手段である。
【００２６】
２つ目は、前記ロジスティック曲線を表すロジスティック方程式を時間ステップで差分化したロジスティック差分方程式のうち厳密解を有するロジスティック差分方程式に前記累積件数集計手段で得られた実績値を与え、回帰分析によってロジスティック方程式の厳密解Ｌ（ｔ）＝ｋ／｛１＋ｍ＊ｅｘｐ（−αｔ）｝のパラメータｋ、ｍを求めるパラメータ推定手段である。
【００２７】
３つ目は、前記ロジスティック差分方程式の厳密解のパラメータｋ、ｍを前記求められた値とし、前記累積件数集計手段で得られた実績値を適用して、パラメータαが時間ステップに対応した確率変数列Ａを求める確率変数列算出手段である。
【００２８】
４つ目は、前記確率変数列Ａが差分ステップごとに独立であるとともに全ての差分ステップで同一な分布である独立同一分布に従うこととして、前記確率変数列Ａに基づいて最尤法を用いることによって前記独立同一分布のパラメータを求める分布係数特定手段である。
【００２９】
５つ目は、パラメータ推定手段で求められた前記パラメータｋ、ｍ、分布係数特定手段で求められた確率変数列Ａの独立同一分布パラメータ、およびロジスティック差分方程式の厳密解を用いて、前記予測する時点の予測値Ｌが前記予測目標値を満たす確率を求める分布推定手段である。
【００３０】
本願の他の装置発明は、前記予測目標値が下限の場合、前記分布推定手段は、前記予測値Ｌが前記下限よりも大きい確率を求める。
【００３１】
本願の他の装置発明は、前記予測目標値が上限の場合、前記分布推定手段は、前記予測値Ｌが前記上限よりも小さい確率を求める。
【００３２】
本願の他の装置発明は、前記予測目標値が下限と上限とからなる範囲の場合、前記分布推定手段は、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める。
【００３３】
なお、本発明の予測装置は、一つのコンピュータによって構築される場合だけでなく、各手段（の有する機能）を適宜分割して二つ以上のコンピュータによって構築される場合も含む。
【００３４】
本願の別の他の装置発明では、前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布を対象とする。
【００３５】
上記予測装置において、前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを要旨とする。
【００３６】
本発明に係わる予測プログラムは、注目する現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測プログラムを提供するものである。
【００３７】
本発明に係わる予測プログラムは、前記予測目標値が下限の場合、前記予測値Ｌが前記下限よりも大きい確率を求め、
前記予測目標値が上限の場合、前記予測値Ｌが前記上限よりも小さい確率を求め、
前記予測目標値が下限と上限とからなる範囲の場合、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める。
【００３８】
本発明に係わる予測プログラムは、前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布を対象とする。
【００３９】
上記予測プログラムにおいて、前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを要旨とする。
【００４０】
本発明に係わる予測プログラムを記録した記録媒体は、上記予測プログラムをそれぞれ記録したコンピュータ読み取り可能な記録媒体を提供するものである。
【００４１】
本発明に係わる予測プログラムを記録した記録媒体は、前記予測目標値が下限の場合、、前記予測値Ｌが前記下限よりも大きい確率を求め、
前記予測目標値が上限の場合、前記予測値Ｌが前記上限よりも小さい確率を求め、
前記予測目標値が下限と上限とからなる範囲の場合、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める。
【００４２】
本発明に係わる予測プログラムを記録した記録媒体は、前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布を対象とする。
【００４３】
上記予測プログラムを記録した記録媒体であって、前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを要旨とする。
【００４４】
請求項１０乃至１２のいずれかに記載した発明によれば、注目する現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測プログラムを、フレキシブルディスク、ＣＤ−ＲＯＭ、ＤＶＤ（Digital Versatile Disk）、光磁気ディスク、ＰＣカード等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることができる。
【００４５】
【発明の実施の形態】
次に、図面を参照して本発明の実施の形態を説明する。
【００４６】
図１は、本発明の一実施形態に係る予測装置１の構成を表すブロック図である。同図に示す予測装置１は、注目する現象で予測対象となる事象が発生した期日やその発生件数等の過去の履歴、および前記予測対象となる事象において予測すべき事象の上限や下限等の値の範囲や予測時点等に代表される予測条件（予測を行うときの条件）を入力する入力手段として、キーボードやマウス等からなる入力部１１を有する。入力部１１で入力される過去の履歴および事象の予測条件は、データベースとしての事象内容記憶部１２および予測条件記憶部１４にそれぞれ記憶される。これらの事象内容記憶部１２および予測条件記憶部１４は、それぞれ事象内容記憶手段および予測条件記憶手段の少なくとも一部を構成するものである。
【００４７】
ここで、「注目する現象」とは、例えばコンピュータウィルスの感染、ソフトウェアの信頼度、人口や生物個体数の変化、並びに新製品、技術、およびサービスの普及等である。各現象で「予測対象とする事象」を具体的に表す値は、コンピュータウィルスの感染を予測する場合にはコンピュータウィルス感染数、ソフトウェア信頼度を予測する場合には故障や欠陥（バグ）の数、人口や生物個体数を予測する場合には人口や生物個体数、新製品、技術、サービスの普及を予測する場合には普及数のことである。
【００４８】
予測装置１には、この他に、事象発生期日や事象発生件数を集計して累積データ（累積事象数）を作成する累積件数集計部１３、パラメータの推定を回帰分析によって行うパラメータ推定手段としての回帰分析部１５、パラメータの推定値を用いて予測すべき事象の統計的分布等を推定する分布推定手段としての分布推定部１６、および推定結果の出力を行うための出力手段としてディスプレイ装置等からなる出力部１７を有している。
【００４９】
予測装置１は、演算機能および制御機能を備えた中央処理装置（ＣＰＵ：Central Processing Unit）、プログラムやデータを格納する機能を備えた主記憶装置、およびハードディスクドライブ、フレキシブルディスクドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤドライブ、光磁気ディスクドライブ、ＰＣカードドライブ等の補助記憶装置を備えたコンピュータによって構成されている。
【００５０】
また、本実施形態に係る各種処理を実行するための予測プログラムは前述した主記憶装置またはハードディスクに格納されている。この予測プログラムは、フレキシブルディスク、ＣＤ−ＲＯＭ、ＤＶＤ、光磁気ディスク、ＰＣカード等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることも勿論可能である。
【００５１】
以上の構成を有する予測装置１が行う予測方法について説明する。まず、予測を行うためのアルゴリズムを説明する。本実施形態においては、ロジスティック曲線モデルを例として説明するが、本発明の予測方法がこのモデルに限定されるわけでないことは勿論である。
【００５２】
ロジスティック曲線を表す微分方程式であるロジスティック方程式
【数１４】
を差分化したロジスティック差分方程式のうち、厳密解を有するロジスティック差分方程式としては、例えば以下に示す式(23), (24), および(25)が知られている（広田良吾著「差分方程式講義−連続より離散へ」（ＳＧＣライブラリ８，サイエンス社）を参照）。
【数１５】
【００５３】
これらのロジスティック差分方程式に対して確率則を導入し、パラメータのうちの一つを確率変数として扱う確率ロジスティック差分方程式について説明する。
【００５４】
以下では、差分化した時間ステップ（差分ステップ）ごとに独立であり、かつ全ての時間ステップが同じ分布に従うような独立同一分布の確率変数列を｛Ａj ｝（ｊ：正の整数）とする。
【００５５】
式(23)に基づいた場合の確率ロジスティック差分方程式は、
【数１６】
と表される。この式(26)からも明らかなように、本実施形態においては、ロジスティック差分方程式におけるパラメータαが時間ステップ数（ｎ）に対応した確率変数Ａn+1 として扱われる。この確率ロジスティック差分方程式(26)の厳密解は、
【数１７】
である。同様に、式(24)に基づいた確率ロジスティック差分方程式は、
【数１８】
と表され、この式の厳密解は、
【数１９】
である。また、式(25)に基づく確率ロジスティック差分方程式は
【数２０】
と表され、この式の厳密解は、
【数２１】
となる。
【００５６】
通常用いる差分方程式は、元の微分方程式を差分間隔の何乗のオーダで近似するかに主眼が置かれる。したがって通常の場合には、差分化がなされた時点で、解の形状や、時間無限大の極限における累積事象数の一定値への収束性などの性質は一般に保存されない。
【００５７】
これに対して本実施形態においては、元の微分方程式が有する性質を保つ差分方程式を用いるので、通常用いられる手法より精度の高いパラメータ推定が可能となる。
【００５８】
確率変数列｛Ａj ｝は、予測対象としている事象に応じて、その事象が従う分布を決定する。一例として、式(27)において、
Ｘj ＝１−δＡj ・・・(32)
とおき、この独立同一確率変数列｛Ｘj ｝がべき関数分布に従う場合について説明する。なお、式(29), (31)を用いる場合には、Ｘj をそれぞれ
【数２２】
とおく。以上の式(32), (33), および(34)の定義によれば、確率ロジスティック差分方程式の厳密解(27), (29), および(31)をすべて同一の式
【数２３】
で表すことができる。したがって、以下の議論は上述した全ての確率ロジスティック差分方程式(26), (28), および(30)に共通である。
【００５９】
注目している事象について、所定の等間隔時間ステップごとに累積件数集計部１３で計測（累積事象数の算出）が行われ、現在時間ステップｉ番目（ｉ：正の整数）までの計測結果（実績値）が得られているものとする。この時点でｎ（＞ｉ）番目の時間ステップの累積事象数として予測される予測値を推定し、その予測値の統計的な分布を求める方法について説明する。
【００６０】
今、ある値Ｌ― を
【数２４】
とする。このとき、任意の時間ステップ数ｎにおける予測値Ｌn の値が、このＬ― よりも大きい確率（補分布）Ｐ｛Ｌn ＞Ｌ― ｝を求める。具体的な計算は次のようになる。
【数２５】
【００６１】
同様に、ある値Ｌ― を
【数２６】
とするとき、任意の時間ステップ数ｎにおける予測値Ｌn の値がＬ― よりも小さい確率（分布）Ｐ｛Ｌn ＜Ｌ― ｝は次のように計算される。
【数２７】
【００６２】
式(38)および式(41)を導出する計算では、Ｘj をべき関数分布としているので、−log Ｘj は指数分布の独立同一確率変数列であり、
【数２８】
はアーラン分布に従うことを用いた。
【００６３】
なお、以上の計算において、Ｌ―（Ｌ― )からｘ―（ｘ― )を求めるためにはｋ、ｍが必要である。そこで、これらのパラメータの推定値ｋ＾、ｍ＾を求めるために、前述した従来技術２の手法を用いる。
【００６４】
また、Ｐ｛Ｌn ＜Ｌ― ｝の計算においては、独立同一確率変数列｛Ｘj ｝がべき関数分布に従うものとしているので、−log Ｘj は指数分布の独立同一確率変数列となり、最尤法を用いることによって、べき数γが、
【数２９】
と求められる。
【００６５】
本実施形態においては、以上説明したように予測値の下限または上限を設定して分布を求める以外にも、予測値の範囲を指定して、その範囲内での分布を求めることも可能である。例えば、任意の時間ステップ数ｎでの予測値がＬ― よりも大きくＬ― よりも小さい確率Ｐ｛Ｌ―＜Ｌn ＜Ｌ― ｝を求めるときには、
【数３０】
から求めることができる。
【００６６】
図２は、予測装置１が以上説明したアルゴリズムにしたがって予測値の分布の推定を行うときの動作手順を説明するフローチャート図である。同図においては、予測装置１の処理の流れに加えて、事象内容記憶部１２ならびに予測条件記憶部１４に対するデータの流れを点線で示している。
【００６７】
まず、入力部１１で対象とする事象の発生期日、発生件数を一定の期間ごとに入力し、事象内容記憶部１２に記憶する（ステップＳ１）。
【００６８】
充分な事象数のデータが記憶された後、予測したい時点（予測値の分布が欲しい時点）と予測値の下限（Ｌ― ）、上限（Ｌ― ）、あるいは範囲（Ｌ―＜Ｌn ＜Ｌ― ）等の予測条件を入力部１１から入力し、予測条件記憶部１４に記憶する（ステップＳ３）。
【００６９】
一方、回帰分析部１５では事象内容記憶部１２で記憶したデータに基づいて累積件数集計部１３で累積された事象数を用いることにより、パラメータｋ、ｍを回帰分析によって推定し、パラメータ推定値ｋ＾、ｍ＾の算出を行う（ステップＳ５）。
【００７０】
なお、ステップＳ３およびＳ５の処理は独立なので、処理の順序を逆にしてもよいし、同時に行ってもよい。
【００７１】
以上の結果に基づいて、分析推定部１６が確率変数列｛Ａj ｝の平均値、予測対象となる事象の予測値の推定、および上記アルゴリズムにしたがって予測値の統計的分布を推定し（ステップＳ７）、この推定結果を出力部１７で表示する（ステップＳ９）。
【００７２】
以上説明した本発明の一実施形態によれば、ロジスティック方程式を差分化し、さらにその差分方程式のパラメータの一つを独立同一分布に従う確率変数として扱う確率ロジスティック差分方程式を用いることによって、注目する現象を記述するモデルのパラメータの推定に加えて、そのパラメータ推定値に基づいて予測対象となる事象の推定値の統計的分布を推定することが可能になる。
【００７３】
本実施形態は、例えばコンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデル等に適用することができ、各モデルに対応する予測値として、コンピュータウィルス感染数、故障・欠陥数（バグ数）、人口、生物個体数、並びに普及数をそれぞれ推定する場合に適用することが可能なものである。
【００７４】
なお、本発明においては、確率変数列が従う独立同一確率分布として、上述したべき関数分布以外に、例えば範囲（０，１）の一様分布を用いても同様の効果を得ることができる。
【００７５】
この場合、式(27)の厳密解において、
δＡj ＝Ｕj ・・・(45)
とおき、独立同一確率変数列｛Ｕj ｝が範囲（０，１）の一様分布に従うものとする。式(29)および式(31)の厳密解を用いるときには、それぞれ
【数３１】
とおく。
【００７６】
予測値Ｌn が式(36)で定義されるＬ― よりも大きい確率Ｐ｛Ｌn ＞Ｌ― ｝は、次のように計算される。
【数３２】
ここでｘ― の定義は式(37)と同じである。
【００７７】
また、予測値Ｌn が式(39)で定義されるＬ― よりも小さい確率Ｐ｛Ｌn ＜Ｌ― ｝は、
【数３３】
と計算される。ｘ― の定義は式(40)と同じである。
【００７８】
式(48)および式(49)を導出する計算においては、−log（１−Ｕj ）は指数分布の独立同一確率変数列であり、
【数３４】
がアーラン分布に従うことを用いた。
【００７９】
なお、時間ステップ数ｎでの予測値がＬ― よりも大きくＬ― よりも小さい確率Ｐ｛Ｌ―＜Ｌn ＜Ｌ― ｝についても、式(43)または式(44)を用いて計算することができるのは勿論である。
【００８０】
このように、本発明は、上記一実施形態同様の効果を奏するさまざまな実施の形態等を含みうるものであることはいうまでもない。
【００８１】
【発明の効果】
以上説明した本発明によれば、ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定するとともに、この推定したパラメータに基づいて、前記現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測方法、予測装置、予測プログラム、および予測プログラムを記録した記録媒体を提供することができる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る予測装置の構成を示すブロック図である。
【図２】本発明の一実施形態に係る予測方法の動作手順を示すフローチャート図である。
【図３】従来法における予測曲線および実績累積件数を表す説明図である。
【符号の説明】
１ 予測装置
１１ 入力部
１２ 事象内容記憶部
１３ 累積件数集計部
１４ 予測条件記憶部
１５ 回帰分析部
１６ 分布推定部
１７ 出力部
【発明の属する技術分野】
本発明は、ロジスティック曲線モデルを用いて推定するパラメータの推定値に基づいて、注目する現象の中で発生する事象の予測を行う予測方法、予測装置、予測プログラム、および予測プログラムを記録した記録媒体に関する。
【０００２】
【従来の技術】
コンピュータウィルスの感染数、ソフトウェアの信頼度、人口や生物個体数の変化、並びに新製品、技術、およびサービスの普及等の現象の予測を行うとき、従来さまざまな手法が用いられている。その中で、ロジスティック曲線等に代表される曲線モデルを用いる統計的手法が知られている。この統計的手法では、曲線モデルのパラメータを推定することにより、注目している現象の中で発生する特定の事象の予測値、発生期間、予測曲線などを得ることができる。
【０００３】
図３は、上述した各種現象を記述するモデルとして好適であると考えられているロジスティック曲線モデルを用いたときに、実績データから得られる予測曲線の例を示した説明図である。同図において、予測曲線は初期時点で求められ、この予測曲線に基づいて予測値の飽和値等の推定が行われる。
【０００４】
ロジスティック曲線モデルを用いたパラメータ推定方法について、より詳細に説明する。ロジスティック曲線モデルは、次の微分方程式（以後、ロジスティック方程式と呼ぶ）で定義される。
【数１】
ロジスティック方程式(1)は、厳密解
【数２】
【０００６】
式(1)より
Ｌ(ｔ) → ｋ （ｔ→∞） …（3）
となるので、このｋが累積事象数の飽和値を意味する。
【０００７】
パラメータｍ、ｋ、αを求めるために、式(1)を次のように書き直す。
【数３】
ここで改めて
【数４】
Ｙ＝Ａ＋ＢＸ …（5）
と表される。ここで現実には、微分値
【数５】
【数６】
Ｙn ＝Ａ＋ＢＬn …（9）
なお、Ｙn として式(8)の代わりに
【数７】
【０００９】
【数８】
式(9)を回帰式として回帰分析を行うことにより、Ａ，Ｂの推定値Ａ＾，Ｂ＾を求め、これらの推定値により、パラメータｍ，ｋ，αの推定値ｍ＾，ｋ＾，α＾が、それぞれ次のように得られる。
【数９】
しかしながら、上述した試験工程初期のデータによるパラメータ推定結果は精度が悪いといわれており、少なくとも適用するロジスティック曲線モデルの変曲点を超えた時点のデータでパラメータ推定を行うことが必要である。三觜武著「ソフトウェアの品質評価法」（日科技連、1981）によれば、
【数１０】
【００１２】
以上説明したパラメータ推定方法（以後「従来技術１」と呼ぶ）の問題点、すなわち、初期のデータによるパラメータ推定結果の精度が悪い、という問題点を解決するために、ロジスティック方程式(1)を通常とは異なる差分方程式に書き直してパラメータ推定を行う方法も提案されている（詳細については、特開２０００−１２２８６０公報を参照。以後、この技術を「従来技術２」と呼ぶ）。
【００１３】
以下、この従来技術２をロジスティック曲線モデルに適用した場合について説明する。ロジスティック方程式(1)を、δを差分間隔として次式のように差分化する。
【数１１】
式(16)の厳密解は、差分間隔０の連続極限（δ→０）でロジスティック方程式(1)の厳密解（式(2)参照）に一致する。さらに、
｜１−δα｜＜１
という条件の下で、
Ｌn →ｋ （ｔn → ∞） …（17）
となり、式(1)が有している性質（式(3)参照）を保存していることがわかる。
【００１５】
パラメータｍ，ｋ，αを求めるために、ｔn ＝ｎδ，δ＝１とおいて式(15)を次のように書き替える。
Ｙn ＝Ａ＋ＢＬn+1 …（18）
ここで、Ｙn ，Ａ，Ｂは、それぞれ
【数１２】
【００１６】
式(18)を回帰式として回帰分析を行うことによって得られるＡ，Ｂの推定値をＡ＾，Ｂ＾とすると、パラメータｍ，ｋ，αの推定値ｍ＾，ｋ＾，α＾はそれぞれ次のように求められる。
【数１３】
この従来技術２の方法によれば、従来技術１よりも早期の段階で正確なパラメータ推定が可能となり、試験工程初期のデータを用いて正確な予測を行うことができる。
【００１８】
【発明が解決しようとする課題】
上述した従来技術１では、初期でのパラメータ推定精度が低く、ある程度の精度を要求するには、少なくとも適用している曲線モデルの変曲点を超えた時点までのデータが必要となる。どの時点までのデータを利用してパラメータ推定を行うかの判定基準の一つとして式(14)があるが、この式におけるｋ＾は経験に頼って予測せざるを得ないため、ｋ＾の精度そのものが低いという問題点があった。
【００１９】
また、従来技術１では、実績データと曲線モデルによる分析が一致しなかった場合、その原因が曲線モデルの選択が不適切であることによるものなのか、あるいはパラメータの推定精度の低さによるものなのかがはっきりしないという問題点があった。
【００２０】
さらに、従来技術１では、使用している方程式が決定論に従うものであるため、確率論的な議論によって分布等の統計的性質を予測することは不可能であった。すなわち、注目する現象で発生する特定の事象の予測値の分布や、予測値がある値を下回る確率（分布）や上回る確率（補分布）を求めることができなかった。
【００２１】
一方、上記従来技術２では、パラメータの推定値は精度よく求められるものの、従来技術１と同様、使用している方程式が決定論に従うものであるため、確率論的な議論によって分布等の統計的性質を予測することは不可能であった。
【００２２】
本発明は、上記に鑑みてなされたものであり、その目的は、ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定するとともに、この推定したパラメータに基づいて、前記現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測方法、予測装置、予測プログラム、および予測プログラムを記録した記録媒体を提供することにある。
【００２３】
【課題を解決するための手段】
本発明は、従来技術２におけるロジスティック差分方程式が厳密解を有するという性質を保持しながら、パラメータの一つを確率変数として扱うことによって、注目する現象を記述するモデルの中で予測対象となる事象の累積事象数の予測値について、その統計的な分布を推定すること主要な特徴とする。
【００２４】
そこで、上記目的を達成するために、本願の装置発明は、
ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定し、この推定したパラメータに基づいて前記現象の中で発生する事象の予測を行う予測装置であって、前記現象の所定の期間ごとの発生件数を記憶する事象内容記憶手段と、予測する時点と予測目標値からなる予測条件を記憶する予測条件記憶手段と、を有し、この２つの記憶手段を活用する下記に示す５つの処理手段を有する。
【００２５】
１つ目は、前記事象内容記憶手段から前記発生件数を読み出し、注目している事象について所定の時間間隔で累積件数を集計し、実績値を得る累積件数集計手段である。
【００２６】
２つ目は、前記ロジスティック曲線を表すロジスティック方程式を時間ステップで差分化したロジスティック差分方程式のうち厳密解を有するロジスティック差分方程式に前記累積件数集計手段で得られた実績値を与え、回帰分析によってロジスティック方程式の厳密解Ｌ（ｔ）＝ｋ／｛１＋ｍ＊ｅｘｐ（−αｔ）｝のパラメータｋ、ｍを求めるパラメータ推定手段である。
【００２７】
３つ目は、前記ロジスティック差分方程式の厳密解のパラメータｋ、ｍを前記求められた値とし、前記累積件数集計手段で得られた実績値を適用して、パラメータαが時間ステップに対応した確率変数列Ａを求める確率変数列算出手段である。
【００２８】
４つ目は、前記確率変数列Ａが差分ステップごとに独立であるとともに全ての差分ステップで同一な分布である独立同一分布に従うこととして、前記確率変数列Ａに基づいて最尤法を用いることによって前記独立同一分布のパラメータを求める分布係数特定手段である。
【００２９】
５つ目は、パラメータ推定手段で求められた前記パラメータｋ、ｍ、分布係数特定手段で求められた確率変数列Ａの独立同一分布パラメータ、およびロジスティック差分方程式の厳密解を用いて、前記予測する時点の予測値Ｌが前記予測目標値を満たす確率を求める分布推定手段である。
【００３０】
本願の他の装置発明は、前記予測目標値が下限の場合、前記分布推定手段は、前記予測値Ｌが前記下限よりも大きい確率を求める。
【００３１】
本願の他の装置発明は、前記予測目標値が上限の場合、前記分布推定手段は、前記予測値Ｌが前記上限よりも小さい確率を求める。
【００３２】
本願の他の装置発明は、前記予測目標値が下限と上限とからなる範囲の場合、前記分布推定手段は、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める。
【００３３】
なお、本発明の予測装置は、一つのコンピュータによって構築される場合だけでなく、各手段（の有する機能）を適宜分割して二つ以上のコンピュータによって構築される場合も含む。
【００３４】
本願の別の他の装置発明では、前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布を対象とする。
【００３５】
上記予測装置において、前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを要旨とする。
【００３６】
本発明に係わる予測プログラムは、注目する現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測プログラムを提供するものである。
【００３７】
本発明に係わる予測プログラムは、前記予測目標値が下限の場合、前記予測値Ｌが前記下限よりも大きい確率を求め、
前記予測目標値が上限の場合、前記予測値Ｌが前記上限よりも小さい確率を求め、
前記予測目標値が下限と上限とからなる範囲の場合、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める。
【００３８】
本発明に係わる予測プログラムは、前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布を対象とする。
【００３９】
上記予測プログラムにおいて、前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを要旨とする。
【００４０】
本発明に係わる予測プログラムを記録した記録媒体は、上記予測プログラムをそれぞれ記録したコンピュータ読み取り可能な記録媒体を提供するものである。
【００４１】
本発明に係わる予測プログラムを記録した記録媒体は、前記予測目標値が下限の場合、、前記予測値Ｌが前記下限よりも大きい確率を求め、
前記予測目標値が上限の場合、前記予測値Ｌが前記上限よりも小さい確率を求め、
前記予測目標値が下限と上限とからなる範囲の場合、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める。
【００４２】
本発明に係わる予測プログラムを記録した記録媒体は、前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布を対象とする。
【００４３】
上記予測プログラムを記録した記録媒体であって、前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを要旨とする。
【００４４】
請求項１０乃至１２のいずれかに記載した発明によれば、注目する現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測プログラムを、フレキシブルディスク、ＣＤ−ＲＯＭ、ＤＶＤ（Digital Versatile Disk）、光磁気ディスク、ＰＣカード等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることができる。
【００４５】
【発明の実施の形態】
次に、図面を参照して本発明の実施の形態を説明する。
【００４６】
図１は、本発明の一実施形態に係る予測装置１の構成を表すブロック図である。同図に示す予測装置１は、注目する現象で予測対象となる事象が発生した期日やその発生件数等の過去の履歴、および前記予測対象となる事象において予測すべき事象の上限や下限等の値の範囲や予測時点等に代表される予測条件（予測を行うときの条件）を入力する入力手段として、キーボードやマウス等からなる入力部１１を有する。入力部１１で入力される過去の履歴および事象の予測条件は、データベースとしての事象内容記憶部１２および予測条件記憶部１４にそれぞれ記憶される。これらの事象内容記憶部１２および予測条件記憶部１４は、それぞれ事象内容記憶手段および予測条件記憶手段の少なくとも一部を構成するものである。
【００４７】
ここで、「注目する現象」とは、例えばコンピュータウィルスの感染、ソフトウェアの信頼度、人口や生物個体数の変化、並びに新製品、技術、およびサービスの普及等である。各現象で「予測対象とする事象」を具体的に表す値は、コンピュータウィルスの感染を予測する場合にはコンピュータウィルス感染数、ソフトウェア信頼度を予測する場合には故障や欠陥（バグ）の数、人口や生物個体数を予測する場合には人口や生物個体数、新製品、技術、サービスの普及を予測する場合には普及数のことである。
【００４８】
予測装置１には、この他に、事象発生期日や事象発生件数を集計して累積データ（累積事象数）を作成する累積件数集計部１３、パラメータの推定を回帰分析によって行うパラメータ推定手段としての回帰分析部１５、パラメータの推定値を用いて予測すべき事象の統計的分布等を推定する分布推定手段としての分布推定部１６、および推定結果の出力を行うための出力手段としてディスプレイ装置等からなる出力部１７を有している。
【００４９】
予測装置１は、演算機能および制御機能を備えた中央処理装置（ＣＰＵ：Central Processing Unit）、プログラムやデータを格納する機能を備えた主記憶装置、およびハードディスクドライブ、フレキシブルディスクドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤドライブ、光磁気ディスクドライブ、ＰＣカードドライブ等の補助記憶装置を備えたコンピュータによって構成されている。
【００５０】
また、本実施形態に係る各種処理を実行するための予測プログラムは前述した主記憶装置またはハードディスクに格納されている。この予測プログラムは、フレキシブルディスク、ＣＤ−ＲＯＭ、ＤＶＤ、光磁気ディスク、ＰＣカード等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることも勿論可能である。
【００５１】
以上の構成を有する予測装置１が行う予測方法について説明する。まず、予測を行うためのアルゴリズムを説明する。本実施形態においては、ロジスティック曲線モデルを例として説明するが、本発明の予測方法がこのモデルに限定されるわけでないことは勿論である。
【００５２】
ロジスティック曲線を表す微分方程式であるロジスティック方程式
【数１４】
【数１５】
これらのロジスティック差分方程式に対して確率則を導入し、パラメータのうちの一つを確率変数として扱う確率ロジスティック差分方程式について説明する。
【００５４】
以下では、差分化した時間ステップ（差分ステップ）ごとに独立であり、かつ全ての時間ステップが同じ分布に従うような独立同一分布の確率変数列を｛Ａj ｝（ｊ：正の整数）とする。
【００５５】
式(23)に基づいた場合の確率ロジスティック差分方程式は、
【数１６】
【数１７】
【数１８】
【数１９】
【数２０】
【数２１】
【００５６】
通常用いる差分方程式は、元の微分方程式を差分間隔の何乗のオーダで近似するかに主眼が置かれる。したがって通常の場合には、差分化がなされた時点で、解の形状や、時間無限大の極限における累積事象数の一定値への収束性などの性質は一般に保存されない。
【００５７】
これに対して本実施形態においては、元の微分方程式が有する性質を保つ差分方程式を用いるので、通常用いられる手法より精度の高いパラメータ推定が可能となる。
【００５８】
確率変数列｛Ａj ｝は、予測対象としている事象に応じて、その事象が従う分布を決定する。一例として、式(27)において、
Ｘj ＝１−δＡj ・・・(32)
とおき、この独立同一確率変数列｛Ｘj ｝がべき関数分布に従う場合について説明する。なお、式(29), (31)を用いる場合には、Ｘj をそれぞれ
【数２２】
【数２３】
【００５９】
注目している事象について、所定の等間隔時間ステップごとに累積件数集計部１３で計測（累積事象数の算出）が行われ、現在時間ステップｉ番目（ｉ：正の整数）までの計測結果（実績値）が得られているものとする。この時点でｎ（＞ｉ）番目の時間ステップの累積事象数として予測される予測値を推定し、その予測値の統計的な分布を求める方法について説明する。
【００６０】
今、ある値Ｌ― を
【数２４】
【数２５】
同様に、ある値Ｌ― を
【数２６】
【数２７】
式(38)および式(41)を導出する計算では、Ｘj をべき関数分布としているので、−log Ｘj は指数分布の独立同一確率変数列であり、
【数２８】
【００６３】
なお、以上の計算において、Ｌ―（Ｌ― )からｘ―（ｘ― )を求めるためにはｋ、ｍが必要である。そこで、これらのパラメータの推定値ｋ＾、ｍ＾を求めるために、前述した従来技術２の手法を用いる。
【００６４】
また、Ｐ｛Ｌn ＜Ｌ― ｝の計算においては、独立同一確率変数列｛Ｘj ｝がべき関数分布に従うものとしているので、−log Ｘj は指数分布の独立同一確率変数列となり、最尤法を用いることによって、べき数γが、
【数２９】
【００６５】
本実施形態においては、以上説明したように予測値の下限または上限を設定して分布を求める以外にも、予測値の範囲を指定して、その範囲内での分布を求めることも可能である。例えば、任意の時間ステップ数ｎでの予測値がＬ― よりも大きくＬ― よりも小さい確率Ｐ｛Ｌ―＜Ｌn ＜Ｌ― ｝を求めるときには、
【数３０】
【００６６】
図２は、予測装置１が以上説明したアルゴリズムにしたがって予測値の分布の推定を行うときの動作手順を説明するフローチャート図である。同図においては、予測装置１の処理の流れに加えて、事象内容記憶部１２ならびに予測条件記憶部１４に対するデータの流れを点線で示している。
【００６７】
まず、入力部１１で対象とする事象の発生期日、発生件数を一定の期間ごとに入力し、事象内容記憶部１２に記憶する（ステップＳ１）。
【００６８】
充分な事象数のデータが記憶された後、予測したい時点（予測値の分布が欲しい時点）と予測値の下限（Ｌ― ）、上限（Ｌ― ）、あるいは範囲（Ｌ―＜Ｌn ＜Ｌ― ）等の予測条件を入力部１１から入力し、予測条件記憶部１４に記憶する（ステップＳ３）。
【００６９】
一方、回帰分析部１５では事象内容記憶部１２で記憶したデータに基づいて累積件数集計部１３で累積された事象数を用いることにより、パラメータｋ、ｍを回帰分析によって推定し、パラメータ推定値ｋ＾、ｍ＾の算出を行う（ステップＳ５）。
【００７０】
なお、ステップＳ３およびＳ５の処理は独立なので、処理の順序を逆にしてもよいし、同時に行ってもよい。
【００７１】
以上の結果に基づいて、分析推定部１６が確率変数列｛Ａj ｝の平均値、予測対象となる事象の予測値の推定、および上記アルゴリズムにしたがって予測値の統計的分布を推定し（ステップＳ７）、この推定結果を出力部１７で表示する（ステップＳ９）。
【００７２】
以上説明した本発明の一実施形態によれば、ロジスティック方程式を差分化し、さらにその差分方程式のパラメータの一つを独立同一分布に従う確率変数として扱う確率ロジスティック差分方程式を用いることによって、注目する現象を記述するモデルのパラメータの推定に加えて、そのパラメータ推定値に基づいて予測対象となる事象の推定値の統計的分布を推定することが可能になる。
【００７３】
本実施形態は、例えばコンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデル等に適用することができ、各モデルに対応する予測値として、コンピュータウィルス感染数、故障・欠陥数（バグ数）、人口、生物個体数、並びに普及数をそれぞれ推定する場合に適用することが可能なものである。
【００７４】
なお、本発明においては、確率変数列が従う独立同一確率分布として、上述したべき関数分布以外に、例えば範囲（０，１）の一様分布を用いても同様の効果を得ることができる。
【００７５】
この場合、式(27)の厳密解において、
δＡj ＝Ｕj ・・・(45)
とおき、独立同一確率変数列｛Ｕj ｝が範囲（０，１）の一様分布に従うものとする。式(29)および式(31)の厳密解を用いるときには、それぞれ
【数３１】
【００７６】
予測値Ｌn が式(36)で定義されるＬ― よりも大きい確率Ｐ｛Ｌn ＞Ｌ― ｝は、次のように計算される。
【数３２】
【００７７】
また、予測値Ｌn が式(39)で定義されるＬ― よりも小さい確率Ｐ｛Ｌn ＜Ｌ― ｝は、
【数３３】
【００７８】
式(48)および式(49)を導出する計算においては、−log（１−Ｕj ）は指数分布の独立同一確率変数列であり、
【数３４】
【００７９】
なお、時間ステップ数ｎでの予測値がＬ― よりも大きくＬ― よりも小さい確率Ｐ｛Ｌ―＜Ｌn ＜Ｌ― ｝についても、式(43)または式(44)を用いて計算することができるのは勿論である。
【００８０】
このように、本発明は、上記一実施形態同様の効果を奏するさまざまな実施の形態等を含みうるものであることはいうまでもない。
【００８１】
【発明の効果】
以上説明した本発明によれば、ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定するとともに、この推定したパラメータに基づいて、前記現象の中で発生する予測対象事象の推定値およびその推定値の統計的分布の推定を可能にする予測方法、予測装置、予測プログラム、および予測プログラムを記録した記録媒体を提供することができる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る予測装置の構成を示すブロック図である。
【図２】本発明の一実施形態に係る予測方法の動作手順を示すフローチャート図である。
【図３】従来法における予測曲線および実績累積件数を表す説明図である。
【符号の説明】
１ 予測装置
１１ 入力部
１２ 事象内容記憶部
１３ 累積件数集計部
１４ 予測条件記憶部
１５ 回帰分析部
１６ 分布推定部
１７ 出力部
Claims (6)
- ロジスティック曲線モデルを用いて注目する現象を記述するモデルのパラメータを推定し、この推定したパラメータに基づいて前記現象の中で発生する事象の予測を行う予測装置であって、
前記現象の所定の期間ごとの発生件数を記憶する事象内容記憶手段と、
予測する時点と予測目標値からなる予測条件を記憶する予測条件記憶手段と、
前記事象内容記憶手段から前記発生件数を読み出し、注目している事象について所定の時間間隔で累積件数を集計し、実績値を得る累積件数集計手段と、
前記ロジスティック曲線を表すロジスティック方程式を時間ステップで差分化したロジスティック差分方程式のうち厳密解を有するロジスティック差分方程式に前記累積件数集計手段で得られた実績値を与え、回帰分析によってロジスティック方程式の厳密解Ｌ（ｔ）＝ｋ／｛１＋ｍ＊ｅｘｐ（−αｔ）｝のパラメータｋ、ｍを求めるパラメータ推定手段と、
前記ロジスティック差分方程式の厳密解のパラメータｋ、ｍを前記求められた値とし、前記累積件数集計手段で得られた実績値を適用して、パラメータαが時間ステップに対応した確率変数列Ａを求める確率変数列算出手段と、
前記確率変数列Ａが差分ステップごとに独立であるとともに全ての差分ステップで同一な分布である独立同一分布に従うこととして、前記確率変数列Ａに基づいて最尤法を用いることによって前記独立同一分布のパラメータを求める分布係数特定手段と、
パラメータ推定手段で求められた前記パラメータｋ、ｍ、分布係数特定手段で求められた確率変数列Ａの独立同一分布パラメータ、およびロジスティック差分方程式の厳密解を用いて、前記予測する時点の予測値Ｌが前記予測目標値を満たす確率を求める分布推定手段と、
を有することを特徴とする予測装置。 - 前記予測目標値が下限の場合、前記分布推定手段は、前記予測値Ｌが前記下限よりも大きい確率を求め、
前記予測目標値が上限の場合、前記分布推定手段は、前記予測値Ｌが前記上限よりも小さい確率を求め、
前記予測目標値が下限と上限とからなる範囲の場合、前記分布推定手段は、前記予測値Ｌが下限よりも大きく、かつ、上限よりも小さい確率を求める、
ことを特徴とする請求項１記載の予測装置。 - 前記独立同一分布は、確率変数の取り得る範囲が０から１までであることを満たす分布である
ことを特徴とする請求項１または請求項２の予測装置。 - 前記現象を記述するモデルは、コンピュータウィルスの感染数を予測するコンピュータウィルス感染数予測モデル、ソフトウェアの信頼度を予測するソフトウェア信頼度成長モデル、人口の変化を予測する人口予測モデル、生物個体数の変化を予測する生物個体数予測モデル、並びに新製品、技術、およびサービスの普及を予測する普及予測モデルのうちの少なくともいずれかであることを特徴とする請求項１ないし請求項３のいずれか１項に記載の予測装置。
- 請求項１ないし４のいずれか１項に記載の予測装置を構成する各手段としてコンピュータを機能させる予測プログラム。
- 請求項５に記載した予測プログラムを記録した、コンピュータ読み取り可能な記録媒体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002241139A JP4498666B2 (ja) | 2002-08-21 | 2002-08-21 | 予測装置、予測プログラム、および記録媒体 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002241139A JP4498666B2 (ja) | 2002-08-21 | 2002-08-21 | 予測装置、予測プログラム、および記録媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004078780A JP2004078780A (ja) | 2004-03-11 |
| JP4498666B2 true JP4498666B2 (ja) | 2010-07-07 |
Family
ID=32023726
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002241139A Expired - Fee Related JP4498666B2 (ja) | 2002-08-21 | 2002-08-21 | 予測装置、予測プログラム、および記録媒体 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4498666B2 (ja) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4770763B2 (ja) * | 2007-03-19 | 2011-09-14 | 日本電信電話株式会社 | 予測モデル選択装置と方法および予測装置と推定値予測方法ならびにプログラム |
| JP4952537B2 (ja) * | 2007-11-27 | 2012-06-13 | 三菱電機株式会社 | ネットワーク性能予測システム、ネットワーク性能予測方法およびプログラム |
| JP5198981B2 (ja) * | 2008-09-03 | 2013-05-15 | インターナショナル・ビジネス・マシーンズ・コーポレーション | 作業評価値予測方法、プログラム及びシステム |
| RU2446454C1 (ru) * | 2011-01-12 | 2012-03-27 | Государственное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (КГЭУ) | Цифровое прогнозирующее устройство |
| RU2551793C1 (ru) * | 2014-07-29 | 2015-05-27 | Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования "ВОЕННАЯ АКАДЕМИЯ СВЯЗИ имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации | Устройство для прогнозирования случайных событий |
| RU2643645C2 (ru) * | 2016-02-09 | 2018-02-02 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Казанский государственный энергетический университет" (ФГБОУ ВПО "КГЭУ") | Цифровое прогнозирующее устройство |
| JP2017157109A (ja) | 2016-03-03 | 2017-09-07 | 富士通株式会社 | 推定装置、推定方法および推定プログラム |
| KR101857271B1 (ko) * | 2016-12-27 | 2018-05-11 | 주식회사 아이디어랩스 | 선택자가 특정 시점에 특정 품목을 선택할지 여부를 예측하기 위한 방법 및 이를 이용한 컴퓨팅 장치 |
| JP6761788B2 (ja) * | 2017-08-22 | 2020-09-30 | 日本電信電話株式会社 | 予測装置、予測方法及び予測プログラム |
| CN108667833B (zh) * | 2018-04-28 | 2020-09-08 | 国网山东省电力公司莱芜供电公司 | 基于耦合的通信系统恶意软件传播建模与最优控制方法 |
-
2002
- 2002-08-21 JP JP2002241139A patent/JP4498666B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2004078780A (ja) | 2004-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6636071B2 (ja) | コンピュータ実装方法、コンピュータシステム及びコンピュータ装置 | |
| US10846643B2 (en) | Method and system for predicting task completion of a time period based on task completion rates and data trend of prior time periods in view of attributes of tasks using machine learning models | |
| JP5982489B2 (ja) | 動的外れ値偏り低減システム及び方法 | |
| JP6571914B2 (ja) | 情報の複数のドメインを組合せることによる仕事の実施データ内の異常の検知 | |
| US20090037402A1 (en) | System and method for predicting clickthrough rates and relevance | |
| US20080178145A1 (en) | Method and System for Generating a Predictive Analysis of the Performance of Peer Reviews | |
| JP2006024017A (ja) | コンピュータ資源のキャパシティを予測するためのシステム、方法およびプログラム | |
| JP6845126B2 (ja) | 故障確率算出装置、故障確率算出方法及びプログラム | |
| JP7188950B2 (ja) | データ処理方法およびデータ処理プログラム | |
| JP4498666B2 (ja) | 予測装置、予測プログラム、および記録媒体 | |
| JP2023052477A (ja) | データ処理方法およびデータ処理プログラム | |
| CN116151353B (zh) | 一种序列推荐模型的训练方法和对象推荐方法 | |
| JP2007323315A (ja) | 協調フィルタリング方法、協調フィルタリング装置、および協調フィルタリングプログラムならびにそのプログラムを記録した記録媒体 | |
| CN111625720B (zh) | 数据决策项目执行策略的确定方法、装置、设备和介质 | |
| JP7231829B2 (ja) | 機械学習プログラム、機械学習方法および機械学習装置 | |
| CN111382909A (zh) | 基于生存分析模型扩展坏样本的拒绝推断方法及相关设备 | |
| JP2021163412A (ja) | 要員計画支援装置および要員計画支援方法 | |
| JP2005063208A (ja) | ソフトウェア信頼度成長モデル選択方法、ソフトウェア信頼度成長モデル選択装置、ソフトウェア信頼度成長モデル選択プログラム、およびプログラム記録媒体 | |
| Singham et al. | Density estimation of simulation output using exponential epi-splines | |
| CN113673811B (zh) | 一种基于session的在线学习绩效评估方法及装置 | |
| US20040039559A1 (en) | Method and apparatus for determining output uncertainty of computer system models | |
| JP2008305229A (ja) | 需要予測方法および装置 | |
| JP3612428B2 (ja) | ソフトウェア信頼性予測方法および装置と記録媒体 | |
| JP4102019B2 (ja) | 製造プロセスにおける解析装置、方法、及びコンピュータ読み取り可能な記憶媒体 | |
| JP2006072670A (ja) | ソフトウェアの信頼性予測方法、そのプログラムを記録した記録媒体および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040729 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061117 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061128 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070129 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20070710 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070906 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20070913 |
|
| A912 | Re-examination (zenchi) completed and case transferred to appeal board |
Free format text: JAPANESE INTERMEDIATE CODE: A912 Effective date: 20071214 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20090709 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20090909 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100414 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |