朱家明，曹綺琦，潘雪航，錢禮會(huì)，李春忠

(1.安徽財(cái)經(jīng)大學(xué)統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院，安徽 蚌埠 233030；2.安徽財(cái)經(jīng)大學(xué)財(cái)政與公共管理學(xué)院，安徽 蚌埠 233030)

根據(jù)中國(guó)互聯(lián)網(wǎng)絡(luò)信息中心發(fā)布的第41次《中國(guó)互聯(lián)網(wǎng)絡(luò)發(fā)展?fàn)顩r統(tǒng)計(jì)報(bào)告》顯示，截至2017年12月，我國(guó)網(wǎng)民規(guī)模達(dá)7.72億，普及率達(dá)55.8%[1]。隨著時(shí)代進(jìn)步，越來(lái)越多的企業(yè)開(kāi)始通過(guò)互聯(lián)網(wǎng)尋求外部資源，一種新的非常有創(chuàng)意的發(fā)展模式正悄然興起，即眾包?！氨姲币辉~于2006年首次提出，指一個(gè)公司或機(jī)構(gòu)把過(guò)去由員工執(zhí)行的工作任務(wù)以自由自愿的形式外包給非特定的(而且通常是大型的)大眾網(wǎng)絡(luò)的做法[2]。而對(duì)于眾包任務(wù)的定價(jià)問(wèn)題，不同的學(xué)者有不同策略。徐芹[3]在壟斷型市場(chǎng)下研究相應(yīng)平臺(tái)對(duì)雙邊用戶的定價(jià)策略，但相關(guān)影響因素的考慮不夠全面；孫信昕[4]則認(rèn)為應(yīng)該根據(jù)每個(gè)工作者提交的成本、希望的任務(wù)數(shù)量以及完成質(zhì)量對(duì)任務(wù)進(jìn)行定價(jià)，但其模型適用范圍小，僅適用于任務(wù)質(zhì)量易測(cè)量的簡(jiǎn)單任務(wù)。為此，本文采用K-means聚類算法定量分析域內(nèi)會(huì)員數(shù)量、任務(wù)點(diǎn)距離、任務(wù)定價(jià)與任務(wù)完成情況之間的相關(guān)性，并為任務(wù)進(jìn)行合理定價(jià)。

1 數(shù)據(jù)的獲取及假設(shè)

數(shù)據(jù)來(lái)源于2017年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽B題。為便于解決問(wèn)題，提出以下研究條件：(1)用戶選擇任務(wù)只與任務(wù)價(jià)格和自己與任務(wù)地點(diǎn)的距離有關(guān)，沒(méi)有其他因素的影響；(2)用戶選擇任務(wù)會(huì)權(quán)衡任務(wù)價(jià)值和自己完成任務(wù)的成本，如果成本高于價(jià)值，則用戶不會(huì)選擇該任務(wù)；(3)APP中的任務(wù)定價(jià)在用戶預(yù)訂前就已確定，與任務(wù)周圍會(huì)員的數(shù)量和會(huì)員距離的遠(yuǎn)近無(wú)關(guān)；(4)未完成任務(wù)的地理位置只考慮廣州市和深圳市(因?yàn)檫@兩個(gè)城市未完成任務(wù)數(shù)量比重最大)，而不考慮其他城市；(5)任務(wù)打包發(fā)布時(shí)，會(huì)員接單會(huì)完成所接的所有任務(wù)。

2 任務(wù)定價(jià)影響因素的相關(guān)性分析

2.1 研究思路

任務(wù)定價(jià)在執(zhí)行之前，故任務(wù)的定價(jià)決定執(zhí)行情況，執(zhí)行情況對(duì)定價(jià)無(wú)影響。根據(jù)已知數(shù)據(jù)，可對(duì)任務(wù)的位置進(jìn)行聚類分析，從任務(wù)與中心點(diǎn)間的距離及任務(wù)的位置等任務(wù)定價(jià)影響因素探究定價(jià)規(guī)律，然后根據(jù)計(jì)算得出各任務(wù)區(qū)中的會(huì)員數(shù)量，與各任務(wù)區(qū)中的價(jià)格進(jìn)行比較，分析定價(jià)規(guī)律[5]。

