预测样品中说受体激动剂含量
[0123] 按照如下方法对实验组和对照组样品中β2受体激动剂的含量进行预测：软件是 SPSS22 · 0，方法为多元统计分析，回归方程(公式7): ρν? = -14 · 102 X xl+4 · 135 X χ2+9 · 903 X χ3+1.398 X χ4+2.477 X χ5-8.862 X χ6+12.486 X χ7+12.914 X χ8-8.496 X χ9-4.082 X χΙΟ +6.775Χχ11-12.455Χχ12-7.357Χχ13-0.457Χχ14-23.512Χχ15-18.022Χχ16+14.833 Χχ17+199.170
[0124] 其中，χ?代表ADCY3基因与内参基因体系的Ct值的差值；
[0125] x2代表ADRB2基因与内参基因体系的Ct值的差值；
[0126] x3代表ATP1A3基因与内参基因体系的Ct值的差值；
[0127] x4代表ATP2A3基因与内参基因体系的Ct值的差值；
[0128] x5代表1L1B基因与内参基因体系的Ct值的差值；
[0129] x6代表MYLK基因与内参基因体系的Ct值的差值；
[0130] x7代表PAKACB基因与内参基因体系的Ct值的差值；
[0131] x8代表PTH基因与内参基因体系的Ct值的差值；
[0132] x9代表TNF-α基因与内参基因体系的Ct值的差值；
[0133] χΙΟ代表CREM基因与内参基因体系的Ct值的差值；
[0134] xll代表FABP5基因与内参基因体系的Ct值的差值；
[0135] xl2代表F0X01基因与内参基因体系的Ct值的差值；
[0136] xl3代表HSL基因与内参基因体系的Ct值的差值；
[0137] xl4代表L0C102185492基因与内参基因体系的Ct值的差值；
[0138] xl 5代表MYH9基因与内参基因体系的Ct值的差值；
[0139] xl6代表NEB基因与内参基因体系的Ct值的差值；
[0140] xl7代表TOE4C基因与内参基因体系的Ct值的差值。
[0141]将以上各值带入所得pvl值为预测值，该值在分析过程中在"保存"选项中选中"预 测值"便可由软件自动给出。与表1对应的实验组和对照组的各山羊肌肉组织样品中β2受体 激动剂含量的预测值如表5所示。
[0142] 表5、样品中β2受体激动剂含量的预测值(g/109g组织）
[0143]
[0144]
[G i%」 壮：衣￡)τ，头版姐忏面屮忏面石孙仪取厄一世个|HJ tfu忏面乃4·日|HJ来兄多亡^女处 理三个山羊的样品，每行的实验组样品与对照组样品（即样品编号相同的实验组样品与对 照组样品）的区别仅在于:实验组样品来源的山羊经莱克多巴胺处理，对照组样品来源的山 羊未经莱克多巴胺处理，实验组样品来源的山羊与对照组样品来源的山羊的其他处理均相 同。
[0146] 实测值与预测值的相关性分析：
[0147] 以实验组样品中莱克多巴胺含量的实测值为因变量，以实验组样品对应17个基因 (ACt值为自变量)进行多元线性回归分析，根据模型得到的莱克多巴胺含量预测值与莱克 多巴胺含量的实测值建立相关关系，pvl = (avl-9.7186)/0.96374，R = 0.9008，pvGPavl* 别表示预测值和实测值。说明相关性较好。根据实测值与预测值的关系可以得出avl = 0.96374 Xpvl+9.7186(公式 8)。
[0148] 3.2莱克多巴胺含量的分析
[0149] 运用EXCEL和SPSS 22.0(IBM)对数据进行相关分析，主成分分析和聚类分析。
[0150] 对计算所得实验组和对照组两组样品的ACt值进行主成分(PCA)分析，结果如图1 所示。由图1可知，利用该组基因对两组样品区分较好。
[0151] 对计算所得实验组和对照组两组样品的ACt值进行聚类分析(HCA)，结果如图2所 示。由图2可知，两组样品基本上各聚为一类，与PCA结果基本一致。
[0152] 根据公式1和公式2确定实验组和对照组两组样品中是否含β2受体激动剂，分类结 果如表6所不。
[0153] 表6、分类结果b'c
[0154]
[0155] a.仅对分析中的案例进行交叉验证。在交叉验证中，每个案例都是按照从该案例 以外的所有其他案例派生的函数来分类的。
[0156] b.已对初始分组案例中的96.3 %个进行了正确分类。
[0157] c.已对交叉验证分组案例中的96.3 %个进行了正确分类。
[0158] 公式5: al = 14 · 924 X xl-6 · 362 X x2+5 · 226 X x3+7 · 777 X x4+2 · 034 X x5+14 · 974 X x6-7.971Xx7+10.516Xx8-4.019Xx9+18.687Xxl0-10.382Xxll+1.690Xxl2+3.882X xl3+5.632 X xl4-9.792 X xl5-18.867 X xl6-〇.499 X xl7-320.475 ；
[0159] 公式6: a2 = 16 · 229 X xl-3 · 494 X x2+6 · 767 X x3+3 · 229 X x4+2 · 444 X x5+13 · 809 X x6-9.465Xx7+12.347Xx8-5.194Xx9+19.640Xxl0-9.493Xxll+2.413Xxl2+3.471 X xl3+5.306Xxl4-10.147Xxl5-16.995Xxl6-0.182Xxl7-321.610；
[0160] 公式5和公式6中，xl-xl7分别代表ADCY3基因、ADRB2基因、ATP1A3基因、ATP2A3基 因、IL1B基因、MYLK基因、PAKACB基因、PTH基因、TNF-α基因、CREM基因、FABP5基因、F0X01基 因、HSL基因、L0C102185492基因、MYH9基因、NEB基因、PDE4C基因与内参基因体系的Ct值的 差值。
[0161] 根据公式5和公式6分别计算al和a2,按照如下方法确定山羊肌肉组织中是否含有 ft?受体激动剂:al 2 a2,山羊肌肉组织含有&受体激动剂，al<a2,动物组织不含&受体激动 剂。
[0162] 3、利用6个基因预测样品中说受体激动剂含量
[0163] 3.1利用6个基因预测样品中说受体激动剂含量
[0164] 按照如下方法对实验组和对照组样品中β2受体激动剂的含量进行预测：软件是 SPSS22.0，方法为多元统计分析，回归方程（&S3):pv2 = 1.284Xyl-6.926Xx2+9.304-17 · 360 X x3-18 · 642 X x4+4 · 348 X x5+209 · 326，其中，
[0165] yl代表ADRB2基因与内参基因体系的Ct值的差值；
[0166] y2代表ATP2A3基因与内参基因体系的Ct值的差值；
[0167] y3代表FABP5基因与内参基因体系的Ct值的差值；
[0168] y4代表F0X01基因与内参基因体系的Ct值的差值；
[0169] y5代表TOE4C基因与内参基因体系的Ct值的差值；
[0170] y6代表NEB基因与内参基因体系的Ct值的差值。
[0171] 将以上各值带入所得pv2值为预测值，该值在分析过程中在"保存"选项中选中"预 测值"便可由软件自动给出。与表1对应的实验组和对照组的各山羊肌肉组织样品中此受体 激动剂含量的预测值如表7所示。
[0172]表7、样品中β2受体激动剂含量的预测值(g/109g组织）
[0173]
[0174] 注:表5中，实验组样品中样品名称仅最后一位不同的样品为相同莱克多巴胺处理 三个山羊的样品，每行的实验组样品与对照组样品（即样品编号相同的实验组样品与对照 组样品）的区别仅在于:实验组样品来源的山羊经莱克多巴胺处理，对照组样品来源的山羊 未经莱克多巴胺处理，实验组样品来源的山羊与对照组样品来源的山羊的其他处理均相 同。
[0175] 实测值与预测值的相关性分析：
[0176] 以实验组样品中莱克多巴胺含量的实测值为因变量，以实验组样品对应6个基因 (ACt值为自变量)进行多元线性回归分析，根据模型得到的莱克多巴胺含量预测值与莱克 多巴胺含量的实测值建立相关关系，。￥2=(&￥2-17.771)/0.68747，1? = 0.8291，？￥2和3￥2分 别表示预测值和实测值。说明相关性较好。根据实测值与预测值的关系可以得出av2 = 0.68747 Xpv2+17.771(公式 4)。
[0177] 3.2莱克多巴胺含量的分析
[0178] 运用EXCEL和SPSS 22.0(IBM)对数据进行相关分析，主成分分析和聚类分析。
[0179] 对计算所得实验组和对照组两组样品的ACt值进行主成分(PCA)分析，结果如图3 所示。由图3可知，利用该组基因对两组样品区分较好。
[0180] 对计算所得实验组和对照组两组样品的ACt值进行聚类分析(HCA)，结果如图4所 示。由图4可知，两组样品基本上各聚为一类，与PCA结果基本一致。
[0181] 根据公式1和公式2确定实验组和对照组两组样品中是否含β2受体激动剂，分类结 果如表8所示。
[0182] 表8、分类结果b，c
[0183]
[0184] a.仅对分析中的案例进行交叉验证。在交叉验证中，每个案例都是按照从该案例 以外的所有其他案例派生的函数来分析的。
[0185] b.已对初始分组案例中的96.3 %个进行了正确分类。
[0186] c.已对交叉验证分组案例中的96.3 %个进行了正确分类。
[0187] 公式1:bl = 0 · 564 X y 1+13 · 288 Xy2-4 · 462 X y3+3 · 554 X y4-10 · 740 Xy5+3 · 724 X y6-101.111；
[0188] &S2:b2=1.603Xyl+9.127Xy2-2.904Xy3+4.800Xy4-9.195Xy5+4.047X y6_86·019;
[0189] 公式1和公式2中，yl-y6分别代表ADRB2基因、ATP2A3基因、FABP5基因、F0X01基因、 PDE4C基因、NEB基因与内参基因体系的Ct值的差值。
[0190] 根据公式1和公式2分别计算bl和b2,按照如下方法确定山羊肌肉组织中是否含有 ft?受体激动剂:bl 2 b2,山羊肌肉组织含有&受体激动剂，bl<b2,动物组织不含&受体激动 剂。
[0191]针对所有实验组样品建立判别模型，，对这6个基因的结果进行验证，判别分析结 果如表9所示，结果显示，直接利用这6个基因对样品中的β2受体激动剂分析是可行的。
[0192] 表9、分类结果b，c
[0193]
[0194] a.仅对分析中的案例进行交叉