R語(yǔ)言編程基礎(chǔ)篇(1)

1. 建立多維數(shù)組

array(1:36, dim = c(2, 2, 3, 3))

2. 使用list.files函數(shù)遍歷文件夾中的文件

比如列出當(dāng)前工作目錄下的所有文件

list.files(getwd())

[1] "1.pdf" "10plots.pdf"

[3] "140408696.txt" "1plots.pdf"

[5] "2plots.pdf" "3plots.pdf"

還可以設(shè)置正則表達(dá)式來(lái)過(guò)濾文件

列出當(dāng)前工作目錄下所有的R語(yǔ)言源代碼文件

list.files(getwd(),pattern = '*.[R|r]$')

[1] "Ask.R" "gg.R"

3. 用sciplot包的畫boxplot圖

#加載數(shù)據(jù)

library(MASS)

cab<-data.frame(cabbages)

cab[1,]

library(sciplot)

bargraph.CI(Cult, HeadWt, group =Date , data =cab,

xlab = NA, ylab = NA, cex.lab = 1.5, x.leg = 1,

col = "black", angle = 45, cex.names = 1.25,

density = c(0,20,100), legend = TRUE,ylim=c(0,5))

box()

bargraph.CI(Date,HeadWt, group =Cult , data =cab,

xlab = NA, ylab = NA, cex.lab = 1.5, x.leg = 1,

col = "black", angle = 45, cex.names = 1.25,

density = c(0,20), legend = TRUE,ylim=c(0,5))

box()

在這個(gè)boxplot中，不用自己計(jì)算均值和誤差，也不用自己調(diào)整數(shù)據(jù)結(jié)構(gòu)，不用轉(zhuǎn)化為matrix，很方便。

4. 多系列圖

x<-seq(from = 1,to = 9,by =2)

y<-seq(from = 2,to = 10,by= 2)

barplot(rbind(x,y),beside=T,col=heat.colors(2))

5. 應(yīng)該注意避免使用的變量名

R語(yǔ)言中預(yù)定義了大量函數(shù)，有些函數(shù)名相當(dāng)簡(jiǎn)單，比如c()函數(shù)，根據(jù)本人的實(shí)踐經(jīng)驗(yàn)，應(yīng)該避免使用過(guò)于簡(jiǎn)單的變量名，以免與R語(yǔ)言已有名稱沖突，而出現(xiàn)意想不到的錯(cuò)誤。自己取變量名字是，最好能加上自己的一些特征，比如公司縮寫，比如本人公司首字母縮寫為MS，則本人使用的變量名都以MS開(kāi)頭，然后接下劃線，比如MS_Alarm，MS_Books等等。

以下是盡量應(yīng)該避免使用的變量名：

單個(gè)字符：a,b,c,d,…,x,y,z,A,B,…,X,Y,Z,

已經(jīng)被R語(yǔ)言使用的名字：data,names,dim,seq,…

另外，命名新變量時(shí)，應(yīng)該先檢查一下變量名是否已經(jīng)存在。

可以使用get()函數(shù)查詢變量名，看是否有返回值。

6. Windows環(huán)境下，R語(yǔ)言調(diào)用C語(yǔ)言庫(kù)

1.安裝Rtools，http://www.murdoch-sutherland.com/Rtools/Rtools.exe這個(gè)在寫R包時(shí)要用，當(dāng)然這里不會(huì)提到。

2.配置環(huán)境變量，“我的電腦”–>“屬性”–>“高級(jí)”–>“環(huán)境變量”–>“系統(tǒng)變量”–>PATH，在后面添加：D:\Rtools\bin; D:\Rtools\perl\bin;D:\Rtools\MinGW\bin;D:\R-2.8.1\bin（前三個(gè)取決于Rtools的安裝位置，最后一個(gè)取決于R安裝的位置）

3.編寫C代碼（命名為fac3.c），下面程序得到i*j*k（1<=i,j, k<=n）的和，注意：參數(shù)必須用指針表示(對(duì)應(yīng)R中的向量)，且編寫的C函數(shù)不能有顯式返回值，即函數(shù)返回要聲明為void；（參見(jiàn)http://www.wentrue.net/blog/?p=72，該文是在linux下的。）

voidfac3(double *n, double *m){

*m=0;

int i,j,k;

for(i=1;i<=*n;i++){

for(j=1;j<=*n;j++){

for(k=1;k<=*n;k++){

*m=*m+i*j*k;

}

}

}

}

4.編譯C，在DOS中（fac3.c所在的文件夾下）輸入RCMD SHLIB fac3.c

即可得到.dll文件

5.R中調(diào)用，并與R做循環(huán)的速度進(jìn)行比較

dyn.load("fac3.dll")

system.time(out<-.C("fac3",a=1000,b=1))#第一個(gè)參數(shù)對(duì)應(yīng)C中的函數(shù)名

我的老機(jī)上運(yùn)行時(shí)間：

user system elapsed

6.67 0.00 6.81

再看看R中的速度

n= 100

m= 0

system.time(for (iin1:n) {

for (jin1:n) {

for (kin1:n) {

m = m + i * j * k

}

}

})

user system elapsed

7.34 0.00 7.44

通過(guò)比較可以發(fā)現(xiàn)，調(diào)用的C做了1000^3次循環(huán)比R中做了100^3循環(huán)的速度還快！