2.2 研究方法：K-means均值聚類[6]

若將n個(gè)樣品分成k類，則先選擇所有樣品中距離最遠(yuǎn)的兩個(gè)樣品xi1、xi2為聚點(diǎn)，使得

d(xi1,xj2)=di1j2=max{dij}。

(1)

然后選擇第3個(gè)聚點(diǎn)xi3，使得xi3與前兩個(gè)聚點(diǎn)的距離最小者等于所有其余的與xi1、xi2的較小距離中最大的，即

min{d(xi3,xir),r=1,2}=max{min{d(xj,xir),r=1,2},j≠i1,i2}。

(2)

最后按相同的原則選取xik，重復(fù)前面的步驟，直至確定k個(gè)聚點(diǎn)xi1,xi2,…,xik。

k-means均值聚類的步驟(樣品之間的距離采用歐氏距離)如下：

設(shè)第k個(gè)初始聚點(diǎn)的集合是

(3)

記

(4)

于是，將樣品分成不相交的k類，得到一個(gè)初始分類

(5)

從初始類G(0)開(kāi)始計(jì)算新的聚點(diǎn)集合L(1)，計(jì)算

(6)

(7)

從L(1)開(kāi)始再進(jìn)行分類，記

(8)

得到一個(gè)新的類

(9)

重復(fù)上述步驟m次，得

(10)

(11)

與

(12)

相同，則結(jié)束計(jì)算。

變量x與變量y間的相關(guān)系數(shù)

(13)

2.3 研究結(jié)果

圖1 任務(wù)結(jié)束地區(qū)位置聚類圖Fig. 1 Map of completed task location cluster

運(yùn)用MATLAB軟件編程，得到已結(jié)束任務(wù)經(jīng)度緯度的位置聚類圖如圖1。

為了研究任務(wù)點(diǎn)位置對(duì)定價(jià)的影響，利用已結(jié)束任務(wù)位置聚成4類的結(jié)果，運(yùn)用EXCEL篩選，分別計(jì)算佛山市、深圳市、廣州市和東莞市任務(wù)定價(jià)的均值及其區(qū)域內(nèi)會(huì)員數(shù)量。結(jié)果見(jiàn)表1。

以4個(gè)聚類中心的位置作為定點(diǎn)，分別計(jì)算4個(gè)區(qū)各任務(wù)點(diǎn)到其中心的距離(單位：m)，分析任務(wù)價(jià)格與任務(wù)距離之間的相關(guān)關(guān)系，結(jié)果見(jiàn)表1。

表1 各地區(qū)任務(wù)價(jià)格和距離的相關(guān)系數(shù)、任務(wù)均價(jià)和會(huì)員數(shù)量Tab. 1 Correlation coefficient of task price and distance,mean task price and number of members by region

從表1可以看出各地區(qū)的任務(wù)價(jià)格和任務(wù)距離存在一定的相關(guān)關(guān)系，并且是正的相關(guān)關(guān)系。也就是說(shuō)各任務(wù)的定價(jià)和距離有關(guān)，任務(wù)點(diǎn)到聚類中心的距離越大，則任務(wù)的定價(jià)越高。

由表1求出各地區(qū)任務(wù)定價(jià)的均值與會(huì)員數(shù)量間的相關(guān)系數(shù)為-0.995 7，表明這兩者顯著相關(guān)，且地區(qū)的會(huì)員數(shù)量越多，其任務(wù)定價(jià)越低。

3 基于最小二乘法對(duì)未完成任務(wù)進(jìn)行重新定價(jià)

3.1 未完成任務(wù)的原因分析

圖2 未完成任務(wù)位置聚類圖Fig. 2 Location clustering of unfinished tasks

3.1.1 研究思路

