注册 登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mobius-Band

 
 
 

日志

 
 

R语言进阶之5:表达式、数学公式与特殊符号  

2015-07-04 15:18:08|  分类: R |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
原文地址:http://developer.51cto.com/art/201306/396747.htm
  1. <- function(x) log(x) + sqrt(x) + x^(1/3) 
  2. lot(y, 1, 1000, main = expression(y == log(x) + sqrt(x) + sqrt(x, 3)), 
  3. lwd = 3col = "blue"

R语言进阶之五:表达式、数学公式与特殊符号 - xxx - xxx的博客

一、R语言的“表达式”

在R语言中,“表达式”的概念有狭义和广义两种意义。狭义的表达式指表达式(expression)类对象,由expression函数产生;而广义的的表达式既包含expression类,也包含R“语言”类(language)。expression和language是R语言中两种特殊数据类:

  1. getClass("expression") 
  2. ## Class "expression" [package "methods"] 
  3. ##  
  4. ## No Slots, prototype of class "expression" 
  5. ##  
  6. ## Extends: "vector" 
  7. getClass("language") 
  8. ## Virtual Class "language" [package "methods"] 
  9. ##  
  10. ## No Slots, prototype of class "name" 
  11. ##  
  12. ## Known Subclasses:  
  13. ## Class "name", directly 
  14. ## Class "call", directly 
  15. ## Class "{", directly 
  16. ## Class "if", directly 
  17. ## Class "<-", directly 
  18. ## Class "for", directly 
  19. ## Class "while", directly 
  20. ## Class "repeat", directly 
  21. ## Class "(", directly 
  22. ## Class ".name", by class "name", distance 2, with explicit coerce 

可以看到expression类由向量派生得到,而language类是虚拟类,它包括我们熟悉的程序控制关键词/符号和name、call 子类。

二、产生“表达式”的函数

虽然我们在R终端键入的任何有效语句都是表达式,但这些表达式在输入后即被求值(evaluate)了,获得未经求值的纯粹“表达式”就要使用函数。下面我们从函数参数和返回值两方面了解expression、quote、bquote和substitute这几个常用函数。

1、expression 函数

expression函数可以有一个或多个参数,它把全部参数当成一个列表,每个参数都被转成一个表达式向量,所以它的返回值是表达式列表,每个元素都是表达式类型对象,返回值的长度等于参数的个数:

  1. (ex <- expression(x = 1, 1 + sqrt(a))) 
  2. ## expression(x = 1, 1 + sqrt(a)) 
  3. length(ex) 
  4. ## [1] 2 
  5. ex[1] 
  6. ## expression(x = 1
  7. mode(ex[1]) 
  8. ## [1] "expression" 
  9. typeof(ex[1]) 
  10. ## [1] "expression" 
  11. ex[2] 
  12. ## expression(1 + sqrt(a)) 
  13. mode(ex[2]) 
  14. ## [1] "expression" 
  15. typeof(ex[2]) 
  16. ## [1] "expression" 

因为expression函数把参数当成列表处理,所以等号‘=’两边的表达式要符合R语言列表元素的书写规则,否则出错,比如:

  1. expression(x+11=1) 

2、quote函数

quote函数只能有一个参数。quote函数的返回值一般情况下是call类型,表达式参数是单个变量的话返回值就是name类型,如果是常量那么返回值的存储模式就和相应常量的模式相同:

  1. (cl <- quote(1 + sqrt(a) + b^c)) 
  2. ## 1 + sqrt(a) + b^c 
  3. mode(cl) 
  4. ## [1] "call" 
  5. typeof(cl) 
  6. ## [1] "language" 
  7. (cl <- quote(a)) 
  8. ## a 
  9. mode(cl) 
  10. ## [1] "name" 
  11. typeof(cl) 
  12. ## [1] "symbol" 
  13. (cl <- quote(1)) 
  14. ## [1] 1 
  15. mode(cl) 
  16. ## [1] "numeric" 
  17. typeof(cl) 
  18. ## [1] "double" 

quote返回值如果是name或常量类型,它的长度就是1;如果是call类型,返回值长度就与函数/运算符的参数个数n对应,长度等于n+1,多出的长度1是函数/符号名。

  1. length(quote(a))  #name或常量类型,返回值长度为1 
  2. ## [1] 1 
  3. length(quote(!a))  #单目运算符,返回值长度为2 
  4. ## [1] 2 
  5. length(quote(-b))  #单目运算符,返回值长度为2 
  6. ## [1] 2 
  7. length(quote(a + b))  #双目运算符,返回值长度为3 
  8. ## [1] 3 
  9. length(quote((a + b) * c))  #多个运算符只算优先级最低的一个 
  10. ## [1] 3 

3、bquote 和 substitute 函数

如果不使用环境变量或环境变量参数,bquote 和 substitute 函数得到的结果与quote函数相同。

  1. bquote(1 + sqrt(a) + b^c) == quote(1 + sqrt(a) + b^c) 
  2. ## [1] TRUE 
  3. substitute(1 + sqrt(a) + b^c) == quote(1 + sqrt(a) + b^c) 
  4. ## [1] TRUE 

但是bquote 和 substitute 函数可以在表达式中使用变量,变量的值随运行进程而被替换。bquote 和 substitute 函数变量替换的方式不一样,bquote函数中需要替换的变量用 .( ) 引用,substitute函数中需要替换的变量用列表参数方式给出。除了这一点,bquote 和 substitute 函数没有差别:

  1. a <- 3 
  2. b <- 2 
  3. (bq <- bquote(y == sqrt(.(a), .(b)))) 
  4. ## y == sqrt(3, 2) 
  5. (ss <- substitute(y == sqrt(a, b), list(a = 3, b = 2))) 
  6. ## y == sqrt(3, 2) 
  7. bq == ss 
  8. ## [1] TRUE 

搞出两个功能完全一样的函数不算很奇怪,R语言里面太多了,可能是照顾不同使用习惯的人们吧。bquote函数的帮助档说这个函数类似于LISP的backquote宏,对于像我这样的LISP盲,使用substitute函数好一些。 substitute函数的典型用途是替换表达式中的变量,如果我们希望在表达式中使用变量并且希望这些变量在运行过程中做出相应改变,就可以使用substitute函数。

  1. par(mar = rep(0.1, 4), cex = 2
  2. plot.new() 
  3. plot.window(c(0, 10), c(0, 1)) 
  4. for (i in 1:9) text(i, 0.5, substitute(sqrt(x, a), list(a = i + 1))) 

R语言进阶之五:表达式、数学公式与特殊符号 - xxx - xxx的博客 

4、parse 函数

parse函数用于从文件读取文本作为表达式,返回的值是expression类型,这函数也很有用。后面有例子。

三、表达式规则与paste函数:

与在R终端直接输入的表达式不一样,expression、quote、bquote和substitute等函数对参数中的(变量)名称都不做任何检查:

  1. <- 1 
  2. x + "x" 
  3. ## Error: 二进列运算符中有非数值参数 
  4. expression(x + "x") 
  5. ## expression(x + "x") 
  6. quote(x + "x") 
  7. ## x + "x" 

但R要检查表达式中的运算符,不符合运算符使用规则的表达式将出错:

  1. expression(x + +++y) 
  2. ## expression(x + +++y) 
  3. expression(x - ---y) 
  4. ## expression(x - ---y) 
  5. # expression(x****y) (Not run) expression(x////y) (Not run) 
  6. # expression(1<=x<=4) (Not run) 
  7. quote(x + +++y) 
  8. ## x + +++y 
  9. quote(x - ---y) 
  10. ## x - ---y 
  11. # quote(x****y) (Not run) quote(x////y) (Not run) quote(1<=x<=4) (Not run) 

+ - 运算连续使用不出错是因为它们还可以当成求正/负值运算的符号。 在表达式产生函数中使用paste函数可以解决这样的问题。在这种条件下,paste对参数的处理方式和表达式产生函数一样,检查运算符但不检查变量名。用NULL作为运算符的参数可以获得意外的效果:

  1. ex <- expression(paste(x, "////", y)) 
  2. cl <- quote(paste(x, "****", y)) 
  3. par(mar = rep(0.1, 4), cex = 2
  4. plot.new() 
  5. plot.window(c(0, 1.2), c(0, 1)) 
  6. text(0.2, 0.5, ex) 
  7. text(0.6, 0.5, cl) 
  8. cl <- quote(paste(1 <= x, NULL <= 4)) 
  9. text(1, 0.5, cl) 

R语言进阶之五:表达式、数学公式与特殊符号 - xxx - xxx的博客 

四、R绘图函数对文本参数中的表达式的处理

quote, bquote 和 substitute 的返回值有三种类型call, name 和 常量,事实上expression 函数的结果最终也是这三种类型。因为expression函数的结果是expression列表,我们取列表元素的值检查看看:

 
  1. (ex <- expression(1 + sqrt(x), x, 1)) 
  2. ## expression(1 + sqrt(x), x, 1) 
  3. ex[[1]] 
  4. ## 1 + sqrt(x) 
  5. mode(ex[[1]]) 
  6. ## [1] "call" 
  7. typeof(ex[[1]]) 
  8. ## [1] "language" 
  9. ex[[2]] 
  10. ## x 
  11. mode(ex[[2]]) 
  12. ## [1] "name" 
  13. typeof(ex[[2]]) 
  14. ## [1] "symbol" 
  15. ex[[3]] 
  16. ## [1] 1 
  17. mode(ex[[3]]) 
  18. ## [1] "numeric" 
  19. typeof(ex[[3]]) 
  20. ## [1] "double" 
确实是这样。所以绘图函数对文本参数中的表达式处理就有三种情况。先看看处理结果:
  1. par(mar = rep(0.1, 4), cex = 2
  2. plot.new() 
  3. plot.window(c(0, 1.2), c(0, 1)) 
  4. text(0.2, 0.5, ex[1]) 
  5. text(0.6, 0.5, ex[2]) 
  6. text(1, 0.5, ex[3]) 
R语言进阶之五:表达式、数学公式与特殊符号 - xxx - xxx的博客

name 和常量类型都很简单,直接输出文本,而call类型就不好判断了。我们前面说过call类型返回值的长度与函数/运算符的参数个数有关。这是怎么体现的呢?由于文本参数最终得到的是文本,我们用as.character函数来看看:

  1. as.character(quote(x - y)) 
  2. ## [1] "-" "x" "y" 
  3. as.character(quote(1 - x + y)) 
  4. ## [1] "+"     "1 - x" "y" 
  5. as.character(quote((1 + x) * y)) 
  6. ## [1] "*"       "(1 + x)" "y" 
  7. as.character(quote(!a)) 
  8. ## [1] "!" "a" 
  9. as.character(quote(sqrt(x))) 
  10. ## [1] "sqrt" "x" 

转换成字符串向量后排在第一位的是运算符或函数名称,后面是参数(如果参数中还有运算符或函数名,R还会对其进行解析)。运算符和函数是相同的处理方式。事实上,在R语言中,所有运算符(包括数学运算符和逻辑运算符)都是函数,你可以用函数的方式使用运算符:


  1. 2 + 4 
  2. ## [1] 6 
  3. 2 - 4 
  4. ## [1] -2 
  5. <= 4 
  6. ## [1] TRUE 
  7. >= 4 
  8. ## [1] FALSE 

R绘图函数对表达式中包含的函数名和它们的参数首先应用Tex文本格式化规则进行处理,这种规则的具体情况可以使用 ?plotmath 进行查看,主要是一些数学公式和符号的表示方法。把这个说明文档中字符串拷贝到maths.txt文件中并保存到当前工作目录后可以用下面的代码做出后面的表格:

  1. ex <- parse("maths.txt") 
  2. labs <- readLines("maths.txt") 
  3. <- length(ex) 
  4. par(mar = rep(0.1, 4), cex = 0.8) 
  5. plot.new() 
  6. plot.window(c(0, 8), c(0, n/4)) 
  7. <- seq(n/4, by = -1, length = n/4) 
  8. <- seq(0.1, by = 2length = 4
  9. xy <- expand.grid(x, y) 
  10. text(xy, labs, adj = c(0, 0.5)) 
  11. xy <- expand.grid(x + 1.3, y) 
  12. text(xy, ex, adj = c(0, 0.5), col = "blue"
  13. box(lwd = 2
  14. abline(v = seq(1.3, by = 2length = 4), lty = 3
  15. abline(v = seq(2, by = 2length = 3), lwd = 1.5) 

R语言进阶之五:表达式、数学公式与特殊符号 - xxx - xxx的博客

 表中奇数列是字符串(表达式),偶数列(蓝色)是Tex格式化的图形。除了上表列出的规则外还有一些拉丁文和希腊文符号,可以在表达式中用 symbol 函数或名称(如alpha)等表示,用到时自己去找吧。 如果函数名(包括运算符)有对应的Tex格式化规则,函数名和参数都按规则进行图形绘制;如果没有,就当成是R语言普通函数:

  1. ex <- expression(sqrt(x), x + y, x^2, x %in% A, x <= y, mean(x, y, z), x | y,  
  2.     x & y) 
  3. <- length(ex) 
  4. par(mar = rep(0.1, 4), cex = 1.5) 
  5. col <- c("red", "blue") 
  6. plot.new() 
  7. plot.window(c(0, n), c(0, 1)) 
  8. for (i in 1:n) text(i - 0.5, 0.5, ex[i], colcol = col[i%%2 + 1]) 

R语言进阶之五:表达式、数学公式与特殊符号 - xxx - xxx的博客

上面例子中前5种运算函数都是有对应数学符号的,所以它出的图(符号和顺序)与数学习惯一致,后三种运算函数没有对应数学符号,所以用普通函数方式(函数名在前,参数在括号内用逗号分隔)出图。其他还有一些琐碎的规则,自己找找吧。

  评论这张
 
阅读(5)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018