上一篇“高级篇:IML(1)”发出来之后,有朋友反映东西东西太简单了,根本不能算“高级”。想想也是,暂时还没有介绍太复杂的SAS程序,于是决定将本篇定为“扩展篇”,SAS Says系列安排如下:
【SAS Says】基础篇
主要介绍SAS软件的一些基本入门知识,包括画图、ODS、导入导出数据数据格式等。已经更新完,可以回复M,进入目录,点击链接查看。
【SAS Says】统计篇
介绍SAS统计方面的一些应用,包括主成分、判别分析、logistic模型、非参数检验等等。统计篇需要晚些才能与大家见面。
【SAS】扩展篇
包括矩阵模块IML、宏语句、sql模块以及贝叶斯等。正在更新中。
本节目录:
1. 元素函数
2. 矩阵函数
(1)矩阵查询函数
(2)矩阵生成函数
【SAS Says】扩展篇:IML(2)
精彩内容
第一部分介绍元素函数,它针对矩阵内元素进行操作,如将矩阵里的数据取绝对值、取余等等;第二部分介绍矩阵的函数操作,它针对整个矩阵操作,如查询矩阵的列数、行数等信息、将矩阵分块、插入一个子矩阵、生成一个随机数矩阵等。
1. 元素函数
元素函数是针对矩阵里的元素进行操作的,比如:
abs( ) | 取绝对值 |
---|---|
exp( ) | e的元素次方 |
int( ) | 取整 |
log( ) | 取ln |
mod( ) | 取余 |
sqrt( ) | 将元素开方 |
举个栗子:
例子
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)矩阵查询函数
例子
proc iml; a={9 -2.38, -1 0, 2 1}; c1=all(a>0); c2=any(a>0); print a,c1,c2; quit;
例子
proc iml; a={9 -1 3, 3 -3 0}; c1=loc(a>0); c2=a[ loc(a>0) ]; print a,c1,c2; quit;
例子
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)矩阵生成函数 通过这类函数得到一些简单、特殊的矩阵。
娘子,快来看大楼~!————→
例子
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;
右边的朋友,我要是能和你一样高就好了。
例子
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;
例子
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;