對(duì)于任務(wù)未完成的原因，可以先利用MATLAB對(duì)未完成任務(wù)位置進(jìn)行可視化處理，得出位置的分布特征，再進(jìn)一步分析造成這種分布可能的原因，然后可以從區(qū)域內(nèi)會(huì)員數(shù)量和任務(wù)價(jià)格等方面探究這些因素是否會(huì)影響任務(wù)完成情況[7-9]。

3.1.2 研究結(jié)果

根據(jù)散點(diǎn)圖可知未完成情況有明顯聚集現(xiàn)象，運(yùn)用MATLAB繪出了未完成任務(wù)位置的聚類圖，結(jié)果如圖2所示。

圖2中聚類中心的位置坐標(biāo)分別是(22.64°N，114.03°E)和(23.12°N，113.25°E)，在地圖上確定分別是深圳市龍華區(qū)和廣州市越秀區(qū)。由于這兩個(gè)地區(qū)內(nèi)會(huì)員數(shù)量較大，經(jīng)濟(jì)也發(fā)達(dá)，所以不存在會(huì)員不夠而任務(wù)過(guò)多的情況。綜合考慮各種影響因素，認(rèn)為價(jià)格的高低直接影響任務(wù)完成情況，并且價(jià)格越高，任務(wù)完成情況越好，反之，則越差。為此，利用EXCEL軟件，計(jì)算每個(gè)價(jià)格的任務(wù)完成率。

在統(tǒng)計(jì)學(xué)中，樣本數(shù)量在30以下的是小樣本， 30以上的是大樣本。小樣本由于樣本數(shù)量過(guò)少，不具有代表性，所以剔除了個(gè)別小樣本數(shù)據(jù)，保留任務(wù)完成數(shù)在30以上的樣本。各價(jià)格下任務(wù)完成情況的統(tǒng)計(jì)結(jié)果如表2所示。

表2 各價(jià)格下任務(wù)完成情況Tab.2 Tasks completed by price

價(jià)格和完成率的相關(guān)系數(shù)高達(dá)0.683，存在明顯的正相關(guān)關(guān)系，擬合方程為y=-0.027x5+9.611x4-1 356x3+95 713x2-3×106x+7×107。擬合曲線見(jiàn)圖3。

圖3 各價(jià)格下任務(wù)完成情況的擬合曲線Fig. 3 Curve of tasks completed by price

圖3中的曲線擬合的可決系數(shù)為0.8544，表明擬合程度較高。由點(diǎn)的分布規(guī)律可知定價(jià)低的任務(wù)完成情況明顯比定價(jià)高的任務(wù)完成情況差，所以價(jià)格的高低能直接影響任務(wù)的完成情況。而任務(wù)的定價(jià)和完成率不存在明顯的線性關(guān)系，可能是因?yàn)橥瓿扇蝿?wù)的難易程度不一樣，包括交通的便利性、完成成本等方面因素，難度大的任務(wù)，定價(jià)再高也必然存在不能完成的現(xiàn)象。

3.2 未完成任務(wù)的重新定價(jià)

3.2.1 研究思路

首先，對(duì)已知任務(wù)的完成情況進(jìn)行篩選，得出已完成任務(wù)與未完成任務(wù)各自的分布情況；然后對(duì)價(jià)格進(jìn)行排序，得到已完成任務(wù)與未完成任務(wù)的價(jià)格走向，已完成任務(wù)的價(jià)格曲線即為任務(wù)定價(jià)模型；最后由模型計(jì)算出未完成任務(wù)的新價(jià)格，再與原方案的未完成任務(wù)的價(jià)格進(jìn)行比較得出結(jié)論。

3.2.2 研究方法——最小二乘法

由給定的一組測(cè)定的離散數(shù)據(jù)(xi,yi)(i=1,2,…,N)，求自變量x和因變量y的近似表達(dá)式y(tǒng)=φ(x)的方法稱為數(shù)據(jù)擬合，φ(x)稱為擬合函數(shù)。記擬合函數(shù)在xi點(diǎn)的偏差為δi=φ(xi)-yi(i=1,2,…,N)，那么選取φ(x)使所有偏差的平方和最小，即

(14)

這一確定擬合函數(shù)的方法稱為最小二乘法[10]。

