R語言基本操作

R語言和Python類似，也是有許許多多的包(Python中是模塊)組成，具體用什么就安裝什么，用的時候再去加載。

1、安裝包：

語法：install.packages("包名")(兩個引號""不能省略)

如:下面就是安裝強大的繪圖包 ggplot2

>install.packages("ggplot2")

2、加載包

語法:library("包名")

>library("ggplot2")

3、查看數(shù)據(jù)類型

語法：class(變量)（數(shù)據(jù)框(data.frame)是R語言特有的數(shù)據(jù)格式）

> class(iris)

[1] "data.frame"

4、取數(shù)據(jù)框中的行和列

語法：<數(shù)據(jù)框>[行，列](如果取從..行（列）到...行（列），中間用:劃分)

還是以iris的數(shù)據(jù)集為例子。

取第十行:

>iris[10,]

取第一行到第100行

>iris[1:100,]

取第三列：

>iris[,3]

取第一列到第三列：

>iris[,1:3]

如果要連續(xù)取多行(列)的話用c來取

iris[c(1:35,50:85,100:135),]

取列：

iris[,c(1:2,4:5)]

5、數(shù)據(jù)框拼接

按行拼接：

rbind

用法：rbind(data.frame1,data.datafram2,........)

> train_data1=iris[1:5,]

> train_data2=iris[7:8,]

> rbind(train_data1,train_data2) 

Sepal.Length Sepal.Width Petal.Length Petal.Width Species

1          5.1         3.5          1.4         0.2  setosa

2          4.9         3.0          1.4         0.2  setosa

3          4.7         3.2          1.3         0.2  setosa

4          4.6         3.1          1.5         0.2  setosa

5          5.0         3.6          1.4         0.2  setosa

7          4.6         3.4          1.4         0.3  setosa

8          5.0         3.4          1.5         0.2  setosa

按列拼接：

cbind

用法：cbind(data.frame1,data.datafram2,........)

> iris[3:4,] 

Sepal.Length Sepal.Width Petal.Length Petal.Width Species

3          4.7         3.2          1.3         0.2  setosa

4          4.6         3.1          1.5         0.2  setosa

> iris[100:101,]   

Sepal.Length Sepal.Width Petal.Length Petal.Width    Species

100          5.7         2.8          4.1         1.3 versicolor

101          6.3         3.3          6.0         2.5  virginica

> train_data1=iris[3:4,1:2]

> train_data2=iris[100:101,4:5]

> cbind(train_data1,train_data2) 

Sepal.Length Sepal.Width Petal.Width    Species

3          4.7         3.2         1.3 versicolor

4          4.6         3.1         2.5  virginica

利用iris的數(shù)據(jù)來比較隨機森林和支持向量機在分類中的作用

#加載支持向量機的包

library("e1071")

library(kernlab)

#加載隨機森林包

library("randomForest")

#加載iris的數(shù)據(jù)集

data_iris=iris

#得到訓練數(shù)據(jù)集

train_data=data_iris[c(1:35,50:85,100:135),]

#得到測試數(shù)據(jù)集

test_data=data_iris[c(35:50,85:100,135:150),]

#隨機森林測試

rFM=randomForest(Species~.,data=train_data,importance=TRUE)

rfm_result=predict(rFM,test_data[,c(1:4)])

#支持向量機測試

svmmodel=svm(Species ~ ., data =train_data)

svm_result=predict(svmmodel,test_data[,c(1:4)])

#隨機森林比較

a=0

b=1

while(b<=48){

if(rfm_result[b]==test_data[,5][b]){

a=a+1

}else{

a=a

}

b=b+1

}

#隨機森林正確率

a/48

#支持向量機比較

j=0

i=1

while(i<=48){

if(svm_result[i]==test_data[,5][i]){

j=j+1

}else{

j=j

}

i=i+1

}

#支持向量機正確率

j/48

結果是：隨機森林正確率百分之百，支持向量機是0.979.

出現(xiàn)這個結果應該是數(shù)據(jù)集太少了吧。目前的目的只是為了完成論文。。。。