R
是非常容易入门的一个统计软件或语言,它是辅助工具,不是主体,应该配合统计相关内容一起学习。
我第一次听说这个软件还是在浙林的时候,邬荣岭老师的讲座或是培训会议上。真正开始接触是2014年,因为要利用asreml
去计算遗传参数。最高效的学习往往发生在你急需这些知识或技巧的时候,所以最好是手上有需要使用R
的任务,当然在这之前可以熟悉一下简单的操作比如+-*/、数据的导入。
我主要用到的参考资料有
-
R语言实战(第2版)
-
R数据可视化手册——这是本工具书,查阅画图参数
-
华南农大林元震老师主编的r与asreml-r统计分析教程——主要是后边asreml-r部分,我觉得林老师参考了一些新西兰Canterbury大学Luis Apiolaza老师的ASReml Cookbook。Luis的个人网页也有一些有用的内容和Twitter。
林老师的这本书出到了第2版,比上一版厚了一倍,如果把程序输出结果拿掉能去掉一半的内容。这本书更像一本讲义或者是笔记,没有很多理论性的内容,通常只有一个简略的概念,一个样例和代码及太多的输出结果。所以这是本工具书,需要相关的代码直接翻书查找(高级代码还得去stackoverflow)。没有统计知识也是它的优点,你看不到容易恐惧的公式和推导,除了LMM那一丁点儿。另外一个优点是第7章——试验设计。《实战》没有这部分内容,也未见到有相关的R
书籍。如果林老师下次能把试验设计交代得更细致一点,比如给个图就更好了。另一个特点还是遗传评估这一章,对于遗传育种专业并且需要进行参数估算的geneticist和breeder可以重点敲这一部分的代码。林老师还有一个非常重要的贡献在于他自己写了一部分批处理的函数AAfun
,在处理数据时,尤其有较多的变量时会比较方便。除了这2章外,其余的内容在其他的R
工具书中一般也都有涉及,所以拿这本书学习R也是合适的(但我首推的还是《实战》,无出其右者)。同样,这本书也事无巨细。达成一个目的途径很多,掌握自己最舒服的一个就可以了,其余的不感兴趣大可不必学习。
所以,自认为一个入门者,以林老师的这本第2版,按照浓缩一下数据管理(导入、处理和导出)这一部分,给还未接触而又有必要学习的同仁铺点路,也给自己一个梳理复习的机会。
21 July 2017 更新
第1章 R简介
- 简介不是重点,不用看。
- Rstudio不是强烈建议要安装,而是必须安装!!!
第2章 基础语法
- 2.1 对象与变量
- 赋值符号
<-
在Rstudio
中的快捷键是'Alt'+'-'。 - 2.1.1 变量创建 比如有树高和胸径要算材积指数,(当然在Excel中算好后导入也可以),学会方法一就够用了。给的例子是直接搬《实战》的,连数据都没改但书的架构为什么不借鉴一下?
- 2.1.2 变量命名用方法二
names()
,矩阵没有变量,行名和列名分别是rownames()
和colnames()
。 - 2.2 运算符为何不放在最前面?
- 2.3 表达式一看就会不是吗?有几个情况会用到表达式:
- 写模型时
- 画图涉及到公式时(需要在图上放拟合模型)
- 构建关联矩阵时(
model.matrix
)
- 2.4 NA是缺失值, NULL是空值。NULL在写函数时会很常用,比如要先建个空的向量或矩阵,然后循环把数值放进去。
- 2.5 控制结构
- 2.5.1
if/else
写函数时比较常用,switch
用的不多。 - 2.2.2
for
循环是批量处理的关键函数,有人诟病因为它比较慢,但它的呈现非常直接,易读。当对象是list
的时候用lapply
则比for要好用很多。while
也常用,比如迭代时。 - 2.6 自编函数是不是放在数据管理后面会更好?因为学习顺序是那样的,数据还不会导进去怎么会写函数?
第3章 数据创建
- 3.1 数据类型以向量、矩阵、数据框和列表最常用,数组用的少(我没用过)。向量、矩阵和数据框不都是矩阵吗?搞得好复杂,对了还有一个叫scalar的东西,不就是维度不一样吗?列表在最开始学的时候比较让我迷糊,没有直观概念。可以这么说,列表是最general的类型,能装前面所有的类型,所以学的多了以后会偏爱列表,尤其是一个函数要返回多个内容,彼此类型又不一致时。
- 3.1.6 因子 学过试验设计、做过ANOVA分析就知道什么是因子,就是变量是类别,没有顺序也不分大小。
- 3.1.7 字符串 就是一种比较general的向量,文本分析或分析序列时会涉及到非常多的字符串操作比如拆分和合并(我在从文献中根据单词提取句子时,数据类型全是字符串)。
- 3.1.8 日期 我处理的数据没有这个类型,所以跳过去吧。
- 3.2 对象的模式和属性 变量应该是因子而不确定是否做过转换时会查看一下,不写高级函数一般不用。需要的话在
Rstudio
中可以直接查看。 - 3.3 数据输入 学会一个
read.csv
就可以了,整那么复杂干啥?read.csv
和read.table
的区别是前者默认表头(header=T),并且不用sep=','
,非常简便。所以请在导入R前把相应的Excel转换成csv
格式。 - 3.4 数据存储 这个分几种情况:
- 处理了一些
csv
要导出来,用write.csv
。 - 不是表格的话可能要导成
txt
,用write.table
,这个和导入数据用的函数是对应的。默认输出会加一列从1到n编码,不需要的话声明中加row.name=F
。字符向量输出会加双引号,不需要的话加quote=F
。 - 同时保存多个文件,下次想直接导入用时可以导出成R自己的格式
.Rdata
,这个格式的文件也比较小,用save()
而非write
。
- 处理了一些
第4章 数据管理(数据操作比较恰当)
这部分是整个数据分析流程中最繁琐与耗时的,要耐心要认真。
- 4.1 数据(类型)转换和4.2 排序(略)
- 可能会遇到这样一种情况,当列表中的每个元素
- 4.3 数据合并 很多时候需要合并数据,我把我遇到的情况列在下边
-
合并2个表格,比如你对同一块林子在2个时间分别做了调查,一般会有出入比如有的死了数据没录入,而数据分析时你得放到一起。放在Excel中肉眼匹配也可以,但
merge()
能节省许多的时间,允许取2者的交集或并集。举个例子:2个数据集df05
和df10
,df_merge <- merge(x = df05, y = df10, by.x = "Treeid05", by.y = "Treeid10")
,取2者并集加all = T
,取x
和2者的交集加all.x = all
,取2这的共同交集加all.x = all, all.y = all
。 -
合并1个字符串,1个字符串向量有多个元素,需要合并这些元素成1个用
paste()
,加collapse = ''
,单引号中加需要的分隔符。 - 一个指标输在多列,比如
-
合并2个表格,比如你对同一块林子在2个时间分别做了调查,一般会有出入比如有的死了数据没录入,而数据分析时你得放到一起。放在Excel中肉眼匹配也可以,但