BoxCox变换及其在STATA软件中实现.docVIP

BoxCox变换及其在STATA软件中实现 　【摘要】 　　 目的：向医疗卫生领域的一线工作者介绍BoxCox变换。 方法 ：利用STATA软件包拟合BoxCox回归模型。结果：通过使用示例数据库进行拟合，成功的进行了BoxCox变换，并进行了回归 分析 。结论：使用BoxCox变换是解决应变量不符合正态分布的重要方法，而STATA软件提供了进行BoxCox变换的可行之道。 【关键词】 BoxCox变换　STATA软件　boxcox命令 　　 　　1 BoxCox变换简介 在统计工作中，特别是在一般线性模型的运用中，经常出现连续型应变量不符合正态分布的情况(在线性模型中对自变量没有类似的要求［1］)，此时如直接使用线性模型，显然是不合适的，而通常的解决办法有两种：一是对应变量y进行变量变换；二是使用广义线性模型。本 研究 讨论前者。 对应变量y进行变量变换有许多方法，如倒数变换、指数变换等，而这些变换都可以通过一个公式而统一起来，这就是Box和Cox在1964年提出的BoxCox变换，其形式为［2］：y(λ)=yλ-1λ, if λ≠0 logy, if λ=0显然，该变换在λ=0时为对数变换，λ=-1时为倒数变换，而在λ=0.5时为平方根变换［3］。BoxCox变换还有一个扩展公式： y(λ)=(y+a)λ-1λ, if λ≠0log(y+a), if λ=0 此时的a是为了使y+agt;0。 论文代写 无论是基本式还是扩展式，真正需要估计的只有一个参数λ(因为a的取值非常明显)，而对λ的估计，在Box和Cox论文中采用了两种方法，其一是最大似然估计，其二是Bayes方法。但对于一般的用户而言，两种方法的 理论 和推证过程都可以不管，最重要的是如何在实践中运用，这也是本研究所要重点关注的。 如何在实践中完成BoxCox变换呢?最直接的当然是使用公式，但最简单的还是使用相关的统计软件。现在我们就用较流行的统计软件STATA来实现BoxCox变换过程，由于STATA7.0以后的版本和以前的版本在变换命令的使用上有很大的不同，故本研究分别进行表述。 毕业论文 　　2 示例一：STATA6.0软件的实现过程 为了说明 问题 ，我们使用一个叫auto的数据库进行示例，这个数据库来自1974年4月美国的消费报告和美国EPA的燃料消耗统计［4］，我们的任务是研究汽车时速(mpg)与车重(wEight)和价格(price)之间的关系。作为应变量的mpg经检验不符合正态分布。 在STATA6.0软件中有直接的命令可供调用，其命令的基本格式为［5］：Boxcox y x1 x2 x3,g(y1) 毕业论文 这里，y为应变量，xi表示的是自变量，y1为进行转换后生成的新变量，其命令名称boxcox，此处STATA只认为最靠近命令的变量为应变量，而其它的变量则被认为是自变量。如果要继续进行BoxCox回归，则在运行boxcox命令后紧接着运行regress命令从而得到回归的结果。下面我们来看实例。 毕业论文 　　2.1 命令 Boxcox mpg wEIght price, g(y1) Regress 　　2.2 结果 论文代写 　　2.2.1 BoxCox变换结果 经过迭代，程序 计算 出λ值(即L)为-0.7618，并对L进行了检验。 论文代写 　　2.2.2 BoxCox回归结果 　　2.2.3 不进行变量变换而直接进行回归分析结果 毕业论文 　　2.2.4 变量在BoxCox变换前后正态性的检验结果 　　从表4可见原来不符合正态分布的mpg现在符合了(见y1的检验)，且对比结果表2和表3就会发现，两者的结果是有差异的，进行了应变量变换后的BoxCox回归的决定系数要高于普通的回归过程，且变量price也变得有意义了。表1 BoxCox变换结果（略）表2 BoxCox回归结果（略）表3 直接进行回归分析结果（略）表4 变量在BoxCox变换前后正态性的检验结果（略） 　　　3 示例二：STATA7.0软件的实现过程 实例同示例一。在STAT7.0软件中也有直接的命令可供调用，其命令的基本格式为［6］：Boxcox y x1 x2 x3, lrtest 注意，此处不能再如6.0版那样生成新变量，另外，此时命令运行的结果已经包含了BoxCox回归，所以在6.0中的regress命令也被取消了。而lrtest的意义是对回归系数进行检验。下面我们来看实例。编辑。 论文代写 　　3.1 命令 Boxcox mpg wEight price, lr