灰色预测模型及SAS实现

中国生统卫计 026 年0 2 月第23 第卷1 　

期·

75

灰色·测模型预 SA及S 现实

南医方科学大珠江医信息院科(5 0182)2 浙万里学院江础基校( 315学01) 0中山大学公共卫学生院统与计流病学行(系 500108)

颜 　杰

　丽驰相

方

积

乾

　　灰色 测是预灰色统系理论的要重组部成分, 由于其 有具所样需本数少 据不、要求型典概的率布分特等点 , 在许多 领域到了广泛的得应 用。灰色 预测模型的基 原本理是对始 数据进原累 加行 , 生成有具一规定的律列 序 然,建立后应的相分方 微程模型, 用 预测来系统动态的发过程 展。但 色灰预测需较复要杂的学数算运 ,在各种 统计软件中都没 提供有成的现程或过块模以直接计可 算 我。们 助 S借SA 供提的 I LM( Int raetivc Mateir L anguxag e, 互交式矩阵言语 )模 块, 制 编SS A程可以较 简序便、 精 地确实现GM ( 1 , 1) 模型 。 .1色 灰M(G1 1) ,型的建模模程

过0 )原始数据设序 列 x(= { x (0) (k), k 1= , 2 … ,,n } 1) 1 )(1 一阶)累加, 生成数据序列 x ( {=x (( k) k=

I,N PU tTyear xt @;@/ 读入原＊数据 始序 列 ＊/y t x t +; /生成＊阶 一加序列 ＊累 /in dx =1 ; zt =e( -yt +LA ( Gyt) )/2; / 为数据矩＊ 阵 B备准数据 /＊ AD T A INLE S;1 　9109 24　39 5 2　 1　991 　22568 3　199 226　091　 4 　19 9 32　3793 5 　1949 27　87 1 6　　 195 　928721 7　196 929　28 　78 　919 　70067 93 19　9 　83791 　100　199 9 31　84 11 　22000 　3376 　121　200 1　 3545 8 P;RC OI LM ;/ ＊用 IM 调 L模块＊/ USE a1; / ＊ 开已打的 有ASS数 集 a据 1 /＊ RADE LL VAA { zt Rndie}Ix N O TBW HRE( Et z.=) ; ＊/ a1 中将量 变tz和 niexd值( 不含 第读入矩阵B ＊ /RED AA L VL A R{x t}I N T O Yn WEREH(z t .=) ;/ ＊将 a1 中量 变xt 矩阵 Yn ＊/ah ta=IN V (` ＊B B)＊B ` ＊ Y n ;/＊ 计算 参数矩 阵 a= a , [u] = [BTB

]1-T

1 2 ,,… ,n } , 其中 x

( )

1(

)k =∑x i 1

=k

(

)

( 0)i

(2 ) 确定G M(1 ,1) 模 : 相应型白化微分方程 (为1 ) xd+ax ( ) 1 u= td 参数列记为 a= [ , au

] T1a =( B BT) BYn T

TB Y ＊/n

,

根 最据二乘法解小:

ah

tat =aaht`; na ={ au ;} /＊将参矩 数转阵 ＊置 /C EATR 2 a FORM haat[t OCL NAM E=na ;] /用＊转置后 的参数 阵数矩建据 立A SS数 集 据2 ＊a/ A PPNEDF OM aRhta t; ＊/将据读入到数据数集 /＊

其中数

据矩阵-[ x B = [ -

x( 1) ( )

1

( 1+ )x (2 )+x 

(

1) (1

) ()2 ]/ 2 3) (]/2

1 1 ,

QU I

T; /＊出 退IML 模块 ＊ /A T D aA3; S ETa2 ; ndiex =1 ;DA T Aa4 ; SE a1 T ;I FDA T Aa5; MERG Ea1 a 34 a ;BYin edx; I

0F)0 N )=1T HNEx p x=t ; ＊/x (( )1=x ( 1(＊)/

n

Y1)

[ x -(1)( n -1)+ x ( (n )] /2 1( )0( 0) = [ (x2 )x , 3)(, …, x (0 ) ( n) ]T

N= 1; x t0=x ;

t(

3 预测模型 ) 上:述一线性阶微方分 的程解为即G M( 1 , 1) 测预

模型x (k +1) =[x ( k )- u ] -ak e+u k, = , 21 ,… n,a a　 　( 4)生成 原始序预测列模型

(1 ) (0

)

EKE Pt0 ixdnx e;

E

SE LO

1)D　 　 yt1=( xt0 -/ u )a＊E PX -(a ＊ t (-1 ))u+/ a; / ＊计算 x (0) k -) (1)=k x[ ((k -) / au e ]-( + a/ ua ＊ 1)/　 　 y0 =t ( x0t- /u )＊aX E( P-a＊ (t - ) 2+u)/ a /; 计＊算x ( 0 k)- 2 ()k 1)-= [ (x (k) - / u]ae -a( + u/ a ＊/0 1)) 1 　)x　p=y t1 - y 0t ;/ ＊计算 测预值 x ((k)= [x ((k) ] x-(

x

x( 0)( 0 )

( 1

=)

( x)0

( 1 ,)( k -)x

( 1)

(

k )=x

(

1 )

(

k- 1)

0

)k ) =(11 - ae) ( x( (k) - u)e a -( a,　k = , 32, …, n

(

-1)＊ /k NED ;er o rr= p xx -t ; / 计算绝对＊差误 /＊ error =errror / x ＊1t0 0;/ 计算＊相对差误 ＊ D/ORP tyi dexn zty t y 10 txt0 ;

.2灰 G色M( 1 ,1)模 型的 ASS 现实 广州市 19以09～ 2001 年 共1 2年 医院位床数标指的动 时间序态列编 SA制S程 如序下:

D T A a1A ;/ ＊建原立数始据集 1 a/

＊( 下第 85转 )

页

中国卫统计生200 6 年 2第 23月 卷 1 期第　

·

8 ·5

显示为组表高 达。 上结以果示我提们这基些可因能与肠癌结的发 有很生大系 。 至 于关些这因是基否真的正导 致机癌体变 或,其诱机理如何 发,都 要进需步的一论和理物生 学验来实证 验 讨 　　。论 于对我们选出筛的来6 4 基组因 它们每,组参都 与同一了个谢通代 ,路 过考通它们察表的 ,达我们 得如 出表 1 示所结果 , 表达的近相组的 α值的

nla. gvo/w i t2)数 库来研究据基之因间功的能相关 。 关于性因基功能相关的性和表达相似的之性间的关 联 ,在还罕有现理上论论证的 大多,工作都的停是留 对现在的考察象 上, 文本也如是此 仅仅,察了考与同一参 谢代路通部分的基因。 基因并不是直接来 挥生物学发功能 。的其间 要经历还D A N转录的 、RN A翻的 、译 白蛋质的饰等修程 过。因此 ,

直接考基因的表达察和 其经多步过化转而的成白质蛋的生物功学能之的关 间联 ,势 必会遗一漏些中步骤间转的规化律 ,削也了结 弱论可信度 。 的此因必有要在续后的究研中合结考基 察在 因cDAN水 平或上蛋质白水上的平达情况表

。 参 考 　　文　

1献.史忠 植. 知发现识. 京北: 清 大华出学版社 ,20 20. . Ro2andl ansJn , Dove reGenabu m , Mrk G aret sine. Re ltian Wgolh-Geenmeo xpEressio n at D w ait h Prto eni-PotrienIn trecations. C lo SpridgnHar bo Laborrt aor yresP s, 202 , 01 :237-46 . .3 A ikihr oNa kyaa , S usuu moGt o, M noir uK nehiasa Ext r.ctioa no forreCltea G end eClsters uy b Mlt ipuelG rah Cpo pmri aon .sG enom eIn fomartcis, 2 010, 1 :2 4-53 . 4. A4l o n U B,raai kN, N o t treamnDA ,t ale Br.oa dapt ernts fog een exper ssoi rnevaeled ybclus etrig analynissof tmou ard norman lolcon i stuessp obedr byol gionucloeidtea rars y Bio.olg y ,919 9, 9( 62)1 :645 -76705 . 5 . 政 , 李霞郭 饶 ,奇 绍 .医 信 息学分析 方 法 .哈滨:尔 哈尔 滨 出 版社 ,0201 6. 孙之荣 . 后.因基组息学信. 北 : 清京华大出版社 , 学200 .27. 沈 王 羽氵 非 方,福德 .真核 基因达表控 调 北京. :等教育出版高 社,1 99 .6[ 7

]

(下转 第57 页

P)RCO PIR TND A AT =a5 ; /＊出计输算果结＊ / U R N /; 运行上述程序 ＊＊ S/A 输S结果出: ye

ra 9190 191 19929 193 9914 1995919 96 1979 91981 999 200 02010 t x429352 528 2669012 739 378271 82721 27298 30067 3091731 24833 7163 548 xp 254935 .00 318.3 3 1errro . 00001. 25892 -. 06276 1-. 03898 07.0 816 -08.0 951 -07 64.629 1 08.04 10 564.7 21 .8066 -61 .97922-1 .40948 error r.0 0000

0

　其　中 p 即为x灰预色值测 ,err o 和rr eror 分r为模 别型回的代绝对差( 残误) 和相差对误差于用衡量型模预测 精度的 。.3论结 灰模色型求要较有的多阵矩运算, 使应用受到 定的限一制。 如果手用工计或算步分计算 ,不仅 杂麻 烦 ,复中间结 果保留的精度也会预给测结果带来定的 一误 差 我们用。序程算了验多篇文的计算结果 献 都存 ,在同程度的误差不。 计在机算普及今天 ,的用 件完软成计算 既简,便 快、 , 又速能得精到确的算计果结 。 在此序程的基上 ,础 做稍改就修可实现以他灰其预色测模 型 。 AS S IM的 L具很强有大矩的阵计功能 ,算其 操作 式方 既可以交 互 进 行 , 也以 可通过 程编 现实 。而 且 ILM 接提供矩直和 SAS 阵数集据间相的互转 , 使用 非常简换 便 是,一非个有常的力计工算 。

具

2346.636 5-4. 5653 2170. 990 -22.90 792 7958 5.72 804.701 2 . 4705 16.- 994

29

35. 58 -7912.12 39390. 148 12723. 62 2173.8 9 332.4 43 8814 .23 689. 9249

3131.1 0 5-640.50 13047.110 - 468 .039