代谢组学数据分析的统计学方法
更新时间:2023-07-21 06:41:01 阅读量: 实用文档 文档下载
20144312·357·
代谢组学数据分析的统计学方法
柯朝甫
1
*
张涛
2
武晓岩
1
李康
1Δ
代谢组学是近年发展快速的一门学科,目前在医学、植物学、微生物学、毒理学、药物研发等诸多领域中
[1-5]
。代谢组学研究产生大量的数得到了广泛的应用
据,这些数据具有高维、小样本、高噪声等复杂特征。
(4)相互作用关系复杂:各种代谢物质可能不仅
具有简单的相加效应,而且可能具有交互作用,从而增加了识别这些具有复杂关系的生物标志物的难度。
(5)相关性和冗余性:各种代谢物并非独立存在,而是相互之间具有不同程度的相关性,同时由于碎片、加合物和同位素的存在使得数据结构存在很大的冗余性,这就需要采用合理的统计分析策略来揭示隐藏其中的复杂数据关系。
(6)分布的不规则和稀疏性:代谢组学数据分布不规则,而且数据具有稀疏性(即有很多值为零),因此,传统的一些线性和参数分析方法此时可能失效。
数据的预处理
代谢组学数据分析的目的是希望从中挖掘出生物相关信息,然而,代谢组学数据的变异来源很多,不仅包括生物变异,还包括环境影响和操作性误差等方面。处理手段主要包括归一化(standardization)、标准化(normalization),即中心化(centering)和尺度化(scal-[10]
。归一化是针ing),以及数据转换(transformation)对样品的操作,由于生物个体间较大的代谢物浓度差
异或样品采集过程中的差异(如取不同时间的尿样),一般使用代谢物的相为了消除或减轻这种不均一性,
对浓度,即每个代谢物除以样品的总浓度,以此来校正个体差异或其他因素对代谢物绝对浓度的影响。标准化是对不同样品代谢物的操作,即统计学意义上的变量标准化。标准化的目的是消除不同代谢物浓度数量级的差别,但同时也可能会过分夸大低浓度组分的重要性,即低浓度代谢物的变异系数可能更大。数据转换是指对数据进行非线性变换,如log转换和power转换等。数据转换的目的是将一些偏态分布的数据转换成对称分布的数据,并消除异方差性的影响,以满足一些线性分析技术的要求。不同的预处理方法会对统计分析结果产生不同的影响(见表1),我在实际应用中,们应该根据具体的研究目的﹑数据类型以及要选用的统计分析方法综合考虑,选择适当的预处理方式。例RobertA.vandenBerg等(2006)通过实际代谢组如,
学数据的分析发现,选用不同预处理方法在很大程度上影响着主成分分析(PCA)的结果,自动尺度化(au-toscaling)和全距尺度化(rangescaling)在对代谢组学数据进行探索性分析时表现更优,其PCA分析后的结
如何从复杂的代谢组学数据中提取出有价值的信息,筛选出潜在的生物标志物成为近年来代谢组学研究的
热点和难点。据此,本文针对目前代谢组学数据分析中的常用统计学方法及其研究进展进行介绍。
代谢组学数据的特点
代谢组学是系统生物学领域中继基因组学和蛋白质组学之后新近发展起来的一门学科,它通过检测生物体在受到外源刺激或基因修饰后其体内代谢物质的
[6]
变化来探索整个生物体的代谢机制。其研究对象为生物体内所有内源性小分子代谢物(分子量<1000Da),研究手段为高通量检测技术和数据处理方法,最终目标是数据建模和生物标志物的筛选。生物
LC/样品如血浆、尿液、组织等,经过GC/MS、NMR、MS等高通量仪器检测后,得到大量的图谱数据,使用XCMS[7]等软件对这些图谱数据进行转换,获得用于统计分析的标准格式的数据。归纳起来,代谢组学数
据具有以下特点:
(1)高噪声:生物体内含有大量维持自身正常功能的内源性小分子,具有特定研究意义的生物标志物只是其中很少一部分,绝大部分代谢物和研究目的无关。
(2)高维、小样本:代谢物的数目远大于样品个数,不适合使用传统的统计学方法进行分析,多变量分
[8]析容易出现过拟合和维数灾难问题。
(3)高变异性:一是不同代谢物质的理化性质差[9]
异巨大,其浓度含量动态范围宽达7~9个数量级,二是生物个体间存在各种来源的变异,如年龄、性别都
可能影响代谢产物的变化,三是仪器测量受各种因素影响,容易出现随机测量误差和系统误差,这使得识别有重要作用的生物标志物可能极其困难。
*国家自然科学基金资助(81172767);高等学校博士学科专项基金(20122307110004)
1.哈尔滨医科大学卫生统计学教研室(150081)2.山东大学卫生统计学教研室
E-mail:likang@Δ通信作者:李康,
·358·ChineseJournalofHealthStatistics,Apr2014,Vol.31,No.2
[11]
果在生物学上能够得到更合理的解释
表1
方法归一化
最大峰归一化总峰和归一化中心化均值中心化尺度化AutoscalingRangescalingParetoscalingVastscalingLevelscaling数据转换Log转换Power转换
珋*:xsi=i=1/J∑xij,
j=1J
。
目的
除了进行传统的单变量假设检验分析,代谢组学分析中通常也计算代谢物浓度在两组间的改变倍数值
(foldchange),如计算某个代谢物浓度在两组中的均值之比,判断该代谢物在两组之间的高低表达。计算ROC曲线下面积(AUC)也是一种经常使用的方法
[14]
常见的数据预处理方法
公式
消除某些混杂
因素对代谢物浓度的干扰
xij'=xij/xij
max
。
多变量分析
xij'=xij/∑xij
消除初始值大小的影响
珋xij'=xij-xi
消除不同代谢
物浓度数量级的差别)
代谢组学产生的是高维的数据,单变量分析不能揭示变量间复杂的相互作用关系,因此多变量统计分析在代谢组学数据分析中具有重要的作用。总体来说,代谢组学数据多变量统计分析方法大致可以分为两类:一类为非监督的学习方法,即在不给定样本标签的情况下对训练样本进行学习,如PCA、非线性映射(NLM)等;另一类为有监督的学习方法,即在给定样本标签的情况下对训练样本进行学习,如偏最小二乘
珋xij'=(xij-xi)/Si
珋xij'=(xij-xi)/(xi
max
-xi
min
珋ixij'=(xij-xi)/珋珋xij'=[(xij-x(xi)/si]i/si)
珋珋xij'=(xij-xi)/xi
线性化数据结
构,消除异方差性的影响
xij'=log2(xij)ijxij='
DA)、判别分析(PLS-基于正交信号校正的偏最小二
DA)、乘判别分析(OPLS-人工神经网络(ANN)、支持PCA、PLS-DA和OPLS-DA向量机(SVM)等。其中,
是目前代谢组学领域中使用最为普遍的多变量统计分
析方法。
PCA是从原始变量之间的相互关系入手,根据变异最大化的原则将其线性变换到几个独立的综合指标上(即主成分),直观地描述不取2~3个主成分作图,同组别之间的代谢模式差别和聚类结果,并通过载荷图寻找对组间分类有贡献的原始变量作为生物标志物。通常情况下,由于代谢组学数据具有高维、小样本PCA的分类结果往的特性,同时有噪声变量的干扰,
PCA作为代谢组学数据的预往不够理想。尽管如此,
分析和质量控制步骤,通常用于观察是否具有组间分[15]
类趋势和数据离群点。在组间分类趋势明显时,说明其中一定有能够分类的标志物。PCA还可以用于分析质控样品是否聚集在一起,如果很分散或具有一定的变化趋势,则说明检测质量存在一定的问题。ZhangZhiyu等(2010)通过PCA成功区分了骨肉瘤患者和正常人,并发现良性骨肿瘤患者中有两例是异常[16]
值。KishoreK.Pasikanti等(2009)利用PCA对尿液膀胱癌代谢组学数据进行分析后观察到质控样品在PCA得分图上紧密聚集,从而验证了仪器检测的稳定
[17]性和代谢组学数据的可靠性。
PLS-DA是目前代谢组学数据分析中最常使用的
(xij∑j=1
2珋-xi)/(J-1)
单变量分析方法
单变量分析方法简便﹑直观和容易理解,在代谢组学研究中通常用来快速考察各个代谢物在不同类别
之间的差异。代谢组学数据在一般情况下难以满足参数检验的条件,使用较多的是非参数检验的方法,如Wilcoxon秩和检验或Kruskal-Wallis检验,t'检验也是一种比较好的统计检验方法。
由于代谢组学数据具有高维的特点,所以在进行
单变量分析时,会面临多重假设检验的问题。如果我们不对每次假设检验的检验水准α进行校正,则总体犯一类错误的概率会明显增加。一种解决方法是采用Bonferion校正,即用原检验水准除以假设检验的次数m作为每次假设检验新的检验水准(α/m)。由于Bonferion校正的方法过于保守,会明显降低检验效能,所以在实际中更为流行的一种做法是使用阳性发
FDR)。这种方法可用现错误率(falsediscoveryrate,
于估计多重假设检验的阳性结果中,可能包含多少假
阳性结果。FDR方法不仅能够将假阳性的比例控制在规定的范围内,而且较之传统的方法在检验效能上也得到显著的提高
[12]
。实际中也可以使用局部FDR
一种分类方法,它在降维的同时结合了回归模型,并利
用一定的判别阈值对回归结果进行判别分析。ZhangTao等(2013)运用PLS-DA技术分析尿液卵巢癌代谢
组学数据,成功将卵巢癌患者和良性卵巢肿瘤患者以及子宫肌瘤患者相互鉴别,并鉴定出组氨酸、色氨酸、
[18]
核苷酸等多种具有判别能力的卵巢癌生物标志物。
(用fdr表示),其结其定义为某一次检验差异显著时,果为假阳性的概率。局部FDR的使用,使得我们能够估计出任意变量为假阳性的概率,通常情况下有FDR≤fdr[13]。
20144312·359·
PLS的思想是,通过最大化自变量数据和应变量数据集之间的协方差来构建正交得分向量(潜变量或主成分),从而拟合自变量数据和应变量数据之间的线性
[19]
关系。PLS的降维方法与PCA的不同之处在于PLS既分解自变量X矩阵也分解应变量Y矩阵,并在分解时利用其协方差信息,从而使降维效果较PCA能
[20]
够更高效地提取组间变异信息。当因变量Y为二
1;分类情况下,通常一类编码为1,另一类编码为0或-当因变量Y为多分类时,则需将其化为哑变量。通
DA模型拟合效果使用R2X、R2Y和Q2Y常,评价PLS-DA模型拟这些指标越接近1表示PLS-这三个指标,
R2X和R2Y分别表示PLS-合数据效果越好。其中,
DA分类模型所能够解释X和Y矩阵信息的百分比,
DAQ2Y则为通过交叉验证计算得出,用以评价PLS-Q2Y越大代表模型预测效果较好。模型的预测能力,
PLS-DA得分图常用来直观地展示模型的分实际中,
类效果,图中两组样品分离程度越大,说明分类效果越
显著。代谢组学数据分析中另一种常用的方法是OPLS-DA,DA的扩展,它是PLS-即首先使用正交信号校正技术,将X矩阵信息分解成与Y相关和不相关的两类信息,然后过滤掉与分类无关的信息,相关的信息主要集中在第一个预测成分。JohanTrygg等认为该方法可以在不降低模型预测能力的前提下,有效减
[21]
少模型的复杂性和增强模型的解释能力。与PLS-2
R2Y、Q2Y和OPLS-DA得DA模型相同,可以用RX、
分图来评价模型的分类效果。CarolynM.Slupsky等(2010)使用OPLS-DA发现卵巢癌患者、乳腺癌患者、
正常人这三者之间的尿液代谢轮廓显著不同,从而推断尿液代谢组学可能为癌症的特异性诊断提供重要
[22]
依据。
使用由于代谢组学数据具有高维、小样本的特性,有监督学习方法进行分析时很容易产生过拟合的现
DA在无差异象。为此,需要使用置换检验考察PLS-[23]
情况下的建模效果。该方法在固定X矩阵的前提下,随机置换Y分类标签n次,每次随机置换后建立新
22
DA模型,与的PLS-并计算相应的RY和QY;然后,
真实标签模型得到的结果进行比较,用图形直观表达是否有过拟合现象。
由于样本量的不足,通常采用上述的交叉验证和置换检验方法作为模型验证方法。而实际中,在样本量允许的情况下,最为有效的模型验证方法即将整个数据集严格按照时间顺序划分为内部训练数据和外部测试数据两部分,利用内部训练数据建立模型,再对外部测试数据进行预测,客观地评价模型的有效性和适用性。
生物标志物的筛选
代谢组学分析的最终目标是希望从中筛选出潜在
的生物相关标志物,从而探索其中的生物代谢机制,因此需要借助一定的特征筛选方法进行变量筛选。对于高维代谢组学数据的特征筛选,研究的目的是从中找出对样本分类能力最强或较强的一个或若干个变量。特征筛选方法主要分为三类:过滤法、封装法和嵌入[24]
法。过滤法主要是采用单变量筛选方法对变量进行筛选,优点是简单而快捷,能够快速的降维,如t'检
Wilcoxon秩和检验、SAM等方法。封装法是一种验、
多变量特征筛选策略,通常是以判别模型分类准确性作为优化函数的前向选择、后向选择和浮动搜索特征变量的算法,它通常是按照“节省原则”进行特征筛选,最终模型可能仅保留其中很少部分的重要变量,如遗传算法等。嵌入法的基本思想是将变量选择与分类模型的建立融合在一起,变量的重要性评价依靠特定分类模型的算法实现,在建立模型的同时,可以给出各
DA方法的VIP统计量变量重要性的得分值,如PLS-等。为了更加客观、全面地评价每个变量的重要性,代谢组学研究中一般采取将上述方法结合起来的方式进
行变量筛选。比较常见的一种策略是先进行单变量分析,再结合多变量模型中变量重要性评分作为筛选标准,如挑选fdr≤0.05和VIP>1.5的变量作为潜在生物标志物。用筛选的潜在生物标志物对外部测试数据集进行预测,评价其预测效果。最后,可以通过研究生物标志物的生物学功能和代谢通路,分析不同生物标志物之间的相互作用和关系,从而为探索生物代谢机制提供重要线索和信息。YangJinglei等(2013)即在代谢组学分析中使用fdr≤0.2和VIP>1.5的双重标准来筛选精神分裂症的特异生物标志物,所筛选出的差异代谢物其AUC在训练数据中达94.5%,外部测试数据中达0.895
[25]
。展
望
由于代谢组学数据变量多、关系复杂的特性,数据分析任务极为艰巨。目前常用的统计学方法在一定程度上为进行代谢组学数据分析提供了有效的工具,但仍然存在诸多不足。如在代谢组学研究中,生物样品DA或之间的变异性往往较大,目前最流行的PLS-OPLS-DA数据分析方法在差异小、噪声大时,模型使PLS-DA和OPLS-DA均是基用效果不够理想。另外,于线性回归的方法,但是代谢组学数据通常不是简单PLS-DA和OPLS-DA模型拟合数的线性关系,因此,
据的结果可能会不够好。基于这些问题,一些学者开始尝试将一些新的高维数据分析方法和思想应用于代谢组学数据分析中,如LinXiaohui等(2011)提出一种将支持向量机、随机森林和遗传算法结合起来进行变量筛选的分析思路,通过比较证实其较单个分析方法
(下转第365页)
ChineseJournalofHealthStatistics,Apr2014,Vol.31,No.2·365·
quit;
参
考
文
献
])+ya[i-1,]);temp=T(hlast*(0:nints[i-1,f=j(nrow(temp),ncol(x),0);dow=1tonrow(temp);dov=1toncol(x);
f[w,v]=last[w,]*pdf('normal',temp[w,],x[,v],stdv[i,]);end;end;
last=t(0.5*hlast*(2*f[+,]-f[1,]-f[nrow(f),]));end;end;end;return(cp);finish;
storemodule=cprob;storemodule=drift;
1.O'BrienPC,FlemingTR.Amultipletestingprocedureforclinicaltrials.
56.Biometrics,1979,35(3):549-2.HwangIK,ShihWJ,DeCaniJS.Groupsequentialdesignsusingafamilyof
typeIerrorprobabilityspendingfunctions.StatMed,1990,9(12):1439-45.3.Reboussin,DM.ProgramsforcomputinggroupsequentialboundariesusingtheLan-DeMetsmethod.1996,GeorgeWashingtonUniversity
Washington,DC.4.Turnbull,BW,JennisonC.Groupsequentialmethodswithapplicationstoclinicaltrials.2000,BocaRaton:Chapman&Hall/CRC.390.
(责任编辑:郭海强)
(上接第359页)
能够发掘出更多的信息,尤其适合分析复杂生物数
[26]
据;ElonCorrea和RoystonGoodacre(2011)提出了BN),一种新型的遗传算法—贝叶斯网络方法(GA-这还能种方法在有效筛选变量并提高分类效果的同时,研究变量间的相互作用和关系
[27]
。毫无疑问,这些新
方法的提出将会为代谢组学数据分析提供新的思路和
契机。随着各种代谢组学检测仪器的快速发展,更有效的代谢组学数据分析技术亟待开发,值得更多的生物统计学者关注和研究。
参
考
文
献
1.DunnWB,EllisDI.Metabolomics:Currentanalyticalplatformsandmethodologies.TrACTrendsinAnalyticalChemistry,2005,24(4):285-294.2.许国旺,路鑫,杨胜利.代谢组学研究进展.中国医学科学院学报,2007,29(6):701-711.3.SpratlinJL,SerkovaNJ,EckhardtSG.Clinicalapplicationsofmetabo-lomicsinoncology:areview.ClinCancerRes,2009,15(2):431-440.4.WishartDS.Applicationsofmetabolomicsindrugdiscoveryanddevel-opment.DrugsRD,2008,9(5):307-322.5.TaylorJ,KingRD,AltmannT,etal.Applicationofmetabolomicstoplantgenotypediscriminationusingstatisticsandmachinelearning.Bioinformatics,2002,18(2):241-248.6.NicholsonJK,LindonJC,HolmesE.'Metabonomics':understandingthemetabolicresponsesoflivingsystemstopathophysiologicalstimuliviamultivariatestatisticalanalysisofbiologicalNMRspectroscopicdata.Xenobiotica,1999,29(11):1181-1189.7.SmithCA,WantEJ,O'MailleG,etal.XCMS:processingmassspectrom-etrydataformetaboliteprofilingusingnonlinearpeakalignment,matc-hing,andidentification.AnalyticalChemistry,2006,78(3):779-787.8.SimaC,DoughertyER.Whatshouldbeexpectedfromfeatureselectioninsmall-samplesettings.Bioinformatics,2006,22(19):2430-2436.
9.DunnWB,EllisDI.Metabolomics:Currentanalyticalplatformsandmethodologies.Trac-TrendAnalChem,2005,24(4):285-294.
10.GoodacreR,BroadhurstD,SmildeA,etal.Proposedminimumrepor-tingstandardsfordataanalysisinmetabolomics.Metabolomics,2007,3
(3):231-241.
11.VandenBergRA,HoefslootHCJ,WesterhuisJA,etal.Centering,scal-ing,andtransformations:improvingthebiologicalinformationcontent
ofmetabolomicsdata.BMCGenomics,2006,7:142-156.
12.BenjaminiY,HochbergY.Controllingthefalsediscoveryrate:apracti-calandpowerfulapproachtomultipletesting.JRStatistSocB,1995,
57(1):289-300.
13.刘晋,张涛,李康.多重假设检验中FDR的控制与估计方法.中国卫
2012,29(2):305-308.生统计,
14.BroadhurstDI,KellaDB.Statisticalstrategiesforavoidingfalsediscov-eriesinmetabolomicsandrelatedexperiments.Metabolomics,2006,2
(4):171-196.
15.TryggJ,HolmesE,LundstedtT.Chemometricsinmetabonomics.JPro-teomeRes,2007,6(2):469-479.16.ZhangZ,QiuY,HuaY,etal.Serumandurinarymetabonomicstudyof
humanosteosarcoma.JProteomeRes,2010,9(9):4861-4868.17.PasikantiKK,EsuvaranathanK,HoPC,etal.Noninvasiveurinary
metabonomicdiagnosisofhumanbladdercancer.JProteomeRes,2009,9(6):2988-2995.18.ZhangT,WuXY,KeCF,etal.IdentificationofPotentialBiomarkers
forOvarianCancerbyUrinaryMetabolomicProfiling.JProteomeRes,2013,12(1):505-516.19.蒋红卫,夏结来.偏最小二乘回归及其应用.第四军医大学学报,
2003,24(3):280-283.20.BoulesteixAL,StrimmerK.Partialleastsquares:aversatiletoolforthe
analysisofhigh-dimensionalgenomicdata.BriefBioinform.2007,8
(1):32-44.
21.TryggJ,WoldS.Orthogonalprojectionstolatentstructures(O-PLS).
JournalofChemometrics,2002,16(3):119-128.22.SlupskyCM,SteedH,WellsTH,etal.Urinemetaboliteanalysisoffers
potentialearlydiagnosisofovarianandbreastcancers.ClinCancerRes,2010,16(23):5835-5841.23.WesterhuisJA,HoefslootHCJ,SmitS,etal.AssessmentofPLSDA
crossvalidation.Metabolomics,2008,4(1):81-89.24.YvanS,IakiI,PedroL.Areviewoffeatureselectiontechniquesin
bioinformatics.Bioinformatics,2007,23(13):273-281.25.YangJ,ChenT,SunL,etal.Potentialmetabolitemarkersofschizo-phrenia.MolecularPsychiatry,2013,18(1):67-78.26.LinXH,WangQC,YinPY,etal.Amethodforhandlingmetabonomics
datafromliquidchromatography/massspectrometry:combinationaluseofsupportvectormachinerecursivefeatureelimination,geneticalgo-rithmandrandomforestforfeatureselection.Metabolomics,2011,7(4):549-558.
27.CorreaE,GoodacreR.Ageneticalgorithm-Bayesiannetworkapproach
fortheanalysisofmetabolomicsandspectroscopicdata:applicationtotherapididentificationofBacillussporesandclassificationofBacillusspecies.BMCBioinformatics,2011,12(1):33-49.
(责任编辑:郭海强)
正在阅读:
代谢组学数据分析的统计学方法07-21
2018-2024年中国积木玩具市场竞争格局报告(目录) - 图文11-11
今年的冬天作文500字07-17
关于31 人们针对某一需要而为计算机编制的指令序列称为_06-24
策划培训课程(040529)游戏的完成度05-22
防水工程施工承包合同格式文本05-22
河南安阳师范学院数字化校园信息化建设方案 - 图文04-18
七年级数学上册 2.6有理数的加减混合运算导学案3(新版)北师大06-06
成交的八大步骤01-19
- 教学能力大赛决赛获奖-教学实施报告-(完整图文版)
- 互联网+数据中心行业分析报告
- 2017上海杨浦区高三一模数学试题及答案
- 招商部差旅接待管理制度(4-25)
- 学生游玩安全注意事项
- 学生信息管理系统(文档模板供参考)
- 叉车门架有限元分析及系统设计
- 2014帮助残疾人志愿者服务情况记录
- 叶绿体中色素的提取和分离实验
- 中国食物成分表2020年最新权威完整改进版
- 推动国土资源领域生态文明建设
- 给水管道冲洗和消毒记录
- 计算机软件专业自我评价
- 高中数学必修1-5知识点归纳
- 2018-2022年中国第五代移动通信技术(5G)产业深度分析及发展前景研究报告发展趋势(目录)
- 生产车间巡查制度
- 2018版中国光热发电行业深度研究报告目录
- (通用)2019年中考数学总复习 第一章 第四节 数的开方与二次根式课件
- 2017_2018学年高中语文第二单元第4课说数课件粤教版
- 上市新药Lumateperone(卢美哌隆)合成检索总结报告
- 统计学
- 数据分析
- 代谢
- 方法
- 郑大内科护理消化系统疾病患者的护理作业满分答案
- 2016年成都理工大学博士研究生招生简章
- 幼儿园小班科学教案范文4篇
- 畜牧业发展不足之处调研报告
- 第三讲诊断试验准确性研究的报告规范_STARD介
- 有关团结合作的成语、谚语和名人名言
- 怎样成为一名优秀英语教师
- 黄冈市2012年中考物理适应性模拟试题及答案(二)
- 公司上市财务培训课程提纲
- 梦幻西游手游降妖秘境副本冲榜攻略
- 第十三章整式的乘除测验题
- 山东省卫生计生系统人员“针对性普法”题库(医护人员)
- 路基附属作业指导书
- 胜利油田工贸公司LFHVI型COD自动分析仪用户手册
- 土壤碱解氮含量的测定
- 大学院校招生代码及专业编码
- 三明市尤溪县英语六年级上册期末模拟复习卷(1)
- 遵守校规校纪主题班会
- 2011年河南省数据结构试题集(10套题并附带答案)知识大全
- 生活老师岗位职责