JPH05314094A - Sales prediction supporting device - Google Patents

Sales prediction supporting device

Info

Publication number
JPH05314094A
JPH05314094A JP11768292A JP11768292A JPH05314094A JP H05314094 A JPH05314094 A JP H05314094A JP 11768292 A JP11768292 A JP 11768292A JP 11768292 A JP11768292 A JP 11768292A JP H05314094 A JPH05314094 A JP H05314094A
Authority
JP
Japan
Prior art keywords
component
value
fluctuation
fluctuation component
trend
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.)
Pending
Application number
JP11768292A
Other languages
Japanese (ja)
Inventor
Mikiaki Sasajima
己喜朗 笹島
Hiroto Yamaguchi
裕人 山口
Masato Baba
正人 馬場
Nobuo Hirai
伸郎 平井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kao Corp
Original Assignee
Kao Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Kao Corp filed Critical Kao Corp
Priority to JP11768292A priority Critical patent/JPH05314094A/en
Publication of JPH05314094A publication Critical patent/JPH05314094A/en
Pending legal-status Critical Current

Links

Abstract

PURPOSE:To provide the supporting device which facilitates the operations, reduces the generation of error as well and is convenient for users to finally judge prediction by simplifying a transient probability distribution function or a membership function into the form of practical use. CONSTITUTION:Sales result data in the past inputted from an input device are stored in a result storage file part 1. Based on these data, a trend component arithmetical prediction part 2 performs the extraction and prediction of trend components. The trend component is the tendency of sale and is smoothed by removing the increase of sales results under any temporary and special conditions, and a predictive value is calculated. Based on this calculated total predictive value, the transient probability distribution function is triangularly expressed for each item of products and displayed at an output device 10. The bottom side of the triangle is defined as the size of fluctuation width, the apex is defined as the total predictive value for each item and when the height is fixed for all the items, the much smoother triangle is judged as a product for which the sales fluctuation in the past is big and the more sale can be estimated corresponding to following sales strategy.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、工場の生産計画や商品
の販売計画を立案するために、商品の販売予測を行うた
めの支援装置として利用する。本発明は、商品の在庫数
量を最小限にコントロールするとともに、商品の注文に
対して欠品や納品遅れをなくし、さらに営業活動による
販売促進を合理的に行うために、過去の販売実績から次
期の販売予測を行うための支援装置として利用する。本
発明の装置はコンピュータ装置を利用するものであっ
て、演算および表示はすべてプログラム制御により実行
される。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention is used as a support device for predicting the sales of goods in order to make a production plan for a factory or a sales plan for goods. The present invention controls the inventory quantity of products to the minimum, eliminates missing items and delays in delivery of product orders, and further promotes sales activities rationally. Used as a support device for making sales forecasts. The device of the present invention utilizes a computer device, and all calculations and displays are executed under program control.

【0002】[0002]

【従来の技術】商品の在庫数量を最小限にして在庫コス
トを小さくすることは商品のコストを小さく抑制するた
めに重要であり、一方、商品の注文に対して在庫がない
状態すなわち欠品があると、その商品の信用を失うとと
もに契約によっては違約金を支払わなければならない。
また、営業活動により商品の売れ行きを促進する場合に
次期の販売目標をどこに設定するかは重要な要素であ
る。
2. Description of the Related Art It is important to minimize the inventory quantity of products and reduce the inventory cost in order to keep the cost of the products small. If you do, you will lose credit for the product and you will have to pay a penalty depending on the contract.
Also, when promoting sales of products through sales activities, where to set the next sales target is an important factor.

【0003】このために従来から、過去の当該商品につ
いて毎月集計される販売実績を長くコンピュータに記録
しておき、それをコンピュータの表示画面にグラフで表
示し、この表示を見ながら担当の管理者が次の月の販売
数量を予測することが一般に行われている。この場合
に、例えば商品が洗剤であるとすると、洗剤は夏期に需
要が大きく冬期には需要が少ない季節変動を有する性質
があるから、次の月の販売数量を予測するには、この季
節変動をトレンド成分として差し引いて考えることが行
われる。このために、コンピュータによりトレンド成分
を抽出して、表示画面に表示されるグラフにはこの抽出
されたトレンド成分を差し引いて表示する技術も知られ
ている。
For this reason, conventionally, the sales record of the past relevant merchandise, which is totaled every month, is recorded in a computer for a long time, which is displayed in a graph on the display screen of the computer, and the administrator in charge while watching this display. It is common practice to forecast sales volume for the next month. In this case, for example, if the product is a detergent, the detergent has a seasonal variation with a large demand in the summer and a small demand in the winter. Is considered as a trend component. Therefore, a technique is known in which a trend component is extracted by a computer, and the extracted trend component is displayed in a graph displayed on a display screen.

【0004】一方、OR(オペレーション・リサーチ)
の研究分野では、販売予測あるいは生産管理のためにこ
れまでさまざまなアルゴリズムが知られている。このア
ルゴリズムによりコンピュータ装置を利用して多数のデ
ータとともに演算を行った報告はいくつか知られてい
る。
On the other hand, OR (Operation Research)
In the research field of, various algorithms have been known so far for sales forecasting or production management. There are some known reports in which a computer device is used to perform calculations with a large number of data using this algorithm.

【0005】[0005]

【発明が解決しようとする課題】販売予測を高級なアル
ゴリズムあるいは高級な装置を用いて実行しても、実際
の販売は自然法則に左右されるだけでなく人間的な要素
により変動する部分があるから、演算の結果をそのまま
うのみにして営業政策に利用することはできない。最終
的には企業の営業あるいは経営の判断として人間的に判
断されることになる。高級なアルゴリズムや高級な装置
はできるだけ正確な予測を行うための研究成果ではある
が、日常的に利用される販売予測のためには、投入する
データの種類や量が多くなり、操作が複雑になり、利用
する上での誤りも発生しやすいなど、必ずしも使用者が
利用し判断するための便利な支援装置とはならない。
Even if the sales forecast is executed by using a high-grade algorithm or a high-grade device, the actual sales are not only influenced by the laws of nature but also fluctuated by human factors. Therefore, the calculation result cannot be used as it is for business policy. Eventually, it will be a human decision as the decision of the business or management of the company. High-end algorithms and high-end devices are the research results for making accurate forecasts as much as possible, but for sales forecasts that are used on a daily basis, the types and amounts of data to be input are large and the operations are complicated. Therefore, it is not always a convenient support device for the user to use and make a judgment such as an error in using the device.

【0006】本発明は、販売予測あるいは生産管理のた
めに、変動の遷移する過程に基づき演算を行う遷移確率
分布関数、あるいはメンバーシップ関数が有用な論理で
あることに着目し、これを実用的な形に単純化して、投
入するデータの種類や量が小さく、操作が簡単であり、
誤りの発生も小さく、最終的に人が予測判断をするため
に便利な支援装置を提供することを目的とする。
The present invention focuses on the fact that a transition probability distribution function or a membership function, which performs an operation based on a transition process of fluctuations, is a useful logic for sales forecasting or production management, and this is put to practical use. Simple, simple to use, the type and amount of data to be input is small, and the operation is easy.
It is an object of the present invention to provide a support device which is small in the occurrence of errors and is convenient for people to finally make a prediction judgment.

【0007】[0007]

【課題を解決するための手段】本発明は、1つまたは複
数の商品について販売実績が蓄積された実績記憶ファイ
ルと、この実績記憶ファイルに蓄積された前記販売実績
を取り出しこの販売実績から周期的な変動であるトレン
ド成分を抽出するトレンド成分演算手段と、前記販売実
績から前記トレンド成分を除去して変動成分を抽出する
変動成分演算手段と、次に到来する期間についてこの変
動成分の予測を行う変動成分予測手段と、前記次に到来
する期間について前記トレンド成分を外挿することによ
りトレンド成分の予測を行うトレンド成分予測手段と、
このトレンド成分予測手段により予測されたトレンド成
分予測値と前記変動成分予測手段により予測された変動
成分予測値との和を前記次に到来する期間についての総
合予測値とする予測演算手段とを備えた販売予測支援装
置において、前記変動成分予測手段は、前記実績記憶フ
ァイルから取り出した実績値を複数個のクラスタに区分
し、最新実績値の変動成分H(k)がi番目のクラスタ
にありi番目のクラスタから変動した実績値をΔaj
(j=1、2、…、N)とするとき、その平均値を
According to the present invention, a sales record file in which sales records of one or a plurality of products are accumulated, and the sales record accumulated in the recording record file are taken out, and the sales records are periodically reproduced. Component calculation means for extracting a trend component that is a large variation, a fluctuation component calculation means for removing the trend component from the sales record to extract a variation component, and a prediction of this variation component for the next coming period. Fluctuating component predicting means, and trend component predicting means for predicting the trend component by extrapolating the trend component for the next coming period,
And a prediction calculation unit that sets the sum of the trend component prediction value predicted by the trend component prediction unit and the fluctuation component prediction value predicted by the fluctuation component prediction unit as a comprehensive prediction value for the next coming period. In the sales forecast support apparatus, the fluctuation component predicting unit divides the performance value extracted from the performance storage file into a plurality of clusters, and the fluctuation component H (k) of the latest performance value is in the i-th cluster i The actual value fluctuated from the second cluster is Δaj
(J = 1, 2, ..., N), the average value is

