R语言门外1公里到大门口

R是非常容易入门的一个统计软件或语言,它是辅助工具,不是主体,应该配合统计相关内容一起学习。

我第一次听说这个软件还是在浙林的时候,邬荣岭老师的讲座或是培训会议上。真正开始接触是2014年,因为要利用asreml去计算遗传参数。最高效的学习往往发生在你急需这些知识或技巧的时候,所以最好是手上有需要使用R的任务,当然在这之前可以熟悉一下简单的操作比如+-*/、数据的导入。

我主要用到的参考资料有

  1. R语言实战(第2版)

    Kabacoff - 2015 - R IN ACTION Data analysis and graphics with R.jpg

  2. R数据可视化手册——这是本工具书,查阅画图参数

    R Graphics Cookbook.jpg

  3. 华南农大林元震老师主编的r与asreml-r统计分析教程——主要是后边asreml-r部分,我觉得林老师参考了一些新西兰Canterbury大学Luis Apiolaza老师的ASReml Cookbook。Luis的个人网页也有一些有用的内容和Twitter

    r与asreml-r统计分析教程(第1版)

林老师的这本书出到了第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 表达式一看就会不是吗?有几个情况会用到表达式:
    1. 写模型时
    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.csvread.table的区别是前者默认表头(header=T),并且不用sep=',',非常简便。所以请在导入R前把相应的Excel转换成csv格式。
  • 3.4 数据存储 这个分几种情况:
    1. 处理了一些csv要导出来,用write.csv
    2. 不是表格的话可能要导成txt,用write.table,这个和导入数据用的函数是对应的。默认输出会加一列从1到n编码,不需要的话声明中加row.name=F。字符向量输出会加双引号,不需要的话加quote=F
    3. 同时保存多个文件,下次想直接导入用时可以导出成R自己的格式.Rdata,这个格式的文件也比较小,用save()而非write

第4章 数据管理(数据操作比较恰当)

这部分是整个数据分析流程中最繁琐与耗时的,要耐心要认真。

  • 4.1 数据(类型)转换和4.2 排序(略)
    • 可能会遇到这样一种情况,当列表中的每个元素
  • 4.3 数据合并 很多时候需要合并数据,我把我遇到的情况列在下边
    1. 合并2个表格,比如你对同一块林子在2个时间分别做了调查,一般会有出入比如有的死了数据没录入,而数据分析时你得放到一起。放在Excel中肉眼匹配也可以,但merge()能节省许多的时间,允许取2者的交集或并集。举个例子:2个数据集df05df10df_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
    2. 合并1个字符串,1个字符串向量有多个元素,需要合并这些元素成1个用paste(),加collapse = '',单引号中加需要的分隔符。
    3. 一个指标输在多列,比如
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 200,738评论 5 472
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,377评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 147,774评论 0 333
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,032评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,015评论 5 361
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,239评论 1 278
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,724评论 3 393
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,374评论 0 255
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,508评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,410评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,457评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,132评论 3 316
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,733评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,804评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,022评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,515评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,116评论 2 341

推荐阅读更多精彩内容