利用subset()函数提取数据
从矩阵中提取某些必要的信息时,可以利用subset()函数轻松实现
subset(data,sunbset,select)#data=要处理的数据,要显示的字符,选择显示哪些字符
#利用鸢尾花(iris)数据集举个例子
> a <- iris[1:150,]
> head(a)
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
6 5.4 3.9 1.7 0.4 setosa
#查看一下Species有多少个分组
> unique(a$Species)
[1] setosa versicolor virginica
Levels: setosa versicolor virginica
#看看setosa组有多少个样本
> length(which(a$Species=="setosa"))
[1] 50
#利用subset()函数将"setosa"提取出来
> a_setsub <- subset(a,Species == "setosa")
> head(a_setsub)
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
6 5.4 3.9 1.7 0.4 setosa
#看一下是不是都提取出来了
> length(a_setsub$Species)
[1] 50
#现在只想显示最后两列
> a_setsub <- subset(a,Species == "setosa",select = c("Petal.Width","Species"))
> head(a_setsub)
Petal.Width Species
1 0.2 setosa
2 0.2 setosa
3 0.2 setosa
4 0.2 setosa
5 0.2 setosa
6 0.4 setosa
#本次的例子中,length(),which()的用法也可以注意一下.