SAS数据展现(tabulate)PPT.ppt

用指数算符产生向量 用指数算符(:)可以方便地产生向量。如果步长不是1，可以用DO函数。 ? 例18.24 用指数算符产生向量。 使用矩阵表达式 矩阵表达式由一系列名字、算符和函数构成。它的作用是完成计算，判断条件和巧妙处理数值。表达式可以出现在赋值语句的任何一边。 算符 矩阵表达式中使用的算符可以分为三种： 前缀算符； 中缀算符； 后缀算符。 算符优先级 复合表达式 a=x+y+z; a=x+y*z; a=(-x)#(y-z); a=x+y*z; a=x/y/z; 二目元素算符 二目元素算符是对两个矩阵的元素分别进行操作。 下标 下标的一般形式： operand[row,column] 其中： operand通常是一个矩阵名，也可以是表达式或矩阵本身； row规定矩阵的行数，可以是表达式，也可以是标量或向量； Column规定矩阵的列数，可以是表达式，也可以是标量或向量。 ? 下标的用途： 选择矩阵的一个元素； 选择矩阵的一整行或列； 选择矩阵的子阵； 简化矩阵的一行或列。 例18.25 选择矩阵的一个元素。 print coffee[rowname=names]; c21=coffee[2,1]; /*矩阵coffee第2行1列元素*/ print c21; 也可以寻找超出行范围的元素，这时得到的是以行为顺序第6个元素。 c6=coffee[6]; print c6;? 例18.26 选择矩阵的一整行或列。 jim=coffee[3,]; /*第3行元素*/ print jim; friday=coffee[,5]; /第5列元素*/ print friday; 例18.27 选择矩阵的子阵。 submat1=coffee[{1 3},{2 3 5}]; /*第一个向量｛1 3｝选择行，第二个向量｛2 3 5｝选择列*/ print submat1; /*可以事先创建一个向量并以他们的名字作为输入*/ rows={1 3}; /*事先创建向量rows*/ cols={2 3 5}; /*事先创建向量cols*/ submat1=coffee[rows,cols];/*以事先创建向量的名字作输入*/ submat2=coffee[1:3,3:5]; /*选择矩阵coffee的第1至3行，3至5列*/ print submat2; 例18.28 根据下标位置进行相应元素的赋值。 coffee[1,2]=4; /*将矩阵coffee的第1行2列处的元素值改为4*/ print coffee; coffee[,5]={0,0,0,0}; /*将矩阵coffee的第5列元素全部设为0*/ print coffee; t={ 3 2 -1, 6 -4 3, 2 2 2 }; print t; i=loc(t0); print i; I 3 5 下标也可以包含运算结果为行数或列数的表达式。例如表达式可以为： t[loc(t0)]=0; t[i]=0; print t; T 3 2 0 6 0 3 2 2 2 牛顿迭代法求f(z)=0一个实根的方法 牛顿迭代法公式： 上例中函数为： 理解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; 创建数值矩阵 含有多个元素的矩阵标识要用括号把元素括起来。用逗号将矩阵的行分开。