资源描述:
《数学模型之软件开发人员的薪金》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、软件开发人员的薪金问题1•基本模型薪金y与资历兀,管理责任兀2,教育程度兀3,之间的多元线性冋归模型为:y=a()+a內+a2x2+a3x3+a4x4+e其屮do,Qr…,04是待估计的回归系数,£是随机误差。利用MATLAB软件可以得到回归系数及其置信区间(置信水平。二0.05),以及检验统计量R,F,p的结果,具体算法如下,相应的结果见表一。1.1具体算法xl=[l11112222333344445556666788881010101011111212131314151616161720]1■fx2=[l01001000011101000010101101100011100101
2、0110000]1•fx3=[l000000100100100000110001010110000100100000101]'•fx4=[0001011001010000110000110100001011010011011010]•■9y=[13876116081870111283117672087211772105351219512313149752137119800114172026313231128841324513677159651236621352138392288416978148031740422184135481446715942231742378025410148
3、611688224170159902633017949256852783718838174831920719346]1;x0=ones(46,1);x=[x0xlx2x3x4];[b,bint,r,rint,stats]=regress(y,x,0.05);vpa(b,8)vpa(bint,8)vpa(stats,8)1・2运算结果112111032.73431
4、10258.255^11807.21312546.127624844486607.80673688253293624810517516.96074・2994.17834-3826.3273-2162.02945147.738
5、5・635.7184931.194412
6、3410.9567226425801057144.8...图1•基本模型运算结果表一Qo11033[1025811807]546[484608]6883[62487517]-2994[-3826-2162]148[-636931]/?2=0.9567F=226p=0参数参数估计值〔信区间2.进一步讨论为了寻找改进的方向,常用残差分析方法,我们将影响因素分成资历与管理-教育组合两类,管理-教育组合的定义如表二。表二组合123456管理010101教育112233下图2给出£与资历兀的关系,图3给出£与管理%?-教育兀3,尢组合间的关系。200
7、010000-1000-2000匚;+丰++打+十十十十十十+++++++101520200010000-1000-2000图2.基本模型£与兀
8、的关系图乳£与兀2一兀3'尢组合的关系附:残差图程序段:x5=[256354315324616533521634232611364413614364313I]1;[b,bint,r,rint,stats]=regress(y,x,0.05);subplot(2,2,1)plot(y)subplot(2,2,2)plot(x5‘r,1b+1)2.更好的模型增加兀2与上,血的交互项后,模型记作利用MATLAB软件求解后的结果如下:1112111
9、203.7538111044.259111363.24852496.86392485.6351508.092837047.999736841.05187254.94774-1726.50424-1938.9889-1514.01955■348.39255■545.2113-151.57386-3070.59626-3371.8328-2769.359571835.967671571.37042100.5648123410.99885544.799030047.0934图4.模型二运行结果表三参数估计值"置信区间.11204*[1104411363]^497^[486508]^(1^7
10、048卩[68417255]“■17273[1939-1514]^/»4*'•34%[-545-152]^“5』卜3372-2769]^畋31836+[15712101]^^2=0.9988Q5544J99p=0a作模型二的两个残差图(图5,图6)如下所示:-500500图5图6附:残差图程序段:xl=[l11112222333344445556666788881010101011111212131314151616161720]1■/x2=[l01001