企业同工同酬案例分析

《企业同工同酬案例分析》由会员分享，可在线阅读，更多相关《企业同工同酬案例分析（22页珍藏版）》请在金锄头文库上搜索。

1、1企业同工同酬案例分析09012320 茅光健09012326 杨俊09012333 桂滕华摘要：资制度的合理性与科学性是激励职工的劳动积极性，提高劳动效率的重要因素。因而通过建立数学模型，分析影响企业员工工资的因素，对于了解以及完善企业工资制度具有建设性意义。我们根据附表所给出的被解释变量平均日工资，以及解释变量性别、工龄等七个因素，进行合理的模型假设，并按照多元线性回归的方法建模。首先用Matlab 软件拟合出平均日工资随工龄的变化曲线，根据图形的变化趋势，选取函数 xy进行第一次建模。通过在交互式画面中对因素的逐条检验，判定工龄和受教育状况为显著性影响因素，并且利用 Matlab 统计功

2、能工具箱中的rcoplot(r,rint)命令剔除数据中的异常点，修正模型。最终得到平均日工资随工龄和受教育状况变化的关系式。 现代社会，女性的社会地位随着自我价值的实现而日益提高。在企业层面上，女性的自我价值是否得到与男性的同等评价仍是众多学者研究的论题。我们在分析了总体的工资模型后，又将数据按性别分为男、女两类，分别按照上述方法进行建模，通过对比进行分析了解是否存在“区别待遇”。以同样的方法再次对已婚与未婚女性进行建模分析，了解婚姻状况对平均日工资的影响，探讨差异存在是否合理。我们在第一次模型的基础上进行第二次修正，即通过剔除异常值，减少非2相关因素，对线性回归相关系数 R进行进一步的修正

3、，使得 R 向 1 趋近，使模型更加完善。最后，我们回归实际，从实际情况中分析模型的合理性，对模型进行第三次修正，即重新分配了受教育情况对工资影响的权重，使得模型更加精确更加符合实际。相信通过我们的研究分析，此次数学建模对于企业的工资制度具有参考性意义。【关键词】：多元线性回归 分类建模 Matlab 统计工具箱 一问题重述职工工资可以说是人们最为关切、议论最多的部分，因此也常常是最受重视的部分。一般说来，现代企业的工资具有补偿职能、激励职能、调节职能、效益职能。科学合理的工资制度，是激励职工的劳动积极性，提高劳动效率的重要手段，正确运用工资的杠杆作用在调动员工积极性方面会起到事半功倍的效果。

4、此外，对于企业中的各种不同的“特殊职务族”，是否要制定和执行专门的倾斜与优惠政策，如对管理干部、高级专家、女工等，也是需要重点考虑的问题。附件（Adata.xls）随机抽取了某企业若干职工的相关数据。请建立适当的数学模型研究下列问题：1.分析平均日工资与其他因素之间的关系，尤其需要说明与哪些因素关系密切；2.考察女工是否受到不公正待遇，以及她们的婚姻状况是否影响其收入；3.继续改进你的模型，并给出模型误差分析。3二问题分析本问题是关于企业职工同工同酬问题的研究分析问题一、我们以日平均工资为被解释变量，以员工性别、工龄、女性婚姻状况、员工受教育状况、工作部门性质、一线工作情况以及培训情况这些因素

5、为解释变量，分析研究出相应的数据模型。从数据模型中观察平均日工资与哪些因素关系密切。问题二、把总体数据分为男女两个类别，根据两个类别的数据分别研究出相应的数据模型，根据两个模型的差异分析，判断女工是否受到不公正待遇。然后把女性类别的数据按照已婚和未婚分为两类，分别研究出两个数据模型，根据它们的差异分析，判断女性婚姻状况是否影响其收入。问题三、通过观察由预测值和实际值得出的残差值，不断修正数据模型，使得相关系数趋向于 1。在已完善的数据模型的基础上，研究讨论误差因素，比如：异常值、无相关变量，从而进行误差分析。三模型假设1.在实际情况中，企业计算工资公式应该对于每位职工都一样，本模型就是为了求解

6、这个计算工资公式。2.在实际情况中，企业职工的工资应该与他的工龄成正比，即工龄越大，其相应的工资越多，而且可以确定工龄就是工资的主要影响因素。3.对于受教育程度的影响，按照实际经验，在其他同等的条件下，本科，硕士，博士，博士后享有的工资呈递增趋势。4.以上两个变量（工龄，受教育程度）可以看成模拟连续变量，而对于性别（男或女）、工作部门（管理或技术）、一线工作情况（有或没有）、培训4（有或没有）这些因素来说，它们不能随机取值，因此可以用 bool 值 0、1 来处理。 5.根据实际情况，工资计算公式应该是各个变量因素的多项式的线性相加的结果，具体模型公式见下面变量符号说明。四变量符号说明 1.

7、性别变量用 1x表示，其对应的 1xf表示该变量产生的工资 1y。2. 工龄变量用 2表示，其对应的 2f表示该变量产生的工资 2。3. 婚姻状况变量用 3x表示，其对应的 3xf表示该变量产生的工资 3y。4. 受教育状况变量用 4表示，其对应的 4f表示该变量产生的工资 4。5. 工作部门变量用 5x表示，其对应的 5xf表示该变量产生的工资 5y。6. 一线工作情况变量用 6表示，其对应的 6f表示该变量产生的工资 67. 培训情况变量用 7x表示，其对应的 7xf表示该变量产生的工资 7y。则职工工资 Y，其计算公式可以表示为： 7654321 yyyy五模型的建立1.首先用 Matl

