pbootcms网站模板|日韩1区2区|织梦模板||网站源码|日韩1区2区|jquery建站特效-html5模板网

詳解R語言中的表達式、數學公式、特殊符號

這篇文章主要介紹了詳解R語言中的表達式、數學公式、特殊符號,本文給大家介紹的非常詳細,對大家的學習或工作具有一定的參考借鑒價值,需要的朋友可以參考下

  在R語言的繪圖函數中,如果文本參數是合法的R語言表達式,那么這個表達式就被用Tex類似的規則進行文本格式化。


y <- function(x) (exp(-(x^2)/2))/sqrt(2*pi)

plot(y, -5, 5, main = expression(f(x) == frac(1,sqrt(2*pi))*e^(-frac(x^2,2))), lwd = 3, col = "blue")

這里寫圖片描述


library(ggplot2)

x <- seq(0, 2*pi, by = 0.01)
y <- sin(x)
data <- data.frame(x, y)
p <- ggplot(data, aes(x, y)) + geom_line()

p + geom_area(fill = 'blue', alpha = 0.3) +
 scale_x_continuous(breaks = c(0, pi, 2*pi), labels = c('0', expression(pi), expression(2*pi))) +
 geom_text(parse = T, aes(x = pi/2,y = 0.3, label = 'integral(sin(x)*dx, 0, pi)'))

這里寫圖片描述

##一、R語言的“表達式”

  在R語言中,“表達式”的概念有狹義和廣義兩種意義。狹義的表達式指表達式(expression)類對象,由expression函數產生;而廣義的的表達式既包含expression類,也包含R“語言”類(language)。expression和language是R語言中兩種特殊數據類:

getClass(“expression”)


# Class "expression" [package "methods"]

# No Slots, prototype of class "expression"

# Extends: "vector"

getClass(“language”)


# Virtual Class "language" [package "methods"]

# No Slots, prototype of class "name"

# Known Subclasses: 
# Class "name", directly
# Class "call", directly
# Class "{", directly
# Class "if", directly
# Class "<-", directly
# Class "for", directly
# Class "while", directly
# Class "repeat", directly
# Class "(", directly
# Class ".name", by class "name", distance 2, with explicit coerce

可以看到expression類由向量派生得到,而language類是虛擬類,它包括我們熟悉的程序控制關鍵詞/符號和name、call 子類。

##二、產生“表達式”的函數

雖然我們在R終端鍵入的任何有效語句都是表達式,但這些表達式在輸入后即被求值(evaluate)了,獲得未經求值的純粹“表達式”就要使用函數。下面我們從函數參數和返回值兩方面了解expression、quote、bquote和substitute這幾個常用函數。

####1、expression 函數

expression函數可以有一個或多個參數,它把全部參數當成一個列表,每個參數都被轉成一個表達式向量,所以它的返回值是表達式列表,每個元素都是表達式類型對象,返回值的長度等于參數的個數:


(ex <- expression(x = 1, 1 + sqrt(a)))

## expression(x = 1, 1 + sqrt(a))

length(ex)

## [1] 2

ex[1]

## expression(x = 1)

mode(ex[1])

## [1] "expression"

typeof(ex[1])

## [1] "expression"

ex[2]

## expression(1 + sqrt(a))

mode(ex[2])

## [1] "expression"

typeof(ex[2])

## [1] "expression" 

因為expression函數把參數當成列表處理,所以等號‘='兩邊的表達式要符合R語言列表元素的書寫規則,否則出錯,比如:


expression(x+11=1)

####2、quote函數

quote函數只能有一個參數。quote函數的返回值一般情況下是call類型,表達式參數是單個變量的話返回值就是name類型,如果是常量那么返回值的存儲模式就和相應常量的模式相同:


(cl <- quote(1 + sqrt(a) + b^c))
## 1 + sqrt(a) + b^c

mode(cl)
## [1] "call"

typeof(cl)
## [1] "language"

(cl <- quote(a))
## a

mode(cl)
## [1] "name"

typeof(cl)
## [1] "symbol"

(cl <- quote(1))
## [1] 1

mode(cl)
## [1] "numeric"

typeof(cl)
## [1] "double" 

quote返回值如果是name或常量類型,它的長度就是1;如果是call類型,返回值長度就與函數/運算符的參數個數n對應,長度等于n+1,多出的長度1是函數/符號名。


length(quote(a)) #name或常量類型,返回值長度為1
## [1] 1

length(quote(!a)) #單目運算符,返回值長度為2
## [1] 2

length(quote(-b)) #單目運算符,返回值長度為2
## [1] 2

