肾炎判断3

肾炎诊断 肾炎诊断 目要用黑体！！！ 【摘要】：肾炎是一种比较常见的疾病，早期症状常不明显，容易被忽略，发展到晚期可引起肾功能衰竭，严重威胁病人的健康和生命，是引起肾功能衰竭最常见的原因，必须早发现早治疗。本文就肾炎病人和健康人的化验结果进行处理，并采用 判别法，马氏距离判别法和支持向量机（ ）模型三种能够准确确诊就诊人员的判别方法，对30名就诊人员进行了确诊并对这三种方法的正确性进行了检验和比较。为了减少化验指标，通过函数对化验结果进行显著性检验，从中找出影响人们患肾炎的关键因素，并用上述三种方法中较为理想的支持向量机（ ）模型对30名就诊人员进行了确诊。最后我们对减少化验指标前后对30名就诊人员进行确诊的结果作了进一步的分析，并对三种方法进行了比较得出了很有实用价值的结论。在此基础上，我们对医生在临床上就肾炎病人确诊方法的科学性进行了进一步证实，并依据人体内各种元素的供需和影响人们患肾炎的关键因素对肾炎病人的饮食给出了合理的建议。摘要格式问题 【关键词】：肾炎诊断 ；判别函数 ；正确性检验；显著性检验。 1.问题重述 人们到医院就诊时，通常要化验一些指标来协助医生的诊断。诊断就诊人员是否患肾炎时通常要化验人体内各种元素含量。表B.1（见附表）是确诊病例的化验结果，其中1－30号病例是已经确诊为肾炎病人的化验结果；31－60号病例是已经确诊为健康人的结果。表B.2是就诊人员的化验结果。 问题一：根据表B.1中的数据，提出一种或多种简便的判别方法，判别属于患者或健康人的方法，并检验你提出方法的正确性。然后按照提出的方法，判断表B.2（附表）中的30名就诊人员的化验结果进行判别，判定他（她）们是肾炎病人还是健康人。 问题二：能否根据表B.1的数据特征，确定哪些指标是影响人们患肾炎的关键或主要因素，以便减少化验的指标。并根据得到的结果，重复判断表B.2（见附表）中的30名就诊人员是肾炎病人还是健康人。 问题三：对2和4的结果作进一步的分析。 2.基本假设 1.通过人体内的7种元素的含量可以判定是否患有肾炎；我们组忽略了 2.不考虑就诊人员的其它疾病因素以及身体素质的差异对各种元素含量的影响；我们组忽略了 3.人处于健康状态时，体内各种元素的百分比含量是相对稳定的。 4，各种元素的含量对肾炎的影响是独立的，不考虑功能的交互作用的影响； 我们组的logistic模型如果加入交叉变量，拟合效果是否会更好呢？？从医学角度来讲，各元素交互作用十分重要。可对此进行模型改进 5.医院所提供的确诊病例的化验结果和就诊人员的化验结果都是准确的； 6.诊断标准符合第二届全国肾脏病学术会议制定的诊断标准。 3.符号约定 ：各种化学元素的编号， ：肾炎病人的病例号， ；那这里的s又表示什么呢？？ ：肾炎病人的化验结果矩阵， 表示病例号为 的肾炎病人体内编号为 的元素的化验结果； ：健康人的病例号， ； ：健康人的化验结果矩阵， 表示病例号为 的健康人体内编号为 的元素的化验结果； ：确诊病人的病例号， ； ：确诊病人的化验结果矩阵， 表示病例号为 的确诊病人体内编号为 的元素的化验结果； ：判别函数； ：判别函数的系数矩阵；其实 4.问题分析 对题中数据进行初步分析，不难发现，无论肾炎患者还是健康人，就不同的人体内相同的元素而言，其含量都是随机的在某个值上下小范围内波动，并使得肾炎患者和健康人之间产生了覆盖现象，比喻10号患者和40号健康人体内的Zn元素相同，都是156，可见，单就一种元素分析是不行的；另一方面，就相同人体内不同的元素而言，由于各种元素同处于人体这个有机整体，它们相互依存，相互制约，即客观上存在一定的联系，却又有明显的不确定性，若能寻找隐藏在随机性后面的统计规律性，确定这种相关关系，问题就可获解。 由上述分析可知，本文的首要任务是对表B.1中的数据进行处理，由于数据量大，可知属于统计分析问题。处理的关键是对表B.1（肾炎病人，视为实验组；健康人，视为对照组）中的数据进行总体比较，从中寻找出一种或多种简便的判别方法并用此法对表B.2中的就诊人员进行确诊。在减少化验指标时，其关键又要确定哪些指标是影响人们患肾炎的关键或主要因素，进而必须确定两组（实验组和对照组）中各元素的显著性大小。最后对结果进行分析，必须能得出能指导实际的正确结论。 5.模型建立及求解（1—2问） 5.1. 模型-的建立及求解 方法 一： 判别法 5.1.1. 模型一的建立 我们先将表B1中的数据分为A和B（A代表健康，B代表患病）。对于p个判别指标，设判别函数y为： 则A组对应判别值为： B组对应判别值为： A，B两组判别值的代表为： 判别的标准就是 的值，为使A组，B组之间有明显的区分，即： 越大越好； 又 属于A组， 属于B组，为了使得它们的代表值能具有很好的代表性，必须使它们与其代表之间的差距即离差平方和： 越小越好。、 综上所述，得fisher最优判别准则为函数： 越大越好。从而解得最优判别函数的系数矩阵 函数 的极大值点，由微分学可知， 为方程组 的解。最终求得判别函数： 其中的临界值为 （1）. 时，若 ， 则判别对象属于A肾炎病人，若 ，则判别对象属于B健康人。 （2）. 时，若 ， 则判别对象属于B肾炎病人，若 ，则判别对象属于A健康人。 在Matlab中编程运行,求得判别函数： 5.1.2.模型一的正确性检验： 现在来对上述方法的正确性作必要的检验，即检验上面求出的最优判别函数是否有效。步骤如下： （1）​ 计算统计量： （2）​ 这里显著性水平 ，从F分布表中查出： （3） 检验结果： ，说明所作出的判别函数是有效的，从而可以用来作判别。 5.1.3.模型一的求解： 将表B1.中的确诊病人进行回代结果见附表A.1其中病例号为32，38，39，60的4个健康人被误判为肾炎病人，正确率为93.3%，再用此法对表B.2中的就诊病人进行判别，结果见附表A.2，患病人数为15人，其病例号为61、62、64、65、66、68、69、71、72、73、75、76、79、83、85。 5.2．模型二的建立： 方法 二：马氏距离判别法 5.2．1.马氏距离判别法模型的建立： 马氏距离判别法仍然是要将表B.1中的化验结果分为肾炎病人A组和健康人B组，化验结果 是一个 维的样品，这里定义 为样品c到A组，B组的距离，这样通过比较 两值的大小就可以判断样品c是属于A组还是B组； 样品c到A组和B组的马氏距离分别为： 其中， 分别为总体A和B的均值和协方差。 则判别函数 为： 当 时，样品c靠近A，c属于A；当 时，样品c靠近B，c属于B。 编程求得： 5.2．2.马氏距离判别法的求解及正确性检验： 将表B.1的确诊病人的化验结果代入马氏距离判别函数中得到结果同附表A.1，说明马氏距离判别法和fisher判别法正确性相当，而将B.2就诊病人的化验结果代入马氏距离判别函数中得到结果同附表A.2，即患病人数为15人。 方法 三：支持向量机（ ）模型 在这里，我们首先把问题限定在线性可分的情况下讨论，然后再进行优化。 模型思想如图1所示：图中实心点和空心点分别表示健康人和肾炎病人两类样本，H为把两类没有错误地分开的分类平面，它是由支持向量支撑起来的，H，H分别为过两类样本中离分类平面最近的点且平行于分类的平面，它们是由支持向量组成的平面，H和H之间的距离叫做两类样本的分类间距。如果我们能够找到这样的一个最优的分类平面：使得分类平面不但能将两类无错误地分开，而且要使得两类的分类间隔最大。 设线性可分样本集为 是类别标号（在本题中， 即 为7维空间实向量，健康人类别编号+1，肾炎病人类别编号为-1）。 维空间中线性判别函数方程的一般形式设为： 则分类平面H的方程为： 分类平面H即为参考平面，由于分类标号的设定，H，H两个分类边界面与分类面H是平行的，可以认为是由H面上下平移一个单位得到的，则H，H的方程为： 若H平面上下不断平移，它会经过两类所有的样本，在这所有平移的平面中，其平移距离满足： 考虑到 与 在正负号上的一致性，上式可以去掉 的绝对值，得: 其中，在分类边界面H，H上的样本，也即是支持向量，能使上式取等号，此式就是使得分类面对所有样本正确分类的约束条件。 分类间距D为分类边界面H，H之间的距离: 要使间隔D最大等价于 或 最小。问题归结为在上述（1）的约束条件下，求函数： 的最小值。可见是一个二元二次线性规划问题，可以定义如下拉格朗日函数： 其中， 为拉格朗日系数，我们的问题是对 和 求拉格朗日函数的极小值。 上述拉格朗日函数分别对 和 求偏微分并令它们等于0，问题最中转化为如下简单的对偶问题：在约束条件： 之下对 求解下列函数的最大值： 设所求的最优解为 ，则 即为最优分类平面H的法向量。 这是一个不等式约束下二次函数极值问题，存在唯一解。并且其解比满足（1）约束取等号的支持向量，它们通常只是全体样本中的很少一部分。 而 是分类的域值，可以由任意一个支持向量用（1）式（此时（1）式取等号）求得，也可以通过分别在健康人和肾炎病人中任意取一对支持向量取平均值求得 。 至此，求得了最优分类函数： 其中, 为符号函数， 时取-1； 时取0； 时取+1。 在Matlab中编程运行，并将确诊病人待入最优分类函数 结果见附表A.3： 从表中可知，A组误判0人，B组误判0人，其正确率为100%， 这个结果是很理想的。 现在运用最优分类函数 对表B2.中的就诊病人进行确诊结果见附表A.4. 结果共有15人患有肾炎，其判别结果与上述两种方法相符。 6.化验指标的减少（3—4问） 6.1．显著性检验： 将各组测得的结果求出 ，经统计学t检验，求得p值，所得数据列于附表A.5。 由附表A.6。知，Zn,Fe,Ca,Mg这几种元素肾炎病人与健康人差异非常显著（其 ），而其他元素差异则没有它们显著（其 ），因此我们有理由认为，Zn,Fe,Ca,Mg这几种元素是影响人们患肾炎的关键或主要因素，其它几种元素则为次要因素。 6.2．减少化验后的模型求解： 在减少指标后，上述三种模型的求解结果见附表A.6，附表A.7，这里马氏距离判别法和fisher判别法的结果仍然是相当的。 其中fisher判别法的判别函数为： 马氏距离判别法的判别函数为： 7.对上述结果的进一步分析 首先，对三种判别方法及其结果进行比较评价，将上述三种方法的判别结果利于表A.7对比可知，在小样本容量（指样本容量 的样本集合）下，马氏距离判别法和fisher判别法的结果仍然是相当的，其正确率也是相当的，相比较而言，支持向量机则能作到更准确的判别。三种方法中fisher判别法，特别是支持向量机，再涉及到高维数据时，都是借助向量这一概念来使问题简化，马氏距离判别法，运用距离这一几何长度来衡量则更易于理解接受。另外，在实际处理中，医生在诊断就诊病人是否患肾炎时，首先是尿常规检查，而化验各种元素含量只是为了协助医生的诊断，从这个意义上讲，上述三种方法在正确性上是能满足要求的。 再分析一下元素，由于Zn,Fe,Ca,Mg这几种元素是影响人们患肾炎的关键或主要因素，因此，肾炎病人必须多进食含这几种元素的食物才能满足机体的正常需求，比如多吃西瓜，牛奶，鸡蛋等等另外，肾炎的根本原因是肾脏病变，其滤过能力改变，不能有效的滤过尿液中的蛋白质等含化学元素的有机物，最终导致体内的化学元素流失，并使得尿液的化学成分，ph值，比重等发生改变，可见实际中诊断肾炎病人首选尿常规检查是科学的。 8．模型的推广 在社会、经济、自然科学领域的研究中，人们经常需要对某一研究对象的一些属性进行观察和研究从而将其进行适当的归类。例如，一个医生要对病人病情进行分析，以便判断到底应该使用何种手段治疗，例如非典型肺炎与典型肺炎治疗的方法就不相同。经济管理人员要对产品进行分类，判别它的销售情形属于“畅销”还是滞销。植物学家要对一类植物进行分析，昆虫学家对一种蠓虫的分类判别研究，地质学家要对某一地层出现的岩芯进行分类分析，在上述事例中经常用到的统计手段就是判别分析法。 判别分析有多种方法和准则，其中常用的判别分析法是主要有：距离判别法、Fisher判别法、Bayes判别法、逐步判别法。除此之外，还有一种重要的的方法就是支持向量机。支持向量机（SVM）是近几年发展的以统计学习理论为基础的新型学习机，它是根据Vapnic提出的统计学习理论发展而来的，该学习机以结构风险为最小化原则取代传统学习机所采用的经验风险最小化原则，使SVM有更强的泛化能力。 9．模型的优缺点 优点： 1．本文建立的统计识别模型能有效的解决实际生活中相类似的样本分类的问题，具有较强的实用性和通用性，而且误差也比较小； 2．本文采用的三种判别分析法成功地为诊断肾炎提供了一个科学、合理的方案，而且准确率达93.33％，具有较高的可信度； 3. 求解的过程中采用matlab数学软件编程求解，计算过程方便快捷，且结果准确。 缺点： 1.在所建立的模型中，对于题目中所提供的少部分就诊人员的化验结果我们并不能作出完全准确的判断；而且采用改进后的判别函数对就诊人员进行诊断时所得的结果中有几个与原先的判定结果不相同； 2.在肾炎病的诊断过程中如果能综合考虑其他因素对人体内元素含量的影响，将使本文所建的模型更加贴近于生活实践，更加具有实用价值；同时也能对就诊人员是否患病提供可靠的依据； 3. 本文所建的费希尔判定模型需要在总体之间差异较大的基础上进行，如果在总体之间差异不是非常明显时就会影响到最后的判定结果；此外，题中所提供的参考数据有限，这在一定程度上也降低了诊断结果的准确性。 十．参考文献及附录 参考文献：【1】宋来忠，王忠明，数学建模与实验，科学出版社，2005.8. 【2】蔡锁章，数学建模原理与方法，海洋出版社，2000.6. 【3】徐全智，杨晋浩，数学建模（第二版），电子科技大学出版社，2008.6. 附录1： 表B.1 确诊病例的化验结果（a1表） 病例号 Zn Cu Fe Ca Mg K Na 1 166 15.8 24.5 700 112 179 513 2 185 15.7 31.5 701 125 184 427 3 193 9.80 25.9 541 163 128 642 4 159 14.2 39.7 896 99.2 239 726 5 226 16.2 23.8 606 152 70.3 218 6 171 9.29 9.29 307 187 45.5 257 7 201 13.3 26.6 551 101 49.4 141 8 147 14.5 30.0 659 102 154 680 9 172 8.85 7.86 551 75.7 98.4 318 10 156 11.5 32.5 639 107 103 552 11 132 15.9 17.7 578 92.4 1314 1372 12 182 11.3 11.3 767 111 264 672 13 186 9.26 37.1 958 233 73.0 347 14 162 8.23 27.1 625 108 62.4 465 15 150 6.63 21.0 627 140 179 639 16 159 10.7 11.7 612 190 98.5 390 17 117 16.1 7.04 988 95.5 136 572 18 181 10.1 4.04 1437 184 101 542 19 146 20.7 23.8 1232 128 150 1092 20 42.3 10.3 9.70 629 93.7 439 888 21 28.2 12.4 53.1 370 44.1 454 852 22 154 13.8 53.3 621 105 160 723 23 179 12.2 17.9 1139 150 45.2 218 24 13.5 3.36 16.8 135 32.6 51.6 182 25 175 5.84 24.9 807 123 55.6 126 26 113 15.8 47.3 626 53.6 168 627 27 50.5 11.6 6.30 608 58.9 58.9 139 28 78.6 14.6 9.70 421 70.8 133 464 29 90.0 3.27 8.17 622 52.3 770 852 30 178 28.8 32.4 992 112 70.2 169 31 213 19.1 36.2 2220 249 40.0 168 32 170 13.9 29.8 1285 226 47.9 330 33 162 13.2 19.8 1521 166 36.2 133 34 203 13.0 90.8 1544 162 98.90 394 35 167 13.1 14.1 2278 212 46.3 134 36 164 12.9 18.6 2993 197 36.3 94.5 37 167 15.0 27.0 2056 260 64.6 237 38 158 14.4 37.0 1025 101 44.6 72.5 39 133 22.8 31.0 1633 401 180 899 40 156 135 322 6747 1090 228 810 41 169 8.00 308 1068 99.1 53.0 289 42 247 17.3 8.65 2554 241 77.9 373 43 166 8.10 62.8 1233 252 134 649 44 209 6.43 86.9 2157 288 74.0 219 45 182 6.49 61.7 3870 432 143 367 46 235 15.6 23.4 1806 166 68.8 188 47 173 19.1 17.0 2497 295 65.8 287 48 151 19.7 64.2 2031 403 182 874 49 191 65.4 35.0 5361 392 137 688 50 223 24.4 86.0 3603 353 97.7 479 51 221 20.1 155 3172 368 150 739 52 217 25.0 28.2 2343 373 110 494 53 164 22.2 35.5 2212 281 153 549 54 173 8.99 36.0 1624 216 103 257 55 202 18.6 17.7 3785 225 31.0 67.3 56 182 17.3 24.8 3073 246 50.7 109 57 211 24.0 17.0 3836 428 73.5 351 58 246 21.5 93.2 2112 354 71.7 195 59 164 16.1 38.0 2135 152 64.3 240 60 179 21.0 35.0 1560 226 47.9 330 表B.2 就诊人员的化验结果（a2表） 病例号 Zn Cu Fe Ca Mg K Na 61 58.2 5.42 29.7 323 138 179 513 62 106 1.87 40.5 542 177 184 427 63 152 0.80 12.5 1332 176 128 646 64 85.5 1.70 3.99 503 62.3 238 762.6 65 144 0.70 15.1 547 79.7 71.0 218.5 66 85.7 1.09 4.2 790 170 45.8 257.9 67 144 0.30 9.11 417 552 49.5 141.5 68 170 4.16 9.32 943 260 155 680.8 69 176 0.57 27.3 318 133 99.4 318.8 70 192 7.06 32.9 1969 343 103 553 71 188 8.28 22.6 1208 231 1314 1372 72 153 5.87 34.8 328 163 264 672.5 73 143 2.84 15.7 265 123 73.0 347.5 74 213 19.1 36.2 2220 249 62.0 465.8 75 192 20.1 23.8 1606 156 40.0 168 76 171 10.5 30.5 672 145 47.0 330.5 77 162 13.2 19.8 1521 166 36.2 133 78 203 13.0 90.8 1544 162 98.9 394.5 79 164 20.1 28.9 1062 161 47.3 134.5 80 167 13.1 14.1 2278 212 36.5 96.5 81 164 12.9 18.6 2993 197 65.5 237.8 82 167 15.0 27.0 2056 260 44.8 72.0 83 158 14.4 37.0 1025 101 180 899.5 84 133 22.8 31.3 1633 401 228 289 85 169 8.0 30.8 1068 99.1 53.0 817 86 247 17.3 8.65 2554 241 77.5 373.5 87 185 3.90 31.3 1211 190 134 649.8 88 209 6.43 86.9 2157 288 74.0 219.8 89 182 6.49 61.7 3870 432 143 367.5 90 235 15.6 23.4 1806 166 68.9 188 附表A.1确诊病人判别结果（ 判别法，马氏距离判别法） 编号（有病组） 判别值（y） 判别类型 正确性 编号（正常组） 判别值（y） 判别类型 正确性 1 -0.0811 有病组 正确 31 -0.1389 正常组 正确 2 -0.0012 有病组 正确 32 -0.0787 有病组 不正确 3 -0.0134 有病组 正确 33 -0.0923 正常组 正确 4 -0.0157 有病组 正确 34 -0.1048 正常组 正确 5 -0.0085 有病组 正确 35 -0.1475 正常组 正确 6 -0.0229 有病组 正确 36 -0.1935 正常组 正确 7 -0.0236 有病组 正确 37 -0.1347 正常组 正确 8 -0.0231 有病组 正确 38 -0.0537 有病组 不正确 9 0.0091 有病组 正确 39 -0.0683 有病组 不正确 10 -0.0147 有病组 正确 40 -0.2426 正常组 正确 11 -0.0098 有病组 正确 41 -0.1567 正常组 正确 12 0.0537 有病组 正确 42 -0.1451 正常组 正确 13 -0.0086 有病组 正确 43 -0.0927 正常组 正确 14 -0.0786 有病组 正确 44 -0.2006 正常组 正确 15 -0.0239 有病组 正确 45 -0.3172 正常组 正确 16 -0.0237 有病组 正确 46 -0.105 正常组 正确 17 -0.0307 有病组 正确 47 -0.149 正常组 正确 18 -0.0013 有病组 正确 48 -0.1183 正常组 正确 19 -0.0725 有病组 正确 49 -0.1701 正常组 正确 20 0.0177 有病组 正确 50 -0.2263 正常组 正确 21 0.0214 有病组 正确 51 -0.2247 正常组 正确 22 0.0385 有病组 正确 52 -0.1289 正常组 正确 23 0.0028 有病组 正确 53 -0.108 正常组 正确 24 -0.064 有病组 正确 54 -0.1225 正常组 正确 25 8.56E-04 有病组 正确 55 -0.2328 正常组 正确 26 -0.0679 有病组 正确 56 -0.1956 正常组 正确 27 0.0183 有病组 正确 57 -0.2407 正常组 正确 28 -0.0072 有病组 正确 58 -0.1649 正常组 正确 29 -0.0038 有病组 正确 59 -0.1177 正常组 正确 30 0.0018 有病组 正确 60 -0.0737 有病组 不正确 附表A.2就诊病人判别结果（ 判别法，马氏距离判别法） 编号 61 62 63 64 65 66 67 68 69 70 判别值 -0.0131 -0.0592 -0.0912 -0.000442 -0.0498 -0.071 -0.1394 -0.0703 -0.047 -0.1576 判别类型 有病组 有病组 正常 有病组 有病组 有病组 正常 有病组 有病组 正常 编号 71 72 73 74 75 76 77 78 79 80 判别值 -0.045 -0.0162 -0.0255 -0.1205 -0.073 -0.0364 -0.0923 -0.1048 -0.0421 -0.1497 判别类型 有病组 有病组 有病组 正常 有病组 有病组 正常 正常 有病组 正常 编号 81 82 83 84 85 86 87 88 89 90 判别值 -0.1849 -0.1448 -0.0035 -0.1075 -0.0301 -0.145 -0.0841 -0.2005 -0.3171 -0.105 判别类型 正常 正常 有病组 正常 有病组 正常 正常 正常 正常 正常 表A.3确诊病人判别结果 患病者 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 健康者 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 表A.4就诊病人判别结果 就诊人员 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 -1 -1 1 -1 -1 -1 1 -1 -1 1 1 -1 -1 1 1 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 -1 1 1 1 1 1 1 -1 1 -1 1 -1 1 1 1 附表A.5 检验结果 元素 有病 正常 P Zn 143.1±5.37 188.6±29.68 0.000334 Cu 12.33±5.05 21.92±23.75 0.038 Fe 23.07±13.75 62.01±75.78 0.0095 Ca 698.17±270.20 2511.13±1271.60 1.29691E-08 Mg 113.39±46.84 295.14±177.57 5.33158E-06 K 201.13±259.78 90.37±51.32 0.029 Na 526.83±300.58 367.21±244.24 0.279 附表A.6就诊病人判别结果（ 判别法，马氏距离判别法） 编号 61 62 63 64 65 66 67 68 69 70 判别值 -0.0258 -0.0482 -0.0907 -0.0431 -0.0632 -0.048 -0.0286 -0.0754 -0.0622 -0.1205 判别类型 有病组 有病组 有病组 有病组 有病组 有病组 有病组 有病组 有病组 正常 编号 71 72 73 74 75 76 77 78 79 80 判别值 -0.0947 -0.0545 -0.0492 -0.143 -0.1165 -0.0743 -0.1027 -0.1246 -0.0865 -0.1308 判别类型 有病组 有病组 有病组 正常 正常 有病组 有病组 正常 有病组 正常 编号 81 82 83 84 85 86 87 88 89 90 判别值 -0.1597 -0.1206 -0.0877 -0.0851 -0.0922 -0.1643 -0.0974 -0.1426 -0.1905 -0.1371 判别类型 正常 正常 有病组 有病组 有病组 正常 有病组 正常 正常 正常 附表A.7三种方法结果比较. 检验指标 判别法 7种微量元素做检验指标 Zn、Fe、Ca、Mg四种微量元素做检验指标 马氏判别法 15 19 费希尔判别法 15 19 支持向量机 15 17 附录2： 文件M1 w0=xlsread('a1');; w1=xlsread('a2'); wz0=[143.1033333;12.33433333; 23.06666667; 698.1666667; 113.3933333; 201.1333333 ;526.8333333]; wp0=[ wz0'; wz0' ;wz0'; wz0'; wz0'; wz0'; wz0'; wz0'; wz0' ; wz0'; wz0'; wz0'; wz0';wz0'; wz0'; wz0' ;wz0' ;wz0' ; wz0' ;wz0'; wz0'; wz0'; wz0'; wz0'; wz0'; wz0'; wz0'; wz0' ;wz0';wz0']; wz1=[186.6; 21.92366667; 62.01166667; 2511.133333; 295.1366667; 90.37; 367.21]; wp1=[ wz1'; wz1' ;wz1'; wz1'; wz1'; wz1'; wz1'; wz1'; wz1' ; wz1'; wz1'; wz1'; wz1';wz1'; wz1'; wz1' ;wz1' ;wz1' ; wz1' ;wz1'; wz1'; wz1'; wz1'; wz1'; wz1'; wz1'; wz1'; wz1' ;wz1';wz1']; a=w0-wp0; b=w1-wp1; s1=a'*a;s2=b'*b;s=s1+s2; q=inv(s)*(wz0-wz1); ya=sum(q.*wz0); yb=sum(q.*wz1); y=(ya+yb)/2 %判别值 A=q'; for n=1:30 y0=sum(A.*w0(n,:));%比较值 if y0>=y u0=1%判断正确 end if y0y u=0%判断正确 end end%以下开始第二问 wx=xlsread('a3'); for n=1:30 yx=sum(A.*wx(n,:));% if yx>=y uy=1%判断有病 end if yx=y u0=1%判断正确 end if y0y u=0%判断正确 end end%以下开始第四问 wx=xlsread('a6'); for n=1:30 yx=sum(A.*wx(n,:));% if yx>=y uy=1%判断有病 end if yx0 n fprintf('第一类判断正确') else if wY<0 n fprintf('第一类判断错误') else n fprintf('第一类无法判断') end end end for n=1:30 Y=w1(n,:)'; wY1=(Y-(wz0'+wz1')/2)'*inv(v)*(wz0'-wz1'); if wY1<0 n fprintf('第二类判断正确') else if wY1>0 n fprintf('第二类判断错误') else n fprintf('第一类无法判断') end end end for n=1:30 Y=wx(n,:)'; wYx=(Y-(wz0'+wz1')/2)'*inv(v)*(wz0'-wz1'); if wYx>0 n fprintf('待测组判断为有病') else if wYx<0 n fprintf('待测组判断为无病') else n fprintf('待测组无法判断') end end end 剔出Cu，K，Na 文件M4： syms x1 x2 x3 x4 ; X=[x1 x2 x3 x4 ]'; w0=xlsread('a4'); w1=xlsread('a5'); wx=xlsread('a6'); wz0=mean(w0); wz1=mean(w1); wp0=[ wz0; wz0 ;wz0; wz0; wz0; wz0; wz0; wz0; wz0 ; wz0; wz0; wz0; wz wz0; wz0 ;wz0 ;wz0; wz0;wz0; wz0; wz0; wz0; wz0; wz0; wz0; wz0; wz0 ;wz0;wz0]; wp1=[ wz1; wz1;wz1; wz1; wz1; wz1; wz1; wz1; wz1 ; wz1; wz1; wz1; wz wz1; wz1 ;wz1 ;wz1 ; wz1 ;wz1; wz1; wz1; wz1; wz1; wz1; wz1; wz wz1 ;wz1;wz1]; s0=(w0-wp0)'*(w0-wp0); s1=(w1-wp1)'*(w1-wp1); v=1/(60-2)*(s0+s1); w=(X-(wz0'+wz1')/2)'*inv(v)*(wz0'-wz1') for n=1:30 Y=w0(n,:)'; wY0=(Y-(wz0'+wz1')/2)'*inv(v)*(wz0'-wz1'); if wY0>0 n fprintf('0ÅжÏÕý³£') else if wY<0 n fprintf('0ÅжϴíÎó') else n fprintf('0ÎÞ·¨ÅжÏ') end end end for n=1:30 Y=w1(n,:)'; wY1=(Y-(wz0'+wz1')/2)'*inv(v)*(wz0'-wz1'); if wY1<0 n fprintf('1ÅжÏÕý³£') else if wY1>0 n fprintf('1ÅжϴíÎó') else n fprintf('1ÎÞ·¨ÅжÏ') end end end for n=1:30 Y=wx(n,:)'; wYx=(Y-(wz0'+wz1')/2)'*inv(v)*(wz0'-wz1'); if wYx>0 n fprintf('xÓв¡') else if wYx<0 n fprintf('xÕý³£') else n fprintf('xÎÞ·¨ÅжÏ') end end end 附录3： x1=xlsread('d4'); x=x1; %y1=-ones(1,30); %y2=ones(1,30); y=xlsread('b5'); for i=1:60 for j=1:60 Y(i,j)=y(i)*y(j); end end G=x*x'; K=G*Y; C=-ones(1,60); Aeq=y';b=0; lb=zeros(1,60); ub=-1.*C; [M,fval]=quadprog(K,C,[],[],Aeq,b,lb,ub); for i=1:60 p(i)=M(i)*y(i); end w=p*x; R=G(:,1); a=-1-p*R; W=[a,w]; f=W(1)+W(2)*x1+W(3)*x2+W(4)*x3+W(5)*x4+W(6)*x5+W(7)*x6+W(8)*x7; sign(f); x2=xlsread('b3'); for i=1:30 F=b+w*x2(i,:)'; p(i)=sign(F); end