8、ab 软件拟合出平均日工资 y随工龄 2x的变化曲线，如图所示：52.观察上图，选取函数 22xky进行建模（令 2x，从而转化为关于 2x的线性运算，为方便起见，以下均用 表示工龄变量）；由于其他几个自变量非随机连续变量，先假设符合线性关系，建立模型，分析并进行修正。而对于 2这一修正也符合实际，因为将 22xky之后，可以判断出方程的斜率越来越小，既满足上图，又符合实际，因为在实际中，随着工龄的增长，不可能工资永远持续增长下去，而且当工龄小的时候它的斜率反而大说明刚刚进入单位的年轻人的工资增长很快，这样就能充分调动年轻人的工作热情，为企业注入活力，这样的 22xky很符合实际，因此很有必要

9、地做这样的模型修正。 ”“,1),( 8x710 11 710 717210观 测 方 程 组这 称 为 模 型 的即 满 足 模 型 关 系那 么 每 组 观 测 值 都 应 该 组 观 测 值模 型 变 量 有因 此 如 果 我 们 所 考 察 的关 系 生 成 的 是 由 模 型 的 变 量机 项 。 因 为 变 量 数 据 都误 差 扰 动 综 合 影 响 的 随 则 中 反 映 各 种个 未 知 参 数是响 的 解 释 变 量 有 影是 七 个 对即 平 均 日 工 资 ，是 被 解 释 变 量其 中 本 变 量 关 系 为多 元 线 性 回 归 模 型 的 基nnnKii XY ，

10、nXY n，YxY LLLL6 XxxY xxxlXlxyY nnn nnniniin 710 7111 110111),( ,LLMM可 写 成 下 列 向 量 方 程那 么 上 述 观 测 方 程 组 就 和 矩 阵 记 法如 果 我 们 引 入 下 列 向 量六数据处理1 1x表示性别，男性用“0”表示，女性用“1”表示；2 2表示工龄变量；73 3x表示婚姻状况，男性用“0”表示，已婚女性用“1”表示，未婚 女性用“2”表示；4 4x表示受教育状况，本科用“0”表示，硕士用“1”，博士用“2” 表示，博士后用“3”表示；5 5x表示工作部门，技术岗位用“0”表示，管理岗位用“1”表示；

11、6 6表示一线工作情况，“0”表示没有参加，“1”表示参加；7 7x表示培训情况，“0”表示未受培训，“1”表示经受培训。8 实验数据组数 n=90；七模型求解1. 利用 Matlab 统计工具箱，在工作区间中导入矩阵 X 和 Y，输入命令：b,bint,r,rint,stats=regress(Y,X,0.05)其中，0.05 为 alpha，输出 b 为 的估计值 ，bint 为 b 的置信区间，r 为残差向量XY，rint 为 r 的置信区间，stats 为回归模型的检验统计量。8stats 有三个值：回归方程的决定系数 2R（ 是相关系数）；F统计量值；与 统计量值相关的概率 P；按照

12、上述步骤得到结果： b = bint =26.0903 19.4684 32.71220.2844 -6.6886 7.25732.3447 -5.4847 3.6484-0.9181 5.1955 13.50089.3482 -3.2266 3.55760.1655 -4.9466 2.7665-1.0901 -4.0906 9.35162.6305 -4.1254 6.9823stats =2RF统计量值 P0.8428 62.7812 09b bint386.8777178 344.5566879 429.19874772.979612886 -77.2875927883.2468185

13、60.887073296 0.754217792 1.0199288-8.385908617 -60.8947945444.1229773107.9310123 60.04799244 155.81403227.438329231 -31.5015540146.37821247-0.201107965 -44.2996389343.8974237.812507558 -69.6452523985.27026751stats=R2F 统计量 P0.7924 44.7022 02.为了检测 7 个变量是否对工资有显著性影响，我们用 stepwise 函数，在 Matlab 工作空间中输入命令：st

14、epwise(X,Y,8,0.05)其中 8 表示矩阵 X 的列数，0.05 为显著性水平 alpha 值；10得到交互式画面：其中蓝色表示较显著性因子，其他表示次要因子，在修正模型时可忽略。上图中的 3x即为表示工龄变量 2x， 5即为受教育情况变量 4x，而且可以从数据中观察出， 7x完全随着 4x的变化而变化， 4x为 0 则 7也为 0， 4x为非 0则 7也为 1，可见 7带来的影响完全可以用 来替代表示，因此在修正模型时只考虑这两个变量，则有11420xxY现在令 nxxlX421142),(M则可以得到 4203.在 Matlab 中利用 rcoplot 函数查找数据之中的异常点

15、，输入命令：rcoplot(r,rint)得到异常值分布图，如下所示：12其中红色表示异常值点，即第 43、47、52、60、61、67、90 组数据，故在进一步修正时可以忽略。将剩余 83 组数据重新带入，进行计算，在 Matlab 中输入命令：b,bint,r,rint,stats=regress(Y,X,0.05)则得到该回归方程的线性回归系数：13b= bint=370.2719 349.6862 390.85760.9284 0.8321 1.0246111.1440 91.2196 131.0683Stats=R2F 统计量 P0.9000 270.2000 0可以看出，修正后回归方程的决定系数 2R为 0.9000，相关系数 R值为0.9487，比修正前模型更准确。可以说明这样的修正是符合实际的，进而得到修正后回归方程为： x421.00.98437.21Y4.下面研究其他次要因素的影响：为分析性别对工资的影响，将数据按性别分为男、女两类，分别按照上述方法进行建模，通过对比进行分析。14通过 Matlab 计算，得到男性工资方程回归系数：b= bint=404.2430 350.3553