length(quote(a + b)) #雙目運算符,返回值長度為3
## [1] 3

length(quote((a + b) * c)) #多個運算符只算優先級最低的一個
## [1] 3 

####3、bquote 和 substitute 函數

如果不使用環境變量或環境變量參數,bquote 和 substitute 函數得到的結果與quote函數相同。


bquote(1 + sqrt(a) + b^c) == quote(1 + sqrt(a) + b^c)
## [1] TRUE

substitute(1 + sqrt(a) + b^c) == quote(1 + sqrt(a) + b^c)
## [1] TRUE 

但是bquote 和 substitute 函數可以在表達式中使用變量,變量的值隨運行進程而被替換。bquote 和 substitute 函數變量替換的方式不一樣,bquote函數中需要替換的變量用 .( ) 引用,substitute函數中需要替換的變量用列表參數方式給出。除了這一點,bquote 和 substitute 函數沒有差別:


a <- 3

b <- 2

(bq <- bquote(y == sqrt(.(a), .(b))))
## y == sqrt(3, 2) 

(ss <- substitute(y == sqrt(a, b), list(a = 3, b = 2)))
## y == sqrt(3, 2) 

bq == ss
## [1] TRUE 

搞出兩個功能完全一樣的函數不算很奇怪,R語言里面太多了,可能是照顧不同使用習慣的人們吧。bquote函數的幫助檔說這個函數類似于LISP的backquote宏,對于像我這樣的LISP盲,使用substitute函數好一些。 substitute函數的典型用途是替換表達式中的變量,如果我們希望在表達式中使用變量并且希望這些變量在運行過程中做出相應改變,就可以使用substitute函數。


par(mar = rep(0.1, 4), cex = 2)

plot.new()

plot.window(c(0, 10), c(0, 1))

for (i in 1:9) text(i, 0.5, substitute(sqrt(x, a), list(a = i + 1))) 

這里寫圖片描述

####4、parse 函數

parse函數用于從文件讀取文本作為表達式,返回的值是expression類型,這函數也很有用。后面有例子。


x <- 1

x + "x"

## Error: 二進列運算符中有非數值參數

expression(x + "x")

## expression(x + "x")

quote(x + "x")

## x + "x" 

但R要檢查表達式中的運算符,不符合運算符使用規則的表達式將出錯:


expression(x + +++y)
## expression(x + +++y) 

expression(x - ---y)
## expression(x - ---y) 
## expression(x****y) (Not run) expression(xy) (Not run) 
## expression(1<=x<=4) (Not run) 

quote(x + +++y)
## x + +++y 

quote(x - ---y)
## x - ---y 
## quote(x****y) (Not run) quote(xy) (Not run) quote(1<=x<=4) (Not run) 

+ - 運算連續使用不出錯是因為它們還可以當成求正/負值運算的符號。 在表達式產生函數中使用paste函數可以解決這樣的問題。在這種條件下,paste對參數的處理方式和表達式產生函數一樣,檢查運算符但不檢查變量名。用NULL作為運算符的參數可以獲得意外的效果:


ex <- expression(paste(x, "", y))

cl <- quote(paste(x, "****", y))

par(mar = rep(0.1, 4), cex = 2)

plot.new()

plot.window(c(0, 1.2), c(0, 1))

text(0.2, 0.5, ex)

text(0.6, 0.5, cl)

cl <- quote(paste(1 <= x, NULL <= 4))

text(1, 0.5, cl) 

這里寫圖片描述

##三、R繪圖函數對文本參數中的表達式的處理

quote, bquote 和 substitute 的返回值有三種類型call, name 和 常量,事實上expression 函數的結果最終也是這三種類型。因為expression函數的結果是expression列表,我們取列表元素的值檢查看看:


(ex <- expression(1 + sqrt(x), x, 1))
## expression(1 + sqrt(x), x, 1)

ex[[1]]
## 1 + sqrt(x)

mode(ex[[1]])
## [1] "call"

typeof(ex[[1]])
## [1] "language"

ex[[2]]
## x

mode(ex[[2]])
## [1] "name"

typeof(ex[[2]])
## [1] "symbol"

ex[[3]]
## [1] 1

mode(ex[[3]])
## [1] "numeric"

typeof(ex[[3]])
## [1] "double" 

確實是這樣。所以繪圖函數對文本參數中的表達式處理就有三種情況。先看看處理結果:


par(mar = rep(0.1, 4), cex = 2)

plot.new()

plot.window(c(0, 1.2), c(0, 1))

