劉忠華，李浩昊，焦忠云，魏 剛

(1.中國煙草總公司 云南省公司，云南 昆明 650011；2.云南省煙草公司 昆明市公司，云南 昆明 650051)

烤煙收購作為烤煙產(chǎn)業(yè)鏈中的關(guān)鍵一環(huán)，上承生產(chǎn)、下接調(diào)撥，重要性不言而喻.長期以來，因烤煙收購管理相對(duì)粗放及有效管理手段缺失，導(dǎo)致烤煙收購執(zhí)行過程難以把控，質(zhì)量水平起伏波動(dòng)問題時(shí)有發(fā)生.為全面做好精益化烤煙收購管理工作，全面提升收購執(zhí)行過程管理水平，有必要在認(rèn)真梳理收購管理流程基礎(chǔ)上，引入現(xiàn)代科學(xué)管理手段，對(duì)烤煙收購工作進(jìn)行精準(zhǔn)“靶向”管控.

統(tǒng)計(jì)預(yù)測模型作為一種數(shù)據(jù)分析處理技術(shù)，現(xiàn)已廣泛應(yīng)用到各行各業(yè).其核心思想是應(yīng)用統(tǒng)計(jì)建模方法探索事物歷史發(fā)展規(guī)律，通過模型反映出的事物歷史發(fā)展規(guī)律來預(yù)測未來事物的發(fā)展趨勢，由此對(duì)未來事物的發(fā)展過程進(jìn)行監(jiān)督指導(dǎo).統(tǒng)計(jì)預(yù)測模型在烤煙收購管理中的應(yīng)用，重點(diǎn)從收購歷史數(shù)據(jù)入手，應(yīng)用logistic曲線擬合[1]、樣條插值方法，對(duì)收購進(jìn)度、收購均價(jià)進(jìn)行合理預(yù)測，從而對(duì)即將開展的烤煙收購工作提供科學(xué)指導(dǎo).

烤煙收購是煙葉產(chǎn)品流通過程的起點(diǎn)，同時(shí)也是卷煙工業(yè)原料供應(yīng)的起點(diǎn).煙葉生產(chǎn)過程中所受的政策、氣候和地域差異等影響因素以及煙葉商品本身的多樣性和特殊性決定了烤煙收購的復(fù)雜性[2].

本文依據(jù)2018—2019年曲靖市烤煙收購業(yè)務(wù)數(shù)據(jù)，擬運(yùn)用統(tǒng)計(jì)分析方法，將烤煙收購過程管理通過量化的形式表現(xiàn)出來.多層次、多維度對(duì)曲靖市麒麟?yún)^(qū)歷年烤煙收購趨勢及慣性合理狀態(tài)進(jìn)行深層次挖掘，發(fā)現(xiàn)烤煙生產(chǎn)中的規(guī)律，結(jié)合當(dāng)前烤煙生產(chǎn)的需求，建立合理有效的烤煙收購量進(jìn)度指標(biāo)和控制限[3]，提高烤煙收購工作的科學(xué)化、標(biāo)準(zhǔn)化的調(diào)控與管理水平，為基層烤煙收購煙站提供可靠的數(shù)據(jù)依據(jù).

1 材料與方法

在對(duì)大量歷史數(shù)據(jù)進(jìn)行分析并經(jīng)過深入研討后，烤煙收購工作的進(jìn)度，即每日的累計(jì)收購量占比，作為本文研究的觀測變量，為完成最終的收購任務(wù)，就要尋找其隨時(shí)間變化的發(fā)展規(guī)律.

1.1 數(shù)據(jù)的選擇與整理

通過對(duì)過去曲靖市麒麟?yún)^(qū)烤煙歷史收購記錄的統(tǒng)計(jì)，選擇2018年和2019年的烤煙收購數(shù)據(jù)作為原始參考數(shù)據(jù).

研究總體分類：將麒麟?yún)^(qū)7個(gè)烤煙收購站的數(shù)據(jù)作為7個(gè)研究總體，把各個(gè)收購站的收購數(shù)據(jù)作為各自總體的樣本.

