健康指數(shù)的數(shù)學(xué)模型

殘差圖5101520253035404550CaseNumber從殘差圖可以看出，除第十四個(gè)數(shù)據(jù)和第三十個(gè)數(shù)據(jù)之外外，其余數(shù)據(jù)的殘差離零點(diǎn)均較近且殘差的置信區(qū)間均包含零點(diǎn)，這說(shuō)明回歸模型y=-16.073+0.7194x能較好的符合原始數(shù)據(jù)而第二個(gè)數(shù)據(jù)可視為異常點(diǎn).x1=[17.7019.9017.1022.3017.0024.6026.4021.2021.70...24.8025.2031.8023.1028.0020.5013.3029.6023.40...20.1018.4012.2034.0016.5021.8018.2010.4011.00...13.9011.9014.1027.0016.1015.8026.5025.9023.60...26.6021.6024.4030.7018.4029.3059.0037.101.9013.90.23.5020.4025.6019.10];x2=[6.807.307.007.107.407.207.007.007.407.007.506.807.707.307.707.107.807.90...7.207.206.307.405.906.106.006.306.205.506.506.007.406.106.406.206.906.50...6.905.905.907.107.706.707.506.907.307.808.008.007.707.80];x3=[413.30396.40553.20585.00540.00522.50516.60565.00534.90..626.80434.80481.30609.50532.60471.50499.50599.30713.30..483.50392.00158.00600.00214.00227.00801.00713.00262.00..673.00445.00470.00113.00377.00128.00217.00499.00446.00..481.00267.0082.00396.00312.00273.001166.00131.0023.00...471.00330.00606.00579.00180.00];x4=[12.0012.2012.5012.0012.5012.2011.9012.3012.10...11.7012.3012.1012.1012.5012.5012.3012.2011.90...12.7012.6012.6013.2012.5011.8011.7012.4012.00...12.9013.5012.3013.5012.4013.6013.3013.2013.10...13.4013.6013.4014.3012.6012.4012.9012.6013.30...12.1013.3012.9013.2012.80];x5=[6.000.000.000.000.000.000.000.000.000.000.006.000.000.000.000.000.000.00...0.000.0042.000.0074.0062.0065.3039.3049.30100.0029.30...70.700.0061.3034.7054.000.0027.304.7076.0071.300.00...

0.0012.700.004.000.000.000.000.000.000.00];x6=[0.000.000.000.000.000.000.000.000.000.000.000.000.000.000.006.000.000.00...0.000.0035.800.0042.200.00100.00100.0090.90100.00100.00...100.000.0065.7067.200.004.1041.200.000.6056.700.0023.70...0.000.000.00100.0067.6014.6080.600.0067.60];x7=[0.000.000.000.000.001.000.000.000.000.000.0049.500.000.000.000.000.00...0.000.000.000.0054.100.000.0011.200.000.000.000.000.0042.300.00...0.006.3048.5032.4016.009.2042.1068.500.0023.80100.00...32.400.000.0029.3015.6021.101.50];x8=[0.000.000.000.000.0064.000.000.00100.000.000.0088.400.003.600.00...0.0031.800.000.000.000.00100.00100.00100.00x8=[0.000.000.000.000.0064.000.000.00100.000.000.0088.400.003.600.00...0.0031.800.000.000.000.00100.00100.00100.00100.000.000.000.00...100.0070.00];49.7095.00100.00100.000.00100.00100.0020.60100.00.x9=[241.30238.70248.60242.90201.7095.20357.30326.5094.40...283.90284.6074.80194.60220.50215.80149.10153.10242.50...187.20230.20147.2049.5032.9064.3027.10120.50216.30...70.5070.90145.6029.30112.9026.8025.8048.7028.60...47.4053.2026.0037.7075.0082.3028.6080.10129.50...53.6050.7065.8046.9057.10];y=[186.90184.90192.30188.20200.2094.30276.80252.4092.50...219.20219.7074.20173.20218.60191.90147.10151.90187.30...144.60177.60145.2048.8031.3063.2026.70107.40167.30...0.0070.40112.508.5012.2025.1023.8048.3028.0046.40...52.5025.4037.2074.2081.6028.2079.00127.7052.30...50.2065.3046.1056.50];0.000.00100.000.00100.00100.00100.00100.00100.00100.00100.00...a=ones(50,1);X=[a,(x1)',x2',x3',x4',x5',x6',x7',x8',x9',(x1.*x6.*x5)',(x2.*x3.*x4)'];Y=[186.90184.90192.30188.20200.2094.30276.80252.4092.50219.20219.7074.20173.20218.60191.90147.10151.90187.30144.60177.60145.2048.8031.3063.2026.70107.40167.300.0070.40112.508.5012.2025.1023.8048.3028.0046.4052.5025.4037.2074.2081.6028.2079.00127.7052.3050.2065.3046.1056.50];alpha=0.05;[b,bint,r,rint,stats]=regress(Y',X,alpha)rcoplot(r,rint)228.7312-0.3041-13.7966-0.1564-8.3400-0.34100.0461-0.15070.05060.7250-0.00020.0018bint=-131.