3.2.3 研究結(jié)果

由之前的分析得知已完成任務(wù)總數(shù)為522個(gè)。4個(gè)地區(qū)未完成任務(wù)的總和為313個(gè)，其中252個(gè)分布在廣州市和深圳市，因此主要針對(duì)廣州、深圳兩市的情況進(jìn)行分析。

在深圳和廣州兩市，未完成任務(wù)為252個(gè)，已完成任務(wù)為229個(gè)，對(duì)兩部分任務(wù)的價(jià)格進(jìn)行升序排列，擬合曲線方程為

y=1.6×10-6x3-1.4×10-4x2+1.2×10-2x+65，

(15)

擬合曲線如圖4。

可以算出可決系數(shù)：

(16)

由可決系數(shù)來(lái)看擬合的效果比較理想。

由于已完成任務(wù)的價(jià)格視為不變，因此可以不用考慮重新調(diào)整它們的價(jià)格。經(jīng)計(jì)算，已完成任務(wù)價(jià)格均值為69.89元，未完成任務(wù)均值為67.93元，因此可以初步判斷未完成任務(wù)價(jià)格偏低。若要改進(jìn)價(jià)格方案，則需適當(dāng)提高其價(jià)格。

從圖4中可以看出，在后半部分已完成任務(wù)的價(jià)格要高于未完成任務(wù)，而已完成任務(wù)價(jià)格不需調(diào)整。因此，以已完成任務(wù)價(jià)格為標(biāo)準(zhǔn)，將擬合曲線方程(15)作為新方案，計(jì)算出未完成任務(wù)的適當(dāng)價(jià)格，將新方案計(jì)算出的價(jià)格和原方案的價(jià)格進(jìn)行了比較，對(duì)比結(jié)果見(jiàn)圖5所示。其中現(xiàn)方案代表計(jì)算給出的任務(wù)價(jià)格。

圖4 任務(wù)價(jià)格曲線圖Fig.4 Task price curve

圖5 兩種方案價(jià)格對(duì)比圖Fig.5 Two pricing options compared

由此，可以看出現(xiàn)方案的價(jià)格相較于原方案的價(jià)格有了一定程度的提高，其原因在于任務(wù)價(jià)格的提高對(duì)于會(huì)員完成任務(wù)的積極性起了很大作用，所以當(dāng)價(jià)格上升，會(huì)員完成任務(wù)的情況也會(huì)提高。

同時(shí)，該任務(wù)價(jià)格的均值為69.10元，相較于原來(lái)的67.93元有適當(dāng)提高，因此推斷，在此種方案下，任務(wù)的完成比例將會(huì)增加。

4 結(jié)語(yǔ)

針對(duì)“拍照賺錢”項(xiàng)目的任務(wù)定價(jià)，基于K-means均值聚類算法分析區(qū)域內(nèi)會(huì)員數(shù)量、任務(wù)點(diǎn)距離、任務(wù)定價(jià)與任務(wù)完成情況之間的相關(guān)性。結(jié)果表明：區(qū)域內(nèi)會(huì)員數(shù)量一定時(shí)，完成任務(wù)的距離越遠(yuǎn)，任務(wù)定價(jià)越高；完成任務(wù)的距離一定時(shí)，區(qū)域內(nèi)會(huì)員數(shù)量越小，任務(wù)定價(jià)越高；在區(qū)域內(nèi)會(huì)員數(shù)量足夠大的前提下，任務(wù)定價(jià)越高任務(wù)完成度越高?；诖诉M(jìn)一步深入分析未完成任務(wù)，利用最小二乘法擬合出新的定價(jià)模型。計(jì)算過(guò)程中綜合利用了MATLAB和EXCEL軟件，使計(jì)算更加簡(jiǎn)便準(zhǔn)確。K-means均值聚類算法不僅可以解決眾包服務(wù)定價(jià)問(wèn)題，而且在商業(yè)、生物、保險(xiǎn)行業(yè)等諸多領(lǐng)域都有廣泛的運(yùn)用，具有一定推廣意義。