研究對(duì)象分類：根據(jù)煙葉商品收購分類標(biāo)準(zhǔn)[4]，把每個(gè)研究總體按總量、部位(上部葉、中部葉、下部葉)、大等級(jí)(上等煙、中等煙、下等煙、低等煙)、42個(gè)主要小等級(jí)劃分作為各自研究對(duì)象.

1.2 定義變量

時(shí)間變量：由于每個(gè)烤煙收購站收購起始時(shí)間不同，所以分別以各收購站收購煙葉的第1天作為起始時(shí)間開始計(jì)數(shù)，用收購天數(shù)作為時(shí)間序列變量，記作ti=1,2,3，…，i代表收購站.

控制變量：計(jì)算各收購站各年度的每日收購量的累計(jì)比例，作為其所屬區(qū)縣公司(各總體)當(dāng)日的樣本數(shù)據(jù).

1.3 Logistic模型原理

Logistic回歸[5]是通過對(duì)數(shù)ln函數(shù)將因變量y來對(duì)應(yīng)一個(gè)概率p，然后將其結(jié)果間接轉(zhuǎn)化成一個(gè)連續(xù)變量.比如我們研究一些現(xiàn)象，其發(fā)生的概率為p，取值范圍為[0，1]，我們就會(huì)很難去用線性模型描述概率p與自變量的關(guān)系，因此我們需要利用Logit變換來進(jìn)行處理，我們通常把出現(xiàn)某種結(jié)果的概率與不出現(xiàn)的概率之比取對(duì)數(shù)處理，變換如下:

(1)

(2)

1.4 樣條插值法

插值思想是在離散-數(shù)據(jù)的基礎(chǔ)上插補(bǔ)連續(xù)函數(shù)，使得連續(xù)函數(shù)所表示曲線經(jīng)過全部的離散點(diǎn)，此事則可利用此連續(xù)函數(shù)計(jì)算出在其他點(diǎn)的估計(jì)值.樣條插值法，就是每2個(gè)點(diǎn)之間確定1個(gè)函數(shù)，這個(gè)函數(shù)就是1個(gè)樣條.函數(shù)不同，樣條也就不同，按照所研究的事項(xiàng)選擇合理的樣條，并按照插值思想進(jìn)行插值即為樣條插值法.

1.5 模型建立

通過各收購站的統(tǒng)計(jì)分析，累計(jì)比例隨時(shí)間的發(fā)展呈現(xiàn)出如圖1中的曲線走勢：麒麟?yún)^(qū)的煙葉收購量累計(jì)比例歷史走勢基本趨近于“S”形分布，前期收購進(jìn)度增長緩慢，隨時(shí)間推移累計(jì)比例增速逐漸加大，到后期又趨于平緩.

圖1 麒麟?yún)^(qū)烤煙收購量累計(jì)比例時(shí)序

為了刻畫累計(jì)比例隨時(shí)間的變化關(guān)系，通過對(duì)圖形的觀察以及用不同曲線擬合作比較后，選擇用logistic函數(shù)模型來擬合曲線.因此建立烤煙收購量累計(jì)比例隨時(shí)間變化的logistic非線性回歸模型.如下：

(3)

公式(3)中，yt為累計(jì)比例均值，K，L，M為估計(jì)參數(shù)，t為天數(shù).

由于烤煙收購受自然因素的影響，每年有不同的計(jì)劃烤煙收購天數(shù).按照麒麟?yún)^(qū)各收購站2020年的收購計(jì)劃天數(shù).用tii表示收購站2020年計(jì)劃收購天數(shù)，i代表收購站.如表1所示.

表1 麒麟?yún)^(qū)各收購站2020年計(jì)劃收購天數(shù)

2 結(jié)果與分析

通過上述數(shù)據(jù)和分析原理，以曲靖市麒麟?yún)^(qū)烤煙收購數(shù)據(jù)為例進(jìn)行建模分析，具體結(jié)果如下：

