企业同工同酬案例分析
09012320茅光健 09012326杨俊 09012333桂滕华
摘要：
资制度的合理性与科学性是激励职工的劳动积极性，提高劳动效率的重要
因素。因而通过建立数学模型，分析影响企业员工工资的因素，对于了解以及完善企业工资制度具有建设性意义。
我们根据附表所给出的被解释变量平均日工资，以及解释变量性别、工龄等
七个因素，进行合理的模型假设，并按照多元线性回归的方法建模。首先用Matlab软件拟合出平均日工资随工龄的变化曲线，根据图形的变化趋势，选取函数y?x进行第一次建模。通过在交互式画面中对因素的逐条检验，判定工龄和受教育状况为显著性影响因素，并且利用Matlab统计功能工具箱中的rcoplot(r,rint)命令剔除数据中的异常点，修正模型。最终得到平均日工资随工龄和受教育状况变化的关系式。
现代社会，女性的社会地位随着自我价值的实现而日益提高。在企业层面上，
女性的自我价值是否得到与男性的同等评价仍是众多学者研究的论题。我们在分析了总体的工资模型后，又将数据按性别分为男、女两类，分别按照上述方法进行建模，通过对比进行分析了解是否存在“区别待遇”。以同样的方法再次对已婚与未婚女性进行建模分析，了解婚姻状况对平均日工资的影响，探讨差异存在是否合理。
我们在第一次模型的基础上进行第二次修正，即通过剔除异常值，减少非相
关因素，对线性回归相关系数R进行进一步的修正，使得R向1趋近，使模型更加完善。最后，我们回归实际，从实际情况中分析模型的合理性，对模型进行第三次修正，即重新分配了受教育情况对工资影响的权重，使得模型更加精确更加符合实际。相信通过我们的研究分析，此次数学建模对于企业的工资制度具有
1
参考性意义。
【关键词】：多元线性回归 分类建模 Matlab统计工具箱
一．问题重述
职工工资可以说是人们最为关切、议论最多的部分，因此也常常是最受重视
的部分。一般说来，现代企业的工资具有补偿职能、激励职能、调节职能、效益职能。科学合理的工资制度，是激励职工的劳动积极性，提高劳动效率的重要 手段，正确运用工资的杠杆作用在调动员工积极性方面会起到事半功倍的效果。此外，对于企业中的各种不同的“特殊职务族”，是否要制定和执行专门的倾斜与优惠政策，如对管理干部、高级专家、女工等，也是需要重点考虑的问题。
附件（Adata.xls）随机抽取了某企业若干职工的相关数据。请建立适当的
数学模型研究下列问题：
1.分析平均日工资与其他因素之间的关系，尤其需要说明与哪些因素关系密
切；
2.考察女工是否受到不公正待遇，以及她们的婚姻状况是否影响其收入； 3.继续改进你的模型，并给出模型误差分析。
二．问题分析
本问题是关于企业职工同工同酬问题的研究分析
问题一、我们以日平均工资为被解释变量，以员工性别、工龄、女性婚姻状
况、员工受教育状况、工作部门性质、一线工作情况以及培训情况这些因素为解释变量，分析研究出相应的数据模型。从数据模型中观察平均日工资与哪些因素关系密切。
问题二、把总体数据分为男女两个类别，根据两个类别的数据分别研究出相
应的数据模型，根据两个模型的差异分析，判断女工是否受到不公正待遇。然后把女性类别的数据按照已婚和未婚分为两类，分别研究出两个数据模型，根据它们的差异分析，判断女性婚姻状况是否影响其收入。
2
问题三、通过观察由预测值和实际值得出的残差值，不断修正数据模型，使
得相关系数趋向于1。在已完善的数据模型的基础上，研究讨论误差因素，比如：异常值、无相关变量，从而进行误差分析。
三．模型假设
1.在实际情况中，企业计算工资公式应该对于每位职工都一样，本模型就是
为了求解这个计算工资公式。
2.在实际情况中，企业职工的工资应该与他的工龄成正比，即工龄越大，其
相应的工资越多，而且可以确定工龄就是工资的主要影响因素。
3.对于受教育程度的影响，按照实际经验，在其他同等的条件下，本科，硕
士，博士，博士后享有的工资呈递增趋势。
4.以上两个变量（工龄，受教育程度）可以看成模拟连续变量，而对于性别
（男或女）、工作部门（管理或技术）、一线工作情况（有或没有）、培训（有或没有）这些因素来说，它们不能随机取值，因此可以用bool值0、1来处理。
5.根据实际情况，工资计算公式应该是各个变量因素的多项式的线性相加的
结果，具体模型公式见下面变量符号说明。
四．变量符号说明
1. 性别变量用x1表示，其对应的f1?x1?表示该变量产生的工资y1。 2. 工龄变量用x2表示，其对应的f2?x2?表示该变量产生的工资
y2。
3. 婚姻状况变量用x3表示，其对应的f3?x3?表示该变量产生的工资y3。
y4。
4. 受教育状况变量用x4表示，其对应的f4?x4?表示该变量产生的工资5. 工作部门变量用x5表示，其对应的f5?x5?表示该变量产生的工资
y5。
y6
6. 一线工作情况变量用x6表示，其对应的f6?x6?表示该变量产生的工资
3
7. 培训情况变量用x7表示，其对应的f7?x7?表示该变量产生的工资y7。
则职工工资 Y，其计算公式可以表示为：
Y?y1?y2?y3?y4?y5?y6?y7
五．模型的建立
1.首先用Matlab软件拟合出平均日工资y随工龄x2的变化曲线，如图所
示
：
2.观察上图，选取函数y2?kx2进行建模（令x‘x2，从而转化为关2?‘于x‘2的线性运算，为方便起见，以下均用x2表示工龄变量）；由于其他几个自
变量非随机连续变量，先假设符合线性关系，建立模型，分析并进行修正。而对于x‘因为将y2?kx2之后，可以判断出方程的斜x2这一修正也符合实际，2?率越来越小，既满足上图，又符合实际，因为在实际中，随着工龄的增长，不可能工资永远持续增长下去，而且当工龄小的时候它的斜率反而大说明刚刚进入单位的年轻人的工资增长很快，这样就能充分调动年轻人的工作热情，为企业注入活力，这样的y2?kx2很符合实际，因此很有必要地做这样的模型修正。
4
本变量关系为 多元线性回归模型的基’Y????x??x01122????7x7??
其中，Y是被解释变量，即平均日工资，x1~x7是七个对Y有影
响的解释变量，?0，?1，?，?7是8个未知参数，?则中反映各种 误差扰动综合影响的随机项。因为变量数据都是由模型的变量模型变量有n组观测值 关系生成的，因此如果我们所考察的满足模型关系， (Yi,X1i,?XKi),i?1,?,n，那么每组观测值都应该即
Y1??0??1X11????7X71??1 ?? Yn??0??1X1n????7X7n??n这称为模型的“观测方程组”
如果我们引入下列向量和矩阵记法??0?xy1?i1??1???1???????,x????,l????,???1?,?????Y???????i????????????n????1??n?1?yn?????xin????n??1x11?x71??X?(l,x1,?xn)????????1x1n?x7n??那么上述观测方程组就可写成下列向量方程Y??0??1x1????7x7???X???
六．数据处理
[1] x1表示性别，男性用“0”表示，女性用“1”表示； [2] x2表示工龄变量；
[3] x3表示婚姻状况，男性用“0”表示，已婚女性用“1”表示，未婚
5
'