如何用R语言编写一些程序

tomorrow

一般情况下, 一个R程序以ASCII 格式保存,扩展名为`.R'. 如果一个工作
要重复好多次, 用R程序是一个不错的选择. 在我们的第一个例子中, 我们想对
三个不同种属的鸟绘制一样的图, 而且数据在三个不同的文件中. 我们将一步
一步的演示看R用不同的方式去完成这个简单的问题.
首先, 我们凭直觉连续键入一系列命令, 而且预先分割图形设备.
layout(matrix(1:3, 3, 1)) #分割图形界面
data <- read.table("Swal.dat") #读入数据
plot(data$V1, data$V2, type="l")
title("swallow") #增加标题
data <- read.table("Wren.dat")
plot(data$V1, data$V2, type="l")
title("wren")
data <- read.table("Dunn.dat")
plot(data$V1, data$V2, type="l")
title("dunnock")
字符`#' 用于在程序中添加注释行: R 会自动跳过注释行.
第一个程序的问题是在我们加入另外一个物种数据时, 它过长. 此外, 一
些命令多次执行, 因此它们可以放在一起,在执行的时候仅仅修改一些参数. 这
里的策略是把参数放到一个字符型的向量中去, 然后用下标去访问这些不同的
值.
layout(matrix(1:3, 3, 1)) # 分割图形界面
species <- c("swallow", "wren", "dunnock")
file <- c("Swal.dat" , "Wren.dat", "Dunn.dat")
for(i in 1:length(species)) {
data <- read.table(file) # 读入数据
plot(data$V1, data$V2, type="l")
title(species) # 增加标题
}
注意代码read.table()里面的参数file 上面没有双引号, 因为这个
参数是字符型.
现在的代码比较紧凑. 它比较容易加入新的物种, 因为设置物种名字和数
据文件的向量都程序的前端.
如果扩展名为`.dat'的数据文件在R的工作目录下面, 程序可以正常运行,
否则用户要设置工作目录, 或者设置绝对路径(例如: file <- "/home/paradis/
data/Swal.dat"). 如果程序保存在文件Mybirds.R 中, 可以通过键入如下命
令执行:
> source("Mybirds.R")
和所有以和所有以文件作为输入对象的函数一样, 如果该文件不在当前工作目录下
面, 用户需要提供该文件的绝对路径.