2.1 利用Logistic模型對(duì)收購站日累計(jì)收購比擬合

選擇2018,2019年度麒麟?yún)^(qū)各收購站分天數(shù)的烤煙收購數(shù)據(jù)，每個(gè)收購站當(dāng)日的收購量累計(jì)比例記作：ytiij，i代表收購站，i=1,…，7,j代表年度，j=2018,2019，ti表示收購站的收購天數(shù).根據(jù)樣本均值得到各收購站當(dāng)日累計(jì)比例均值：

(4)

由公式(4)得到麒麟?yún)^(qū)烤煙各收購站每日收購量的累計(jì)比例均值，對(duì)麒麟?yún)^(qū)烤煙各收購站每日收購量的累計(jì)比例均值做logistic非線性回歸函數(shù)擬合.

運(yùn)行SPSS 26.0的非線性回歸程序[6]，得到茨營收購站的輸出結(jié)果，表2和表3是關(guān)于茨營收購站累積比例的logistic非線性回歸函數(shù)的擬合結(jié)果.從回歸效果上看，擬合優(yōu)度R2=0.988，參數(shù)系數(shù)估計(jì)在95%的置信水平下均顯著.

表2 茨營收購站方差分析

表3 茨營收購站參數(shù)估算

把參數(shù)估計(jì)值代入(1)式，得到茨營收購站累計(jì)比例均值的預(yù)測方程：

為了保證累計(jì)比例最后一天是等于1，需要對(duì)參數(shù)K進(jìn)行修正，得到最終茨營收購站累計(jì)比例均值的預(yù)測方程：

圖2是茨營收購站2020年實(shí)際數(shù)據(jù)和預(yù)測方程擬合數(shù)據(jù)的對(duì)比圖.

同理可以得到其他6個(gè)收購站的累計(jì)比例均值的預(yù)測方程，方程系數(shù)如下表4所示.

表4 麒麟?yún)^(qū)各收購站累積比例均值的預(yù)測方程系數(shù)

圖2 茨營收購站累積比例均值原始值和預(yù)測值曲線

2.2 計(jì)算各收購站小等級(jí)日收購量占比

利用歷史收購數(shù)據(jù)，各收購站各小等級(jí)日收購量每天的占比，由于預(yù)測收購天數(shù)與歷史收購天數(shù)并不一致，這里利用下面公式(5)運(yùn)用一階樣條插值方法，對(duì)各小等級(jí)日收購量每天的占比做插值，使得該占比的天數(shù)與收購進(jìn)度預(yù)測的天數(shù)一致.對(duì)各小等級(jí)日收購量每天的占比插值后會(huì)出現(xiàn)每天各小等級(jí)占比之和不為1的情況，這里需要對(duì)其做修正(用各小等級(jí)占比除以當(dāng)天所有等級(jí)占比總和).

(5)

運(yùn)用上述方法，得到各收購站小等級(jí)分日預(yù)計(jì)收購占比.第一天茨營的小等級(jí)分日預(yù)計(jì)收購量占比數(shù)據(jù)如表5所示.

2.3 收購進(jìn)度預(yù)測應(yīng)用

結(jié)合收購站最新年度煙葉收購計(jì)劃、logistic回歸模型擬合得出的分日收購占比估計(jì)值和樣條插值法得出的小等級(jí)煙葉分日收購占比估計(jì)值，即可計(jì)算得出收購站各小等級(jí)煙葉的分日預(yù)計(jì)收購數(shù)量和累計(jì)預(yù)計(jì)收購數(shù)量.第一天茨營各小等級(jí)煙葉的分日預(yù)計(jì)收購數(shù)量數(shù)據(jù)如表6所示.

表5 麒麟?yún)^(qū)茨營收購站小等級(jí)分日預(yù)計(jì)收購量占比

表6 麒麟?yún)^(qū)茨營各小等級(jí)煙葉第一天的預(yù)計(jì)收購數(shù)量

