r - 计算一个符号 Ryacas 表达式

Question

这是一个可重现的例子：

a <- 0.05
za.2 <- qnorm(1-a/2)
b <- 0.20
zb <- qnorm(1-b)

lambda12 <- -log(1/2)/12
lambda18 <- -log(1/2)/18
theta <- lambda18/lambda12
(d = round(4*(za.2+zb)^2/log(theta)^2))   

Tf<-36
library(Ryacas)
n <- Sym("n")

Solve(n/2*(2-exp(-lambda12*Tf)-exp(-lambda18*Tf))==d , n)

最后一行返回

expression(list(n == 382/1.625))

有没有办法提取商并将其分配给另一个变量（235.0769）？

Answer 1

G.Grothendieck 在评论中指出，您需要首先捕获要在下面操作的表达式：

soln <- Solve(n/2*(2-exp(-lambda12*Tf)-exp(-lambda18*Tf))==d , n) 
X <- yacas(soln)$text

然后，为了提取商，您可以利用许多 R 语言对象是或可以被强制为列表的事实。

 X <- expression(list(n == 382/1.625))
 res <- eval(X[[1]][[2]][[3]])
 res
 [1] 235.0769

以下只是说明了为什么该索引序列提取了表达式的正确部分：

as.list(X)
# [[1]]
# list(n == 382/1.625)

as.list(X[[1]])
# [[1]]
# list
# 
# [[2]]
# n == 382/1.625

as.list(X[[1]][[2]])
# [[1]]
# `==`
# 
# [[2]]
# n
# 
# [[3]]
# 382/1.625