【0008】[0008]

【数4】 とし、変動成分予測値を H' (k+1) = H(k) + Δi として演算する手段を備えたことを特徴とする。[Equation 4] And a means for calculating the fluctuation component predicted value as H ′ (k + 1) = H (k) + Δi.

【0009】また、前記実績値、前記トレンド成分、お
よび前記変動成分を時間軸上に表示するとともに、前記
総合予測値を頂点位置とし前記変動成分演算手段により
抽出された変動成分の振幅を底辺長さとし次に到来する
期間の長さを高さとする三角形で表示する手段を備える
ことが望ましい。さらに、前記最新実績値の変動成分H
(k)が変動成分が零であるクラスタ内にあるときに
は、変動成分予測値を零として、前記トレンド成分予測
値を前記総合予測値とすることが望ましい。
Further, the actual value, the trend component, and the fluctuation component are displayed on a time axis, and the amplitude of the fluctuation component extracted by the fluctuation component computing means is set as a vertex position and the base length is set as the base length. It is desirable to provide means for displaying a triangle whose height is the length of the next period. Furthermore, the fluctuation component H of the latest actual value
When (k) is in a cluster in which the fluctuation component is zero, it is desirable that the fluctuation component prediction value be zero and the trend component prediction value be the comprehensive prediction value.

【0010】本発明のもう一つの観点として、1つまた
は複数の商品について販売実績が蓄積された実績記憶フ
ァイルと、この実績記憶ファイルに蓄積された前記販売
実績から季節変動に伴うトレンド成分を抽出するトレン
ド成分演算手段と、前記販売実績から前記トレンド成分
を除去して変動成分を抽出する変動成分演算手段と、次
に到来する期間についてこの変動成分の予測を行う変動
成分予測手段と、前記次に到来する期間について前記ト
レンド成分を外挿することによりトレンド成分の予測を
行うトレンド成分予測手段と、このトレンド成分予測手
段により予測されたトレンド成分予測値と前記変動成分
予測手段により予測された変動成分予測値との和を前記
次に到来する期間についての総合予測値とする予測演算
手段とを備えた販売予測支援装置において、前記変動成
分予測手段は、前記実績記憶ファイルから取り出した実
績値を複数個のクラスタに区分し、最新実績値の変動成
分H(k)がi番目のクラスタにあり、その変動成分H
(k)がその1期前の変動成分実績値H(k−1)に対
し増加傾向にあるとき、i番目のクラスタから増加変動
した実績値をΔcj、(j=1、2、…、N)として、
その増加変動した実績値の平均値を
As another aspect of the present invention, a history storage file in which sales records of one or a plurality of products are stored, and trend components associated with seasonal fluctuations are extracted from the sales records stored in the sales record file. Trend component calculating means, a fluctuation component calculating means for removing the trend component from the sales record to extract a fluctuation component, a fluctuation component predicting means for predicting the fluctuation component for the next coming period, and Trend component predicting means for predicting the trend component by extrapolating the trend component for the period of arrival, and the trend component prediction value predicted by the trend component predicting means and the fluctuation predicted by the fluctuation component predicting means. Sales including a prediction calculation means for making the sum of the component predicted values into a total predicted value for the next coming period In the prediction support apparatus, the fluctuation component prediction means divides the performance value extracted from the performance storage file into a plurality of clusters, and the fluctuation component H (k) of the latest performance value is in the i-th cluster, and the fluctuation thereof Ingredient H
When (k) has an increasing tendency with respect to the fluctuation component actual value H (k−1) of the previous period, the actual value increased and varied from the i-th cluster is Δcj, (j = 1, 2, ..., N). ) As
The average value of the actual values that increased and changed

【0011】[0011]

【数5】 とし、その変動成分H(k)がその1期前の変動成分実
績値H(k−1)に対し減少傾向にあるとき、i番目の
クラスタから減少変動した実績値をejとして、その減
少変動した実績値の平均値を
[Equation 5] When the fluctuation component H (k) has a decreasing tendency with respect to the fluctuation component actual value H (k−1) of the previous period, the actual decreasing value of the i-th cluster is set as ej, and the decreasing fluctuation The average value of the actual results

【0012】[0012]

【数6】 とし、変動成分予測値を H' (k+1) = H(k) + Δi として演算する手段を備えたことを特徴とする。[Equation 6] And a means for calculating the fluctuation component predicted value as H ′ (k + 1) = H (k) + Δi.

【0013】また、前記実績値、前記トレンド成分、お
よび前記変動成分を時間軸上に表示するとともに、前記
総合予測値を頂点位置とし前記変動成分演算手段により
抽出された変動成分の振幅を底辺長さとし次に到来する
期間の長さを高さとする三角形で表示する手段を備える
ことが望ましい。前記最新実績値の変動成分H(k)が
変動成分が零であるクラスタ内にあるときには、 H' (k+1) = 0 として、前記トレンド成分予測値を前記総合予測値とす
ることが望ましい。さらに、複数の商品品目について前
記総合予測値、前記変動成分予測値、および前記三角形
を演算する手段を備えることが望ましい。
The actual value, the trend component, and the fluctuation component are displayed on the time axis, and the amplitude of the fluctuation component extracted by the fluctuation component computing means is the base length with the total predicted value as the vertex position. It is desirable to provide means for displaying a triangle whose height is the length of the next period. When the fluctuation component H (k) of the latest actual value is in the cluster in which the fluctuation component is zero, it is preferable to set H ′ (k + 1) = 0 and use the trend component prediction value as the comprehensive prediction value. Further, it is desirable to include means for calculating the total predicted value, the fluctuation component predicted value, and the triangle for a plurality of product items.

【0014】[0014]

【作用】入力装置(キーボードなど)から入力された過
去の販売実績データが実績記憶ファイル部に蓄積されて
いる。このデータを元にしてトレンド成分演算予測部で
は、トレンド成分の抽出および予測が行われる。トレン
ド成分とは売行きの傾向のことで、一時的で特異な条件
下での販売実績の伸びを除去して平滑化され、予測値が
算出される。これは特異なデータが予測値の精度に影響
するのを防ぐためである。
[Operation] Past sales record data input from the input device (keyboard or the like) is accumulated in the record storage file section. Based on this data, the trend component calculation / prediction unit extracts and predicts the trend component. The trend component is the tendency of sales, which is smoothed by removing the growth of the sales performance under a temporary and unique condition to calculate the forecast value. This is to prevent peculiar data from affecting the accuracy of the predicted value.

【0015】次に、変動成分演算予測部において変動成
分の抽出および予測が行われる。変動成分とは売行きの
増減およびその変動幅のことである。変動幅は数段階の
クラスタに分解され、それぞれのクラスタごとに変動の
予測値が算出されるが増加傾向のときと、減少傾向のと
きでさらに分類され、増加傾向のときは過去の増加傾向
のデータを元に予測値を算出し、減少傾向のときは過去
の減少傾向のデータを元に予測値を算出する。
Next, the fluctuation component calculation / prediction unit extracts and predicts the fluctuation component. The fluctuation component is the fluctuation in sales and the fluctuation range. The fluctuation range is decomposed into clusters of several stages, and the predicted value of fluctuation is calculated for each cluster, but it is further classified according to the increasing trend and the decreasing trend, and the increasing trend of the past increasing trend. The predicted value is calculated based on the data, and when the tendency is decreasing, the predicted value is calculated based on the data of the past decreasing tendency.

【0016】このように算出されたトレンド成分および
変動成分の予測値を元にして、総合予測演算部は次月の
売上予測値を算出する。その結果は総合予測値データ出
力インタフェースから出力される。
Based on the predicted values of the trend component and the fluctuation component calculated in this way, the comprehensive prediction calculation unit calculates the sales predicted value for the next month. The result is output from the comprehensive predicted value data output interface.

【0017】また、この算出された総合予測値を元にし
て、製品の品目ごとの遷移確率分布関数を三角形で表わ
し、出力装置(CRTなど)に表示することもできる。
三角形の底辺を変動幅の大きさで表わし、頂点を品目ご
との総合予測値で表わし、高さは全品目とも一定とすれ
ば、よりなだらかな三角形ほど過去に売上の変動が大き
く、今後の販売戦略によってはまだまだ売れる見込みが
あると判断される。
Based on the calculated total predicted value, the transition probability distribution function for each product item can be represented by a triangle and displayed on an output device (CRT or the like).
If the bottom of the triangle is represented by the size of the fluctuation range, the apex is represented by the total forecast value for each item, and the height is constant for all items, the smoother the triangle, the greater the fluctuation in sales in the past and the future sales. Depending on the strategy, it is judged that there is still a possibility of selling.