各收購站小等級(jí)煙葉的分日預(yù)計(jì)收購數(shù)量和累計(jì)預(yù)計(jì)收購數(shù)量后，通過逐層匯總的方式即可計(jì)算得出麒麟?yún)^(qū)的分日預(yù)計(jì)收購數(shù)量和累計(jì)預(yù)計(jì)收購數(shù)量，見表7.

麒麟?yún)^(qū)分公司以下屬收購站每日的收購量為樣本，計(jì)算得出麒麟?yún)^(qū)分公司每日收購量的標(biāo)準(zhǔn)差，應(yīng)用置信區(qū)間估計(jì)法可計(jì)算得出各管理層級(jí)的累計(jì)收購進(jìn)度上下控制限，從而實(shí)現(xiàn)指導(dǎo)收購過程、跟蹤管控的目的.以2020年麒麟?yún)^(qū)煙草公司收購進(jìn)度預(yù)測為例，如圖3所示.

表7 麒麟?yún)^(qū)預(yù)計(jì)每天的收購數(shù)量

圖3 2020年麒麟?yún)^(qū)煙草公司收購進(jìn)度預(yù)測

在圖3中，實(shí)線為2020年實(shí)際收購進(jìn)度，虛線為2020年預(yù)計(jì)收購進(jìn)度，菱形線和方形線為預(yù)計(jì)收購進(jìn)度控制上下限.實(shí)線高于菱形線時(shí)表示實(shí)際收購進(jìn)度超出預(yù)計(jì)收購進(jìn)度上限，說明收購進(jìn)度存在非理性過快問題，需要調(diào)整烤煙收購進(jìn)度；同理虛線低于方形線表示實(shí)際收購進(jìn)度低于預(yù)計(jì)收購進(jìn)度下限，說明收購進(jìn)度極度過慢問題，需要調(diào)整收購管理要求(避免收購成本過大).可以看出總體范圍在可控之間，前 3 d 的實(shí)際收購進(jìn)度低于預(yù)期收購進(jìn)度下限，應(yīng)該調(diào)整收購管理要求.

2.4 收購均價(jià)預(yù)測應(yīng)用

根據(jù)上文的預(yù)測及小等級(jí)每日的收購量預(yù)測、結(jié)合2020年度煙葉收購的品種等級(jí)單價(jià)，可計(jì)算得出收購站每日的預(yù)計(jì)收購金額.預(yù)計(jì)收購金額除以預(yù)計(jì)收購量，即可得到收購站每日預(yù)計(jì)收購均價(jià).第1天各收購站預(yù)計(jì)收購均價(jià)數(shù)據(jù)如表8所示：

表8 麒麟?yún)^(qū)各收購站第一天預(yù)計(jì)收購均價(jià)

得出收購站分日預(yù)計(jì)收購均價(jià)和累計(jì)預(yù)計(jì)收購均價(jià)后,通過逐層匯總可計(jì)算得出麒麟?yún)^(qū)的分日預(yù)計(jì)收購均價(jià)和累計(jì)預(yù)計(jì)收購均價(jià).

麒麟?yún)^(qū)分公司以下屬收購站每日的收購均價(jià)為樣本，計(jì)算得出麒麟?yún)^(qū)分公司每日收購均價(jià)的標(biāo)準(zhǔn)差，應(yīng)用置信區(qū)間估計(jì)法可計(jì)算得出各管理層級(jí)的累計(jì)收購均價(jià)上下控制限，從而實(shí)現(xiàn)對(duì)收購均價(jià)的宏觀把握調(diào)控.