text(0.2, 0.5, ex[1])

text(0.6, 0.5, ex[2])

text(1, 0.5, ex[3])

這里寫圖片描述

name 和常量類型都很簡單,直接輸出文本,而call類型就不好判斷了。我們前面說過call類型返回值的長度與函數/運算符的參數個數有關。這是怎么體現的呢?由于文本參數最終得到的是文本,我們用as.character函數來看看:


as.character(quote(x - y))
## [1] "-" "x" "y"

as.character(quote(1 - x + y))
## [1] "+" "1 - x" "y"

as.character(quote((1 + x) * y))
## [1] "*" "(1 + x)" "y"

as.character(quote(!a))
## [1] "!" "a"

as.character(quote(sqrt(x)))
## [1] "sqrt" "x" 

轉換成字符串向量后排在第一位的是運算符或函數名稱,后面是參數(如果參數中還有運算符或函數名,R還會對其進行解析)。運算符和函數是相同的處理方式。事實上,在R語言中,所有運算符(包括數學運算符和邏輯運算符)都是函數,你可以用函數的方式使用運算符:


2 + 4
## [1] 6

2 - 4
## [1] -2

2 <= 4
## [1] TRUE

2 >= 4
## [1] FALSE 

R繪圖函數對表達式中包含的函數名和它們的參數首先應用Tex文本格式化規則進行處理,這種規則的具體情況可以使用 ?plotmath 進行查看,主要是一些數學公式和符號的表示方法。把這個說明文檔中字符串拷貝到maths.txt文件中并保存到當前工作目錄后可以用下面的代碼做出后面的表格:


ex <- parse("maths.txt")

labs <- readLines("maths.txt")

n <- length(ex)

par(mar = rep(0.1, 4), cex = 0.8)

plot.new()

plot.window(c(0, 8), c(0, n/4))

y <- seq(n/4, by = -1, length = n/4)

x <- seq(0.1, by = 2, length = 4)

xy <- expand.grid(x, y)

text(xy, labs, adj = c(0, 0.5))

xy <- expand.grid(x + 1.3, y)

text(xy, ex, adj = c(0, 0.5), col = "blue")

box(lwd = 2)

abline(v = seq(1.3, by = 2, length = 4), lty = 3)

abline(v = seq(2, by = 2, length = 3), lwd = 1.5) 

這里寫圖片描述

右鍵查看圖片,瀏覽大圖

表中奇數列是字符串(表達式),偶數列(藍色)是Tex格式化的圖形。除了上表列出的規則外還有一些拉丁文和希臘文符號,可以在表達式中用 symbol 函數或名稱(如alpha)等表示,用到時自己去找吧。 如果函數名(包括運算符)有對應的Tex格式化規則,函數名和參數都按規則進行圖形繪制;如果沒有,就當成是R語言普通函數:


ex <- expression(sqrt(x), x + y, x^2, x %in% A, x <= y, mean(x, y, z), x | y, x & y)

n <- length(ex)

par(mar = rep(0.1, 4), cex = 1.5)

col <- c("red", "blue")

plot.new()

plot.window(c(0, n), c(0, 1))

for (i in 1:n) text(i - 0.5, 0.5, ex[i], col = col[i%%2 + 1]) 

這里寫圖片描述

上面例子中前5種運算函數都是有對應數學符號的,所以它出的圖(符號和順序)與數學習慣一致,后三種運算函數沒有對應數學符號,所以用普通函數方式(函數名在前,參數在括號內用逗號分隔)出圖。其他還有一些瑣碎的規則,自己找找吧。

到此這篇關于詳解R語言中的表達式、數學公式、特殊符號的文章就介紹到這了,更多相關R語言表達式、數學公式、特殊符號內容請搜索html5模板網以前的文章希望大家以后多多支持html5模板網!

【網站聲明】本站部分內容來源于互聯網,旨在幫助大家更快的解決問題,如果有圖片或者內容侵犯了您的權益,請聯系我們刪除處理,感謝您的支持!

相關文檔推薦

