《R语言初学指南》一2.3 找到R脚本中的错误

简介:

本节书摘来自异步社区《R语言初学指南》一书中的第2章,第2.3节,作者【美】Brian Dennis(布莱恩·丹尼斯),更多章节内容可以访问云栖社区“异步社区”公众号查看

2.3 找到R脚本中的错误

R语言初学指南
复杂项目中的R脚本会非常长。即使是R专家,也很少能一次性将其编写正确。脚本中的所有错误都是通过调试来修改的。

在脚本中称追踪错误或“bug”为“调试”。调试包括一些适用性检测工作。

下面在脚本中故意制造一个bug,由此来练习如何调试。首先来做一些清理工作。使控制台成为活动窗口,并输入下列命令:

>objects()

按照惯例,这里在展示控制台中的命令时,依然在其前面显示R提示符(“>”),但要记住,在控制台中不需要输入该提示符。不要忘记命令末尾的括号。在运行该命令后,就可在工作区内看到一串已用R命令建立好的“对象”:

[1] "a"       "b"       "k"       "kill.rate"  
[5] "m"       "moose.density"

这些对象已存储在R中,可供日后调用。现在输入下列命令:

>rm(a,b,k,kill.rate,m,moose.density)
>a

哎?找不到对象a了,其他几个对象也找不到。命令rm()表示“移除”,它会删除括号中的所有对象。由于用户可能会忘记之前在R会话中建立的对象,而导致后面的计算出现错误,这时命令objects()和rm()就会非常有用。

现在开始通过制造bug来搞破坏!回到R编辑器中,将脚本的第一个命令moose.density的首字母改为大写字母M,然后重新运行脚本。这时会弹出图形窗口,但其中没有图形显示——出现错误了。看一下控制台,找到与下列语句相似的内容(在不同版本的R中,出现的错误信息可能会稍有不同):

> Moose.density=c(.17,.23,.23,.26,.37,.42,.66,.80,1.11,1.30,1.37,
+ 1.41,1.73,2.49)
> kill.rate=c(.37,.47,1.90,2.04,1.12,1.74,2.78,1.85,1.88,1.96,
+ 1.80,2.44,2.81,3.75)
> plot(moose.density,kill.rate,type="p")
Error in plot(moose.density, kill.rate, type = "p") :
 object 'moose.density' not found
> m=2.5*(0:100)/100
> a=3.37
> b=0.47
> k=a*m/(b+m)
> points(m,k,type="l")
Error in plot.xy(xy.coords(x, y), type = type, ...) :
 plot.new has not been called yet

控制台中重现了脚本,其中还包含两条错误信息。这里有两点非常重要,需要注意。

首先,尽管没有修改plot()语句和points()语句,但它们下面出现了错误信息,而实际上有错误(大写字母M)的语句却没有任何提示。在脚本中,错误可以影响下面一系列命令出错。这里,定义向量Moose.density的命令是一个在R中完全正确的语句,它在R的内存中建立了对象Moose.density以便日后调用,在控制台中使用object()命令即可看到这一对象。问题出在绘图语句上。plot()语句使用的向量moose.density并不存在,故plot()语句无法进行绘图。接着由于没有元素可以添加到图形中,points()语句也就无法执行。

其次,在通常情况下,若可以定位并修正某个错误,则后面的一些错误也会被修正。因此,建议用户在策略上,要从上到下一个一个地修改错误。此外,在每修改完一个语句或一小部分R脚本后,就对它们进行测试,以防修改本身包含错误并对后续的语句产生影响。不过,即使再小心谨慎,也几乎没有很长的R脚本可在第一次就能完美地运行出来。

若要在脚本中检测并定位不明显的错误,需要对脚本从头开始进行一系列搜索及实验。在空白的工作区中,可选中代码并运行这部分脚本,从而很方便地完成这一任务。读者可以试一下:首先像刚才一样,在控制台中删除工作区内的所有对象,然后在R编辑器中选中第一个命令,使其高亮。在“Edit”的下拉菜单中(在微软的Windows中,若对高亮区点击右键,会弹出菜单),找到并点击“Run line selection”选项,即可运行脚本中的第一行命令。在控制台中检查一下,这条命令建立了哪些对象。

若该命令的运行情况与预期一致,则继续运行下一条语句。单独选中并运行下一条语句,然后再连同脚本中前面的语句一起运行一遍。

当检查到绘图语句时,错误信息就会出现。绘图命令在寻找一个不存在的对象。若在之前的检查中没有发现Moose.density与moose.density是两个不同的对象,这时也许能清晰地发现这一点。R的用户应当时刻注意变量名的大小写以及拼写错误问题,这是错误的一个主要来源。

找到并改正大写的M,然后运行脚本,即可成功绘图。但是等一下:过早地认为这是一个完美的脚本可能会导致失败。有时即使脚本可以很好地运行,其运行结果也可能错得离谱。当计算本身有概念性的错误时,这种情况就会发生。R会精确完成用户的指令,且当用户的指令有错误时,R也会按错误指令忠实地运行完成。

数据输入是否正确?作者在输入本章数据时,是否与第1章保持一致了呢?要想有效地使用R,仔细一点儿总不会错。

相关文章
|
8月前
|
存储 Go
速绘丨GO富集气泡图绘制方法,利用R语言ggplot2包快速绘制,完整脚本可重复绘图
速绘丨GO富集气泡图绘制方法,利用R语言ggplot2包快速绘制,完整脚本可重复绘图
|
8月前
|
Linux Python
彩蛋丨利用R语言脚本实现批量合并Excel表格,再也不用手动点来点去了!
彩蛋丨利用R语言脚本实现批量合并Excel表格,再也不用手动点来点去了!