模拟研究中随机数据的产生及SAS实现.pdfVIP

JournalofMathematicalMedicine Vo1．21 No．6 2008 文章编号：1004—4337(2008)06—0704—03 中图分类号 ：R195．1 文献标识码：A · 统计分析 · 模拟研究中随机数据的产生及 SAS实现 钱 俊 张 超 陈平雁 (广州南方医科大学公卫学院生物统计系 广州 510515) 摘 要： 随机数据的产生是随机模拟研究中重要的一环，通过数学推导和软件应用可以产生各种分布类型的随机数据，仅以单 变量、多变量以及生存时间随机数据的生成及其实例介绍其 SAS实现。 关键词： 随机模拟； 随机数据； 分布类型； 生存时间 随机模拟是一种广义的数值计算方法 ，它利用计算机模 分布的随机数据，则需要运用概率论的知识，变换数据从 “标 拟实际过程，然后加以统计处理并求得实际问题 的解 。由于 准”到 “一般”。 模拟的过程可以反复进行，系统结构和参数的改变也都 比较 产生单变量的随机数据的步骤分三步：确定分布类型和 容易，因此应用广泛，适用性强，特别在传统理论难以解决的 参数取值，利用随机数发生器产生相应 “标准”分布的随机数 问题上不失为一种有效可行的方法[1j。近年在医药学领域， 据 ，变换数据从 “标准”到 “一般”。下举例说明： 随机模拟的应用研究不断深入l2]。 例 1 生成 1000个服从参数为a，b的均匀分布U(a，6) 在随机模拟的实现中，随机数据的产生是最基本而重要 的随机数据(a≠O，6≠1)。 的一个环节。随机数据应根据问题实际背景和数据的数字特 方法 1：调用RAND()函数，产生均匀分布U(O，1)的随 征产生，具有可信性。较普遍的方法是以实际观测的数据集 机数据，再进行变量变换，转换为U(a，6)上的随机数据。 为例，借鉴其数据结构，协变量的取值产生随机数据_5]。 dataExamplel； doi一1tO1000 1 单变量随机数据的生成 -z—RAND(uniform)； ／*产生服从 “标准”分布 U SAS中提供了产生常见分布随机数的函数 RAND()。 (O，1)的随机数据*／ 它易于使用，基本上实现了 “标准分布”数据的随机生成。语 一(6一口)* + ； ／*转换为 “一般”分布U(a，6)的 法为：RAND(分‘布类型’，参数 1，参数2……)，涵盖了常用的 随机数据 *／ 分布类型：二项分布、正态分布、指数分布、对数正态分布、泊 output； 松分布、均匀分布、威布尔分布等。但它产生的随机数据一般 end； 都服从 “标准分布”：O～1上的均匀分布 u(0，1)，参数值为 1 方法2：调用RANUNI(seed)函数或UNIFORM(seed)函 的指数分布 E(1)等。如果模拟研究 中需要获得一组 “一般” 数。 著性差异前提下，男婴和女婴出生时间高峰无显著性差异 (P O．05)，与黄细铭等报道相同l4]，可以合并。初产妇和经产 参 考 文 献 妇分娩时间，以及其年龄对分娩时间影响的圆形分布分析，到 目前未见报导。本资料结果显示：初产妇和经产妇分娩时间 1 金丕焕，主编．医用统计方法 ．第 1版．上海医科大学 出版社， 高峰期无显著性差异(PO．05)，可以合并。本资料结果提 1993，2O8～ 214． 示，孕妇的分娩时间可能与孕妇本身人体生物节律分泌的缩 2 黄志勇，陈峰，刘关键 ．《中国医学百科全书 ·医学统计学》统计软 宫素有关，与婴儿的性别、孕妇年龄、孕妇初产与否无关。根