主站蜘蛛池模板: 加中寰球移民官网-美国移民公司,移民机构,移民中介,移民咨询,投资移民 | 志高装潢官网-苏州老房旧房装修改造-二手房装修翻新 | 电脑刺绣_绣花厂家_绣花章仔_织唛厂家-[源欣刺绣]潮牌刺绣打版定制绣花加工厂家 | 郑州爱婴幼师学校_专业幼师培训_托育师培训_幼儿教育培训学校 | 杭州门窗厂家_阳光房_包阳台安装电话-杭州窗猫铝合金门窗 | 耐酸碱泵-自吸耐酸碱泵型号「品牌厂家」立式耐酸碱泵价格-昆山国宝过滤机有限公司首页 | 防爆电机生产厂家,YBK3电动机,YBX3系列防爆电机,YBX4节防爆电机--河南省南洋防爆电机有限公司 | 酒水灌装机-白酒灌装机-酒精果酒酱油醋灌装设备_青州惠联灌装机械 | 超声波破碎仪-均质乳化机(供应杭州,上海,北京,广州,深圳,成都等地)-上海沪析实业有限公司 | 我车网|我关心的汽车资讯_汽车图片_汽车生活! | 新疆十佳旅行社_新疆旅游报价_新疆自驾跟团游-新疆中西部国际旅行社 | 天助网 - 中小企业全网推广平台_生态整合营销知名服务商_天助网采购优选 | 不锈钢螺丝,不锈钢螺栓,不锈钢标准件-江苏百德特种合金有限公司 交变/复合盐雾试验箱-高低温冲击试验箱_安奈设备产品供应杭州/江苏南京/安徽马鞍山合肥等全国各地 | 知名电动蝶阀,电动球阀,气动蝶阀,气动球阀生产厂家|价格透明-【固菲阀门官网】 | 德州万泰装饰 - 万泰装饰装修设计软装家居馆 | 医用酒精_84消毒液_碘伏消毒液等医用消毒液-漓峰消毒官网 | 制丸机,小型中药制丸机,全自动制丸机价格-甘肃恒跃制药设备有限公司 | 全自动实验室洗瓶机,移液管|培养皿|进样瓶清洗机,清洗剂-广州摩特伟希尔机械设备有限责任公司 | 耐高温硅酸铝板-硅酸铝棉保温施工|亿欧建设工程 | 国资灵活用工平台_全国灵活用工平台前十名-灵活用工结算小帮手 | 南京泽朗生物科技有限公司-液体饮料代加工_果汁饮料代加工_固体饮料代加工 | 浩方智通 - 防关联浏览器 - 跨境电商浏览器 - 云雀浏览器 | 旋转/数显粘度计-运动粘度测定仪-上海平轩科学仪器 | 不锈钢电动球阀_气动高压闸阀_旋塞疏水调节阀_全立阀门-来自温州工业阀门巨头企业 | 家庭教育吧-在线家庭教育平台,专注青少年家庭教育 | 中国品牌排名投票_十大品牌榜单_中国著名品牌【中国品牌榜】 | 上海公司注册-代理记账-招投标审计-上海昆仑扇财税咨询有限公司 上海冠顶工业设备有限公司-隧道炉,烘箱,UV固化机,涂装设备,高温炉,工业机器人生产厂家 | 沉降天平_沉降粒度仪_液体比重仪-上海方瑞仪器有限公司 | 艺术生文化课培训|艺术生文化课辅导冲刺-济南启迪学校 | 耐酸碱胶管_耐腐蚀软管总成_化学品输送软管_漯河利通液压科技耐油耐磨喷砂软管|耐腐蚀化学软管 | 青岛侦探调查_青岛侦探事务所_青岛调查事务所_青岛婚外情取证-青岛狄仁杰国际侦探公司 | 湖南教师资格网-湖南教师资格证考试网 | 防爆电机_ybx3系列电机_河南省南洋防爆电机有限公司 | 不锈钢复合板厂家_钛钢复合板批发_铜铝复合板供应-威海泓方金属复合材料股份有限公司 | 耐磨焊丝,堆焊焊丝,耐磨药芯焊丝,碳化钨焊丝-北京耐默公司 | 杰福伦_磁致伸缩位移传感器_线性位移传感器-意大利GEFRAN杰福伦-河南赉威液压科技有限公司 | 电缆故障测试仪_电缆故障定位仪_探测仪_检测仪器_陕西意联电气厂家 | 碳化硅,氮化硅,冰晶石,绢云母,氟化铝,白刚玉,棕刚玉,石墨,铝粉,铁粉,金属硅粉,金属铝粉,氧化铝粉,硅微粉,蓝晶石,红柱石,莫来石,粉煤灰,三聚磷酸钠,六偏磷酸钠,硫酸镁-皓泉新材料 | 真空冷冻干燥机_国产冻干机_冷冻干燥机_北京四环冻干 | 不锈钢丸厂家,铝丸,铸钢丸-淄博智源铸造材料有限公司 | 不发火防静电金属骨料_无机磨石_水泥自流平_修补砂浆厂家「圣威特」 |