知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023059706
(43)【公開日】2023-04-27
(54)【発明の名称】泌乳量算出システム
(51)【国際特許分類】
A01K 29/00 20060101AFI20230420BHJP
【ＦＩ】
A01K29/00 A
A01K29/00 C
【審査請求】未請求
【請求項の数】16
【出願形態】ＯＬ
(21)【出願番号】P 2021169871
(22)【出願日】2021-10-15
(71)【出願人】
【識別番号】314012076
【氏名又は名称】パナソニックＩＰマネジメント株式会社
(74)【代理人】
【識別番号】100118924
【弁理士】
【氏名又は名称】廣幸 正樹
(72)【発明者】
【氏名】藤原 誠二
(72)【発明者】
【氏名】蓑島 国彦
(72)【発明者】
【氏名】藤山 広光
(72)【発明者】
【氏名】柴田 泰匡
(57)【要約】
【課題】予測起算日以前のデータに基づいて、前記予測起算日以後の泌乳量を算出する泌乳量算出システム。
【解決手段】前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測日における前記分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部を含む牛群推定泌乳量算出部を有する制御装置を有する泌乳量算出システム。
【選択図】図１
【解決手段】前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測日における前記分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部を含む牛群推定泌乳量算出部を有する制御装置を有する泌乳量算出システム。
【選択図】図１
【特許請求の範囲】
【請求項1】
予測起算日以前の牛群を構成する個体牛ごとの分娩後日数と実乳量、および気候データを含む牧場データに基づいて、前記予測起算日以後の予測日における前記個体牛の泌乳量を牛群推定泌乳量として算出する泌乳量算出システムであって、
前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に前記牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測日における前記分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部を含む牛群推定泌乳量算出部を有する制御装置を有する泌乳量算出システム。
前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に前記牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測日における前記分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部を含む牛群推定泌乳量算出部を有する制御装置を有する泌乳量算出システム。
【請求項2】
前記牛群推定泌乳量と、前記予測日において想定される牛乳単価に基づいて売上金を算出する牛群予測売上金算出部をさらに有する請求項１に記載された泌乳量算出システム。
【請求項3】
前記牛群予測売上金算出部は、前記予測日の牛群推定泌乳量に前記牛乳単価を乗じて得る牛群推定個体売上金を算出する牛群推定個体売上金算出部を有することを特徴とする請求項２に記載された泌乳量算出システム。
【請求項4】
前記牛群予測売上金算出部は、前記予測日における前記牛群に属する個体牛の牛群推定泌乳量の総和である牛群推定総泌乳量に前記牛乳単価を乗じて得られる牛群推定総売上金を算出する牛群推定総売上金算出部を有することを特徴とする請求項２または３に記載された泌乳量算出システム。
【請求項5】
前記牛群予測売上金算出部は、前記予測日における前記牛群に属する個体牛の前記牛群推定泌乳量に前記牛乳単価を乗じて得られる牛群推定個体売上金の総和を牛群推定総売上金として算出する牛群推定総売上金算出部を有することを特徴とする請求項２または３に記載された泌乳量算出システム。
【請求項6】
前記牛群推定泌乳量と、乾乳基準値に基づいて、前記予測日において前記牛群に属する乳牛のうち、乾乳牛候補を示す牛群乾乳牛判定部をさらに有する請求項１に記載された泌乳量算出システム。
【請求項7】
前記牛群乾乳牛判定部は、前記予測日における前記牛群推定泌乳量が前記乾乳基準値より低い前記個体牛を前記牛群乾乳牛候補として表示することを特徴とする請求項６に記載された泌乳量算出システム。
【請求項8】
前記牛群乾乳牛判定部は、前記予測日における前記牛群推定泌乳量が前記乾乳基準値となる牛群乾乳日より前記分娩後日数が大きくなる前記個体牛を前記牛群乾乳牛候補として表示することを特徴とする請求項６に記載された泌乳量算出システム。
【請求項9】
予測起算日以前の牛群を構成する個体牛ごとの分娩後日数と実乳量、および気候データを含む牧場データに基づいて、前記予測起算日以後の予測日における前記個体牛の推定泌乳量を個体推定泌乳量として算出する泌乳量算出システムであって、
前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測起算日における前記個体牛の前記分娩後日数に対応する前記回帰式に、
前記予測起算日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部と、
前記予測起算日の前記牛群推定泌乳量と、前記予測起算日の特定個体牛の実乳量に基づいて、前記特定個体牛の前記回帰式に対する割合を乳量係数として求める乳量係数算出部と、
前記予測日における前記特定個体牛の分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、得られた結果にさらに前記乳量係数を乗じた値を個体推定泌乳量として算出する個体回帰式演算部を含む個体推定泌乳量算出部を有する制御装置を有する泌乳量算出システム。
前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測起算日における前記個体牛の前記分娩後日数に対応する前記回帰式に、
前記予測起算日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部と、
前記予測起算日の前記牛群推定泌乳量と、前記予測起算日の特定個体牛の実乳量に基づいて、前記特定個体牛の前記回帰式に対する割合を乳量係数として求める乳量係数算出部と、
前記予測日における前記特定個体牛の分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、得られた結果にさらに前記乳量係数を乗じた値を個体推定泌乳量として算出する個体回帰式演算部を含む個体推定泌乳量算出部を有する制御装置を有する泌乳量算出システム。
【請求項10】
前記個体推定泌乳量と、前記予測日において想定される牛乳単価に基づいて売上金を算出する個体予測売上金算出部をさらに有する請求項９に記載された泌乳量算出システム。
【請求項11】
前記個体予測売上金算出部は、前記個体推定泌乳量に、前記牛乳単価を乗じて前記特定個体牛の前記予測日における個体推定売上金を算出する個体推定売上金算出部をさらに有する請求項１０に記載された泌乳量算出システム。
【請求項12】
前記個体予測売上金算出部は、前記牛群に属する全ての前記個体牛の前記予定日における前記個体推定泌乳量の合計である個体推定総泌乳量に前記牛乳単価を乗じた額を個体推定売上金として算出する個体推定総売上金算出部をさらに有する請求項１０または１１に記載された泌乳量算出システム。
【請求項13】
前記個体予測売上金算出部は、前記牛群中の全ての前記個体牛について、前記個体推定泌乳量に前記牛乳単価を乗じて前記個体毎の売上金を算出し、それらの総和を、個体推定総売上金として算出する個体推定総売上金算出部をさらに有する請求項１０または１１に記載された泌乳量算出システム。
【請求項14】
前記個体推定泌乳量と、乾乳基準値に基づいて、前記予測日において前記牛群に属する乳牛のうち、乾乳牛候補を示す個体乾乳牛判定部をさらに有する請求項９に記載された泌乳量算出システム。
【請求項15】
前記個体乾乳牛判定部は、前記予測日における前記個体推定泌乳量が前記乾乳基準値より低くなる前記個体牛を前記乾乳牛候補として表示することを特徴とする請求項１４に記載された泌乳量算出システム。
【請求項16】
前記個体乾乳牛判定部は、前記予測日における前記個体推定泌乳量が前記乾乳基準値となる個体乾乳日より前記分娩後日数が大きくなる前記個体牛を個体乾乳牛候補として表示することを特徴とする請求項１４に記載された泌乳量算出システム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、分娩後日数に応じた乳牛の搾乳量を、過去のデータに基づいて推定する泌乳量算出システムに関するものである。
【背景技術】
【0002】
酪農業にとっては、泌乳量の増減は、牧場経営に直接影響する大きな問題である。より多くの乳量を得るために、飼料の配合を工夫したり、飼育環境を改善し、乳牛のストレスを少なくするといった工夫が行われている。
【0003】
一方、このような改善は、個々の牧場に依存する部分が多いため、他の牧場との比較は行いにくい。一般に乳牛は暑さに弱いため、緯度の高い地方の方が泌乳量は多いと考えられるが、気候の違いがどの程度泌乳量に影響しているのかは明らかになっているわけではない。
【0004】
まして、泌乳量に影響を及ぼす要因は非常に多く考えられるため、どのようなケースの時に、どのような要因が泌乳量に大きく影響を及ぼしているかについては、経験に頼っているのが現状である。
【0005】
したがって、考えられる要因を入力すると泌乳量を推定できるシステムは、酪農経営に大きな指針を与えるものである。
【0006】
特許文献１は、情報端末を用いて、乳牛を泌乳ステージ別若しくは繁殖ステージ別に群で把握するシステムを提供している。
【0007】
また、特許文献２では、牛個体ごとに、過去乳量から回帰式により泌乳曲線を算出するものが開示されている。
【先行技術文献】
【特許文献】
【0008】
【特許文献1】特開２０２０－１５６３５９号公報
【特許文献2】特開２０２０―０２０７０７号公報
【発明の概要】
【発明が解決しようとする課題】
【0009】
特許文献１は、乳牛の現在の状況を把握しようとするものであり、特許文献２は、過去の泌乳量から未来の泌乳量を推定しようとするものである。しかし、乳牛の泌乳量を増加させるための要因を調べることができるものではない。
【0010】
高温やストレスによるオキシトシンの減少とアドレナリンの分泌が泌乳量に影響があるという生理的な理屈は理解されているものの、どういった飼育環境や乳牛自身の状態が、乳牛にとって快適な状況であるのかは、人間には理解しにくい。そのため、多くの事象（情報）を集めて要因分析し、泌乳量を高めようとする方法は有用であると考えられる。
【0011】
ここで、多くの情報とは、できるだけ多くの牧場の、飼育環境や飼料配合、乳牛自身の個別情報が該当する。これらの情報は日々更新されている。また、情報の多様性を確保するには、異なる地域の情報もあれば好ましい。したがって、多くの牧場からのデータを利用できるようにするのが好ましい。
【0012】
しかしながら、規模や事情の異なる複数の牧場からのデータの定期的な収集は困難を極める。牧場によっては、毎日各種データが数値化される場合もあれば、数日毎でなければ、データがそろわない場合もあるからである。
【0013】
特に所謂泌乳曲線は分娩後日数と乳量の関係を表すものであるが、分娩後日数が１日違いでそろっている乳牛がいる牧場はよほど大規模でないかぎり、想定しにくい。したがって、実際に収集できる泌乳量のデータでさえ、抜けが生じたり、情報の範囲（分娩後日数の範囲）が異なり、画一的に分析できないという課題があった。
【0014】
また、酪農経営には、泌乳量の将来予測や、予測売上高が重要となる。しかし、分娩後日数が日々変わっていく複数の乳牛の将来的な泌乳量の推定は、従来あまり行われてこなかった。
【課題を解決するための手段】
【0015】
本発明に係る泌乳量算出システムは、上記の課題に鑑みて想到されたものであり、収集した泌乳量に係る情報が、同一範囲にそろわない場合であっても、要因分析を可能し、泌乳量の将来予測値を提示するシステムを提供するものである。
【0016】
より具体的に本発明に係る泌乳量算出システムは、
予測起算日以前の牛群を構成する個体牛ごとの分娩後日数と実乳量、および気候データを含む牧場データに基づいて、前記予測起算日以後の予測日における前記個体牛の泌乳量を牛群推定泌乳量として算出する泌乳量算出システムであって、
前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に前記牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測日における前記分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部を含む牛群推定泌乳量算出部を有する制御装置を有することを特徴とする。
予測起算日以前の牛群を構成する個体牛ごとの分娩後日数と実乳量、および気候データを含む牧場データに基づいて、前記予測起算日以後の予測日における前記個体牛の泌乳量を牛群推定泌乳量として算出する泌乳量算出システムであって、
前記牧場データを蓄積するメモリと、
前記予測起算日以前の前記牧場データ中の一定期間の前記分娩後日数と前記実乳量の関係を補間する補間式を作成する補間式作成部と、
前記補間式から算出される特定の前記分娩後日数における補間乳量を目的変数とし、前記牧場データから選ばれた少なくとも１つの要因パラメータを説明変数として前記分娩後日数毎に前記牛群推定泌乳量を算出する回帰式を作成する回帰式作成部と、
前記予測日における前記分娩後日数に対応する前記回帰式に、
前記予測日において想定される前記説明変数の値を代入し、前記牛群推定泌乳量を算出する牛群回帰式演算部を含む牛群推定泌乳量算出部を有する制御装置を有することを特徴とする。
【発明の効果】
【0017】
本発明に係る泌乳量算出システムは、収集した泌乳量に関するデータを、補間することで、一定の範囲の泌乳量を１日単位でさえ得ることができるので、泌乳量に影響を及ぼす要因の推定に役立つという効果を与える。
【0018】
また、気温といった予測可能な要因が泌乳量に及ぼす影響が判れば、将来の泌乳量の予測を行うことができ、酪農経営における重要な指針を得ることができるという効果も奏する。
【図面の簡単な説明】
【0019】
【図1】本発明に係る泌乳量算出システムの構成を示す図である。
【図2】泌乳量算出システムの全体（メイン）フローを示す図である。
【図3】泌乳量算出システムの算出する出力を例示する図である。
【図4】推定泌乳量を算出する工程のフローを示す図である。
【図5】補間式の作成を説明する図である。
【図6】補間式を作成するフローを示す図である。
【図7】回帰式の作成を説明する図である。
【図8】回帰式を作成する工程のフローを示す図である。
【図11】要因分析において、一定期間の乳量の実測値と回帰式から求めた推定泌乳量ＥＹを比較する様子を示す図である。
【図12】牛群回帰式演算のフローを示す図である。
【図13】牛群推定総泌乳量を説明する図である。
【図14】牛群総回帰式演算のフローを示す図である。
【図15】牛群推定総泌乳量を用いた牛群推定総売上金算出部を有する牛群予測売上金を算出するフローを示す図である。
【図16】個体牛毎の牛群推定泌乳量を用いた牛群推定総売上金算出部を有する牛群予測売上金を算出するフローを示す図である。
【図17】牛群乾乳日を説明する図である。
【図18】個体牛毎の牛群推定泌乳量と乾乳基準値を比較する牛群乾乳牛判定のフローを示す図である。
【図19】牛群乾乳日と個体牛毎の分娩後日数を比較する牛群乾乳牛判定のフローを示す図である。
【図20】個体推定泌乳量の算出手順を説明する図である。
【図21】個体推定総泌乳量を説明する図である。
【図22】個体推定泌乳量を算出するフローを示す図である。
【図23】個体回帰式演算のフローを示す図である。
【図24】個体総回帰式演算のフローを示す図である。
【図25】個体推定総泌乳量を用いた個体推定総売上金算出部を有する個体予測売上金を算出するフローを示す図である。
【図26】個体牛毎の個体推定泌乳量を用いた個体推定総売上金算出部を有する牛群予測売上金を算出するフローを示す図である。
【図27】個体乾乳日を説明する図である。
【図28】個体牛毎の個体推定泌乳量と乾乳基準値を比較する個体乾乳牛判定のフローを示す図である。
【図29】個体牛毎の個体乾乳日と個体牛毎の分娩後日数を比較する個体乾乳牛判定のフローを示す図である。
【図30】泌乳量算出システムを実際の牧場に使用し、推定泌乳量は、ＴＨＩと日照時間を要因パラメータとして重回帰できることを示す図である。
【発明を実施するための形態】
【0020】
以下に本発明に係る泌乳量算出システムについて図面を示し説明を行う。なお、以下の説明は、本発明の一実施形態および一実施例を例示するものであり、本発明が以下の説明に限定されるものではない。以下の説明は本発明の趣旨を逸脱しない範囲で改変することができる。
【0021】
図１に本発明に係る泌乳量算出システムの構成を示す。本発明に係る泌乳量算出システム１は、端末１０と本体１２およびメモリ１４で構成される。端末１０は、それ自身がＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｏｒ Ｕｎｉｔ）を有するコンピュータとメモリおよび表示画面で構成することができ、移動通信端末（所謂スマートフォン等）を含めることができる。
【0022】
本体１２はＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｏｒ Ｕｎｉｔ）で構成される。１台のＣＰＵであってもよいし、複数のＣＰＵを接続したものであってもよい。本体１２はサービスを提供する装置なので、サーバーと呼んでよい。メモリ１４は本体１２が使用するメモリである。ＣＰＵは制御装置と呼んでもよい。したがって、本体１２は制御装置である。また、本体１２には、端末１０からの信号を受けとる入力部１２ａがあるとしてもよい。すなわち、本体１２は入力部１２ａで指定されたコマンドや数値を受ける。
【0023】
本体１２と端末１０は、通信可能に接続されている。通信は、有線であっても、無線であってもよい。端末１０は、クライアントと呼ぶこともできる。端末１０は、実際に乳牛を飼育している牧場１６（ここでは牧場１６Ａ、１６Ｂ、１６Ｃを例示した。）に設置されるのが好ましい。ここで、「設置」とは、端末１０を牧場１６に物理的に設置する意味の他、牧場１６の責任者若しくは担当者が移動通信端末を保有することであってもよい。もちろん、端末１０は、牧場１６以外の場所に設置されてもよい。また、牧場１６は、乳牛が飼育されていればよい。
【0024】
本体１２には、複数の牧場１６が接続されていてもよい。泌乳量算出システム１は、泌乳量に係るデータを多量に集め、これらを分析することで、泌乳量を算出するので、多くのデータが集まるのが望ましいからである。また、本体１２は、外部情報１８とインターネットを介して接続されていてもよい。牧場１６から得られるデータだけでなく、全国的な天気に関する情報をネット経由で取得して利用することもできるからである。
【0025】
このように、泌乳量算出システム１は、多量のデータを様々な場所から収集し、これらを分析することで、泌乳量を算出するので、ネットワーク経由でサービスを提供する形態が好適な実施形態としてあり得る。すなわち、泌乳量算出システム１はクラウドで構成されてもよい。
【0026】
したがって、泌乳量算出システムを実施するクラウドの形態としては、代表的なＳａａＳ（Ｓｏｆｔｗａｒｅ ａｓ ａ Ｓｅｒｖｉｃｅ）が好適であるが、ＰａａＳ（Ｐｌａｔｆｏｒｍ ａｓ ａ Ｓｅｒｖｉｃｅ）、ＨａａＳ（Ｈａｒｄｗａｒｅ ａｓ ａ Ｓｅｒｖｉｃｅ）若しくはＩａａＳ（Ｉｎｆｒａｓｔｒｕｃｔｕｒｅ ａｓ ａ Ｓｅｒｖｉｃｅ）といった形態であってもよい。
【0027】
本発明に係る泌乳量算出システム１は、複数頭の乳牛の集合を１つの牛群とし、牛群での乳量の予測を行う。また、牛群に属する個々の乳牛を個体牛（以下単に「個体」とも呼ぶ。）とし、その個体牛の泌乳量予測も行う。さらに、これらの予測される泌乳量から、売上金の推定や、乾乳日および乾乳牛の推定も行うことができる。乳牛は出産によって乳の分泌が始まるが、分娩後日数が相当日経過すると、泌乳量が低下する。このような牛は次の出産までの間休ませる。乾乳日以降その乳牛は、次の出産までの間、泌乳量が期待できなくなるからである。これを「乾乳牛」と呼ぶ。また、「乾乳日」とは、乳牛を乾乳牛扱いとする日のことである。
【0028】
より具体的には、泌乳量算出システム１は、端末１０から予測起算日より将来の予測日に想定される要因パラメータＰ（説明変数）の値を入力すると、予測日若しくは予測期間における分娩後日数Ｎ日の乳牛から搾乳できる牛群推定泌乳量ＮＥＹ（ｋｇ）または個体推定泌乳量ＮＤＥＹ（ｋｇ）を得ることができる。また、これらは、分娩後日数Ｎを省略し、単に牛群推定泌乳量ＥＹ、個体推定泌乳量ＤＥＹとも呼ぶ。牛群推定泌乳量ＥＹまたは個体推定泌乳量ＤＥＹを総称する際には、「牛群推定泌乳量ＮＥＹ（ｋｇ）等」若しくは「牛群推定泌乳量ＥＹ等」と呼ぶ。
【0029】
なお、牛群推定泌乳量ＥＹは、牛群と冠しているものの１頭の乳牛の推定泌乳量である。牛群に属する個体牛の全ての泌乳量（後述する「牛群推定総泌乳量」）は、各個体牛毎の分娩後日数における牛群推定泌乳量ＥＹの総和を求める必要がある。
【0030】
ここで予測起算日とは、予測を行う日である。また、予測起算日は週や月であることを排除しない。通常予測起算日は予測を行う日であってよいが、過去の日や期間であってもよい。また予測日とは、予測起算日より将来の日であり、牛群推定泌乳量ＮＥＹ（ｋｇ）等を予測したい日である。また、予測期間とは、予測起算日より将来の期間であり、牛群推定泌乳量ＮＥＹ（ｋｇ）等を予測したい期間である。予測日と言った場合は予測期間を含めてもよい。
【0031】
また、要因パラメータＰとは、気候データ、乳牛自体の体重、飼料摂取量といった個別データ、摂取させている飼料の種類といった飼料データ等を含む牧場データの内、泌乳量を最もよく説明することのできる項目である。この要因パラメータＰは、直接測定可能な項目以外に、直接測定可能な項目のデータを加工したデータであってもよい。この要因パラメータＰは、泌乳量算出システム１内で、探すことができる。
【0032】
気候データは、日の出時刻、日の入り時刻、気温、湿度、日照時間、日射量、風向、風速等である。また、牛舎内温度、牛舎内湿度を含んでもよい。また、ＴＨＩ（Ｔｅｍｐｅｒａｔｕｒｅ Ｈｕｍｉｄｉｔｙ Ｉｎｄｅｘ：温度湿度指数）等、これらの加工データを含んでもよい。なお、ＴＨＩとは一例として（１）式で表される指数である。
【0033】
【数1】
【0034】
また、気候データは、牛舎内、牛舎外の気候だけでなく、牧場の１ｋｍ四方の気候、地域の気候また全国の天気予報サービスからの気候データであってもよい。
【0035】
個別データは、乳牛自体のデータである。具体的には、登録番号、血統、産次数、分娩後日数、乳量、体重、病歴、所有者、飼育場所等である。
【0036】
飼料データは、過去に摂取させた飼料の種類や配合、回数などが含まれる。
【0037】
［基本データ］
このように要因パラメータＰは、予測起算日以前のデータを含む概念である。推定のために用いる予測起算日以前のデータを「基本データ」と呼ぶ。「基本データ」は予測のための回帰式を作成するためのデータである。また牧場データの一部といってもよい。本発明に係る泌乳量算出システム１は、予測起算日以前の要因パラメータＰに基づいて一定の乳牛群若しくは個体牛の乳牛の泌乳量を目的変数として説明できる回帰式を求め、その回帰式に予測起算日以降の予測日に想定される要因パラメータＰ（説明変数）を当てはめることで、予測起算日以降の予測日若しくは予測期間における分娩後日数Ｎ日の乳牛から搾乳できる牛群推定泌乳量ＮＥＹ（ｋｇ）等を得ることができるものと言い換えることもできる。
このように要因パラメータＰは、予測起算日以前のデータを含む概念である。推定のために用いる予測起算日以前のデータを「基本データ」と呼ぶ。「基本データ」は予測のための回帰式を作成するためのデータである。また牧場データの一部といってもよい。本発明に係る泌乳量算出システム１は、予測起算日以前の要因パラメータＰに基づいて一定の乳牛群若しくは個体牛の乳牛の泌乳量を目的変数として説明できる回帰式を求め、その回帰式に予測起算日以降の予測日に想定される要因パラメータＰ（説明変数）を当てはめることで、予測起算日以降の予測日若しくは予測期間における分娩後日数Ｎ日の乳牛から搾乳できる牛群推定泌乳量ＮＥＹ（ｋｇ）等を得ることができるものと言い換えることもできる。
【0038】
牧場１６に設置された端末１０は、牧場の乳牛に係る個別データ、実際に摂取させた飼料データ、日々の気候データを牧場データＦＤとして本体１２に送信する。本体１２はこれらの牧場データをメモリ１４に蓄積しておく。そして、予測起算日が端末１０から入力されると、予測起算日より以前の要因パラメータＰに基づいて牛群推定泌乳量等を算出する式を求める。続いて入力される予測日若しくは予測期間において想定される要因パラメータＰに従って、分娩後日数Ｎに応じた牛群推定泌乳量ＮＥＹ等を算出し、端末１０に表示させる。
【0039】
なお、分娩後日数Ｎに応じた牛群推定泌乳量ＮＥＹ等の算出の際に入力される予測起算日以降の予測日において想定される要因パラメータＰ（説明変数）は、上記の回帰式に入力される要因パラメータＰの値であり、回帰式の説明変数の値である。この要因パラメータＰ（説明変数）の値は、天気予報等で公表される値だけでなく、例年若しくは過去の該当する日若しくは期間の値および、それらに基づいて算出した何らかの値であってもよい。
【0040】
つまり、説明変数となる要因パラメータＰは過去の値を使い得る。毎年の気候は四季を通じてそれほど大きくは変化しないと考えられるので、予測起算日以降の予測日や予測期間は、昨年若しくはそれ以前の予測日に相当する日や期間と大きくは変わらないからである。
【0041】
なお、本体１２は、気候データを外部情報１８から取得してもよい。クローズタイプの牛舎であれば、牛舎内の環境は管理されており、データも送信しやすい。しかし、オープンタイプの牛舎の場合、その日の牛舎内の風速は取得しにくい。そのような場合は、外部の気候データを参考にできるからである。
【0042】
気候に関するデータは、さまざまな物が提供されており、それらも有効に使用することができる。また、牧場データＦＤは端末１０以外の方法で本体１２に送信されてもよい。もちろん、牧場データＦＤを本体１２に送信する端末１０と、要因パラメータＰを入力し、泌乳量を算出させる端末１０が別々であってもよい。
【0043】
図２に泌乳量算出システム１の全体（メイン）フローを示す。なお、牧場データＦＤを本体１２に送信するのは、各牧場が適宜行うとする。また、扱いの単位となる牛群に属する乳牛は予め登録されているとする。
【0044】
泌乳量算出システム１は、計算条件や要因パラメータＰ（説明変数を含む）といった端末１０等から入力された入力に対して、牛群推定泌乳量ＥＹを算出する工程、牛群予測売上金Ｅ＄を算出する工程、牛群乾乳牛ＤＣを判定する工程、個体推定泌乳量ＤＥＹを算出する工程、個体予測売上金ＤＥ＄を算出する工程、個体乾乳牛ＤＤＣを判定する工程による結果を返すことができる。
【0045】
また、牛群推定総泌乳量ＳＵＭＥＹや個体推定総泌乳量ＳＵＭＤＥＹを算出する工程があってもよい。また、牛群予測売上金Ｅ＄を算出する工程や個体予測売上金ＤＥ＄を算出する工程には、それぞれ牛群推定個体売上金ＥＩ＄を算出する工程、牛群推定総売上金ΣＥ＄を算出する工程や、個体推定売上金ＤＥＩ＄を算出する工程、個体推定総売上金ΣＤＥ＄を算出する工程などを含むことができる。
【0046】
これらの結果は端末１０に返してよい。また、これらの結果は、ネット上の特定のメモリ上に置いてもよい。端末１０はこれらの結果を表示することができる。これらを算出する工程は、メインフローのメニューと呼ぶことができる。すなわち、端末１０からはこれらの算出を指示することができる。
【0047】
なお、本発明は構成の大部分がソフト的に実施される。そこで、処理の「工程」とは、処理のまとまりをいい、本体１２には、「工程」を実行する「部」が存在すると考えてよい。具体的には、本体１２（制御装置）には、牛群推定泌乳量算出部、牛群予測売上金算出部、牛群乾乳牛判定部、個体推定泌乳量算出部、個体予測売上金算出部、個体乾乳牛判定部が存在すると言える。また、牛群推定総泌乳量算出部、個体推定総泌乳量算出部、牛群推定個体売上金算出部、牛群推定総売上金算出部、個体推定売上金算出部、個体推定総売上金算出部が存在すると言える。
【0048】
また、これらの処理工程の中に下位の処理のまとまりがある場合は、それも「部」と呼んでよい。例えば、牛群推定泌乳量算出部には、補間式を作成する工程や回帰式を作成する工程が含まれ、それぞれ「補間式作成部」、「回帰式作成部」と呼ぶことができる。また、得られた回帰式に要因パラメータＰ（説明変数）の値が代入された際に、牛群推定泌乳量を具体的に算出する工程（ステップＳ２１０）があり、牛群回帰式演算部と呼ぶことができる。なお、泌乳量算出システム１は少なくとも牛群推定泌乳量ＥＹを算出する工程（牛群推定泌乳量算出部）を有していればよい。
【0049】
＜泌乳量算出システム＞
図２を参照して、泌乳量算出システム１がスタートすると（ステップＳ１００）、終了判断を行う（ステップＳ１０２）。終了判断は、端末１０からの終了指示や、本体１２と端末１０の通信の切断であってよい。終了する場合（ステップＳ１０２のＹ分岐）は、終了する（ステップＳ１０４）。終了しない場合（ステップＳ１０２のＮ分岐）は、次に処理を進める。
図２を参照して、泌乳量算出システム１がスタートすると（ステップＳ１００）、終了判断を行う（ステップＳ１０２）。終了判断は、端末１０からの終了指示や、本体１２と端末１０の通信の切断であってよい。終了する場合（ステップＳ１０２のＹ分岐）は、終了する（ステップＳ１０４）。終了しない場合（ステップＳ１０２のＮ分岐）は、次に処理を進める。
【0050】
次に牛群推定泌乳量ＥＹを算出するか、個体推定泌乳量ＤＥＹを算出するかを選択する（ステップＳ１０６、ステップＳ１０８）。それぞれの選択肢では、選択する場合（Ｙ分岐）は、それぞれの処理に移り、選択しない場合（Ｎ分岐）は処理を次に移す。ステップＳ１０８の個体推定泌乳量ＤＥＹの算出を行わない場合（ステップＳ１０８のＮ分岐）は、終了判断（ステップＳ１０２）に戻る。
【0051】
牛群推定泌乳量ＥＹを選択した場合（ステップＳ１０６のＹ分岐）は、牛群推定泌乳量ＥＹが算出される（ステップＳ１１０）。そして、その後牛群予測売上金Ｅ＄の算出を行うか（ステップＳ１１２）、牛群乾乳牛ＤＣを判定するか（ステップＳ１１４）が選択される。それぞれの処理が選択された場合（ステップＳ１１２およびステップＳ１１４でのＹ分岐）は、それぞれ牛群予測売上金Ｅ＄が算出され、牛群乾乳牛ＤＣが判定される（ステップＳ１１６およびステップＳ１１８）。その後処理は、個体推定泌乳量ＤＥＹの選択（ステップＳ１０８）に処理が移される。
【0052】
個体推定泌乳量ＤＥＹを選択した場合（ステップＳ１０８のＹ分岐）は、個体推定泌乳量ＤＥＹが算出される（ステップＳ１２０）。そして、その後個体予測売上金ＤＥ＄の算出を行うか（ステップＳ１２２）、個体乾乳牛ＤＤＣを判定するか（ステップＳ１２４）が選択される。それぞれの処理が選択された場合（ステップＳ１２２およびステップＳ１２４でのＹ分岐）は、それぞれ個体予測売上金ＤＥ＄が算出され、個体乾乳牛ＤＤＣが判定される（ステップＳ１２６およびステップＳ１２８）。その後処理は、終了判定（ステップＳ１０２）に処理が移される。
【0053】
図３には、それぞれの出力例を示す。図３（ａ）は牛群推定泌乳量ＥＹの出力例である。横軸は分娩後日数（日）であり、縦軸は乳量（ここでは単位を「ｋｇ」とする）である。特定の要因パラメータＰ（例えば温度）の下で、分娩後日数に対する牛群推定泌乳量ＥＹを示すことができる。図では温度がＴ１℃の時と、Ｔ２℃の時の牛群推定泌乳曲線Ｍが示されている。なお、本発明の泌乳量算出システム１が提供する牛群推定泌乳曲線Ｍは、最小で１日単位の算出値（牛群推定泌乳量ＥＹ）の集合である。これらの算出値は直線もしくは曲線で結ばれていてもよい。牛群推定泌乳量ＥＹで形成した牛群推定泌乳曲線を符号Ｍで表す。言い換えると、牛群推定泌乳曲線Ｍは牛群推定泌乳量ＥＹを適当な分娩後日数Ｎ毎にプロットし、直線若しくは曲線で繋げたものである。
【0054】
図３（ｂ）は、牛群予測売上金Ｅ＄の出力例である。横軸は予測日であり、縦軸は売上金（ここでは単位を「円」とする。）である。牛群推定泌乳量ＥＹが求められるので、牛乳の単価（円／ｋｇ）を乗ずることで、予測日における牛群予測売上金Ｅ＄を算出することができる。なお、牛群予測売上金Ｅ＄を算出する工程には、個体牛毎の売上金である牛群推定個体売上金ＥＩ＄と牛群を構成する全乳牛による売上金である牛群推定総売上金ΣＥ＄を求める様にしてもよい。
【0055】
図３（ｃ）は、牛群乾乳牛ＤＣの判定結果を示すものである。横軸は分娩後日数（日）であり、縦軸は乳量（ここでは単位を「ｋｇ」とする）である。ある予測日ηにおける牛群推定泌乳量ＥＹ（牛群推定泌乳曲線Ｍ）を算出し、乾乳牛と判定する乾乳基準値ＤＹＴＨ以下となるのは、分娩後日数が何日目以降となるかを求める。この日を牛群乾乳日ＤＹｄａｙと呼ぶ。牛群乾乳日ＤＹｄａｙは、減少する牛群推定泌乳量が乾乳基準値となる日といってもよい。
【0056】
したがって、分娩後日数が牛群乾乳日ＤＹｄａｙより大きい個体牛は牛群乾乳牛ＤＣと判定される。図３（ｃ）で、符号ＤＣ（黒三角）は、ＤＹｄａｙより分娩後日数が大きいので牛群乾乳牛ＤＣである。一方、ＮＯＴＤＣ（黒丸）は、ＤＹｄａｙより分娩後日数が小さいので牛群乾乳牛ＤＣではない。
【0057】
図３（ｄ）は、個体推定泌乳量ＤＥＹの出力例である。横軸は分娩後日数（日）であり、縦軸は乳量（ここでは単位を「ｋｇ」とする）である。個体推定泌乳量ＤＥＹは、牛群推定泌乳量ＥＹに対する割合として求められる。個体推定泌乳量ＤＥＹを表した曲線を個体推定泌乳曲線ＤＭとする。図３（ｄ）では、牛群推定泌乳量ＥＹを表す牛群推定泌乳曲線Ｍより下方に個体推定泌乳量ＤＥＹを表す曲線ＤＭがある。これは、この特定個体牛は牛群として推定された泌乳量よりも少ない乳を分泌する乳牛であることを表している。したがって、他の個体牛では、曲線ＤＭが曲線Ｍより上に描かれる場合もある。
【0058】
図３（ｅ）は、個体予測売上金ＤＥ＄を表すグラフである。横軸は予測日であり、縦軸は売上金（例えば円である。）ここで例示した個体牛は、図３（ｄ）で示したように牛群推定泌乳量ＥＹより産出量が少ないので、個体予測売上金ＤＥ＄も牛群予測売上金Ｅ＄よりも少ない。個体予測売上金ＤＥ＄も、個体牛毎の売上金である個体推定売上金ＤＥＩ＄と個体推定総売上金ΣＤＥ＄を算出するようにしてもよい。
【0059】
図３（ｆ）は、個体乾乳牛ＤＤＣを表すグラフである。横軸は分娩後日数（日）であり、縦軸は乳量（ここでは単位を「ｋｇ」とする）である。例示した個体牛は、牛群として求めた推定泌乳量ＥＹより乳量が少ないので、個体乾乳日ＤＤＹｄａｙも牛群乾乳日ＤＹｄａｙよりも早まることとなる。図３（ｆ）では、予測日ηにおいて、この個体牛の分娩後日数が個体乾乳日ＤＤＹｄａｙ以上であれば個体乾乳牛ＤＤＣと判定され、個体乾乳日ＤＤＹｄａｙより小さければ個体乾乳牛ＤＤＣと判定されない（ＮＯＴＤＤＣ）。産次数が少ないうちは乳量が少ない個体牛も牛群として扱うことができるが、産次数が多くなり、除籍を検討する場合は、個体乾乳牛ＤＤＣの判定は有用である。
【0060】
図２を再度参照する。以上の全体フローは、本体１２と端末１０の間で分担して行ってよい。例えば、各処理の選択および最終結果の表示は端末１０で行い、各処理は本体１２で行うなどである。以下、各工程を詳説する。
【0061】
＜牛群推定泌乳量ＥＹ＞
図４に牛群推定泌乳量ＥＹの算出工程のフローを示す。牛群推定泌乳量ＥＹの算出工程が始まったら（ステップＳ１１０）、母集団データの選択と作成条件を決める（ステップＳ２００）。図１で示したメモリ１４中には、多数の牧場データが蓄積されている。その中で、必要な場所（牛群を構成する個体牛の特定を含む。）および基本データの期間を選択する。また、牛群推定泌乳量ＥＹを算出する条件もここで入力する。
図４に牛群推定泌乳量ＥＹの算出工程のフローを示す。牛群推定泌乳量ＥＹの算出工程が始まったら（ステップＳ１１０）、母集団データの選択と作成条件を決める（ステップＳ２００）。図１で示したメモリ１４中には、多数の牧場データが蓄積されている。その中で、必要な場所（牛群を構成する個体牛の特定を含む。）および基本データの期間を選択する。また、牛群推定泌乳量ＥＹを算出する条件もここで入力する。
【0062】
例えば、場所は、自分の牧場で、牛群は全ての飼育乳牛とし、基本データ期間は昨年度１年等である。また、条件とは、牛群推定泌乳量ＥＹを算出するための条件すべてを含んでよい。例えば、牧場データを参照する際にまとめる期間の長さを決める等である。より具体的にいうと、実乳量は１日毎の乳量にするか、３日毎、１週間毎、１か月毎の平均乳量であるとか、１か月の最大乳量値である等である。また、これに伴い、その他の変数も同様に一定期間毎に代表値を決めてよい。
【0063】
また、予測起算日と予測日を入力する。予測起算日は通常「予測操作を行おうとしている本日」であるが、本日より過去の日であってもよい。また、予測日とは、実質的には予測日に想定される要因パラメータＰである。すでに述べたように、過去の予測日に相当する日（若しくは期間）の実績値を予測日に想定される要因パラメータＰとして用いてもよい。したがって、要因パラメータＰを気候データ由来のものにしておけば、長期にわたる予測をすることができる。
【0064】
次に補間式を作成する（ステップＳ２０２）。図５に補間式について説明する。図５は横軸が分娩後日数Ｎ（日）であり、縦軸は実乳量ＲＹ（ｋｇ）である。例えば、ある牧場の複数の乳牛の、ある期間（例えば２０２０年８月１日）のデータである。ここでプロットされた各点は、個体牛毎の実乳量ＲＹなので、散布図となる。しかし、実乳量であるために、プロットできない範囲も生じる。該当する分娩後日数の乳牛がいない場合である。図５では、無データ領域ＶＲで表した。
【0065】
このように、牧場データに抜けがあると後の回帰式を作成する際のデータに不足が生じ、回帰の精度が低下する。そこで、この散布図を適当な関数で近似する。好適に利用できるのはＷＯＯＤ曲線である。ＷＯＯＤ曲線は、（２）式で表される曲線で、分娩後日数Ｎに対する泌乳量Ｙを表す曲線としてよく知られている。なお、実際の散布図を補間するには、ＷＯＯＤ曲線に限定される必要はなく他の関数であってもよい。
【0066】
【数2】
【0067】
なお、ここでＡ、Ｂ、Ｃは定数であり、Ｙは乳量、Ｎは分娩後日数、ｅはネイピア数である。このＷＯＯＤ曲線を最小二乗法で散布図にフィットするように定数Ａ、Ｂ、Ｃを決めることができる。このように、分娩後日数Ｎと乳量Ｙの関係を連続関数でフィッテングさせることを「補間式を作成する」という。散布図を補間した式をより一般化して（３）式と表す。
【0068】
【数3】
【0069】
つまり（３）式の補間式は、乳量Ｙは分娩後日数Ｎの関数として表される。関数の形としては、ＷＯＯＤ曲線は好適に利用できるが、これに限定されるものではない。図４の補間式作成工程（ステップＳ２０２）は、上記の様に実乳量ＲＹの散布図から補間式（３）を決める工程である。
【0070】
［補間式作成］
図６に図４の補間式を作成する工程（ステップＳ２０２）の処理を示す。補間式を作成する工程が始まると（ステップＳ２０２）、図４（ＥＹ算出フロー）のステップＳ２００で決めた母集団と条件より分娩後日数Ｎと実乳量ＲＹのデータを取り出す（ステップＳ２３０）。これらは実データと呼ぶことができる。具体的に散布図を描いてもよい。実データとは、各乳牛毎の分娩後日数Ｎと実乳量ＲＹの組である。
図６に図４の補間式を作成する工程（ステップＳ２０２）の処理を示す。補間式を作成する工程が始まると（ステップＳ２０２）、図４（ＥＹ算出フロー）のステップＳ２００で決めた母集団と条件より分娩後日数Ｎと実乳量ＲＹのデータを取り出す（ステップＳ２３０）。これらは実データと呼ぶことができる。具体的に散布図を描いてもよい。実データとは、各乳牛毎の分娩後日数Ｎと実乳量ＲＹの組である。
【0071】
次に分娩後日数Ｎと乳量Ｙの実データを最もよく反映する補間式を求める。例えば、上記に示した、ＷＯＯＤ曲線を使った最小二乗法を適用して、定数Ａ、Ｂ、Ｃを決定する（ステップＳ２３２）（図５も参照）工程である。これは、実乳量ＲＹにフィッテングさせる関数を求めると言える。そして、（３）式の補間式を得る（ステップＳ２３４）。その後ＥＹ算出ルーチン（図４）に戻る（ステップＳ２３６）。
【0072】
補間式は、毎日分作成することができる。搾乳はほぼ毎日行われるからである。しかし、環境や各個体牛について、大きな変化がない場合は、３日毎や、１週間毎の実データの平均を実データとみなして補間式（３）を作成してもよい。
【0073】
また、時間軸を大きくとらえ、１か月の実データの平均値を実データとみなして補間式を作成してもよい。より具体的には、ある個体牛の１か月の乳量の平均をその乳牛の実乳量ＲＹとする。またその乳牛の分娩後日数はその１か月の平均（つまり、月初めにＮ＝１０日であれば、その月の分娩後日数を２５日とする）とする等である。これらは、図４の「母集団データ、作成条件」（ステップS２００）で決められていても良い。
【0074】
引き続き図４を参照する。補間式を作成したら要因分析を行うか否か判断する（ステップＳ２０４）。この判断は泌乳量算出システム１の利用者が端末１０から指示することができる。要因分析とは、牧場データの中から（２）式の補間式から得られる乳量Ｙを最もよく説明できる要因パラメータＰを調べるか否かの判断である。例えば、泌乳量算出システム１を初めて使用する場合や、母集団を大きく変える場合等に利用する。
【0075】
要因分析を行う場合（ステップＳ２０４のＹ分岐）は、要因分析の工程を行う（ステップＳ２０６）。要因分析工程の詳細は後述する。行わない場合（ステップＳ２０４のＮ分岐）は、処理を次に移す。要因分析については、図１０で詳説する。
【0076】
次に回帰式を作成する（ステップＳ２０８）。図７に回帰式の作成について説明する。図５、図６で示した補間式（３）（図７（ａ））によって、例えば各月毎の分娩後日数Ｎの乳牛の補間乳量が算出できる。ここでは補間式を一般形Ｆｗ（）で表し、乳量Ｙは月ごとの平均泌乳量とする。したがって、Ｙ１月は１月の補間乳量を表し、分娩後日数Ｎの乳牛の補間式（３）を求める式をＦｗ１月（Ｎ）と表す。
【0077】
ここで、分娩後日数が５０日後をＮ５０と表す。１月から１２月までの各月における分娩後日数が５０日後の乳牛の補間乳量は、Ｆｗ１月（Ｎ５０）、Ｆｗ２月（Ｎ５０）、・・・、Ｆｗ１２月（Ｎ５０）で算出される。もちろんＦｗｍ（Ｎ）は、各月毎に補間式が作成されている（「ｍ」は「月」を表す。）。
【0078】
この各月の補間乳量を最もよく説明することのできる因子を要因パラメータＰとして決める。どのような要因パラメータＰがよいのかは、ステップＳ２０６の要因分析の工程を行って知ることになる。たとえば、温度、湿度、日照時間等である。すでに牛群推定泌乳量ＥＹを求めるために使用する要因パラメータＰが決まっている場合は、それを使用する。要因パラメータＰは１つでなくてもよい。すなわち、これは各月ごとの補間乳量Ｙを１以上の因子によって最小二乗法等により回帰式を求めることである。
【0079】
よく知られているように回帰式は（４）式のように表される。
【0080】
【数4】
【0081】
ここで、ＥＹは牛群推定泌乳量であり、ｘ１、ｘ２、・・・、ｘｋは因子（要因パラメータＰ）であり、ａ１、ａ２、・・・、ａｋ、ｃは定数である。
【0082】
図７（ｃ）には、因子が１つの場合の回帰式の例を示す。図７（ａ）で示す補間式（３）から、分娩後日数５０日の乳牛の各月毎の補間乳量Ｙｍを算出する（図７（ｂ）参照。）。ここで、ｍは月を表す。各月毎の補間乳量とは、月初めの３日間の平均とするなどであってもよい。特定の１日だけのデータにその月の補間乳量を代表させてもよい。
【0083】
次に、これら月ごとの補間乳量を要因パラメータＰ（ここでは、例えば正午の温度）で整理すると図７（ｃ）のグラフを得ることができる。なおここでは、各月の分娩後５０日の乳牛の牛群推定泌乳量ＥＹは、正午の気温でよく説明できるとした場合の結果を例示する。回帰式は、（５）式で表される。
【0084】
【数5】
【0085】
ここで、ｘ１は、正午の温度であり、ＥＹは牛群推定泌乳量である。なお、このような回帰式は分娩後日数分だけ作成される。つまり、分娩後日数Ｎの最終日を３００日とすると、１日から３００日後までの回帰式を得ることができる（図７（ｅ）参照。）。
【0086】
因子が複数ある場合は、図７（ｃ）のように２次元では記載できなくなるが、複数の要因パラメータＰによる重回帰式を得れば、分娩後日数Ｎと補間乳量Ｙの関係をよく近似することができる場合がある。このようにして求めた回帰式を（６）式と表す（図７（ｄ）参照。）。
【0087】
【数6】
【0088】
ここで、ＮＥＹは、分娩後日数Ｎの乳牛が要因パラメータＰの時の牛群推定泌乳量ＥＹと呼ぶ。また、（６）式に要因パラメータＰを代入し結果を求めることを「牛群推定泌乳量を算出する」と言ってもよい。算出した推定泌乳量は牛群推定泌乳量ＥＹである。なお、要因パラメータＰの時の分娩後日数Ｎの牛群推定泌乳量をＮＥＹ（Ｐ）と表す（Ｎは省略する場合もある。）。また、牛群に属するｉ番目の個体牛ＣＷｉの牛群推定泌乳量を表す場合は、ＮＥＹ[ＣＷｉ]等と[]（カギカッコ）を用いて表す。
【0089】
［回帰式作成］
図８に回帰式を作成する工程のフローを示す。図４で示したＥＹ算出フローで回帰式を作成する工程（ステップＳ２０８）が開始されると、図８に飛び、まず要因パラメータＰが入力される（ステップＳ２５０）。要因パラメータＰは、端末１０から入力することができる。つまり、泌乳量算出システム１の利用者が入力する。ここでは仮に要因パラメータＰがｐ１（温度）であったとする。
図８に回帰式を作成する工程のフローを示す。図４で示したＥＹ算出フローで回帰式を作成する工程（ステップＳ２０８）が開始されると、図８に飛び、まず要因パラメータＰが入力される（ステップＳ２５０）。要因パラメータＰは、端末１０から入力することができる。つまり、泌乳量算出システム１の利用者が入力する。ここでは仮に要因パラメータＰがｐ１（温度）であったとする。
【0090】
次に要因パラメータｐ１について指定された一定期間ｍ（これはステップＳ２００で定められる）毎に、分娩後日数Ｎ毎の推定泌乳量Ｙｍを、補間式（３）によって求める（ステップＳ２５２）。つまり、図７（ｂ）の補間乳量Ｙｍを求める。これによって、図７（ｃ）の散布図を描くことができる。
【0091】
次に要因パラメータｐ１を説明変数、補間乳量Ｙｍを目標変数として、回帰式を求める（ステップＳ２５４）。すなわち、この時要因パラメータｐ１について、分娩後日数Ｎが同じ場合の乳量の関係（６）式（図７（ｄ））を得ることができる。回帰式の求め方は公知の方法でよい。
【0092】
このように本発明では、要因パラメータＰを使って牛群推定泌乳量を算出する際に、実乳量ＲＹを直接使用するのではなく、実乳量ＲＹを補間する補間式から算出した補間乳量を利用する。そのため、実乳量ＲＹにデータの抜けがあっても、妥当な値を平均的な乳量として算出することができる。
【0093】
（６）式は、分娩後日数Ｎが１から最終日（ここで最終日は、分娩後日数の最長日の意味）までの数を作ることができる（図７（ｅ）参照。）。回帰式（６）は、図７（ｅ）をまとめて表している。回帰式が求まったら、図４の牛群推定泌乳量ＥＹを求めるルーチンに戻る（ステップＳ２５６）。
【0094】
図４を再び参照し、牛群推定泌乳量ＥＹを求める回帰式（６）が求められれば、回帰式（６）に要因パラメータＰを入力すれば牛群推定泌乳量ＥＹを求めることができる（ステップＳ２１０）。
【0095】
［牛群回帰式演算ＥＹＩ］
図１２には回帰式（６）を実際に演算し、牛群推定泌乳量ＥＹを求めるフローを示す。これは牛群回帰式演算（ＥＹＩ）工程である。牛群回帰式演算が開始されたら（ステップＳ２１０）、分娩後日数Ｎと要因パラメータＰの入力を行う（ステップＳ２８０）。次に入力された分娩後日数Ｎと要因パラメータＰと（６）式を用いて要因パラメータＰの時の分娩後日数Ｎにおける牛群推定泌乳量ＮＥＹを算出する（ステップＳ２８２）。算出後表示してもよい。
図１２には回帰式（６）を実際に演算し、牛群推定泌乳量ＥＹを求めるフローを示す。これは牛群回帰式演算（ＥＹＩ）工程である。牛群回帰式演算が開始されたら（ステップＳ２１０）、分娩後日数Ｎと要因パラメータＰの入力を行う（ステップＳ２８０）。次に入力された分娩後日数Ｎと要因パラメータＰと（６）式を用いて要因パラメータＰの時の分娩後日数Ｎにおける牛群推定泌乳量ＮＥＹを算出する（ステップＳ２８２）。算出後表示してもよい。
【0096】
次に、再入力を確認する（ステップＳ２８４）。再入力をして他の分娩後日数Ｎや要因パラメータＰの時の牛群推定泌乳量ＮＥＹを求めたい場合（ステップＳ２８４のＹ分岐）は、入力（ステップＳ２８０）に戻る。そうでない場合は、ＥＹルーチンに戻る（ステップＳ２８６）。
【0097】
分娩後日数Ｎに対応する回帰式（６）は要因パラメータＰ（説明変数）の値を入力すると、牛群推定泌乳量ＥＹを算出する。したがって、分娩後日数Ｎ毎を表す回帰式（６）に順次同一の要因パラメータＰの値を入力することで、同一の要因パラメータＰの値における分娩後日数Ｎと牛群推定泌乳量ＥＹのデータセットを得ることができる。図３（ａ）の牛群推定泌乳量に係る曲線は、このようにして得られた分娩後日数Ｎと牛群推定泌乳量ＥＹのデータセット（データの対）をプロットし、若しくはこの点同士を直線でつないだものである。
【0098】
図９に図３（ａ）の牛群推定泌乳量ＥＹをプロットした牛群推定泌乳曲線Ｍを拡大した図を示す。分娩後日数毎に牛群推定泌乳量ＥＹをプロットしたものである。要因パラメータＰは、ｐ１とｐ２の２つの要因で重回帰式を求めた場合を示している。したがって、説明変数もｐ１とｐ２の２つの要因パラメータＰである。分娩後日数Ｎは、飛び飛びの値（１０日、２０日、３０日、５０日、８０日、１００日、１５０日、２００日、２５０日）で求めたが、１日単位で求めることもできる。各牛群推定泌乳量ＥＹは、（６）式が表す別々の回帰式（分娩後日数違い。図７（ｅ）参照。）から求められたものである。すなわち、これらの式の定数（例えば（５）式におけるａ１やｃ等）は異なる。
【0099】
牛群回帰式を演算するステップＳ２１０（図４参照）は、クライアント側が端末１０から予測日に想定される要因パラメータＰ（説明変数）の値と、所望する分娩後日数Ｎを入力し、本体１２が牛群推定泌乳量ＥＹを返す。クライアント側は、少なくとも１つの牛群推定泌乳量ＥＹを返せば本発明に係る泌乳量算出システム１としての要件を満たす。複数の分娩後日数Ｎに対して牛群推定泌乳量ＥＹを求め、図９のようにグラフにしてもよい。
【0100】
次に牛群総回帰式演算ＥＹＳを実行し牛群推定総泌乳量ＳＵＭＥＹを求める（ステップＳ２１２）。図１３に牛群推定総泌乳量ＳＵＭＥＹの算出の概念図を示す。横軸は分娩後日数（日）であり、縦軸は牛群推定泌乳量（ｋｇ）である。牛群推定泌乳量ＥＹが求まるので、牛群推定泌乳曲線Ｍを描くことができる。具体的には予測日に想定される要因パラメータＰを（６）式の全ての回帰式に代入することで、全ての分娩後日数Ｎに対する牛群推定泌乳量ＥＹをプロットすることができる。ここで、牛群は９頭の乳牛から構成されているとする。
【0101】
そして、予測日の要因パラメータＰ（説明変数）がｐ１であるとし、予測日における各乳牛の分娩後日数がＮ１、Ｎ２、Ｎ３、Ｎ４、Ｎ５日であるとする。また、各分娩後日数における乳牛の頭数は、それぞれ１頭、２頭、４頭、１頭、１頭であるとする。図１３では、横軸の分娩後日数に括弧をつけてその中に頭数を表示した。
【0102】
すると牛群推定総泌乳量ＳＵＭＥＹは、各乳牛の各分娩後日数Ｎから求められる牛群推定泌乳量ＥＹの総和として求められる。
【0103】
すなわち、一般的に牛群推定総泌乳量ＳＵＭＥＹは、各分娩後日数ｋ日の乳牛がそれぞれＣｋ頭いるとすると（７）式のように表される。
【0104】
【数7】
【0105】
また、各乳牛の分娩後日数Ｎから牛群推定泌乳量ＮＥＹを求め、その牛群推定泌乳量ＮＥＹをその乳牛の乳量とし、全乳牛について加算することで牛群推定総泌乳量ＳＵＭＥＹとしてもよい。（７）式と同じ結果である。
【0106】
なお、牛群推定泌乳量ＥＹが一定期間（例えば１週間、１か月等）の単位で求められている場合は、（７）式にさらに期間（日数）を乗じてもよい。
【0107】
［牛群総回帰式演算］
図１４に牛群推定総泌乳量を求めるフローを示す。これは牛群総回帰式演算工程である。牛群総回帰式演算が開始されたら（ステップＳ２１２）、初期設定が行われる（ステップＳ３００）。初期設定としては、牛群推定総泌乳量ＳＵＭＥＹと、個体牛の識別指標ｉを初期化し、識別指標の最後の値（ＥＮＤ）をｎ＋１にセットする。ここでｎは牛群に属する個体牛の頭数である。
図１４に牛群推定総泌乳量を求めるフローを示す。これは牛群総回帰式演算工程である。牛群総回帰式演算が開始されたら（ステップＳ２１２）、初期設定が行われる（ステップＳ３００）。初期設定としては、牛群推定総泌乳量ＳＵＭＥＹと、個体牛の識別指標ｉを初期化し、識別指標の最後の値（ＥＮＤ）をｎ＋１にセットする。ここでｎは牛群に属する個体牛の頭数である。
【0108】
次にｉ番目の個体牛ＣＷｉの分娩後日数Ｎｉを取得する（ステップＳ３０２）。そして、（６）式の回帰式を用いて牛群推定泌乳量ＮｉＥＹを算出し、牛群推定総泌乳量ＳＵＭＥＹに加える（ステップＳ３０４）。なお、要因パラメータＰ（説明変数）は、ステップＳ２１０（牛群回帰式演算工程：図１２参照）の際に用いた値を使ってもよいし、初期設定（ステップＳ３００）の前後に利用者に問い合わせてもよい。
【0109】
全ての個体牛についての演算が終了したかを判断しテップＳ３０６）、まだ残っていれば（ステップＳ３０６のＮ分岐）、識別指標ｉをインクリメントし（ステップＳ３１０）、個体牛ＣＷｉの分娩後日数Ｎｉを取得するステップ（ステップＳ３０２）に戻る。全ての個体牛についての牛群推定泌乳量の総和が終了したら（ステップＳ３０６のＹ分岐）、結果を表示し（ステップＳ３０８）、ＥＹルーチンに戻る（ステップＳ３１２）。
【0110】
図４を再度参照し、牛群回帰式演算若しくは牛群総回帰式演算を繰り返す場合（ステップＳ２１４のＹ分岐）は、再びステップＳ２１０を行う。繰り返しを行わない場合（ステップＳ２１４のＮ分岐）は、終了判定を行う（ステップＳ２１６）。終了する場合（ステップＳ２１６のＹ分岐）は、メインルーチンへ戻る（ステップＳ２１８）。終了しない場合（ステップＳ２１６のＮ分岐）は、ステップＳ２００に戻り、ＥＹ算出工程を再実行する。
【0111】
［要因分析］
図１０には、図４のステップＳ２０６の要因分析の工程の処理を示す。要因分析の工程は、回帰式の作成（ステップＳ２０８）と途中までは同じである。具体的には、図８で示すステップＳ２５０、ステップＳ２５２、ステップＳ２５４は、図８の場合と同じである。したがって、ステップ番号も同じ番号を用いる。
図１０には、図４のステップＳ２０６の要因分析の工程の処理を示す。要因分析の工程は、回帰式の作成（ステップＳ２０８）と途中までは同じである。具体的には、図８で示すステップＳ２５０、ステップＳ２５２、ステップＳ２５４は、図８の場合と同じである。したがって、ステップ番号も同じ番号を用いる。
【0112】
したがって、要因パラメータＰを入力し（ステップＳ２５０）、要因パラメータＰについて指定された一定期間ｍ（これはステップＳ２００で定められる）毎に、分娩後日数Ｎ毎の乳量Ｙｍを、補間式（３）によって求め（ステップＳ２５２）、要因パラメータＰを説明変数、乳量Ｙｍを目標変数として、回帰式を求める（ステップＳ２５４）工程が行われる。
【0113】
要因分析では、一定期間ｍの乳量の実乳量ＲＹと回帰式から求めた牛群推定泌乳量ＥＹを比較する（ステップＳ２６０）。図１１はこの様子を図にしたものである。要因パラメータとしては、ｐ１とｐ２の２つを選んだ場合を示す。例えば、温度と風速等である。一定期間は異なる４日の場合を示した。例えば、春夏秋冬の代表的な1日としてよい。具体的には、ｍ１月ｎ１日、ｍ２月ｎ２日、ｍ３月ｎ３日およびｍ４月ｎ４日とする。これらの日の具体的な要因パラメータｐ１、ｐ２は記録として残っている。要因パラメータの選択は、別途主要因分析を行ってもよいし、試行錯誤によって選択してもよい。要因パラメータｐ１、ｐ２を入力すれば回帰式（６）から牛群推定泌乳量ＥＹを算出することができる。
【0114】
分娩後日数Ｎは１から最大日まで求めることができるが、ここでは、３０日、５０日、１００日、１５０日の４種について求めた状態を示している。これらの値をそれぞれの日時の実乳量ＲＹとの差の絶対値ＥＲＲとする。なお、図１１で３０ＲＹは、分娩後日数が３０日の乳牛の実乳量である。また、３０ＥＲＲは、分娩後日数が３０日の場合の牛群推定泌乳量ＥＹと実乳量ＲＹの差の絶対値を示す。このＥＲＲは、分娩後日数が５０日、１００日、１５０日の場合も算出している。なお、後述する実施例では、分娩後日数が５０日の場合について、牛群推定泌乳量ＥＹと実乳量ＲＹの比較を行った例を示した。
【0115】
４日分について、牛群推定泌乳量ＥＹと実乳量ＲＹの差の絶対値ＥＲＲを合計したものを総誤差ＴΣで表す。このようにステップＳ２６０では、牛群推定泌乳量ＥＹと実乳量ＲＹを比較する。
【0116】
再度図１０を参照して、牛群推定泌乳量ＥＹと実乳量ＲＹを比較した後は、終了判定を行う（ステップＳ２６２）。終了判定は、ステップＳ２５０で入力した要因パラメータＰが牛群推定泌乳量ＥＹと実乳量ＲＹの差の絶対値ＥＲＲが妥当であるか否かで決める。ＥＲＲが妥当であるか否かについては、予め値を決めておいてもよいし、利用者が判断してもよい。また、最適値を選択するようにプログラム（所謂「ＡＩ」を含んでよい。）を設定してもよい。また、牛群推定泌乳量ＥＹと実乳量ＲＹの差は絶対値として説明したが、２乗値を用いてもよい。
【0117】
終了する場合（ステップＳ２６２のＹ分岐）は、ＥＹ算出フローに戻る（ステップＳ２６４）。継続する場合（ステップＳ２６２のＮ分岐）は、再びステップＳ２５０に戻り、要因パラメータＰを入力しなおす。
【0118】
再度図４を参照する。要因分析（ステップＳ２０６）を抜けたら、牛群回帰式演算の工程（ステップＳ２１０）に処理を移す。要因分析の工程（ステップＳ２０６）を抜けたら、回帰式および好適な要因パラメータＰを見つけていると考えられるので、牛群推定泌乳量ＥＹを算出する牛群回帰式演算の工程（ステップＳ２１０）を行ってよいと考えられるからである。牛群回帰式演算の工程（ステップＳ２１０）では、すでに述べたように、予想日に想定される要因パラメータＰの値（説明変数）と所望の分娩後日数Ｎから（６）式の回帰式を使って牛群推定泌乳量ＥＹを演算し求める。これ以降の工程は上記に説明した通りである。
【0119】
［回帰式構成工程］
なお、本発明に係る泌乳量算出システム１においては、回帰式（６）を求める工程は、いずれのメニューでも利用される。したがって、ステップＳ２００、ステップＳ２０２、ステップＳ２０４、ステップＳ２０６、ステップＳ２０８までをまとめてステップＳ２９０とする。ステップＳ２９０は、回帰式構成工程と呼んでもよい。
なお、本発明に係る泌乳量算出システム１においては、回帰式（６）を求める工程は、いずれのメニューでも利用される。したがって、ステップＳ２００、ステップＳ２０２、ステップＳ２０４、ステップＳ２０６、ステップＳ２０８までをまとめてステップＳ２９０とする。ステップＳ２９０は、回帰式構成工程と呼んでもよい。
【0120】
牛群推定泌乳量ＥＹの算出（ステップＳ１１０）は、要因パラメータＰ（説明変数）の値を入力すると、牛群推定泌乳量ＥＹを算出する。特に指定されなければ、全ての分娩後日数について牛群推定泌乳量ＥＹを算出する。もちろん、特定の分娩後日数についてだけ算出してもよい。この場合の要因パラメータＰの値は、予測日において想定される要因パラメータＰである。予測日において想定される要因パラメータＰは、母集団データ作成条件（ステップＳ２００）の設定の際に入力されていてもよい。また、牛群推定泌乳量ＥＹを算出する工程であるステップＳ２１０で入力されてもよい。また、要因パラメータＰの再入力を行う際（ステップＳ２１４のＹ分岐でステップＳ２１０に戻った際）にも、再度入力してよい。
【0121】
以上のように牛群推定泌乳量算出は、特定の母集団に対して、要因パラメータＰが、ある値の時の牛群に属する乳牛の平均的な乳量を推定した値を返す。言い返すと、牛群推定泌乳量値とは、要因パラメータＰが、ある値の時のその母集団に属する乳牛から見込める予測乳量である。
【0122】
＜牛群予測売上金Ｅ＄＞
図２を再度参照する。牛群推定泌乳量ＥＹを算出していると（ステップＳ１１０）、そこから牛群予測売上金Ｅ＄を求めることができる。牛群予測売上金Ｅ＄を算出すると選択された場合（ステップS１１２のＹ分岐）は、牛群予測売上金Ｅ＄を算出する（ステップＳ１１６）工程を実施する。牛群予測売上金Ｅ＄を算出する工程では、牛群推定個体売上金ＥＩ＄と、牛群推定総売上金ΣＥ＄が算出可能である。牛群推定個体売上金ＥＩ＄は、牛群推定泌乳量ＥＹに基づいて算出される個体牛の予測売上金である。また、牛群推定総売上金ΣＥ＄は、牛群推定泌乳量ＥＹに基づいて算出される牛群に属する全ての個体牛の総売上金である。
図２を再度参照する。牛群推定泌乳量ＥＹを算出していると（ステップＳ１１０）、そこから牛群予測売上金Ｅ＄を求めることができる。牛群予測売上金Ｅ＄を算出すると選択された場合（ステップS１１２のＹ分岐）は、牛群予測売上金Ｅ＄を算出する（ステップＳ１１６）工程を実施する。牛群予測売上金Ｅ＄を算出する工程では、牛群推定個体売上金ＥＩ＄と、牛群推定総売上金ΣＥ＄が算出可能である。牛群推定個体売上金ＥＩ＄は、牛群推定泌乳量ＥＹに基づいて算出される個体牛の予測売上金である。また、牛群推定総売上金ΣＥ＄は、牛群推定泌乳量ＥＹに基づいて算出される牛群に属する全ての個体牛の総売上金である。
【0123】
牛群推定個体売上金ＥＩ＄は、特定の個体牛ＣＷｉの牛群推定泌乳量ＥＹ[ＣＷｉ]を（６）式から算出し牛乳単価ｑ（円／ｋｇ）を乗ずることで得ることができる。また、牛群推定総売上金ΣＥ＄は個体牛ＣＷｉの牛群推定泌乳量ＥＹ[ＣＷｉ]を全個体牛分加算するか、牛群推定総泌乳量ＳＵＭＥＹに牛乳単価ｑ（円／ｋｇ）を乗ずることで得ることができる。なお、牛乳単価ｑは予測日において想定される金額であるが、特に根拠は必要とせず、季節ごとの過去の実績であってもよい。
【0124】
［牛群予測売上金算出１］
図１５には、牛群推定総売上金ΣＥ＄を牛群推定総泌乳量ＳＵＭＥＹに牛乳単価ｑ（円／ｋｇ）を乗ずることで得る場合のフローを示す。
図１５には、牛群推定総売上金ΣＥ＄を牛群推定総泌乳量ＳＵＭＥＹに牛乳単価ｑ（円／ｋｇ）を乗ずることで得る場合のフローを示す。
【0125】
牛群予測売上金Ｅ＄の工程がスタートしたら（ステップＳ１１６Ａ）、牛群推定個体売上金ＥＩ＄の演算を行うか（ステップＳ３３０）、牛群推定総売上金ΣＥ＄の演算を行うか（ステップＳ３３２）を選択する。
【0126】
牛群推定個体売上金ＥＩ＄の演算を行う場合（ステップＳ３３０のＹ分岐）は、特定の個体牛ＣＷｉの分娩後日数Ｎｉと要因パラメータＰ（説明変数）を入力し、牛群推定泌乳量ＮｉＥＹを算出し、牛乳単価ｑ（円／ｋｇ）を乗じ、牛群推定個体売上金ＥＩ＄を演算する（ステップＳ３４０）。結果は表示してもよい。
【0127】
そして繰り返す場合（ステップＳ３４２のＹ分岐）は、再びステップＳ３４０を繰り返し、そうでない場合（ステップＳ３４２のＮ分岐）は、処理を牛群推定総売上金ΣＥ＄の算出工程（ステップＳ３３２）に移す。
【0128】
牛群推定総売上金ΣＥ＄の演算を行う場合（ステップＳ３３２のＹ分岐）は、まず、牛群推定総泌乳量ＳＵＭＥＹを算出する（ステップＳ２１２）。なお、牛群推定泌乳量ＥＹを算出する工程（ステップＳ１１０：図２参照）において、すでに牛群推定総泌乳量ＳＵＭＥＹを算出していて、その値を使う場合は、このステップは飛ばしてもよい。
【0129】
次に牛群推定総泌乳量ＳＵＭＥＹに牛乳単価ｑ（円／ｋｇ）を乗じ、牛群推定総売上金ΣＥ＄を演算する（ステップＳ３５０）。結果は表示してもよい。
【0130】
繰り返しを行う場合（ステップS３５２のＹ分岐）は、ステップＳ２１２に戻る。そうでない場合（ステップＳ３５２のＮ分岐）は、牛群予測売上金Ｅ＄の工程に戻る（ステップＳ３５４）。
【0131】
以上の工程により、予測日における牛群に属する個体牛の牛群推定泌乳量ＥＹの総和である牛群推定総泌乳量ＳＵＭＥＹに牛乳単価ｑを乗じて得られる牛群推定総売上金ΣＥ＄を算出することができる。
【0132】
［牛群予測売上金算出２］
図１６には、牛群推定総売上金ΣＥ＄を個体牛ＣＷｉの牛群推定泌乳量ＥＹ[ＣＷｉ]に牛乳単価ｑを乗じたものを全個体牛分加算することで得る場合のフローを示す。
図１６には、牛群推定総売上金ΣＥ＄を個体牛ＣＷｉの牛群推定泌乳量ＥＹ[ＣＷｉ]に牛乳単価ｑを乗じたものを全個体牛分加算することで得る場合のフローを示す。
【0133】
牛群予測売上金Ｅ＄の工程がスタートしたら（ステップＳ１１６Ｂ）、牛群推定個体売上金ＥＩ＄の演算を行うか（ステップＳ３３０）、牛群推定総売上金ΣＥ＄の演算を行うか（ステップＳ３３４）を選択する。
【0134】
牛群推定個体売上金ＥＩ＄の演算については（ステップＳ３３０のＹ分岐）、図１５の場合と同じであるので省略する。牛群推定総売上金ΣＥ＄の演算を行う場合（ステップＳ３３４のＹ分岐）は、牛群推定総売上金ΣＥ＄および個体牛の識別指標ｉを初期化し、識別指標の最後の値（ＥＮＤ）をｎ＋１にセットする（ステップＳ３６０）。ここでｎは牛群に属する個体牛の頭数である。また、予測日の要因パラメータＰ（説明変数）の値を取得する。これは別途入力され、若しくは牛群推定泌乳量を求めた際（ステップＳ１１０）に用いた値を用いてもよい。
【0135】
次に予測日のｉ番目の個体牛ＣＷｉの分娩後日数Ｎｉを取得する（ステップＳ３６２）。これは予測起算日における分娩後日数に予測日までの日数を加えることで得ることができる。
【0136】
そして、個体牛の牛群推定泌乳量ＮｉＥＹを算出し、牛乳単価ｑ（円／ｋｇ）を乗じ、牛群推定総売上金ΣＥ＄に加える（ステップＳ３６４）。つまり個体牛ＣＷｉの牛群推定個体売上金ＥＩ＄を求めて牛群推定総売上金ΣＥ＄に加える。
【0137】
全ての個体牛についての加算が終了していなければ（ステップＳ３６６のＮ分岐）、識別指数ｉをインクリメントし（ステップＳ３７０）、予測日の個体牛ＣＷｉの分娩後日数Ｎｉを取得する工程（ステップＳ３６２）に戻る。
【0138】
全ての個体牛について加算が終了したら（ステップＳ３６６のＹ分岐）、結果を表示し、再入力を問い合わせる（ステップＳ３６８）。要因パラメータＰの値を変更したい（予測日を変更する）場合は再入力を行い（ステップＳ３６８のＹ分岐）、初期設定（ステップＳ３６０）に戻る。繰り返さない場合（ステップＳ３６８のＮ分岐）は、牛群予測売上金Ｅ＄に戻る（ステップＳ３７２）。
【0139】
以上の工程により、予測日における牛群に属する個体牛の牛群推定泌乳量ＥＹに牛乳単価ｑを乗じて得られる牛群推定個体売上金ＥＩ＄の総和として牛群推定総売上金ΣＥ＄を算出することができる。
【0140】
なお、上記のステップＳ１１６のＡおよびＢはどちらも図２の牛群予測売上金を算出する工程（ステップＳ１１６）とすることができる。
【0141】
＜牛群乾乳牛ＤＣの判定＞
図２を参照する。牛群推定泌乳量ＥＹを算出していると（ステップＳ１１０）、牛群乾乳牛ＤＣを判定することができる。図１７を用いて牛群乾乳牛ＤＣの判定を説明する。図１７は、横軸が分娩後日数（日）であり、縦軸は牛群推定泌乳量ＥＹ（ｋｇ）である。牛群推定泌乳量ＥＹを算出しているので、牛群推定泌乳曲線Ｍを描くことができる。
図２を参照する。牛群推定泌乳量ＥＹを算出していると（ステップＳ１１０）、牛群乾乳牛ＤＣを判定することができる。図１７を用いて牛群乾乳牛ＤＣの判定を説明する。図１７は、横軸が分娩後日数（日）であり、縦軸は牛群推定泌乳量ＥＹ（ｋｇ）である。牛群推定泌乳量ＥＹを算出しているので、牛群推定泌乳曲線Ｍを描くことができる。
【0142】
牛群乾乳牛ＤＣは、牛群推定泌乳曲線Ｍに対して、乾乳牛候補と判定できる乾乳基準値ＤＹＴＨを適用することで決めることができる。乾乳基準値ＤＹＴＨは、分娩後日数が一定期間を過ぎた後の乳牛を乾乳牛候補と判断する泌乳量である。図１７では、分娩後日数Ｎ５の個体牛は、牛群推定泌乳量ＥＹが乾乳基準値ＤＹＴＨより低いので、牛群乾乳牛ＤＣの候補と決めることができる。一方、分娩後日数Ｎ４より少ない乳牛は、牛群推定泌乳量ＥＹが乾乳基準値ＤＹＴＨより高いので牛群乾乳牛ＤＣの候補にはならない（ＮＯＴＤＣ）。
【0143】
なお、牛群推定泌乳曲線Ｍと乾乳基準値ＤＹＴＨとの交点の日は、牛群乾乳日ＤＹｄａｙである。また、乾乳基準値ＤＹＴＨは、経験値的に求められるもので、図４のステップＳ２００の際に入力しておくなどする。また、出産後直後で泌乳量が少ない個体牛は乾乳牛の候補からは外される。
【0144】
［牛群乾乳牛判定１］
図１８に牛群推定泌乳量ＥＹと乾乳基準値ＤＹＴＨを直接比較して牛群乾乳牛ＤＣの候補を表示するフローを示す。牛群乾乳牛ＤＣ判定の工程が開始されると（ステップＳ１１８Ａ）、初期設定を行う（ステップＳ４００）。初期設定としては、個体識別指標ｉを初期化（ｉ＝１）し、識別指標の最終の値（ＥＮＤ）をｎ＋１にセットする。ここでｎは牛群に属する乳牛の頭数である。
図１８に牛群推定泌乳量ＥＹと乾乳基準値ＤＹＴＨを直接比較して牛群乾乳牛ＤＣの候補を表示するフローを示す。牛群乾乳牛ＤＣ判定の工程が開始されると（ステップＳ１１８Ａ）、初期設定を行う（ステップＳ４００）。初期設定としては、個体識別指標ｉを初期化（ｉ＝１）し、識別指標の最終の値（ＥＮＤ）をｎ＋１にセットする。ここでｎは牛群に属する乳牛の頭数である。
【0145】
次にｉ番目の個体牛ＣＷｉの予測日における分娩後日数Ｎｉを取得する（ステップＳ４０２）。これは、個体牛ＣＷｉの予測起算日における分別後日数に予測日までの日数を加算することで得られる。そして、分娩後日数Ｎｉの時の牛群推定泌乳量ＮｉＥＹを求め、乾乳基準値ＤＹＴＨと比較する（ステップＳ４０４）。
【0146】
牛群泌乳量ＮｉＥＹが乾乳基準値ＤＹＴＨより少なければ（ステップＳ４０４のＹ分岐）、個体牛ＣＷｉを牛群乾乳牛ＤＣの候補として表示し（ステップＳ４０６）、ＣＷｉが最後の個体牛か否かを判断する（ステップＳ４０８）。牛群推定泌乳量ＮｉＥＹが乾乳基準値ＤＹＴＨより少なくなければ（ステップＳ４０４のＮ分岐）、ステップＳ４０６をスキップする。
【0147】
ＣＷｉが最後の個体牛でなければ（ステップＳ４０８のＮ分岐）識別指標ｉをインクリメントし（ステップＳ４１０）、ｉ番目の個体牛ＣＷｉの予測日における分娩後日数Ｎｉを取得するステップ（ステップＳ４０２）に戻る。最後の個体牛まで判定が終了したら（ステップＳ４０８のＹ分岐）、ＤＣ判定に戻る（ステップＳ４１２）。
【0148】
以上の工程により、予測日における牛群推定泌乳量ＥＹが乾乳基準値ＤＹＴＨより低い個体牛を牛群乾乳牛候補として表示することができる。
【0149】
［牛群乾乳牛判定２］
また、牛群乾乳牛の判定は、予測日の牛群推定泌乳曲線Ｍにおいて、牛群乾乳日ＤＹｄａｙを求め、分娩後日数がその牛群乾乳日ＤＹｄａｙ以上の個体牛を表示することでも得ることができる。また、牛群に属する全ての個体牛の予測日における牛群推定泌乳量が乾乳基準値ＤＹＴＨより少ない個体牛を表示することでも得ることができる。
また、牛群乾乳牛の判定は、予測日の牛群推定泌乳曲線Ｍにおいて、牛群乾乳日ＤＹｄａｙを求め、分娩後日数がその牛群乾乳日ＤＹｄａｙ以上の個体牛を表示することでも得ることができる。また、牛群に属する全ての個体牛の予測日における牛群推定泌乳量が乾乳基準値ＤＹＴＨより少ない個体牛を表示することでも得ることができる。
【0150】
図１９には、牛群乾乳牛ＤＣを牛群乾乳日から判定するフローを示す。牛群乾乳牛ＤＣ判定が実行されると（ステップＳ１１８Ｂ）、分娩後日数Ｎにａ日を入れる（ステップＳ４３０）。ａ日は、牛群推定泌乳量ＥＹが最初に乾乳基準値ＤＹＴＨを超える日（ＤＣ０：図１７参照）以降の日を入れる。ａ日は、通常３０日～５０日を設定すればよい。「ａ日」を早期日数と呼ぶ。
【0151】
次に分娩後日数Ｎは最終日か否かを判断する（ステップＳ４３２）。最終日であれば（ステップＳ４３２のＹ分岐）、「ＤＣは該当なし」として（ステップＳ４４２）、ＤＣ判定ステップ（Ｓ１１８Ｂ）に戻る（ステップＳ４４４）。
【0152】
分娩後日数Ｎが最終日でない場合（ステップＳ４３２のＮ分岐）は、分娩後日数がＮの時の牛群推定泌乳量ＮＥＹが乾乳基準値ＤＹＴＨ以下であるか否かを判断する（ステップＳ４３４）。
【0153】
分娩後日数がＮの時の牛群推定泌乳量ＮＥＹが乾乳基準値ＤＹＴＨ以下であった場合（ステップＳ４３４のＹ分岐）、牛群乾乳日ＤＹｄａｙをＮ日とし（ステップＳ４３６）、牛群乾乳日ＤＹｄａｙより分娩後日数Ｎｉが大きくなる個体牛ｉを牛群乾乳牛ＤＣ候補として表示する（ステップＳ４３８）。その後ＤＣ算出ステップ（Ｓ１１８）に戻る（ステップＳ４４４）。
【0154】
分娩後日数がＮの時の牛群推定泌乳量ＮＥＹが乾乳基準値ＤＹＴＨ以下でなかった場合（ステップＳ４３４のＮ分岐）は、分娩後日数Ｎをインクリメントし（ステップＳ４４０）、ステップＳ４３２に処理を戻す。
【0155】
以上の工程によって、記予測日における牛群推定泌乳量ＥＹが乾乳基準値ＤＹＴＨとなる牛群乾乳日ＤＹｄａｙより分娩後日数が多くなる個体牛を牛群乾乳牛ＤＣ候補として表示することができる。
【0156】
なお、上記のステップＳ１１８のＡおよびＢはどちらも図２の牛群乾乳牛を判定する工程（ステップＳ１１８）とすることができる。
【0157】
＜個体推定泌乳量ＤＥＹ＞
再度図２を参照する。次に個体推定泌乳量ＤＥＹについて説明する。牛群推定泌乳量ＥＹでは、全ての乳牛が同一の牛群推定泌乳量ＥＹに従うとした。しかし、個体推定泌乳量ＤＥＹは、個体牛毎に泌乳量を推定するものである。個体推定泌乳量ＤＥＹは、個体牛毎に算出された乳量係数を牛群推定泌乳量ＥＹに乗じたものとして得られる。
再度図２を参照する。次に個体推定泌乳量ＤＥＹについて説明する。牛群推定泌乳量ＥＹでは、全ての乳牛が同一の牛群推定泌乳量ＥＹに従うとした。しかし、個体推定泌乳量ＤＥＹは、個体牛毎に泌乳量を推定するものである。個体推定泌乳量ＤＥＹは、個体牛毎に算出された乳量係数を牛群推定泌乳量ＥＹに乗じたものとして得られる。
【0158】
図２０に個体推定泌乳量ＤＥＹを求める様子を示す。牛群中のｉ番目の個体牛を個体牛ＣＷiとする。個体牛ＣＷiが牛群推定泌乳量ＥＹに対してどの程度の乳量係数を持つかを算出する日を「個体特定日」と呼ぶ。個体特定日は、予測起算日以前の日であればよく、個体牛ＣＷiの実乳量ＲＹ[ＣＷi]が決定している最新の日でよい。なお、ＲＹ[ＣＷi]は、個体牛ＣＷiの実乳量である。
【0159】
図２０（ａ）は個体特定日（この日の要因パラメータはｐ１）における牛群推定泌乳量ＥＹを表す。牛群推定泌乳量ＥＹは牛群推定泌乳曲線Ｍで表している。横軸は分娩後日数であり、縦軸は牛群推定泌乳量ＥＹである。個体特定日の個体牛ＣＷiは分娩後日数がＮｒで、実乳量がＲＹ[ＣＷi]であったとする。なお、分娩後日数がＮｒの個体牛ＣＷiの実乳量はＮｒＲＹ[ＣＷi]と表す。また、分娩後日数Ｎｒの時の牛群推定泌乳量は、ＮｒＥＹで表される。
【0160】
なお、ここでは個体牛ＣＷiの実乳量ＮｒＲＹ[ＣＷi]が牛群推定泌乳量ＮｒＥＹより大きいとした。もちろん逆であってもよい。そして、個体牛ＣＷiの乳量係数αｉを（８）式のように求める。このようなステップは乳量係数算出工程と言える。
【0161】
【数8】
【0162】
αｉは、各個体牛毎に決まる値となる。αｉを算出する対象となる個体牛ＣＷｉを特定個体牛とも呼ぶ。なお、乳量係数αｉは、実乳量ＮｒＲＹ[ＣＷi]および牛群推定泌乳量ＮｒＥＹともに、一定期間の平均値であってもよい。また、連続する一定期間の泌乳量係数を平均してもよい。例えば、予測起案日の直近の３か月において、１か月毎の実乳量ＮｒＲＹ[ＣＷi]および牛群推定泌乳量ＮｒＥＹから各月の乳量係数αｉを求め、３か月分の乳量係数を平均することで乳量係数αｉとするなどである。なお、予測起算日は当日だけでなく、予測起算日以前の一定期間を含めることができるので、このような算出をしても乳量係数αｉは予測起算日における値と言ってもよい。
【0163】
図２０（ｂ）を参照して、次に予測日（この日の要因パラメータはｐｋ）の牛群推定泌乳量ＥＹを算出する。図２０（ｂ）では、牛群推定泌乳曲線Ｍとして表されている。図２０（ａ）の牛群推定泌乳曲線Ｍと図２０（ｂ）の牛群推定泌乳曲線Ｍは要因パラメータＰの値が異なっていれば、異なる曲線となる。個体牛ＣＷｉの個体推定泌乳量ＤＥＹ[ＣＷi]は、牛群推定泌乳量ＥＹのαｉ倍として求める。
【0164】
牛群推定泌乳量ＥＹは（６）式の回帰式で表される。したがって、個体推定泌乳量ＤＥＹを求める工程は、（６）式の回帰式をαｉ倍して求められる個体回帰式を算出する工程と言ってもよい。なお、個体推定泌乳量ＤＥＹ[ＣＷi]は個体推定泌乳曲線ＤＭとして描くこともできる。したがって、個体推定泌乳曲線ＤＭは、牛群推定泌乳曲線Ｍのαｉ倍である（図２０（ｂ）参照）。
【0165】
次に予測日の個体牛ＣＷｉの分娩後日数Ｎｒ＋Ｘを求める。これは予測日において個体牛ＣＷｉの分娩後日数が何日になるかで求めることができる。ここでＸは個体特定日から予測日までの日数である。そして、個体推定泌乳量Ｎｒ＋ＸＤＥＹ[ＣＷｉ]を求める。これはαｉ（Ｎｒ＋ＸＥＹ）を求めることに等しい。なお、Ｎｒ＋ＸＤＥＹ[ＣＷｉ]は、分娩後日数がＮｒ＋Ｘ（日）の個体牛ＣＷｉの個体推定泌乳量である。また、Ｎｒ＋ＸＥＹは、分娩後日数がＮｒ＋Ｘ（日）の牛群推定泌乳量である。
【0166】
つまり、予測日における個体推定泌乳量ＤＥＹは、予測日における分娩後日数がＮｒ＋Ｘ（日）の牛群泌乳量ＥＹのαｉ倍である。この工程は個体推定泌乳量を算出する工程である。
【0167】
図２１を参照して、個体推定総泌乳量ＳＵＭＤＥＹについて説明する。個体推定泌乳量ＤＥＹは、図２０のようにして求めることとした。したがって、牛群に属する全ての乳牛について個々に個体推定泌乳量ＤＥＹが求められる。予測日における個体推定泌乳量ＤＥＹの総和が個体推定総泌乳量ＳＵＭＤＥＹである。
【0168】
図２１は横軸が分娩後日数（日）であり、縦軸は乳量（ｋｇ）である。簡単な説明のため、今牛群は個体牛Ａおよび個体牛Ｂで構成されているとする。グラフ中には、予測日における牛群推定泌乳量ＥＹ（牛群推定泌乳曲線Ｍ）と個体牛Ａおよび個体牛Ｂの個体推定泌乳量ＮＡＤＥＹ[Ａ]、ＮＢＤＥＹ[Ｂ]を示した。ＮＡ、ＮＢは個体牛Ａおよび個体牛Ｂの予測日における分娩後日数である。なお、個体推定泌乳量ＤＥＹ[Ａ]、ＤＥＹ[Ｂ]は、個体推定泌乳曲線ＤＭＡ、ＤＭＢとして表した。また、個体牛Ａおよび個体牛Ｂの乳量係数をそれぞれαＡ、αＢとした。
【0169】
また、個体牛Ａの個体推定泌乳量ＤＥＹ[Ａ]は牛群推定泌乳量ＥＹより多く、個体牛Ｂの個体推定泌乳量ＤＥＹ[Ｂ]は牛群推定泌乳量ＥＹより少ないとした。また、個体牛Ａおよび個体牛Ｂの予測日の分娩後日数は予測日においてそれぞれＮＡ（日）とＮＢ（日）とした。
【0170】
個体推定総泌乳量ＳＵＭＤＥＹは、個体牛Ａの分娩後日数ＮＡでの個体推定泌乳量ＮＡＤＥＹ[Ａ]と個体牛Ｂの分娩後日数ＮＢでの個体推定泌乳量ＮＢＤＥＹ[Ｂ]の和として求めることができる。この工程は個体推定泌乳量を算出する工程である。
【0171】
図２２に個体推定泌乳量ＤＥＹ算出のフローを示す。図２の個体推定泌乳量ＤＥＹ算出の工程がスタートしたら（ステップＳ１２０）、個体特定日の回帰式構成の工程を行う（ステップＳ５００）。これは、図４のステップＳ２９０を実行し、個体特定日（乳量係数αを決める日）における（６）式の回帰式を得ることである。個体特定日における牛群泌乳量ＥＹを求めると言ってもよい。
【0172】
次に個体回帰式演算ＤＥＹＩの工程を行う（ステップＳ５０２）。ステップＳ５００で、（６）式の回帰式が得られたので、個体特定日の実乳量ＲＹを用いることで乳量係数αを決定することができる。その結果、予測日における個体牛ＣＷｉの個体推定泌乳量ＤＥＹを算出することができる。個体回帰式演算ＤＥＹＩの工程は結果としてＣＷｉの個体推定泌乳量ＤＥＹを算出する。
【0173】
次に個体総回帰式演算ＤＥＹＳの工程（ステップＳ５０４）を行う。予測日における個体推定泌乳量ＤＥＹの算出が可能となったので、牛群に属する全ての個体牛についての個体推定泌乳量ＤＥＹを加算することで、個体推定総泌乳量ＳＵＭＤＥＹを求める。個体総回帰式演算ＤＥＹＳの工程は結果として個体推定総泌乳量ＳＵＭＤＥＹを算出する。
【0174】
そして、再入力の有無を問う（ステップＳ５０６）。再入力を行う場合（ステップＳ５０６のＹ分岐）は、個体特定日の回帰式構成の工程（ステップＳ５００）に戻り、個体特定日の変更などで繰り返し処理を実行できる。また、再入力を行わない場合（ステップＳ５０６のＮ分岐）は、処理をＤＥＹ算出工程（ステップＳ１２０）に戻す（ステップＳ５０８）。
【0175】
［個体回帰式演算］
図２３に個体回帰式演算の工程のフローを示す。個体回帰式演算のフローがスタートしたら（ステップS５０２）、個体特定日の個体牛ＣＷｉの分娩後日数Ｎｒ（日）と個体特定日の実乳量ＮｒＲＹ[ＣＷｉ]を取得する（ステップＳ５２０）。これは基本データ中に記録されている。
図２３に個体回帰式演算の工程のフローを示す。個体回帰式演算のフローがスタートしたら（ステップS５０２）、個体特定日の個体牛ＣＷｉの分娩後日数Ｎｒ（日）と個体特定日の実乳量ＮｒＲＹ[ＣＷｉ]を取得する（ステップＳ５２０）。これは基本データ中に記録されている。
【0176】
次に個体特定日の牛群推定泌乳量ＮｒＥＹを算出する（ステップＳ５２２）。具体的には（６）式が求められる。これは図２２のステップＳ５００で求めた結果をそのまま用いてもよい。
【0177】
次に実乳量ＮｒＲＹ[ＣＷｉ]と牛群推定泌乳量ＥＹから（８）式に基づいて乳量係数αｉを算出する（ステップＳ５２４）。そして、予測日における牛群推定泌乳量Ｎｒ＋ＸＥＹを算出し（ステップＳ５２６）、それをαｉ倍することで、個体推定泌乳量ＤＥＹを得る（ステップＳ５２８）。Ｘは個体特定日から予測日までの日数である。その後個体回帰式演算の工程に戻る（ステップＳ５３０）。ステップＳ５２８は（９）式の演算が行われる。
【0178】
【数9】
【0179】
なお、ここで得られる個体推定泌乳量ＤＥＹは、露わに書くとＮｒ＋ＸＤＥＹ[ＣＷｉ]である。これは牛群中のｉ番目の個体牛ＣＷｉの分娩後日数がＮｒ＋Ｘ（日）の推定泌乳量を意味する。算出された個体推定泌乳量ＤＥＹは端末等に表示することができる。また、Ｎｒ（日）は個体特定日の個体牛ＣＷｉの分娩後日数であり、Ｘ（日）は個体特定日から予測日までの日数である。
【0180】
［個体総回帰式演算］
図２４に個体総回帰式演算の工程のフローを示す。個体総回帰式演算の工程がスタートすると（ステップＳ５０４）、まず初期設定が行われる（ステップＳ５４０）。初期設定は、識別指標ｉと、個体推定総泌乳量ＳＵＭＤＥＹの初期化（ｉ＝１、ＳＵＭＤＥＹ＝０）および個体数の最終の値（ＥＮＤ）の設定である。最終の値（ＥＮＤ）はｎ＋１に設定し、「ｎ」は牛群中の乳牛の個体数である。
図２４に個体総回帰式演算の工程のフローを示す。個体総回帰式演算の工程がスタートすると（ステップＳ５０４）、まず初期設定が行われる（ステップＳ５４０）。初期設定は、識別指標ｉと、個体推定総泌乳量ＳＵＭＤＥＹの初期化（ｉ＝１、ＳＵＭＤＥＹ＝０）および個体数の最終の値（ＥＮＤ）の設定である。最終の値（ＥＮＤ）はｎ＋１に設定し、「ｎ」は牛群中の乳牛の個体数である。
【0181】
次に予測日の個体推定泌乳量ＤＥＹ（より正確には「Ｎｒ＋ＸＤＥＹ[ＣＷｉ]」）を算出する（ステップＳ５４２）。Ｎｒ（日）は個体特定日の分娩後日数であり、Ｎｒ＋Ｘ（日）で予測日の分娩後日数となる。そして、個体推定泌乳量ＤＥＹと個体推定総泌乳量ＳＵＭＤＥＹを加算し、新たな個体推定総泌乳量ＳＵＭＤＥＹとする（ステップＳ５４４）。
【0182】
全ての個体牛についての加算が終了したかを判定し（ステップＳ５４６）、終了してなければ（ステップＳ５４６のＮ分岐）、識別指標ｉをインクリメントし（ステップＳ５４８）、ステップＳ５４２に処理を戻す。終了していれば（ステップＳ５４６のＹ分岐）、処理を個体総回帰式演算の工程に戻す（ステップＳ５５０）。
【0183】
なお、ここで得られる個体推定総泌乳量ＳＵＭＤＥＹは、予測日における個体推定泌乳量の総和である。算出された個体推定総泌乳量ＳＵＭＤＥＹは、端末等に表示することができる。
【0184】
＜個体予測売上金ＤＥ＄＞
図２を参照して、個体推定泌乳量ＤＥＹが求まったら、個体予測売上金ＤＥ＄を算出することができる。個体予測売上金ＤＥ＄を算出すると選択された場合（ステップS１２２のＹ分岐）は、個体予測売上金ＤＥ＄を算出する（ステップＳ１２６）。個体予測売上金ＤＥ＄は、各個体牛の推定売上金である個体推定売上金ＤＥＩ＄と、牛群に属する全ての個体牛による推定売上金である個体推定総売上金ΣＤＥ＄を求めることができる。
図２を参照して、個体推定泌乳量ＤＥＹが求まったら、個体予測売上金ＤＥ＄を算出することができる。個体予測売上金ＤＥ＄を算出すると選択された場合（ステップS１２２のＹ分岐）は、個体予測売上金ＤＥ＄を算出する（ステップＳ１２６）。個体予測売上金ＤＥ＄は、各個体牛の推定売上金である個体推定売上金ＤＥＩ＄と、牛群に属する全ての個体牛による推定売上金である個体推定総売上金ΣＤＥ＄を求めることができる。
【0185】
個体推定売上金ＤＥＩ＄は、ステップＳ１２０（図２参照）で求めた個体推定泌乳量ＤＥＹに牛乳の単価ｑ（円／ｋｇ）を乗じることで得ることができる。一方、個体推定総売上金ΣＤＥ＄は、個体推定泌乳量ＤＥＹに牛乳単価ｑを乗じて、個体牛毎の個体推定売上金を算出し、その総和を求める方法と、個体推定総泌乳量ＳＵＭＤＥＹに牛乳単価ｑを乗じて求める方法がある。個体推定総売上金ΣＤＥ＄を求める工程としてはどちらの工程であってもよい。
【0186】
【0187】
［個体推定総売上金算出１］
図２５には、個体予測売上金ＤＥ＄算出工程のフローを示す。なおこのフローは個体推定総売上金ΣＤＥ＄を求めるのに、個体推定総泌乳量ＳＵＭＤＥＹを用いる。個体予測売上金ＤＥ＄を算出する工程がスタートしたら（ステップＳ１２６Ａ）、個体推定売上金ＤＥＩ＄の算出か、個体推定総売上金ΣＤＥ＄の算出かを選択する（ステップＳ５６０、５６２）。それぞれ選択しない場合は処理を次に移す。つまりステップＳ５６０のＮ分岐では、ステップＳ５６２に処理が移り、ステップＳ５６２のＮ分岐では、処理を個体予測売上金ＤＥ＄算出工程に戻す。
図２５には、個体予測売上金ＤＥ＄算出工程のフローを示す。なおこのフローは個体推定総売上金ΣＤＥ＄を求めるのに、個体推定総泌乳量ＳＵＭＤＥＹを用いる。個体予測売上金ＤＥ＄を算出する工程がスタートしたら（ステップＳ１２６Ａ）、個体推定売上金ＤＥＩ＄の算出か、個体推定総売上金ΣＤＥ＄の算出かを選択する（ステップＳ５６０、５６２）。それぞれ選択しない場合は処理を次に移す。つまりステップＳ５６０のＮ分岐では、ステップＳ５６２に処理が移り、ステップＳ５６２のＮ分岐では、処理を個体予測売上金ＤＥ＄算出工程に戻す。
【0188】
個体推定売上金ＤＥＩ＄を選択すると（ステップＳ５６０のＹ分岐）、予測日における個体牛ＣＷｉの個体推定泌乳量を算出し、牛乳単価ｑを乗じて個体推定泌乳量ＤＥＩ＄とする（ステップＳ５７０）。どの個体牛にするかはこのステップの始めに入力されているとする。次に再入力の確認を行い（ステップＳ５７２）、再入力を行う場合（ステップＳ５７２のＹ分岐）は、個体牛の変更などを行い、再度ステップＳ５７０を実行する。再入力を行わない場合（ステップＳ５７２のＮ分岐）は、処理を次に移す（ステップＳ５６２）。
【0189】
個体推定総売上金ΣＤＥ＄を選択すると（ステップＳ５６２のＹ分岐）、まず。個体推定総泌乳量ＳＵＭＤＥＹを算出する（ステップＳ５０４）。これは図２４のステップＳ５０４を実行することで求められる。しかし、個体推定総泌乳量ＳＵＭＤＥＹがすでに算出されている場合は、その値を用いてもよい。
【0190】
そして、個体推定総泌乳量ＳＵＭＤＥＹに牛乳単価ｑを乗じた値を個体推定総売上金ΣＤＥ＄とする（ステップＳ５８０）。算出された個体推定総売上金ΣＤＥ＄は端末等に表示してもよい。
【0191】
次に再入力の確認を行い（ステップＳ５８２）、再入力を行う場合（ステップＳ５８２のＹ分岐）は、予測日等の変更などを行い、再度ステップＳ５０４を実行する。再入力を行わない場合（ステップＳ５８２のＮ分岐）は、処理を個体予測売上金ＤＥ＄の工程に戻す（ステップS５８４）。
【0192】
以上の工程により牛群に属する全ての個体牛の予定日における個体推定泌乳量ＤＥＹの合計である個体推定総泌乳量ＳＵＭＤＥＹに牛乳単価ｑを乗じた額を個体推定総売上金ΣＤＥ＄として算出することができる。
【0193】
［個体推定総売上金算出２］
図２６には、個体推定総売上金ΣＤＥ＄を、各個体牛の推定売上金の総和として求める方法の場合のフローを示す。個体予測売上金ＤＥ＄を算出する工程がスタートしたら（ステップＳ１２６Ｂ）、個体推定売上金ＤＥＩ＄の算出か、個体推定総売上金ΣＤＥ＄の算出かを選択する（ステップＳ５６０、５６４）。それぞれ選択しない場合は処理を次に移す。つまりステップＳ５６０のＮ分岐では、ステップＳ５６４に処理が移り、ステップＳ５６４のＮ分岐では、処理を個体予測売上金ＤＥ＄算出工程に戻す（ステップＳ６００）。
図２６には、個体推定総売上金ΣＤＥ＄を、各個体牛の推定売上金の総和として求める方法の場合のフローを示す。個体予測売上金ＤＥ＄を算出する工程がスタートしたら（ステップＳ１２６Ｂ）、個体推定売上金ＤＥＩ＄の算出か、個体推定総売上金ΣＤＥ＄の算出かを選択する（ステップＳ５６０、５６４）。それぞれ選択しない場合は処理を次に移す。つまりステップＳ５６０のＮ分岐では、ステップＳ５６４に処理が移り、ステップＳ５６４のＮ分岐では、処理を個体予測売上金ＤＥ＄算出工程に戻す（ステップＳ６００）。
【0194】
個体推定売上金ＤＥＩ＄の算出フローは、図２５の場合と同様であるので省略する。個体推定総売上金ΣＤＥ＄の算出が選択されたら（ステップＳ５６４のＹ分岐）、初期設定を行う（ステップＳ５９０）。初期設定は、識別指標ｉと、個体推定総売上金ΣＤＥ＄の初期化（ｉ＝１、ΣＤＥ＄＝０）および個体数の最終の値（ＥＮＤ）の設定である。最終値ＥＮＤはｎ＋１に設定し、「ｎ」は牛群中の乳牛の個体数である。
【0195】
次に予測日の個体推定泌乳量ＤＥＹを求め、牛乳単価ｑを乗じて、個体推定総売上金ΣＤＥ＄に加え、新たな個体推定総売上金ΣＤＥ＄とする（ステップＳ５９２）。最後の個体牛までの加算が終了したかを確認し（ステップＳ５９４）、終わってなければ（ステップＳ５９４のＮ分岐）、識別指標ｉをインクリメントし（ステップＳ５９８）、個体推定総売上金ΣＤＥ＄の加算（ステップＳ５９２）を繰り返す。
【0196】
最後の個体牛まで算出が終了したら（ステップＳ５９４のＹ分岐）、再入力の有無を確認する（ステップＳ５９６）。この際、得られた個体推定総売上金ΣＤＥ＄を端末などに表示してもよい。再入力を行う場合（ステップＳ５９６のＹ分岐）は、予測日などを変更し、初期化工程（ステップＳ５９０）から処理を繰り返す。再入力しない場合（ステップＳ５９６のＮ分岐）は、処理を個体予測売上金ＤＥ＄に戻す（ステップＳ６００）。
【0197】
以上の工程によって、牛群中の全ての個体牛について、個体推定泌乳量に牛乳単価ｑを乗じて個体毎の売上金を算出し、それらの総和を、個体推定総売上金ＤＥ＄として算出することができる。
【0198】
なお、上記のステップＳ１２６のＡおよびＢはどちらも図２の個体予測売上金ＤＥ＄を算出する工程（ステップＳ１２６）とすることができる。
【0199】
＜個体乾乳牛ＤＤＣの判定＞
図２を参照する。個体推定泌乳量ＤＥＹを算出していると（ステップＳ１２０）、個体乾乳牛ＤＤＣも判定することができる。図２７を用いて個体乾乳牛ＤＤＣの判定を説明する。図２７は、横軸が分娩後日数（日）であり、縦軸は個体推定泌乳量ＤＥＹ（ｋｇ）である。なお、牛群推定泌乳量ＥＹも、牛群推定泌乳曲線Ｍとして描いてある。
図２を参照する。個体推定泌乳量ＤＥＹを算出していると（ステップＳ１２０）、個体乾乳牛ＤＤＣも判定することができる。図２７を用いて個体乾乳牛ＤＤＣの判定を説明する。図２７は、横軸が分娩後日数（日）であり、縦軸は個体推定泌乳量ＤＥＹ（ｋｇ）である。なお、牛群推定泌乳量ＥＹも、牛群推定泌乳曲線Ｍとして描いてある。
【0200】
個体乾乳牛ＤＤＣは、個体推定泌乳曲線ＤＭに対して、乾乳牛と判断できる乾乳基準値ＤＹＴＨを適用することで求めることができる。図２７は、この個体牛の予測日における分娩後日数がＮ６より多ければ、この個体牛は個体乾乳牛ＤＤＣの候補となる。分娩後日数がＮ６より少なければ、この個体牛は個体乾乳牛ではない（ＮＯＴＤＤＣ）。なお、個体推定泌乳曲線ＤＭと乾乳基準値ＤＹＴＨの交点の日（Ｎ６）は個体乾乳日ＤＤＹｄａｙである。
【0201】
【0202】
個体乾牛の判別は２通りの方法が考えられる。個体牛毎に個体推定泌乳量ＤＥＹを算出し、それを乾乳基準値ＤＹＴＨと比較する方法と、各個体牛毎に個体推定泌乳量ＤＥＹと乾乳基準値ＤＹＴＨから個体乾乳日ＤＤＹｄａｙを求め、予測日におけるその個体の分娩後日数と個体乾乳日ＤＤＹｄａｙを比較する方法である。いずれの方法も図２の個体乾乳牛判別の工程（ステップＳ１２８）となり得る。
【0203】
［個体乾乳牛判定１］
図２８には、個体牛毎に個体推定泌乳量ＤＥＹを算出し、それを乾乳基準値ＤＹＴＨと比較する方法の場合のフローを示す。個体乾乳牛判断の工程がスタートしたら（ステップＳ１２８Ａ）、初期設定を行う（ステップＳ６１０）。初期設定は、識別指標ｉの初期化および個体数の最終の値（ＥＮＤ）の設定である。識別指標ｉはｉ＝１に設定し、最終の値（ＥＮＤ）はｎ＋１に設定する。「ｎ」は牛群中の乳牛の個体数である。
図２８には、個体牛毎に個体推定泌乳量ＤＥＹを算出し、それを乾乳基準値ＤＹＴＨと比較する方法の場合のフローを示す。個体乾乳牛判断の工程がスタートしたら（ステップＳ１２８Ａ）、初期設定を行う（ステップＳ６１０）。初期設定は、識別指標ｉの初期化および個体数の最終の値（ＥＮＤ）の設定である。識別指標ｉはｉ＝１に設定し、最終の値（ＥＮＤ）はｎ＋１に設定する。「ｎ」は牛群中の乳牛の個体数である。
【0204】
次に、予測日におけるｉ番目の個体牛ＣＷｉの分娩後日数Ｎｒ＋Ｘを取得する（ステップＳ６１２）。ここでＮｒは個体特定日の分娩後日数であり、Ｘは個体特定日から予測日までの日数である。
【0205】
そして、予測日における分娩後日数Ｎｒ＋Ｘが個体早期日数Ｄａ以下であるか否かを判断する（ステップＳ６１４）。個体早期日数Ｄａは、牛群推定泌乳量ＥＹが最初に乾乳基準値ＤＹＴＨを超える日（ＤＣ０：図１３参照）である早期日数ａに対応する日であり、個体牛ＣＷｉの個体推定泌乳量ＤＥＹが最初に乾乳基準値ＤＹＴＨを超える日である。通常３０日～５０日を設定すればよい。
【0206】
分娩後日数Ｎｒ＋Ｘが個体早期日数Ｄａ以下であれば、識別指標ｉをインクリメントし別の個体牛で分娩後日数Ｎｒ＋Ｘを取得する工程を繰り返す（ステップＳ６１２）。分娩後間もない個体牛は、乾乳牛の判断から外すためである。
【0207】
そして、個体牛ＣＷｉの個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨ以下か否かを判断する（ステップＳ６１６）。個体牛ＣＷｉの個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨ以下であったら（ステップＳ６１６のＹ分岐）、個体牛ＣＷｉを個体乾乳牛の候補として表示し（ステップＳ６１８）、全ての個体牛についての判断が終了したか否かを判断する（ステップＳ６２０）。個体牛ＣＷｉの個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨ以下でなければ（ステップＳ６１６のＮ分岐）、ステップＳ６１８をスキップし、ステップＳ６２０に処理を移す。
【0208】
最後の個体牛まで判断が終了してなかった(ステップＳ６２０のＮ分岐）、識別指標ｉをインクリメントし、予測日におけるｉ番目の個体牛ＣＷｉの分娩後日数Ｎｒ＋Ｘを取得するステップ（ステップＳ６１２）まで戻る。最後の個体牛まで判断が終了していたら(ステップＳ６２０のＹ分岐）、個体乾乳牛判断の工程に戻る（ステップＳ６２４）。
【0209】
以上の工程によって、予測日における個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨより低くなる個体牛を乾乳牛候補として表示することができる。
【0210】
［個体乾乳牛判定２］
図２９には、個体乾乳日ＤＤＹｄａｙを用いて個体乾乳牛を判断する場合のフローを示す。個体乾乳牛判断の工程がスタートしたら（ステップＳ１２８Ｂ）、初期設定を行う（ステップＳ６４０）。初期設定は、識別指標ｉの初期化および個体数の最終の値（ＥＮＤ）の設定である。識別指標ｉはｉ＝１に設定し、最終の値（ＥＮＤ）はｎ＋１に設定する。「ｎ」は牛群中の乳牛の個体数である。
図２９には、個体乾乳日ＤＤＹｄａｙを用いて個体乾乳牛を判断する場合のフローを示す。個体乾乳牛判断の工程がスタートしたら（ステップＳ１２８Ｂ）、初期設定を行う（ステップＳ６４０）。初期設定は、識別指標ｉの初期化および個体数の最終の値（ＥＮＤ）の設定である。識別指標ｉはｉ＝１に設定し、最終の値（ＥＮＤ）はｎ＋１に設定する。「ｎ」は牛群中の乳牛の個体数である。
【0211】
次に分娩後日数Ｎを個体早期日数Ｄａに設定する（ステップＳ６４２）。そして牛群に属する全ての個体牛についての判定が終了したかを判断する（ステップＳ６４４）。終了した場合（ステップＳ６４４のＹ分岐）は、処理を個体乾乳牛判断工程に戻す（ステップＳ６６０）。まだ判定する個体牛が残っている場合（ステップＳ６４４のＮ分岐）は処理を次に移す。
【0212】
次の処理では、分娩後日数Ｎが最終か否かを判断する（ステップＳ６４６）。分娩後日数Ｎの最終値は、予め設定されているとする。通常は２６０日以降の日数で設定される。Ｎが最終値でなければ（ステップＳ６４６のＮ分岐）、処理を次の処理（ステップＳ６４８）に移し、Ｎが最終値であった場合（ステップＳ６４６のＹ分岐）は、識別指標ｉをインクリメントし（ステップＳ６５８）、分娩後日数Ｎの設定ステップ（ステップＳ６４２）に戻る。
【0213】
Ｎが最終値でない場合（ステップＳ６４６のＮ分岐）は、ｉ番目の個体の予測日の個体推定泌乳量ＤＥＹにおいて、分娩後日数がＮの時の個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨ以下であるか否かを判断する（ステップＳ６４８）。これは、個体牛ＣＷｉの予測日における個体回帰式から求めることができる。個体牛ＣＷｉの予測日における個体回帰式は（６）式に個体牛ＣＷｉの乳量係数αｉを乗じて得られ、その個体回帰式の分娩後日数がＮの場合を求めればよいからである。
【0214】
個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨ以下とならなかったら（ステップＳ６４８のＮ分岐）、Ｎをインクリメントし（ステップＳ６５６）、Ｎが最終日か否かの判断に戻る（ステップＳ６４６）。個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨ以下となったら、その分別後日数Ｎがこの個体牛にとっての個体乾乳日ＤＤＹｄａｙであるとする（ステップＳ６５０）。
【0215】
そして、予測日における個体牛ＣＷｉの分娩後日数が個体乾乳日ＤＤＹｄａｙより大きいかを判断する（ステップＳ６５２）。なお、Ｎｒ[ＣＷｉ]は個体特定日における個体牛ＣＷｉの分娩後日数であり、Ｘは個体特定日から予測日までの日数である。予測日における個体牛ＣＷｉの分娩後日数が個体乾乳日ＤＤＹｄａｙより大きければ（ステップＳ６５２のＹ分岐）、個体牛ＣＷｉを個体乾乳牛ＤＤＣとして表示する（ステップＳ６５４）。
【0216】
予測日における個体牛ＣＷｉの分娩後日数が個体乾乳日ＤＤＹｄａｙより大きくなければ（ステップＳ６５２のＮ分岐）、識別指標ｉをインクリメントし、分娩後日数Ｎの設定ステップ（ステップＳ６４２）に戻る。このように、牛群に属する全ての個体牛について、それぞれの個体乾乳日ＤＤＹｄａｙを求め、予測日におけるその個体牛の分娩後日数と比較することで個体乾乳牛ＤＤＣか否かを判断することができる。
【0217】
以上の工程により、予測日における個体推定泌乳量ＤＥＹが乾乳基準値ＤＹＴＨとなる個体乾乳日ＤＤＹｄａｙより分娩後日数Ｎｒ＋Ｘ（日）が多くなる個体牛を個体乾乳牛候補として表示することができる。
【0218】
なお、上記のステップＳ１２８のＡおよびＢはどちらも図２の個体乾乳牛を判定する工程（ステップＳ１２８）とすることができる。
【0219】
以上のように本発明に係る泌乳量算出システムは、過去の基本データに基づいて将来の乳量を予測する。将来の乳量が予測できれば、その量に基づいて売上高を予測することができる。
【0220】
また、売上高が予測できると、飼料費用（これが経費となる。）との比較が可能となり、生産性を見積もることができる。例えば、飼料費用をユーザーが入力することで、予測した売上金と比較し、飼料費用が予測売上金の３割となるなどの表示をすることも可能となる。その結果、飼料のロスや効率的な飼料配合を検討することができる。特に濃厚飼料の配合の検討なども可能となる。
【0221】
また、脂肪やタンパク質などの含有量を表す乳成分についても同様に予測し表示させてもよい。所定濃度以上時の乳量単価増加、所定濃度未満時の乳量単価減少の取引を行う場合は、所定濃度に応じた乳量単価増減を考慮した売上金を表示させてもよい。
【実施例0222】
日本国内の牧場で得たデータを用いて本発明に係る泌乳量算出システムを実施した。この牧場は、ある地方の複数の牧場で乳牛数は、およそ５０００頭ほどのデータである。母集団は、２０１９年の月毎のデータである。要因パラメータＰは、ＴＨＩと日照時間である。ＴＨＩは月の平均温度と平均湿度の値を用いた。日照時間は月合計値である。結果を図３０に示す。
【0223】
図３０（ａ）は、２０１９年の２月、５月、８月、１１月の５か月のＴＨＩと日照時間の値である。対象の牧場の平均実乳量（ｋｇ／日）と分娩後日数ＮをＷＯＯＤ曲線で近似し、そこから得た値を用いて推定泌乳量ＥＹを要因パラメータＰとしてＴＨＩだけを用いた単回帰で求めた結果を図３０（ｂ）に示す。また、要因パラメータＰとして、ＴＨＩに日照時間（月合計）の２つを用いた重回帰によって推定泌乳量ＥＹを求めた結果を図３０（ｃ）に示す。なお、いずれも分娩後日数が５０日の場合である。
【0224】
要因パラメータＰとしてＴＨＩだけを用いると、標準誤差が０．６５とかなりばらついた。しかし、要因パラメータＰとしてＴＨＩと日照時間を用いると、標準誤差が０．１７
となり、かなり実乳量を推定できていることが分かった。
となり、かなり実乳量を推定できていることが分かった。
【0225】
泌乳量がＴＨＩと日照時間で推定できる点は、これまで知られることがなく、今後の泌乳量の増加に対する対応に大きなインパクトを与えるものである。
【産業上の利用可能性】
【0226】
本発明は、酪農の情報化のために、好適に利用することができる。
【符号の説明】
【0227】
１ 泌乳量算出システム
１０ 端末
１２ 本体
１４ メモリ
１６ 牧場
１８ 外部情報
１０ 端末
１２ 本体
１４ メモリ
１６ 牧場
１８ 外部情報
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
