化验结果诊断模型参考答案.docx

1、化验结果诊断模型参考答案化验结果诊断模型问题重述与分析人们到医院就诊时，通常要化验一些指标来协助医生的诊断。本题给出了人们是否患某种疾病时通常要化验的几种指标以及其检验值。表1是确诊病例的化验结果，其中130号病例是已经确诊为患该种疾病的化验结果；3160号病例是已经确诊为健康人的结果。表2是某些就诊人员的化验结果，但未确诊其是否患有该种疾病。根据已知数据，需要解答如下问题：1) 问题：根据表1中的数据，提出一种简便的判别方法，判别属于患者或健康人的方法，并检验你提出方法的正确性。分析：根据表1当中60个化验结果，将Zn、Cu、Fe、Ca、Mg、K、Na看成是七个指标，则前30个为该疾病患者的

2、指标值，后30个为健康人的指标值， 可以将这些数据进行标准化处理，再采用主成分分析方法，将多个指标转化为几个综合指标，当给定一个患者的各指标值时，可以算出各综合指标的得分，当这些得分满足一定条件时，如根据正负值可以判定为健康或疾病。2) 问题：按照（1）提出的方法，对表2中的15名就诊人员的化验结果进行判别，判定他们是患该种疾病的病人还是健康人。分析： 由（1）中已有的综合指标，根据给定的15名就诊人员的指标值计算出综合指标的得分，以此判断他们的健康状况。3) 问题：能否根据表1的数据特征，确定哪些指标是影响人们患该疾病的关键或主要因素，以便减少化验的指标。并根据你给出的结果，重复2的工作。分

3、析： 为了确定哪些指标是影响该疾病的主要因素，则需要确定出哪些因素在判别中起的权重最大，可以考虑采取回归模型，通过去除一些变量，然后比较各组的显著性与正确率，正确率最高的那组中的变量即为影响该疾病的主要因素。一、 模型假设1) 假设医院化验设备先进，化验过程科学可靠，化验结果真实可信，确诊情况（有病/健康人）符合实际。2) 在解决本题过程中，所有的化验结果只是针对该类疾病检验，并不考虑其他疾病的影响。3) 本文所建模型的检验结果只是作为医生为病人诊断的一个参考，医生为问诊人员作出最终判定还需考虑其他因素，但与本题求解无关。二、 符号说明三、 模型建立与求解（一） 问题一的求解： 模型一：1、

4、数据“标准化”题目已给出了60为确诊病例的化验结果以及诊断结果，但是60个病例中各元素的含量的呈无规律性。所以我们需要对原始数据进行处理，首先对其进行标准化分析：用向量 =表示每个就诊人员的化验结果，则 =表示第病人的化验结果。将每个指标“标准化”，即做如下变换：其中，标准化的数据见附录一。2、 主成分分析对标准化的数据运用SPSS软件进行主成分分析，结果如表1、表2：表1解释的总方差成份初始特征值提取平方和载入合计方差的 %累积 %合计方差的 %累积 %13.12944.70244.7023.12944.70244.70221.97328.19272.8941.97328.19272.894

5、3.72310.32783.2214.5708.14791.3685.2844.05295.4206.2042.91298.3327.1171.668100.000由表1可以看出，前两个主成分，的方差和占全部方差的比例为72.894%，我们就选取为第一主成分，为第二主成分，基本上保留了原来7个指标的信息，这样得到了2个新指标。SPSS软件得到的这成分系数矩阵如表2：表2：成份矩阵a成份12x1.453-.538x2.852.293x3.682.195x4.898-.051x5.941.094x6-.206.856x7-.005.904由表2得到前2个主成分，的线性组合为：=（4.1）3、 模型

6、验证将60个就诊人员的化验结果带入（4.1）式得到结果如表3，我们的判别标准为：第一主成分为正值表示健康，为负值表示患病。表3：病例号第一主成分第二主成分正误判标志（正=0；误=1）1-1.491420.2399702-1.14428-0.1898303-1.281340.0753804-1.496641.2768505-0.5735-1.8200906-1.58142-1.3215907-1.22966-1.9461608-1.761880.8638609-2.16782-0.99101010-1.700220.09606011-3.453138.3533012-1.790410.81834

7、013-0.42078-0.98245014-1.80823-0.49929015-2.032760.68222016-1.44607-0.50165017-2.014460.70872018-0.85694-0.3205019-1.061552.25198020-3.601883.81686021-3.60084.08631022-1.470191.01608023-0.89605-1.51287024-4.295350.05951025-1.54492-1.80111026-2.12941.19027-3.17737-0.35639028-2.924070.77313029-4.09657

8、4.48592030-0.27663-1.237290311.32902-1.89114032-0.22283-0.92755133-0.59422-1.621041340.49476-0.728640350.16255-1.646280360.60017-1.801560370.49891-1.12276038-1.10205-1.664681390.904982.1258604017.600424.835980411.86782-0.258560421.34229-1.57748043-0.188020.541531441.42776-1.53630452.81954-0.52070460

9、.40085-2.066580471.1396-0.992880481.622951.903360496.182691.173070503.73368-0.482520514.016640.792220522.24093-0.429760531.013090.4541054-0.23626-1.035621551.9522-2.265040561.389-1.78010573.45006-1.066840582.94115-1.719420590.06593-1.093790600.45845-0.900860由表3可以看出，前30个就诊人员的第一主成分均为负值，判定为患病，后30个就诊人员的

