SAS Says 擴(kuò)展篇IML：函數(shù)玩一玩

1. 元素函數(shù)

元素函數(shù)是針對(duì)矩陣?yán)锏脑剡M(jìn)行操作的，比如：

舉個(gè)栗子：

例子
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. 矩陣函數(shù)

矩陣函數(shù)將對(duì)整個(gè)矩陣進(jìn)行操作：

矩陣查詢(xún)函數(shù)可以查詢(xún)矩陣的一些基本信息，如行數(shù)、列數(shù)、行最大值、列最大值等等；

矩陣生成函數(shù)可以對(duì)矩陣進(jìn)行分塊、取對(duì)角線(xiàn)等，生成一個(gè)新的矩陣。
（1）矩陣查詢(xún)函數(shù)

all (條件)：當(dāng)矩陣所有元素都滿(mǎn)足條件時(shí)，返回1，否則返回0；

any (條件)：只要矩陣中有一個(gè)元素滿(mǎn)足條件，就返回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(條件)：返回滿(mǎn)足條件的元素的標(biāo)號(hào)，還可以與截取運(yùn)算符[]搭配使用，獲取標(biāo)對(duì)應(yīng)的數(shù)值。

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

nrow(矩陣)：求矩陣行數(shù)；

ncol(矩陣)：求矩陣列數(shù)；

type(矩陣)：得到矩陣的類(lèi)型，數(shù)值型返回N、字符型返回C，如果矩陣不含任何值，返回U。

length(矩陣)：求矩陣每個(gè)元素的長(zhǎng)度，只能是字符矩陣。

例子
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）矩陣生成函數(shù)

通過(guò)這類(lèi)函數(shù)得到一些簡(jiǎn)單、特殊的矩陣。

block(M1,M2...)：創(chuàng)建分塊對(duì)角陣；

diag(M)：只保留矩陣M的對(duì)角線(xiàn)，其余元素均賦值為0；

vecdiag(M)：將矩陣M的對(duì)角線(xiàn)元素變成列向量。

娘子，快來(lái)看大樓~！————→

例子

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)：創(chuàng)建n×n的單位陣；

j(row,col,value)：創(chuàng)建row×col的矩陣，矩陣的數(shù)據(jù)均為value；

repeat(value,row,col)：和上面的一樣，生成row×col、元素均為value的矩陣；

insert(M1,M2,插入第n行,或插入第m列)：將矩陣M2插入矩陣M1中，插入第n行，或者第m列，注意m和n必須有一個(gè)為0，因?yàn)檫@兩個(gè)只能指定一個(gè)。

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

例子

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）均勻分布的偽隨機(jī)數(shù)；

normal(seed)：生成均值為0，方差為1的偽隨機(jī)數(shù)；

例子
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;