【0018】[0018]

【実施例】本発明実施例装置の構成を図1を参照して説
明する。図1は本発明実施例装置のブロック図である。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS The configuration of the device of the embodiment of the present invention will be described with reference to FIG. FIG. 1 is a block diagram of an apparatus according to an embodiment of the present invention.

【0019】本発明は、1つまたは複数の商品について
それぞれ十分に長い期間にわたる販売実績が蓄積された
実績記憶ファイル部1と、この実績記憶ファイル部1に
蓄積された前記販売実績を取り出しこの販売実績から周
期的な変動や季節変動に伴うトレンド成分を抽出して前
記販売実績から前記トレンド成分を除去して変動成分を
抽出し、次に到来する期間についてこの変動成分の予測
を行う変動成分演算予測部3と、前記次に到来する期間
について前記トレンド成分を外挿することによりトレン
ド成分の予測を行うトレンド成分演算予測部2と、この
トレンド成分演算予測部2により予測されたトレンド成
分予測値と変動成分演算予測部3により予測された変動
成分予測値との和を前記次に到来する期間についての総
合予測値とする総合予測演算部4とを備えた販売予測支
援装置において、変動成分演算予測部3は、実績記憶フ
ァイル部1から取り出した長い期間にわたる実績値を複
数個のクラスタに区分し、最新実績値の変動成分H
(k)がi番目のクラスタにありi番目のクラスタから
増加傾向に変動した実績値をΔaj、(j=1、2、
…、N)とするとき、その変動成分H(k)がその1期
間前の変動成分実績値H(k−1)に対し増加傾向にあ
るとき、その平均値を
The present invention takes out the sales record file unit 1 in which the sales records of one or a plurality of products are accumulated for a sufficiently long period of time, and the sales records accumulated in the record storage unit 1 are taken out and sold. A trend component calculation that extracts a trend component associated with periodic fluctuations and seasonal variations from the actual results, removes the trend component from the sales results, extracts the fluctuation component, and predicts the fluctuation component for the next coming period. A prediction unit 3, a trend component calculation prediction unit 2 that predicts a trend component by extrapolating the trend component for the next coming period, and a trend component prediction value predicted by the trend component calculation prediction unit 2. And the variation component prediction value predicted by the variation component calculation prediction unit 3 as a total prediction value for the next coming period. In the sales prediction support device including the prediction calculation unit 4, the fluctuation component calculation prediction unit 3 divides the long-term result values extracted from the result storage file unit 1 into a plurality of clusters, and changes the latest result value fluctuation component. H
(K) is in the i-th cluster, and the actual value that has changed from the i-th cluster in an increasing tendency is Δaj, (j = 1, 2,
, N), when the fluctuation component H (k) tends to increase with respect to the fluctuation component actual value H (k−1) one period before, the average value is calculated.

【0020】[0020]

【数7】 とし、その変動成分H(k)がその1期間前の変動成分
実績値H(k−1)に対し減少傾向にあるとき、i番目
のクラスタから減少変動した実績値をej(j=1、
2、…、N)とするとき、その平均値を
[Equation 7] Then, when the fluctuation component H (k) has a decreasing tendency with respect to the fluctuation component actual value H (k−1) one period before, the actual value decreased from the i-th cluster is changed to ej (j = 1,
2, ..., N), the average value is

【0021】[0021]

【数8】 とし、変動成分予測値を H' (k+1) = H(k) + Δi として演算する手段を備えたことを特徴とする。[Equation 8] And a means for calculating the fluctuation component predicted value as H ′ (k + 1) = H (k) + Δi.

【0022】また、前記実績値、前記トレンド成分、お
よび前記変動成分を時間軸上に表示するとともに、前記
総合予測値を頂点位置とし変動成分演算予測部3により
抽出された変動成分の振幅を底辺長さとし次に到来する
期間の長さを高さとする三角形でCRT10に表示する
手段を備える構成である。さらに、前記最新実績値の変
動成分H(k)が変動成分が零であるクラスタ内にある
ときには、 H' (k+1) = 0 として、前記トレンド成分予測値を前記総合予測値とす
る。また、複数の商品品目について前記総合予測値、前
記変動成分予測値、および前記三角形を演算する手段を
備える構成である。
The actual value, the trend component, and the fluctuation component are displayed on the time axis, and the amplitude of the fluctuation component extracted by the fluctuation component calculation / prediction unit 3 is set at the apex position. The CRT 10 is provided with a means for displaying a triangle with the length as the length and the length of the next period as the height. Further, when the fluctuation component H (k) of the latest actual value is in the cluster in which the fluctuation component is zero, H ′ (k + 1) = 0 is set and the trend component prediction value is set as the comprehensive prediction value. Further, it is configured to include means for calculating the comprehensive predicted value, the fluctuation component predicted value, and the triangle for a plurality of product items.

【0023】次に、本発明実施例装置の動作を図2を参
照して説明する。図2は本発明実施例装置の動作を示す
フローチャートである。十分に長い期間にわたる販売実
績を入力して、まずトレンド予測を行う。トレンド成分
の計算アルゴリズムは、3つのステップからなるモデル
フィッティングとAIC(赤池の情報量基準)によるモ
デル選択からなる。モデルは以下の4つの式の中からA
ICの値が最小になるものが選ばれる、 モデル1:y=a+b・k+c・sin(πk/6)+
d・cos(πk/6) モデル2:y=a+b・sin(πk/6)+c・co
s(πk/6) モデル3:y=a+b・k モデル4:y=a (ただし a、b、c、d:パラメータ,k:通算月) モデル1は、季節変動(1年周期の正弦波)とリニア
(直線的)なトレンドの上昇または下降を考慮したモデ
ルである。モデル2は季節変動のみを考慮したモデルで
あり、モデル3はリニアなトレンドの上昇または下降の
みを考慮したモデルである。また、モデル4はほぼ一定
の販売量があるというモデルである。モデルは各月ごと
に機械的に選択されるので、データセット(品目)によ
っては月ごとに選ばれるモデルが異なる可能性がある。
しかし、特異なデータ以外の販売量の変化をこれら4つ
のモデルのいずれかで表現することが可能な商品につい
ては、選択されるモデルが変化することがトレンド成分
の予測に重大な影響を与えることはないと考えられる。
モデルのフィッティングは3つのステップにより行われ
る。第一のステップは3ヵ月移動メディアンが用いら
れ、第二のステップは通常の最小2乗法(OLS:Or
dinary Least Squares)によって
モデルフィッティングが行われる。そして、第三のステ
ップは反復重み付き最小2乗法(IWLS:Itera
ted Weighted LeastSquare
s)によって最終的なパラメータの計算がなされる。
Next, the operation of the apparatus of the present invention will be described with reference to FIG. FIG. 2 is a flow chart showing the operation of the apparatus of the present invention. Enter the sales results over a sufficiently long period and first make a trend forecast. The trend component calculation algorithm consists of model fitting consisting of three steps and model selection by AIC (Akaike's information criterion). The model is A from the following four formulas
The one with the smallest IC value is selected, model 1: y = a + bk + csin (πk / 6) +
d · cos (πk / 6) model 2: y = a + b · sin (πk / 6) + c · co
s (πk / 6) Model 3: y = a + b · k Model 4: y = a (however, a, b, c, d: parameter, k: total month) Model 1 is a seasonal fluctuation (a sine wave with a one-year cycle). ) And a linear (linear) trend up or down trend is considered. The model 2 is a model that considers only seasonal fluctuations, and the model 3 is a model that considers only rising or falling of a linear trend. Further, the model 4 is a model in which there is an almost constant sales volume. Since the model is mechanically selected for each month, the model selected for each month may differ depending on the data set (item).
However, for products for which changes in sales volume other than peculiar data can be expressed by any of these four models, changes in the selected model have a significant impact on the forecast of trend components. Not considered.
Model fitting is done in three steps. The first step uses a 3-month mobile median, and the second step is the usual least squares method (OLS: Or).
The model fitting is carried out by means of the directory least squares. The third step is the iterative weighted least squares method (IWLS: Itera).
ted Weighted LeastSquare
The final parameter calculation is performed according to s).

