函数玩一玩 | 【SAS Says·扩展篇】IML：2.函数

【SAS Says·扩展篇】IML

分6集，回复【SASIML】查看全部：

• 入门 | SAS里的平行世界
• 函数 | 函数玩一玩
• 编程 | IML的条件与循环
• 模块 | 5分钟懂模块
• 穿越 | 矩阵与数据集的穿越
• 作业 | 编一个SAS回归软件

本节目录：

1. 元素函数

2. 矩阵函数

（1）矩阵查询函数

（2）矩阵生成函数

【SAS Says · 扩展篇】IML：函数玩一玩

精彩内容

第一部分介绍元素函数，它针对矩阵内元素进行操作，如将矩阵里的数据取绝对值、取余等等；第二部分介绍矩阵的函数操作，它针对整个矩阵操作，如查询矩阵的列数、行数等信息、将矩阵分块、插入一个子矩阵、生成一个随机数矩阵等。

1. 元素函数

元素函数是针对矩阵里的元素进行操作的，比如：

举个栗子：

例子

proc iml; a={9 -2.38, -1 0, 2 1}; c1=abs(a); c2=exp(a); c3=int(a); c4=mod(a); print a,c1,c2,c3,c4; quit;

2. 矩阵函数

矩阵函数将对整个矩阵进行操作：

矩阵查询函数可以查询矩阵的一些基本信息，如行数、列数、行最大值、列最大值等等；

矩阵生成函数可以对矩阵进行分块、取对角线等，生成一个新的矩阵。 （1）矩阵查询函数

• all (条件)：当矩阵所有元素都满足条件时，返回1，否则返回0；
• any (条件)：只要矩阵中有一个元素满足条件，就返回1，否则返回0。

例子

proc iml; a={9 -2.38, -1 0, 2 1}; c1=all(a>0); c2=any(a>0); print a,c1,c2; quit;

• loc(条件)：返回满足条件的元素的标号，还可以与截取运算符[]搭配使用，获取标对应的数值。

例子

proc iml; a={9 -1 3, 3 -3 0}; c1=loc(a>0); c2=a[ loc(a>0) ]; print a,c1,c2; quit;

• nrow(矩阵)：求矩阵行数；
• ncol(矩阵)：求矩阵列数；
• type(矩阵)：得到矩阵的类型，数值型返回N、字符型返回C，如果矩阵不含任何值，返回U。
• length(矩阵)：求矩阵每个元素的长度，只能是字符矩阵。

例子

proc iml; a={9 -1 3, 3 -3 0}; b={"hello","world"}; c1=nrow(a); c2=ncol(a); c3=type(a); c4=length(b); print a,c1,c2,c3,c4; quit;

（2）矩阵生成函数 通过这类函数得到一些简单、特殊的矩阵。

• block(M1,M2...)：创建分块对角阵；
• diag(M)：只保留矩阵M的对角线，其余元素均赋值为0；
• vecdiag(M)：将矩阵M的对角线元素变成列向量。

娘子，快来看大楼~！————→

例子

proc iml; a={9 -1 , 3 -3}; b={1 2, 4 5}; c=block(a,b); c1=diag(c); c2=vecdiag(c); print a,c1,c2; quit;

• i(n)：创建n×n的单位阵；
• j(row,col,value)：创建row×col的矩阵，矩阵的数据均为value；
• repeat(value,row,col)：和上面的一样，生成row×col、元素均为value的矩阵；
• insert(M1,M2,插入第n行,或插入第m列)：将矩阵M2插入矩阵M1中，插入第n行，或者第m列，注意m和n必须有一个为0，因为这两个只能指定一个。

右边的朋友，我要是能和你一样高就好了。

例子

proc iml; a=i(3); b=j(2,4,7); c=repeat(7,2,4); d=insert(a,{1 -1 0},2,0); print a,b,c,d; quit;

• uniform(seed)：生成（0,1）均匀分布的伪随机数；
• normal(seed)：生成均值为0，方差为1的伪随机数；

例子

proc iml; a=normal(1); b=normal({1 1 1,1 1 1, 1 1 1}); c=normal(repeat(1,3,3)); d=uniform({2 2 2}); print a,b,c,d quit;