18IML编程技术讲述.pptVIP

第18章 IML编程技术 本章内容包括： 概述； 理解IML语言； 线型回归IML模块。 矩阵操作； IML编程语句； SAS 数据集操作； 访问外部文件。 概述 SAS/IML(Interactive Matrix Language)是在一种在动态和交互式环境下实现的强大而灵活的过程语言（交互式矩阵语言）。IML最基本的数据对象是矩阵。在SAS/IML模式下，可以立刻看到语句运行的结果，也可以把语句贮存在一个模块中为以后使用。 SAS/IML软件特点 利用SAS/IML软件可以轻易地进行高效编程和使用一系列子程序。因为SAS/IML软件是SAS系统的一个模块，所以，可以使用SAS数据集和外部文件，也可以对一个已有的SAS数据集进行编辑或者创建一个新的SAS数据集。 SAS/IML软件是一种编程语言； SAS/IML软件对矩阵作运算； SAS/IML软件有强大的算符功能，有许多函数和CALL子程序； SAS/IML软件将算符应用到整个矩阵； SAS/IML软件是交互式的； SAS/IML软件是动态的； SAS/IML软件可以处理数据； SAS/IML软件可以作图。 一个简单IML交互程序 例18.1 求平方根IML程序。 proc iml; /* 激活IML */ reset deflib=ResDat; /*指定ResDat为缺省的SAS逻辑库*/ start approx(x); /* 定义模块开始 */ y=1; /* 初始化y */ do until(w1e-3); /* 循环开始, le-3=0.001 */ z=y; y=.5#(z+x/z); /* 估计开方根，牛顿迭代 */ w=abs(y-z); /* 计算估计值之间的差异 */ end; /* 循环结束 */ return(y); /* 返回近似值 */ finish approx; /* 模块结束 */ t=approx({3,5,7,9}); /* 调用函数APPROX */ print t; /* 打印矩阵 */ quit; run; 理解IML语言 定义矩阵 SAS/IM软件对矩阵作运算，这里矩阵的定义同矩阵代数中矩阵的定义一样，是一个两维的（行×列）数组。SAS/IM软件可以定义数值矩阵也可以定义字符矩阵。 ? 特殊矩阵定义： 1×N矩阵称作行向量; M×1矩阵称作列向量; 1×1矩阵称作标量. 矩阵名字与标识 矩阵名字是一个合法的SAS名字，不能超过8字节。 矩阵标识是由它的值表示的矩阵。 ? 例18.2 矩阵名举例。 A, XX, MATRIX1. 例18.3矩阵标识举例。 100, ., HI, {1 2 2 8 8 }, {3, 3, 6,6}. 由矩阵标识创建矩阵 创建标量矩阵 例18.4创建标量矩阵。 proc iml; reset deflib=ResDat; a=100; b=. ; c=IML; 创建数值矩阵 含有多个元素的矩阵标识要用括号把元素括起来。用逗号将矩阵的行分开。 ? 例18.5 由矩阵标识创建的数值阵。 X={1 2 3 4 5 6}; /*定义一个行向量X*/ Y={1,2,3,4,5}; /*定义一个列向量Y*/ Z={1 2, 3 4, 5 6}; /*定义一个3×2矩阵Z*/ W=3#Z; /*创建一个新矩阵W，它的每个元素是矩阵Z相应元素的3倍*/ 创建字符矩阵 例18.6创建字符矩阵时，字符长度取最长字符的长度。 a={IML Modual}; /*字符长度以较长的字符为准，这里为6*/ B={abc DEFG}; 重复因子 重复因子用括弧括起来放在要重复的元素前面。 ? 例18.7 用与不用重复因子的等价语句。 ANS={[2] Yes, [2] No}; ANS={Yes Yes, No No}; 赋值语句 结果＝表达式： w=3#z; ? 函数作为表达式： a=sqrt(b); y=inv(x); r=rank(x); 表达式中的三种运算符 表达式中的算符 语句类型