【0024】次に、第一のステップについて説明する。
第一のステップの主な目的は、ピーク除去にある。「3
ヵ月移動メディアン」とは各月とその前後1ヵ月を合わ
せた3ヵ月のデータを用い、それらのメディアン(中央
値)を値とする方法である。3ヵ月移動メディアンを図
3および図4を参照して説明する。図3および図4は3
ヵ月移動メディアンを説明する図である。図3におい
て、1989年7月の場合は実績が 89年8月<89年7月<89年6月 の順番なので、移動メディアン値は7月のデータ値その
ものになる。また、1989年12月の場合は実績が 90年1月<89年11月<89年12月 の順番なので、移動メディアン値は11月のデータ値に
なる。3ヵ月移動メディアンを用いると、3ヵ月に1回
以下の大きな変動は除去することができる。移動メディ
アンの特徴は、図4に示したようにデータの値が単調に
増加または減少している場合には両端のデータ以外はも
とのデータ値そのものになるという点である。図4の9
0年1月〜7月では、移動メディアンのデータと実績の
データはほとんど一致している。また、移動メディアン
を用いた場合は両端の点(図3および図4では89年1
月と90年8月)の値は定義されないのでデータ数はも
とのデータ数より2つ少なくなることになる。
Next, the first step will be described.
The main purpose of the first step is peak removal. "3
"Monthly mobile median" is a method that uses the median (median) of each month and the data for three months including the month before and after that. A 3-month mobile median will be described with reference to FIGS. 3 and 4. 3 and 4 are 3
It is a figure explaining a mobile median for six months. In FIG. 3, in the case of July 1989, the actual results are in the order of August 1989 <July 1989 <June 1989, so the mobile median value is the July data value itself. In the case of December 1989, the actual results are January 1990 <November 1989 <December 1989, so the mobile median value becomes the data value for November. With a three-month mobile median, large fluctuations less than once every three months can be eliminated. The characteristic of the mobile median is that, as shown in FIG. 4, when the data value is monotonically increasing or decreasing, the original data values are the same except for the data at both ends. 9 in FIG.
From January to July 0, the data of the mobile median and the actual data are almost the same. In addition, when the mobile median is used, dots at both ends (in Fig. 3 and Fig. 4, 89
Since the values for month and August 1990) are not defined, the number of data will be two less than the original number of data.

【0025】次に、第二のステップについて図5を参照
して説明する。図5は最小2乗法を説明する図である。
第二のステップの目的は、第三のステップの初期値を与
えることにある。もとのデータ数をK個とすると、第一
のステップの移動メディアンによって両端の2ヵ月を除
いたK−2個のデータが得られている。これらK−2個
のデータに対して、通常の最小2乗法を用いた回帰によ
ってモデルのあてはめを行う。初めに示した4つのモデ
ルについてそれぞれあてはめを行うので、最小2乗法の
計算は4回行われることになる。最小2乗法は残差の2
乗和の最小化であり、
Next, the second step will be described with reference to FIG. FIG. 5 is a diagram for explaining the least squares method.
The purpose of the second step is to give an initial value for the third step. Assuming that the original number of data is K, K-2 data excluding two months at both ends are obtained by the moving median in the first step. A model is fitted to these K−2 pieces of data by regression using the usual least squares method. Since the fitting is performed for each of the four models shown at the beginning, the least squares calculation is performed four times. The least squares method is the residual 2
Is the minimization of the sum of sums,

【0026】[0026]

【数9】 (ただしr:残差)となり、図5に示すような回帰が得
られる。
[Equation 9] (However, r: residual), and a regression as shown in FIG. 5 is obtained.

【0027】次に、第三のステップについて図6を参照
して説明する。図6は反復重み付け最小2乗法を説明す
るための図である。移動メディアンでは3ヵ月に一度の
ピークは除去できるものの、2ヵ月続く特異なデータに
対しては効果が少ない。図6の販売実績を例にとると1
1月と12月の販売量が他の月と比較して特に大きくな
っている。そのため第二のステップで得られる回帰式
(初期解)はそれらの2ヵ月のデータの影響を大きく受
けていることがわかる。この影響を取り除くために反復
重み付き最小2乗法を用いている。反復重み付き最小2
乗法とは、通常の最小2乗法が残差の2乗和の最小化で
あり、
Next, the third step will be described with reference to FIG. FIG. 6 is a diagram for explaining the iterative weighted least squares method. Although the mobile median can remove the peak once every three months, it has little effect on the unique data that continues for two months. Taking the sales record of Figure 6 as an example, 1
The sales volume in January and December is particularly large compared to other months. Therefore, it can be seen that the regression equation (initial solution) obtained in the second step is greatly affected by those two-month data. In order to remove this effect, the iterative weighted least squares method is used. Repeated weighted minimum 2
The multiplicative method is a method of minimizing the sum of squares of residuals by using the ordinary least square method.

【0028】[0028]

【数10】 (ただしr:残差)をもとにしているのに対して、(基
準化)残差の適当な偶関数ρの和の最小化
[Equation 10] (However, r: residual) based on the other hand, (normalized) minimize the sum of appropriate even function ρ of residuals

【0029】[0029]

【数11】 (ただしu:基準化残差uk =rk /σ,σ:標準偏差
の推定値)をもとにしている。この式は便宜的に、
[Equation 11] (However, u: standardized residual u k = r k / σ, σ: estimated value of standard deviation). For convenience, this formula

【0030】[0030]

