热线电话:13121318867

登录
首页精彩阅读SAS Says 扩展篇IML:函数玩一玩
SAS Says 扩展篇IML:函数玩一玩
2016-02-03
收藏



SAS Says 扩展篇IML:函数玩一玩

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)矩阵查询函数

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;

D3

数据分析咨询请扫描二维码

最新资讯
更多
客服在线
立即咨询