圖4為2020年麒麟?yún)^(qū)煙草公司收購均價(jià)預(yù)測，實(shí)線為2020年實(shí)際累計(jì)收購均價(jià)，虛線為2020年預(yù)測累計(jì)收購均價(jià)，菱形線和方形線為預(yù)計(jì)累計(jì)收購均價(jià)控制上下限.虛線高于菱形線表示實(shí)際累計(jì)收購均價(jià)超出預(yù)計(jì)累計(jì)收購均價(jià)上限，說明收購價(jià)格非理性虛高，收購質(zhì)量管控相對(duì)松懈，需要調(diào)整收購質(zhì)量控制眼光；同理虛線低于方形線時(shí)表示實(shí)際累計(jì)收購均價(jià)低于預(yù)計(jì)累計(jì)收購均價(jià)下限，說明收購價(jià)格明顯偏低，不利于煙農(nóng)增收指標(biāo)、收購稅利指標(biāo)的完成，需要調(diào)整收購管理要求.

圖4 2020年麒麟?yún)^(qū)煙草公司收購均價(jià)預(yù)測

從圖4中可以看出開收后的 4 d 之內(nèi)，均價(jià)都低于預(yù)測值，且這 4 d 的每日均價(jià)都在大幅度上漲.通過對(duì)收購數(shù)據(jù)的分析，造成這樣的原因是前期收購進(jìn)度過低.

3 結(jié)論

此次研究在大量調(diào)研的基礎(chǔ)上，依據(jù)麒麟?yún)^(qū)2年烤煙收購業(yè)務(wù)數(shù)據(jù)，運(yùn)用統(tǒng)計(jì)學(xué)模型，從宏觀著眼、以微觀入手，多層次、多維度對(duì)曲靖烤煙歷年收購趨勢及當(dāng)前理性狀態(tài)進(jìn)行深入透徹分析，有效地解決了烤煙收購政策把握和過程管理控制之間的矛盾問題.

在分析研究麒麟?yún)^(qū)歷年烤煙收購數(shù)據(jù)資料基礎(chǔ)上，加大研究力度，將統(tǒng)計(jì)學(xué)理論引入收購管理，能夠在確保市局(公司)年初制定的嚴(yán)控指標(biāo)不破及全市煙農(nóng)收入不減的前提下，理性推斷各等級(jí)預(yù)期收購?fù)瓿闪?，便于科學(xué)指導(dǎo)年度收購工作的開展，這對(duì)今后麒麟?yún)^(qū)烤煙收購精細(xì)化管理和烤煙生產(chǎn)的長遠(yuǎn)發(fā)展意義深遠(yuǎn).

本研究解決了烤煙生產(chǎn)經(jīng)營過程中較為重要的生產(chǎn)收購控制管理的突出問題，確保了年度收購目標(biāo)任務(wù)的圓滿實(shí)現(xiàn)，為科學(xué)、綜合、系統(tǒng)地提升曲靖烤煙生產(chǎn)收購管理能力產(chǎn)生了積極貢獻(xiàn).

煙草行業(yè)發(fā)展的基礎(chǔ)是烤煙，而烤煙生產(chǎn)價(jià)值的最終實(shí)現(xiàn)還將依托于收購，而收購過程中最關(guān)鍵的問題就是如何實(shí)現(xiàn)等級(jí)的合理把握和均價(jià)利稅的穩(wěn)步保增，巧妙運(yùn)用非線性規(guī)劃模型理論，創(chuàng)新應(yīng)用思路，在收購前期合理設(shè)定等級(jí)配比和預(yù)期均價(jià)目標(biāo)，并賦予其他方法有效指導(dǎo)、控制烤煙收購全過程，及時(shí)把握收購趨勢，及時(shí)調(diào)整等級(jí)偏誤.

將研究得到的烤煙收購進(jìn)程控制指標(biāo)體系用于實(shí)際指導(dǎo)烤煙收購過程，通過對(duì)收購進(jìn)程的控制，并定期分析控制情況，對(duì)超出控制限的要研究狀況，確?？緹熒a(chǎn)規(guī)模和經(jīng)濟(jì)效益的總體平衡，并力爭年度間實(shí)現(xiàn)煙農(nóng)戶煙葉收入穩(wěn)步增長，減少可能帶來的對(duì)烤煙生產(chǎn)未來可持續(xù)健康發(fā)展的不定性影響，為保障烤煙種植和后續(xù)收購管理工作提供科學(xué)依據(jù).