【数12】 (ただしw:重み関数)という重み付き最小2乗法に書
き直すことができる。今回ρにはbiweight(ま
たはbisquare)とよばれている、 ρ(u)=(B2 /6)(1−(1−(u/B)2 3 ) if|u|<B または、 ρ(u)=(B2 /6) otherwise (ただしB:定数(本発明実施例では6.0))を用い
ており、この場合重み関数wはρの導関数ρ′から w(u)= ρ′(u)/u=(1−(u/B)2 2 if|u|<B または、 ρ′(u)/u=0 otherwise となる。今回の場合は定数Bを6.0にしているので、
回帰直線から標準偏差の6倍以上離れたデータについて
は重みが0、つまり実質的にはそのデータが除去された
ことと同じことになる。重みを決定するためには基準化
残差が必要であり、基準化残差を求めるためには標準偏
差の推定が必要になる。標準偏差の推定は初期値(移動
メディアン値に対する回帰)の場合には、
[Equation 12] (However, w: weight function) can be rewritten as a weighted least squares method. The time [rho is called Biweight (or bisquare), ρ (u) = (B 2/6) (1- (1- (u / B) 2) 3) if | u | <B or, [rho (u) = (B 2/ 6) otherwise ( except B: constant (6.0) in the present invention embodiment) is used, in which case the weight function w is from the derivative [rho 'of [rho w (u) = ρ ′ (u) / u = (1- (u / B) 2 ) 2 if | u | <B or ρ ′ (u) / u = 0 otherwise. In this case, the constant B is set to 6.0, so
For data that is more than 6 times the standard deviation from the regression line, the weight is 0, which means that the data is substantially removed. A standardized residual is required to determine the weight, and a standard deviation must be estimated to obtain the standardized residual. If the standard deviation is estimated as an initial value (regression with respect to the moving median value),

【0031】[0031]

【数13】 (ただしr′:移動メディアンの値と回帰式の値の差,
K:データ数(月数))を用いている。第二のステップ
の最小2乗法は元データではなく移動メディアン値に対
して行われるので、標準偏差の推定値は便宜的に通常の
場合の式の2倍にしてある。回帰式が決定した後では、
[Equation 13] (However, r ': difference between the value of the moving median and the value of the regression equation,
K: number of data (months)) is used. Since the second-step least-squares method is performed on the moving median value, not the original data, the estimated standard deviation is expediently doubled from the usual case. After the regression equation is determined,

【0032】[0032]

【数14】 を標準偏差の推定値として用いている。標準偏差の推定
値は各データの重みが変われば変化する。また逆に、各
データの重みは標準偏差の推定値が変われば変化する。
そこで、これらの変化が充分小さくなるまで反復計算を
行う。今回の場合には収束判定条件として「データの重
みの変化が充分小さい」という意味で
[Equation 14] Is used as an estimate of standard deviation. The estimated value of the standard deviation changes as the weight of each data changes. On the contrary, the weight of each data changes if the estimated value of the standard deviation changes.
Therefore, iterative calculation is performed until these changes are sufficiently small. In this case, it means that the change of the weight of the data is small enough as the convergence judgment condition.

【0033】[0033]

【数15】 (ただし wk :回帰式を求めるために用いた重み,
w′k :得られた回帰式から計算される重み)を用いて
いる。
[Equation 15] (However, w k : the weight used to obtain the regression equation,
w ′ k : weight calculated from the obtained regression formula) is used.

【0034】次に、図6の販売実績の重みの変化を図7
を参照して説明する。図7は重みの変化を説明する図で
ある。図6にモデル1をあてはめた場合の重み付けの変
遷を図7に示している。図7(a)〜(d)には実績、
トレンド成分(つまり重み1の線)、および重み0.
9、0.5、0の等高線が示してある。図7(a)か
ら、移動メディアン値に対する回帰式から得られた初期
重みでは、12月のデータの重みは0にはなっていない
ことがわかる。また、11月のデータの重みは0.9以
上である。その「初期重み」を重みとした重み付き最小
2乗法の結果が図7(b)の「反復1」のトレンド成分
である。この結果から重みを計算し直すと、12月のデ
ータでは重みが0になっており、また11月のデータの
重みも低くなっていることがわかる。反復計算の過程を
追っていくと、11月のデータは図7(d)の3回目の
反復計算の結果から得られる重みが0になっており、そ
れ以降の計算には影響していないことがわかる。この場
合は図7(e)の5回目の反復で収束条件が成立し、反
復計算が終了する。移動メディアン値に対して第二のス
テップまでの回帰を行った場合と比較して、11月、1
2月といった特異なデータの影響が充分小さくなってい
ることがわかる。第二のステップが各モデルに対して計
算を行ったのと同様に、第三のステップも各モデルに対
して計算が行われる。当然ながら各モデルに対する初期
解は、それぞれのモデルについて第二のステップで計算
された回帰式を用いる。
Next, the change in the weight of the sales record of FIG. 6 is shown in FIG.
Will be described. FIG. 7 is a diagram illustrating a change in weight. FIG. 7 shows changes in weighting when the model 1 is applied to FIG. 7 (a)-(d),
Trend components (that is, lines with weight 1) and weights 0.
The contour lines at 9, 0.5, 0 are shown. From FIG. 7A, it can be seen that the initial weight obtained from the regression equation for the moving median value does not have a weight of 0 in December. Further, the weight of the data in November is 0.9 or more. The result of the weighted least squares method using the “initial weight” as a weight is the trend component of “iteration 1” in FIG. 7B. When the weights are recalculated from this result, it can be seen that the weight of the December data is 0 and the weight of the November data is also low. As the iterative calculation process is followed, the data for November has a weight of 0 obtained from the result of the third iterative calculation in FIG. 7 (d), which indicates that it does not affect subsequent calculations. Recognize. In this case, the convergence condition is satisfied in the fifth iteration of FIG. 7E, and the iterative calculation is completed. Compared with the case where the regression up to the second step was performed on the moving median value, November, 1
It can be seen that the influence of unique data such as February is sufficiently small. Just as the second step performed calculations for each model, the third step also performed calculations for each model. Of course, the initial solution for each model uses the regression equation calculated in the second step for each model.

【0035】次に、モデル選択について説明する。それ
ぞれフィッティングを行った4つのモデルの中からどの
モデルを選択するのかという問題にはAIC(赤池の情
報量基準)を用いている。今回の場合、AICは、
Next, model selection will be described. AIC (Akaike's information criterion) is used for the problem of which model is selected from the four fitted models. In this case, the AIC

【0036】[0036]

【数16】 (ただしp:モデルのパラメータの数)でモデル1の場
合4個、モデル2の場合3個、モデル3の場合2個、モ
デル4の場合1個という式で表される。4つのモデルの
それぞれの収束値について上式の値を計算し、その値が
最小になるモデルを選択する。この選択の結果により選
ばれたモデルの回帰式がトレンド成分として用いられ
る。さらに選択された回帰式にK+1を代入した値が次
月の予想値として計算される。
[Equation 16] (Where p is the number of parameters of the model), it is represented by an equation of 4 in the case of model 1, 3 in the case of model 2, 2 in the case of model 3, and 1 in the case of model 4. The value of the above equation is calculated for each convergence value of the four models, and the model having the smallest value is selected. The regression equation of the model selected by the result of this selection is used as the trend component. Further, a value obtained by substituting K + 1 into the selected regression equation is calculated as the expected value for the next month.

【0037】次に、このトレンド成分を用いて変動成分
を計算し、変動成分予測値を算出する手順を図8を参照
して説明する。図8は変動成分予測値算出を説明するた
めの図である。図8(a)にトレンド成分Tと売上実績
データの関係を示した(T'(k+1):トレンド成分
Tの変動成分予測値)。変動成分Hは、 変動成分H=(売上実績データ)−(トレンド成分T) という関係を用いて計算し、図8(b)に示すようなグ
ラフを得る。さらに変動成分予測を算出する手段は、図
8(c)に示すように変動成分データを等間隔の9個の
クラスタに分け、今月の変動成分データがあるクラスタ
内における過去の変動成分データの次月への変動量の平
均を求める。そして、その値を今月の変動成分データに
加え、次月の変動成分の予測値とする。図8(c)の場
合には今月の変動成分H(k)はクラスタ「8」に属す
るため、今月の変動成分H(k)から変動成分予測値H
(k+1)への変動量Δ8は、過去にクラスタ「8」に
属する変動成分の変動量(Δa、Δb)の平均値とな
り、 Δ8=(Δa+Δb)/2 である。従って、変動成分予測値H' (k+1)は、 H' (k+1)=H(k)+Δ8 となる。各クラスタごとに変動量の平均Δ1〜Δ9の値
が定まるが今月の変動成分が0の近傍にある場合、変動
成分を予測すると逆に予測精度が劣化するというシミュ
レーション結果が得られた。従って、変動成分0のある
クラスタとその近傍のクラスタ、図8(c)では「6」
とその近傍のクラスタ「5」および「7」に今月の変動
成分があるときは変動成分予測値は0とし、 次月の売上予測値=トレンド成分予測値 とすることにする。クラスタiの変動量Δiが図8
(d)に示すような場合、 Δi=(Δa+Δb+Δc+Δd)/4 としていたが、これでは売上が増加傾向の場合も減少傾
向の場合もいずれも同じ予測値を与えることになる。そ
こでクラスタiにある今月の変動成分が前月に比して増
加傾向の場合と減少傾向の場合の予測値を別々に与える
ことにする。つまり、今月の変動成分H(k)が前月の
変動成分H(k−1)から増加傾向にあるときは、増加
傾向でクラスタiに存在する変動成分を次月への変動量
の平均値とする。図8(d)のクラスタiを見ると、前
月から今月へは増加傾向にあるので、過去の同クラスタ
内の増加傾向の平均値Δi′を求めると、 Δi′=(Δb+Δd)/2 となる。従って、変動成分予測値は、 H' (k+1)=H(k)+Δi′ となる。このようにして、次月の売上予測値はトレンド
成分予測値と変動成分予測値の和として与えられる。
Next, the procedure for calculating the fluctuation component using this trend component and calculating the fluctuation component predicted value will be described with reference to FIG. FIG. 8 is a diagram for explaining the calculation of the fluctuation component predicted value. FIG. 8A shows the relationship between the trend component T and the sales result data (T ′ (k + 1): the predicted variation component of the trend component T). The fluctuation component H is calculated using the relationship of fluctuation component H = (sales performance data) − (trend component T), and a graph as shown in FIG. 8B is obtained. Further, the means for calculating the variation component prediction divides the variation component data into nine clusters at equal intervals as shown in FIG. 8C, and next to the past variation component data in the cluster having the variation component data of this month. Calculate the average of fluctuations over the month. Then, the value is added to the fluctuation component data of this month to be the predicted value of the fluctuation component of the next month. In the case of FIG. 8C, since the fluctuation component H (k) of this month belongs to the cluster “8”, the fluctuation component prediction value H is calculated from the fluctuation component H (k) of this month.
The fluctuation amount Δ8 to (k + 1) is the average value of the fluctuation amounts (Δa, Δb) of the fluctuation components belonging to the cluster “8” in the past, and Δ8 = (Δa + Δb) / 2. Therefore, the fluctuation component predicted value H ′ (k + 1) is H ′ (k + 1) = H (k) + Δ8. The values of averages Δ1 to Δ9 of the fluctuation amount are determined for each cluster, but when the fluctuation component of this month is in the vicinity of 0, the prediction result of the fluctuation component conversely deteriorates the prediction result. Therefore, the cluster having the fluctuation component 0 and the clusters in the vicinity thereof, “6” in FIG.
If there are fluctuation components of this month in the clusters “5” and “7” in the vicinity of and, the fluctuation component prediction value is set to 0, and the sales prediction value of the next month = trend component prediction value. The fluctuation amount Δi of the cluster i is shown in FIG.
In the case shown in (d), Δi = (Δa + Δb + Δc + Δd) / 4 was set, but this gives the same predicted value regardless of whether the sales tend to increase or decrease. Therefore, the prediction values for the case where the fluctuation component of this month in cluster i has an increasing tendency and a decreasing tendency as compared with the previous month will be given separately. That is, when the fluctuation component H (k) of this month has an increasing tendency from the fluctuation component H (k−1) of the previous month, the fluctuation component existing in the cluster i in the increasing tendency is set as the average value of the fluctuation amounts to the next month. To do. Looking at the cluster i in FIG. 8D, there is an increasing tendency from the previous month to this month. Therefore, if the average value Δi ′ of the increasing tendency in the same cluster in the past is obtained, Δi ′ = (Δb + Δd) / 2. .. Therefore, the fluctuation component predicted value is H ′ (k + 1) = H (k) + Δi ′. In this way, the sales forecast value for the next month is given as the sum of the trend component forecast value and the fluctuation component forecast value.

【0038】ここで、図9を参照して総生産量と遷移確
率の関係を説明する。図9は変動成分と遷移確率を示す
図である。これは図9の左側に変動成分の推移を示し、
右側に第k期時点での遷移確率分布を模式的に示す。こ
の図では、縦軸に変動成分の変位量を取り、横軸に遷移
確率を示しており、構成要素m、n(品目に対応)の確
率分布をそれぞれρm 、ρn とした時の最大確率を与え
るときの変動成分の変位量をXm 、Xn として示す。つ
まり、このXm 、Xn は変動成分の予測値そのものであ
る。通常、総生産量Gは予測値Xiの合計と等しくなる
とは限らず、
Here, the relationship between the total production amount and the transition probability will be described with reference to FIG. FIG. 9 is a diagram showing variation components and transition probabilities. This shows the transition of the fluctuation component on the left side of FIG.
The transition probability distribution at the time point k is schematically shown on the right side. In this figure, the vertical axis is the displacement of the fluctuation component, and the horizontal axis is the transition probability. The maximum when the probability distributions of the constituent elements m and n (corresponding to the item) are ρ m and ρ n , respectively. The displacement amounts of the fluctuation components when the probability is given are shown as X m and X n . That is, these X m and X n are the predicted values of the fluctuation component itself. Normally, the total production G is not always equal to the sum of the predicted values Xi,

【0039】[0039]

【数17】 となる。[Equation 17] Becomes

【0040】次に、図10を参照して遷移確率分布関数
を三角形で表す手順を説明する。図10は遷移確率分布
関数を三角形で表す手順を説明する図である。図9に示
したような遷移確率と変位の関係は図10(a)に示す
ような三角形であらわされる。次に、図10(b)に示
すように総生産量Gと予測値合計との差分を各構成要素
i(i∈N)に最適配分する。
Next, the procedure for representing the transition probability distribution function by a triangle will be described with reference to FIG. FIG. 10 is a diagram illustrating a procedure of representing a transition probability distribution function by a triangle. The relationship between the transition probability and the displacement as shown in FIG. 9 is represented by a triangle as shown in FIG. Next, as shown in FIG. 10B, the difference between the total production amount G and the total predicted value is optimally distributed to each component i (iεN).

【0041】次に、図11を参照して遷移確率分布関数
から最適配分を決定する手順を説明する。図11は遷移
確率分布関数から最適配分を決定する手順を説明する図
である。ここで、 ρi(ΔXi)=yi とおくことにより、遷移確率の増減を方程式で表す。こ
れらの遷移確率分布関数を用いて総生産量1000単位
を達成するように線形最適問題を解く、 ただし、G−ΣXi>0ならばai=−1/ΔXi ma
x G−ΣXi<0ならばai=−1/ΔXi min (ここでΔXi maxは変動成分の最大値、ΔXi mi
nは変動成分の最小値を表す)として、図11(a)の
ような形で表わせば遷移確率分布関数は、 となる。ここで図11(b)に示すように、 y1=y2=…=yn=α (0<α≦1) として、これらの式を満たすΔXi を T−ΣXi =ΣΔXi の制約の中で求めればよい。
Next, the procedure for determining the optimum distribution from the transition probability distribution function will be described with reference to FIG. FIG. 11 is a diagram illustrating a procedure for determining the optimum distribution from the transition probability distribution function. Here, by setting ρi (ΔXi) = yi, an increase / decrease in transition probability is expressed by an equation. A linear optimal problem is solved using these transition probability distribution functions so as to achieve a total production amount of 1000 units, provided that if G-ΣXi> 0, then ai = -1 / ΔX i ma
If x G-ΣXi <0, then ai = -1 / ΔX i min (where ΔX i max is the maximum value of the fluctuation component, ΔX i mi
(n represents the minimum value of the fluctuation component), the transition probability distribution function is expressed as shown in FIG. Becomes Here, as shown in FIG. 11B, assuming that y1 = y2 = ... = yn = α (0 <α ≦ 1), ΔX i satisfying these equations is set within the constraint of T−ΣX i = ΣΔX i. Just ask.

【0042】以上のように、変動成分の系列のうち、変
動成分の最大値および最小値を底辺とし、図8に示す過
程により算出した変動成分予測値を頂点とする三角形を
構成する。この三角形が、遷移確率ρを与える遷移確率
分布関数を近似し、構成要素毎(本実施例では品目に相
当する)にこれを決定する。このようにして決定した遷
移確率分布関数をもとに、前述の最適配分法によって予
測値の合計と売上目標金額の差分を図10(b)に示す
ように構成要素に振り分ける。この場合の遷移確率分布
関数は売行きの感度を表わす関数ということになる。な
お、遷移確率分布関数を三角形で定義することにより、
誰にもわかりやすい形つまりデータメンテナンスに有利
な形にすることが可能であるばかりでなく、これにより
演算の高速化を実現できる。
As described above, a triangle having the maximum and minimum values of the fluctuation component as the bases and the fluctuation component predicted value calculated by the process shown in FIG. This triangle approximates the transition probability distribution function that gives the transition probability ρ, and determines this for each component (corresponding to the item in this embodiment). Based on the transition probability distribution function determined in this way, the difference between the total of the predicted values and the sales target amount is distributed to the constituent elements as shown in FIG. In this case, the transition probability distribution function is a function that represents sales sensitivity. By defining the transition probability distribution function as a triangle,
Not only can it be made in a form that is easy for anyone to understand, that is, a form that is advantageous for data maintenance, but this also makes it possible to speed up computation.

【0043】次に、4品目の製品の変動成分データに対
する、三角形で表した遷移確率分布関数を図12を参照
して説明する。図12は4品目の製品の変動成分データ
と遷移確率分布関数の三角形の関係を示す図である。図
12に示すように、変動成分データから三角形で表され
た遷移確率分布関数が得られるが、品目1、2、4と品
目3の三角形を比較したとき、品目3の三角形は底辺が
長く、変動の振幅が大きいことがうかがえる。すなわ
ち、販売戦略によっては予測値以上に売れる可能性がま
だかなりあるということが推測できる。従って品目3に
予測値と売上目標との差分が大きな割合で割り振れるこ
とが推測できる。一方、何らかの理由(例えば、生産力
の向上、新規市場の開拓、他社の動向等)によって、遷
移確率分布が今後の状況にそぐわなくなった場合、それ
に対応する手段が必要となる。図12(a)は本発明実
施例装置にて、各品目の変動成分および遷移確率分布関
数を表示したときのイメージである。図12(b)に遷
移確率分布関数を90度回転した図を示す。ここで、直
接的に遷移確率分布関数(ピークの位置や変位幅)を変
えることにより、前述の対応手段とすることが可能とな
る。また、この機能は最適化を図る上でのシミュレーシ
ョン機能でもあり、データ量の少ない系列(新製品、季
節限定商品)の予測に対しても有効である。このように
して遷移過程の情報を有することから、種々の特性を持
った構成要素への最適配分が可能になった。図13は新
手法による販売予定予測を示す図であるがその結果、図
13に示すように販売予定の予測精度が向上し、販売予
定策定までの工数が短縮できる。
Next, the transition probability distribution function represented by a triangle with respect to the fluctuation component data of the four products will be described with reference to FIG. FIG. 12 is a diagram showing the relationship between the variation component data of four products and the triangle of the transition probability distribution function. As shown in FIG. 12, a transition probability distribution function represented by a triangle is obtained from the variation component data. When comparing the triangles of item 1, 2, 4 and item 3, the triangle of item 3 has a long base, It can be seen that the amplitude of fluctuation is large. That is, it can be inferred that there is still a considerable possibility that the sales will exceed the predicted value depending on the sales strategy. Therefore, it can be estimated that the difference between the predicted value and the sales target is assigned to item 3 at a large rate. On the other hand, if for some reason (for example, productivity improvement, new market development, trends of other companies, etc.), the transition probability distribution does not fit the future situation, a means for responding to it is necessary. FIG. 12A is an image when the fluctuation component and the transition probability distribution function of each item are displayed on the device of the present invention. FIG. 12B shows a diagram obtained by rotating the transition probability distribution function by 90 degrees. Here, by directly changing the transition probability distribution function (peak position and displacement width), it becomes possible to use the above-mentioned coping means. This function is also a simulation function for optimization, and is also effective for forecasting series (new products, seasonal products) with a small amount of data. In this way, since the information of the transition process is included, it becomes possible to optimally distribute the components to various characteristics. FIG. 13 is a diagram showing a sales schedule forecast by the new method. As a result, as shown in FIG. 13, the forecast accuracy of the sales schedule is improved, and the man-hour until the sales schedule is formulated can be shortened.

【0044】[0044]

【発明の効果】以上説明したように、本発明は販売予測
あるいは生産管理のために、変動の遷移する過程に基づ
き演算を行う遷移確率分布関数、あるいはメンバーシッ
プ関数が有用な論理であることに着目し、これを実用的
な形に単純化して、投入するデータの種類や量が小さ
く、操作が簡単であり、誤りの発生も小さく、最終的に
人が予測判断をするために便利である。遷移確率分布関
数を実用的な形に単純化して、投入するデータの種類や
量が小さく、操作も簡単であり、誤りの発生も少なく、
種々の特性を持った構成要素の販売予測に基づいて、目
標値の品目ごとの最適配分が可能となり、その結果、販
売予定の予測精度が向上し、販売予定策定までの工数が
短縮できる。
As described above, according to the present invention, the transition probability distribution function or the membership function, which is operated based on the transition process of fluctuation, is a useful logic for sales forecasting or production management. Focusing on this, simplifying this to a practical form, the type and amount of data to be input is small, the operation is easy, the occurrence of errors is small, and it is convenient for people to finally make prediction judgments. .. The transition probability distribution function is simplified to a practical form, the type and amount of data to be input is small, the operation is easy, and the number of errors is small,
Based on the sales forecast of components having various characteristics, it becomes possible to optimally distribute the target value for each item. As a result, the forecast accuracy of the sales schedule is improved, and the man-hours until the formulation of the sales schedule can be shortened.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明実施例装置のブロック図。FIG. 1 is a block diagram of an apparatus according to an embodiment of the present invention.

【図2】本発明実施例装置の動作を示すフローチャー
ト。
FIG. 2 is a flowchart showing the operation of the apparatus according to the present invention.

【図3】移動メディアンを説明する図。FIG. 3 is a diagram illustrating a mobile median.

【図4】移動メディアンを説明する図。FIG. 4 is a diagram illustrating a mobile median.

【図5】最小2乗法を説明する図。FIG. 5 is a diagram illustrating a least squares method.

【図6】反復重み付け最小2乗法を説明する図。FIG. 6 is a diagram illustrating an iterative weighted least squares method.

【図7】重み付けの変遷を説明する図。FIG. 7 is a diagram for explaining changes in weighting.

【図8】変動成分予測値算出を説明する図。FIG. 8 is a diagram illustrating calculation of a fluctuation component predicted value.

【図9】変動成分と遷移確率を示す図。FIG. 9 is a diagram showing fluctuation components and transition probabilities.

【図10】遷移確率分布関数を三角形で表わす手順を説
明する図。
FIG. 10 is a diagram illustrating a procedure of representing a transition probability distribution function by a triangle.

【図11】遷移確率分布関数から最適配分を決定する手
順を説明する図。
FIG. 11 is a diagram illustrating a procedure for determining an optimum distribution from a transition probability distribution function.

【図12】4品目の製品の変動成分データと遷移確率分
布関数の三角形の関係を示す図。
FIG. 12 is a diagram showing a relationship between variation component data of four products and a triangle of a transition probability distribution function.

【図13】新手法による販売予定予測を示す図。FIG. 13 is a diagram showing a forecast of sales schedule by a new method.

【符号の説明】[Explanation of symbols]

1 実績記憶ファイル部 2 トレンド成分演算予測部 3 変動成分演算予測部 4 総合予測演算部 10 CRT 1 Actual result storage file part 2 Trend component calculation prediction part 3 Fluctuation component calculation prediction part 4 Total prediction calculation part 10 CRT

Claims (7)

【特許請求の範囲】[Claims] 【請求項1】 1つまたは複数の商品について販売実績
が蓄積された実績記憶ファイルと、 この実績記憶ファイルに蓄積された前記販売実績を取り
出しこの販売実績から周期的な変動であるトレンド成分
を抽出するトレンド成分演算手段と、 前記販売実績から前記トレンド成分を除去して変動成分
を抽出する変動成分演算手段と、 次に到来する期間についてこの変動成分の予測を行う変
動成分予測手段と、 前記次に到来する期間について前記トレンド成分を外挿
することによりトレンド成分の予測を行うトレンド成分
予測手段と、 このトレンド成分予測手段により予測されたトレンド成
分予測値と前記変動成分予測手段により予測された変動
成分予測値との和を前記次に到来する期間についての総
合予測値とする予測演算手段とを備えた販売予測支援装
置において、 前記変動成分予測手段は、 前記実績記憶ファイルから取り出した実績値を複数個の
クラスタに区分し、最新実績値の変動成分H(k)がi
番目のクラスタにあり、i番目のクラスタから変動した
実績値をΔaj(j=1、2、…、N)とするとき、そ
の平均値を 【数1】 とし、変動成分予測値を H' (k+1) = H(k) + Δi として演算する手段を備えたことを特徴とする販売予測
支援装置。
1. A result storage file in which sales results of one or a plurality of products are accumulated, and the sales results accumulated in the result storage file are taken out, and a trend component that is a periodic fluctuation is extracted from the sales results. A trend component calculating means, a fluctuation component calculating means for removing the trend component from the sales record to extract a fluctuation component, a fluctuation component predicting means for predicting the fluctuation component for the next coming period, Trend component predicting means for predicting the trend component by extrapolating the trend component for the period of arrival, and the trend component predicted value predicted by the trend component predicting means and the fluctuation predicted by the fluctuation component predicting means. Sales including a prediction calculation means for making the sum of the component predicted values into a total predicted value for the next coming period In predicting support device, the fluctuation component predicting means, dividing the actual value extracted from the actual storage file into a plurality of clusters, the variation component of the latest actual value H (k) is i
When the actual value that is in the i-th cluster and fluctuates from the i-th cluster is Δaj (j = 1, 2, ..., N), the average value is And a means for calculating the fluctuation component predicted value as H '(k + 1) = H (k) + Δi.
【請求項2】 前記実績値、前記トレンド成分、および
前記変動成分を時間軸上に表示するとともに、前記総合
予測値を頂点位置とし前記変動成分演算手段により抽出
された変動成分の振幅を底辺長さとし次に到来する期間
の長さを高さとする三角形で表示する手段を備えた請求
項1記載の販売予測支援装置。
2. The actual value, the trend component, and the fluctuation component are displayed on a time axis, and the amplitude of the fluctuation component extracted by the fluctuation component computing means is the base length with the comprehensive predicted value as a vertex position. 2. The sales forecasting support device according to claim 1, further comprising means for displaying a triangle having a height of the next coming period as a height.
【請求項3】 前記最新実績値の変動成分H(k)が変
動成分が零であるクラスタ内にあるときには変動成分予
測値を零として、前記トレンド成分予測値を前記総合予
測値とする請求項1記載の販売予測支援装置。
3. The fluctuation component prediction value is set to zero and the trend component prediction value is set to the total prediction value when the fluctuation component H (k) of the latest actual value is in a cluster in which the fluctuation component is zero. 1. The sales forecast support device described in 1.
【請求項4】 1つまたは複数の商品について販売実績
が蓄積された実績記憶ファイルと、 この実績記憶ファイルに蓄積された前記販売実績から季
節変動に伴うトレンド成分を抽出するトレンド成分演算
手段と、 前記販売実績から前記トレンド成分を除去して変動成分
を抽出する変動成分演算手段と、 次に到来する期間についてこの変動成分の予測を行う変
動成分予測手段と、 前記次に到来する期間について前記トレンド成分を外挿
することによりトレンド成分の予測を行うトレンド成分
予測手段と、 このトレンド成分予測手段により予測されたトレンド成
分予測値と前記変動成分予測手段により予測された変動
成分予測値との和を前記次に到来する期間についての総
合予測値とする予測演算手段とを備えた販売予測支援装
置において、 前記変動成分予測手段は、 前記実績記憶ファイルから取り出した実績値を複数個の
クラスタに区分し、最新実績値の変動成分H(k)がi
番目のクラスタにあり、その変動成分H(k)がその1
期前の変動成分実績値H(k−1)に対し増加傾向にあ
るとき、 i番目のクラスタから増加変動した実績値をΔcj(j
=1、2、…、N)として、その増加変動した実績値の
平均値を 【数2】 とし、その変動成分H(k)がその1期前の変動成分実
績値H(k−1)に対し減少傾向にあるとき、 i番目のクラスタから減少変動した実績値をΔej(j
=1、2、…、N)として、その減少変動した実績値の
平均値を 【数3】 とし、変動成分予測値を H' (k+1) = H(k) + Δi として演算する手段を備えたことを特徴とする販売予測
支援装置。
4. A record storage file in which sales records of one or a plurality of products are stored, and trend component calculation means for extracting trend components associated with seasonal fluctuations from the sales records stored in the record storage file. A fluctuation component calculation means for removing the trend component from the sales record to extract a fluctuation component, a fluctuation component prediction means for predicting the fluctuation component for the next coming period, and a trend for the next coming period. The trend component predicting means for predicting the trend component by extrapolating the component, and the sum of the trend component predicting value predicted by the trend component predicting means and the fluctuation component predicting value predicted by the fluctuation component predicting means. In a sales forecast support device including a forecast calculation unit that provides a comprehensive forecast value for the next coming period, The variation component predicting unit divides the performance value extracted from the performance storage file into a plurality of clusters, and the fluctuation component H (k) of the latest performance value is i.
In the second cluster, whose variation component H (k) is 1
When the fluctuation component actual value H (k-1) before the period is on the increase, the actual value increased from the i-th cluster is changed by Δcj (j
= 1, 2, ..., N), the average value of the increased and changed actual values is When the fluctuation component H (k) has a decreasing tendency with respect to the fluctuation component actual value H (k−1) of the preceding period, the actual value decreased from the i-th cluster is Δej (j).
= 1, 2, ..., N), the average value of the reduced and changed actual values is given by And a means for calculating the fluctuation component predicted value as H '(k + 1) = H (k) + Δi.
【請求項5】 前記実績値、前記トレンド成分、および
前記変動成分を時間軸上に表示するとともに、前記総合
予測値を頂点位置とし前記変動成分演算手段により抽出
された変動成分の振幅を底辺長さとし次に到来する期間
の長さを高さとする三角形で表示する手段を備えた請求
項4記載の販売予測支援装置。
5. The actual value, the trend component, and the fluctuation component are displayed on a time axis, and the amplitude of the fluctuation component extracted by the fluctuation component calculation means is the base length with the comprehensive predicted value as the vertex position. 5. The sales forecasting support device according to claim 4, further comprising means for displaying a triangle whose height is the length of the next coming period.
【請求項6】 前記最新実績値の変動成分H(k)が変
動成分が零であるクラスタ内にあるときには H' (k+1) = 0 として、前記トレンド成分予測値を前記総合予測値とす
る請求項4記載の販売予測支援装置。
6. When the fluctuation component H (k) of the latest actual value is in a cluster in which the fluctuation component is zero, H ′ (k + 1) = 0 is set, and the trend component prediction value is set as the comprehensive prediction value. Item 4. A sales forecast support device according to item 4.
【請求項7】 複数の商品品目について前記総合予測
値、前記変動成分予測値、および前記三角形を演算する
手段を備えた請求項1ないし2のいずれかに記載の販売
予測支援装置。
7. The sales forecast support apparatus according to claim 1, further comprising means for calculating the total forecast value, the fluctuation component forecast value, and the triangle for a plurality of product items.
JP11768292A 1992-05-11 1992-05-11 Sales prediction supporting device Pending JPH05314094A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP11768292A JPH05314094A (en) 1992-05-11 1992-05-11 Sales prediction supporting device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP11768292A JPH05314094A (en) 1992-05-11 1992-05-11 Sales prediction supporting device

Publications (1)

Publication Number Publication Date
JPH05314094A true JPH05314094A (en) 1993-11-26

Family

ID=14717684

Family Applications (1)

Application Number Title Priority Date Filing Date
JP11768292A Pending JPH05314094A (en) 1992-05-11 1992-05-11 Sales prediction supporting device

Country Status (1)

Country Link
JP (1) JPH05314094A (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08166465A (en) * 1994-12-12 1996-06-25 East Japan Railway Co Wind velocity prediction method and traffic means operating method based thereon
JP2003203138A (en) * 2001-12-28 2003-07-18 Toshiba Corp Environmental influence assessment method and program
JP2004234471A (en) * 2003-01-31 2004-08-19 Ricoh Co Ltd Demand prediction method and demand prediction program
US6820060B1 (en) * 1996-06-24 2004-11-16 Jack Eisner Apparatus for generating sales probability
JP2006011542A (en) * 2004-06-22 2006-01-12 Ns Solutions Corp Demand predicting device, demand predicting method and program
US7054797B2 (en) 2000-06-30 2006-05-30 Kabushiki Kaisha Toshiba Apparatus for environmental impact estimation and method and program stored in a computer readable medium for executing the same
JP2007072844A (en) * 2005-09-08 2007-03-22 Osaka Gas Co Ltd Display system of forecast information
JP2007293624A (en) * 2006-04-25 2007-11-08 Ricoh Co Ltd Demand prediction method and demand prediction program
JP2008234094A (en) * 2007-03-19 2008-10-02 Nippon Telegr & Teleph Corp <Ntt> Prediction model selection method, prediction method, parameter estimation device, and its program
US7848943B2 (en) 2005-12-22 2010-12-07 International Business Machines Corporation System and method for supporting purchase or production of products by potential demand prediction
JP2012150537A (en) * 2011-01-17 2012-08-09 Shimizu Corp Rent calculation system and rent calculation program
WO2016072474A1 (en) * 2014-11-06 2016-05-12 エコノミックインデックス株式会社 Information processing device, method, and program
WO2020240994A1 (en) * 2019-05-24 2020-12-03 株式会社日立製作所 Reservation contracting system and reservation contracting method
JP2021135586A (en) * 2020-02-25 2021-09-13 株式会社日立製作所 Time series data prediction device and time series data prediction method
US11182808B2 (en) 2019-02-05 2021-11-23 Target Brands, Inc. Method and system for attributes based forecasting

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08166465A (en) * 1994-12-12 1996-06-25 East Japan Railway Co Wind velocity prediction method and traffic means operating method based thereon
US6820060B1 (en) * 1996-06-24 2004-11-16 Jack Eisner Apparatus for generating sales probability
US7054797B2 (en) 2000-06-30 2006-05-30 Kabushiki Kaisha Toshiba Apparatus for environmental impact estimation and method and program stored in a computer readable medium for executing the same
US7219042B2 (en) 2000-06-30 2007-05-15 Kabushiki Kaisha Toshiba Apparatus for environmental impact estimation and method and program stored in a computer readable medium for executing the same
JP2003203138A (en) * 2001-12-28 2003-07-18 Toshiba Corp Environmental influence assessment method and program
US7412365B2 (en) 2001-12-28 2008-08-12 Kabushiki Kaisha Toshiba Environmental impact estimation method and apparatus
JP2004234471A (en) * 2003-01-31 2004-08-19 Ricoh Co Ltd Demand prediction method and demand prediction program
JP4559776B2 (en) * 2004-06-22 2010-10-13 新日鉄ソリューションズ株式会社 Demand forecasting device, demand forecasting method and program
JP2006011542A (en) * 2004-06-22 2006-01-12 Ns Solutions Corp Demand predicting device, demand predicting method and program
JP2007072844A (en) * 2005-09-08 2007-03-22 Osaka Gas Co Ltd Display system of forecast information
US7848943B2 (en) 2005-12-22 2010-12-07 International Business Machines Corporation System and method for supporting purchase or production of products by potential demand prediction
JP2007293624A (en) * 2006-04-25 2007-11-08 Ricoh Co Ltd Demand prediction method and demand prediction program
JP2008234094A (en) * 2007-03-19 2008-10-02 Nippon Telegr & Teleph Corp <Ntt> Prediction model selection method, prediction method, parameter estimation device, and its program
JP2012150537A (en) * 2011-01-17 2012-08-09 Shimizu Corp Rent calculation system and rent calculation program
WO2016072474A1 (en) * 2014-11-06 2016-05-12 エコノミックインデックス株式会社 Information processing device, method, and program
US11182808B2 (en) 2019-02-05 2021-11-23 Target Brands, Inc. Method and system for attributes based forecasting
WO2020240994A1 (en) * 2019-05-24 2020-12-03 株式会社日立製作所 Reservation contracting system and reservation contracting method
JP2021135586A (en) * 2020-02-25 2021-09-13 株式会社日立製作所 Time series data prediction device and time series data prediction method

Similar Documents

Publication Publication Date Title
US11398000B2 (en) Methods and systems for machine-learning for prediction of grid carbon emissions
JP6034890B2 (en) Specification, estimation, causal driver discovery and market response elasticity or lift coefficient automation
CA2153976C (en) System and method for the advanced prediction of weather impact on managerial planning applications
US8014983B2 (en) Computer-implemented system and method for storing data analysis models
Cheng et al. Estimate at completion for construction projects using evolutionary support vector machine inference model
US7251589B1 (en) Computer-implemented system and method for generating forecasts
US8082175B2 (en) System and method for optimization of a promotion plan
JPH05314094A (en) Sales prediction supporting device
US20030110066A1 (en) Generating an optimized pricing plan
CN109858706A (en) A kind of project cost data management system and method
US20030101087A1 (en) Lease rent optimizer revenue management system
US20140310059A1 (en) System , method and computer program forecasting energy price
US20040117236A1 (en) Automated optimization tool for electric utility sypply services
US11182808B2 (en) Method and system for attributes based forecasting
WO2010019897A1 (en) Automatically prescribing total budget for marketing and sales resources and allocation across spending categories
EP4092592A1 (en) Generating demand forecasts based on trend features derived from historical data
US6546303B1 (en) Computation of supply chain planning process efficiency
CN112613637A (en) Method and device for processing charging load
Shi Application of the model combining demand forecasting and inventory decision in feature based newsvendor problem
Fu et al. Bike-sharing inventory management for market expansion
CN101324939A (en) System and method for forecasting new business market based on data development
US8473331B2 (en) Computer-implemented systems and methods for determining future profitability
CN113469598A (en) Intelligent supply chain system and server platform
CN115860786A (en) Big data based sales management method and device
KR102638339B1 (en) System and method for managing bid based on intelligent demand response