10、第一主成分大致上为正值，判定为健康，正确率为91.6667%。 （二） 问题二的求解由模型一得到前两个主成分的线性组合为：=将15名待诊人员的化验结果带入上式得：表4：病例号第一主成分第二主成分61-5.200578411.4897649562-3.70917326-1.1761822263-3.33417236-0.1852121964-2.46849776-1.17181594653.44950972-0.3016820166-0.010136984.8819197667-1.538760260.882981368-0.4575774-0.70249665690.7083415-1.078

11、35557703.5182755-0.25750581711.40770589-0.63732864722.04054597-1.29664076732.69524135-0.85109768743.11945951-1.0433454975-0.220201161.44700283用第一主成分来判定化验结果，由表4可知，15名待诊人员中有8名患有该疾病，7名健康。（三）模型一的改进：模型二：Logistic回归模型问题一的模型的正确率为91.6667%，因此考虑正确率更高的其他模型，且模型一中忽略了第二主成分的作用，故解释时有较大误差。以Y=0表示健康，Y=1表示不健康，考虑的因变量为一个二

12、元变量，且只取0与1两个值，因变量取1的概率为要研究的对象，且是的线性函数，故考虑采用Logistic线性回归模型。对附录一中的数据运用SPSS进行Logistic回归分析得表5：表5：方程中的变量BS.E,WalsdfSig.Exp (B)步骤 1ax1.48948.943.0001.9921.630x2.347276.987.0001.9991.415x3-1.479160.310.0001.993.228x4-.0887.972.0001.991.916x5.02162.318.00011.0001.021x6.234109.431.0001.9981.264x7.01532.972.0

13、0011.0001.016常量33.4707350.783.0001.9963.435E14由表5可以看出这7个变量都是显著的，因而最终的回归方程为：根据以上公式，我们可以将这个模型计算出来的应用于实际病例的判别。只要给出某一个受检者的化验结果，就能应用此计算公式算出其患病几率，我们以0.5为参照，当0.5时表示该受检者患病， 当0.5时表示该受检者健康。具体数据如下：病号x1x2x3x4x5x6x7是否健康判别结果116615.824.570011217951311218515.731.57011251844271131939.825.954116312864211415914.239.78

14、9699.223972611522616.223.860615270.32181161719.299.2930718745.525711720113.326.655110149.414111814714.5306591021546801191728.857.8655175.798.4318111015611.532.5639107103552111113215.917.757892.413141372111218211.311.376711126467211131869.2637.19582337334710.999999839141628.2327.162510862.44651115150

15、6.6321627140179639111615910.711.761219098.5390111711716.17.0498895.5136572111818110.14.041437184101542111914620.723.812321281501092112042.310.39.762993.7439888112128.212.453.137044.1454852112215413.853.3621105160723112317912.217.9113915045.221810.9999998192413.53.3616.813532.651.618210.9999999982517

16、55.8424.980712355.6126112611315.847.362653.616862710.9999999942750.511.66.360858.958.913910.999999932878.614.69.742170.81334641129903.278.1762252.3770852113017828.832.499211270.216910.9999999383121319.136.2222024940168003217013.929.8128522647.9330003316213.219.8152116636.213300342031390.8154416298.9

17、394003516713.114.1227821246.3134003616412.918.6299319736.394.500371671527205626064.6237003815814.437102510144.672.5003913322.831163340118089900401561353226747109022881000411698308106899.153289004224717.38.65255424177.937300431668.162.8123325213464900442096.4386.921572887421900451826.4961.73870432143

18、367004623515.623.4180616668.8188004717319.117249729565.8287004815119.764.22031403182874004919165.4355361392137688005022324.486360335397.7479005122120.1155317236815073900522172528.22343373110494005316422.235.5221228115354900541738.99361624216103257005520218.617.737852253167.3005618217.324.8307324650.

19、710900572112417383642873.5351005824621.593.2211235471.7195005916416.138213515264.324000601792135156022647.9330006185.51.73.9950362.3238762.61621440.715.154779.771218.516385.71.094.279017045.8257.90.999999987641760.5727.331813399.4318.81651927.0632.919693431035530661888.2822.61208231131413721671535.8

20、734.8328163264672.516817110.530.567214547330.516916213.219.8152116636.21330702031390.8154416298.9394.507116420.128.9106216147.3134.50.7126521647216713.114.1227821236.596.507316412.918.6299319765.5237.80741671527205626044.87207515814.4371025101180899.51由上表可以看出，改进后的模型正确率非常高。且15名就诊人员中有9人患病，6人健康。（四）问题三的

21、模型建立与求解 由模型二中的表5知，的各显著性水平分别为：X1X2X3X4X5X6X70.9920.9990.9930.9911.0000.9980.996因此我们将各种元素进行组合排除，将分别去掉，比较去掉后该回归的各个参量的值，以标准误差和正确率作为评判假设是否合理的依据。再将去除，重新建立回归模型，以相同标准进行评判。同理，分别从去除23个量作为一组，进行评判。最后依照相同的方法每两个、每三个进行分组，均依上述标准评判，具体分组不再赘述。下面仅就回归后正确率比较大的两组组给出分析结果。1、 去除三个变量，由SPSS软件得到的结果见附录二回归方程为：病号x1x2x3x4x5x6x7是否健康

22、判别结果116615.824.570011217951311.00 218515.731.570112518442711.00 31939.825.954116312864211.00 415914.239.789699.223972611.00 522616.223.860615270.321810.71 61719.299.2930718745.525710.87 720113.326.655110149.414110.95 814714.53065910215468011.00 91728.857.8655175.798.431811.00 1015611.532.563910710355211.00 1113215.917.757892.41314137211.00 1218211.311.376711126467211.00 131869.2637.19582337334710.04 141628.2327.162510862.446511.00 151506.632162714017963911.00 1615910.711.761219098.539010.98 1711716.17.0498895.513657211.00 181811