1

我想将公式传递给函数,但我无法从公式的字符串中删除双引号。我怎么能做到?谢谢。

clust_start=c(1,100)
i <- 1 
Lambexp <- 1:99
rank <- 1:99
alpha <- 1:99
d <- 1
b <- 1:2
ww=as.matrix(c(1:99), nrow=99)
expression.logdens <- paste(Lambexp[rank[(clust_start[i]:(clust_start[i + 1] - 1))]],         
"*exp(", alpha[(clust_start[i]:(clust_start[i + 1] - 1))], " + bb*", ww[(clust_start[i]:    
(clust_start[i + 1] - 1)), d], " - ", b[d]*ww[(clust_start[i]:(clust_start[i + 1] - 1)),     
d], ")", sep="")  

expr.formula <- paste(expression.logdens, sep="", collapse=" + ")

我尝试了 expr.formula <- as.expression(expr.formula) 和 as.name(expr.formula),它们没有用。

logdens <- function(bb){
    y <- expr.formula    
        return(y)
    }
logdens(1)
4

2 回答 2

4

你要这个:

logdens <- function(bb){
  eval(parse(text=expr.formula))
}

输出:

> logdens(1)
[1] 1.539673e+45
于 2012-02-09T00:40:16.290 回答
0

它就像as.formula(expr.formula)......一旦你修复你的公式实际上是有效的;)

于 2012-02-09T02:50:13.477 回答