Listing 1.0 补充一些常用操作
# 代码内容(如注释)包含中文时出现乱码,可使用UTF-8 或者 GB2312(File → Reopen with Encoding → UTF-8)
getwd() # 查看当前工作路径
ls() # 查看环境中的变量
rm() # 删除变量 rm(list=ls()) 为删除所有的变量。慎重!
# 符号 # 后面的内容R不会识别,用于注释。
q() # q()用于关闭软件
save(variable1,variable2,file="保存名.Rdata") # 保持数据
Listing 1.1 A Sample R session
# 赋值方法
age <- c(1,3,5,2,11,9,3,9,12,3) # 将向量c()的内容赋值给age; <-符号为右向左赋值
c(1,3,5,2,11,9,3,9,12,3) ->age # 也可以用 -> 符号左向右赋值。
weight <- c(4.4,5.3,7.2,5.2,8.5,7.3,6.0,10.4,10.2,6.1) # 将向量c()的内容赋值给weight
# 几个基本统计量
mean(weight) # mean()求平均
sd(weight) # sd()求方差
cor(age,weight) # cor()求相关性
plot(age,weight) #做一个散点图
Listing 1.2 An example of commands used to manage the R Workspace
setwd("Path/") # 设置工作路径
# 另外 getwd(),为显示路径,括号内无内容时为显示当前工作路径。
options() # 显示环境变量,即包含程序设定的一般参数的集合。
options(digits=3) #修改环境变量中的数值参数,所有数值呈现为有效数为后3位(会对第4位进行“四舍五入”)
a<-1.346
a # 0.1235,显示有效位数3位,第4位为6,“四舍五入”后进一位。
# runif(n,min=0,max=1);n表示生成的随机数数量,min表示均匀分布的下限,max表示均匀分布的上限;若省略参数min、max,则默认生成[0,1]上的均匀分布随机数。
x <- runif(20) #生成20个区间为[0,1]的均匀分布随机数
x # 这里x的值是20个均匀分布的随机数,此处x不执行(options(digits=3))有效位数为3。
# summary(),简单统计:从最小值到最大值,将极差区间一分为4,以次展示。
summary(x) # Min(最小值),1st Qu(25%值),Median(中位数),Mean 3rd Qu.(75%),Max(最大值)
hist(x) # 柱形图(频数分布图)
savehistory() # 保持执行历史操作代码,可点击保存按钮。
save.image() # 生成一个“.RData”的文件,可点击查看
Listing 1.3 - Working with a new package
help.start() # 打开查看 R 中文件的小窗口,
# install.packages("vcd") # 已安装
help(package="vcd") # 查看vcd包的帮助,查看包不用 ?vcd
library(vcd)
help(Arthritis) # 或者 ?Arthritis
Arthritis # 为 vcd 包中的数据,"关节炎"治疗研究的数据,也可用 View(Arthritis) 查看
example(Arthritis)
# q() # 关闭软件