R在数据框中按照要求提取指定行

我们这里的用到的是R自带的数据集iris

gerp函数的使用如下：

grep(pattern, x, ignore.case = FALSE, perl = FALSE, value = FALSE, fixed = FALSE, useBytes = FALSE, invert = FALSE)

各参数的含义如下：

subset函数可以从向量、矩阵、数据框中提取子集

使用如下：

subset(x, subset, select, drop = FALSE, ...)

有时候我们需要提取非表格文件中含有某些符号的特定的行，例如，我们需要从下面fasta文件中提取所有包含">"的行，

基本思路是逐行读入，定位，按照定位提取每一行

可以看到，使用grep后，我们得到包含“>”的行号，我们利用这个行号进行定位

即可得到如下结果

首先，既然你要画残差(residual)的qqplot 那一定要先做回归。

你没说明回归自变量和因变量都是什么，假设第二列是因变量pressure[2]，第三列是自变量pressure[3]

1.回归

>lm.fit<-lm(pressure[2]~1+pressure[3])

2.做qq图

这里不用qqplot这个指令，用qqnorm和qqline就好

>qqnorm(lm.fit$res)

>qqliine(lm.fit$res)

提取？就直接write.table()

比如你要第二列和第三列

>write.table(pressure[c(2,3)],file="c:\test.txt",quote=F,row.names=F,sep="\t")

第一个参数pressure[c(2,3)]是你要提取的变量

第二个参数file= 是你要存放的地址和文件名

第三个参数quote=F 是你打出来的数字都没有引号（默认是字符型）

第四个参数row.names=F 是不要行名称（为什么不要？自己试试就知道了）

第五个参数sep="\t" 每列数据之间空一个tab的距离，也就是8个字节，清晰明了。

---