近年来DNA芯片技术日益成为研究大量基因表达连续变化的实验室工具。芯片技术的发展使得同时获得成千上万个基因的表达谱成为可能。DNA芯片在产生的短短几年时间已经显现出在基因序列分析、基因诊断、基因表达研究、基因组研究、发现新基因及各种病原体的诊断等生物医学领域中的应用价值。利用芯片数据，“癌变基因”的发现以及对飞速增长的基因组数据库增加功能注释等难题将迎刃而解。DNA芯片数据具有高维（成千上万个基因）和样本小（通常小于30）的特点，为了避免对伪结果进行分析，准确估计抽样方差很重要。在微阵列试验中既要包括真实的变化，又需要随机的变异。大量研究表明，聚类分析及其相关技术对于挖掘基因表达的相关模式非常有用。但是仅用这些方法不能对分析结果进行统计推断，难以得到具有生物学意义的结论，尤其是不适合分析前后时间点数据高度相关的动态基因表达数据。 本文描述的统计框架包含了基因表达分析的众多目标，与现有的分析方法完全一致，同时提高了这些方法的效用。本文着重研究差别表达基因的鉴定。本研究提出了基于混合线性模型的分析微阵列数据的方法，并将其应用于差别表达基因的鉴定、在动态或静态过程中估算基因主效应以及预测基因与环境的互作效应。用蒙特卡罗模拟对该方法的有效性和可靠性进行了比较系统的研究。这种方法可以有效地将基因表达水平根据变异来源的不同剖分为几个组成部分。主要研究内容和结论概述如下： 1．提出了分析芯片数据的一般模型，其中包括了基因、阵列效应、染料、处理效应以及基因×阵列、基因×染料、基因×处理互作效应。根据不同的试验设计，该模型可以做适当的调整。本文提出的方法主要分为两步来进行：首先，将芯片数据通过噪音过滤消除大的试验系统误差，然后在一个比较宽松的标准下通过单基因模型初步判断差异表达基因；其次，用多基因模型分析这些初定的差异表达基因以便在较严的标准下控制假阳性。用MINQUE法估计各项效应的方差和协方差分量，用AUP法预测随机效应。基因和处理的互作效应作为鉴定差异表达基因的具体指标。 2．对新提出的基于混合线性模型分析DNA芯片数据的方法用蒙特卡罗模拟进行了验证。模拟结果表明该方法在绝人多数情况下忧于传统的t检验和 WOlfinger提出的混合模型方法。验证了基因和处理的互作效应可以作为鉴定差异表达基因的更为恰当的指标。 3．研究表明我们提出的基于混合线性模型的方法可以无偏或近无偏地估算固定效应和预测随机效应。对基因主效应的无偏估计值和基因与处理互作效应的无偏预测值进行聚类可以获得具有统计学和生物学意义的结果。 4．将我们提出的混合线性模型进行拓展，可以用来分析动态的基因表达数据。我们定义了一个新变量度量给定卜1时刻的基因表达量来确定1时刻的基因表达情况，用条件变量的方法来估计条件方差、预测条件遗传效应，可以揭示在特定时间段基因表达的变异情况。 5．对新提出的基于条件变量的分析芯片数据的方法进行了蒙特卡罗模拟研究。结果表明基于条件变量的分析方法在大多数情况下表现得比差值法更有效。同时结果还进一步显示了将基囚和环境的互作效应作为鉴定差异表达基因的指标是非常有效的。 6．为了适应实际分析的需要，用C／C＋＋语言编写了软件，可以用于分析基因芯片的表达数据，估算基因表达变异来源的方差组成和预测遗传效应，同时寻找差异表达基因。 7．以几种药物处理特异癌症细胞系的实际芯片实验数据的分析为例，说明了本研究所提方法的分析过程及分析所得结果的生物学意义。