如何从 R 中的函数输出创建绘图

Question

我有一个函数可以计算基因 (x) 中表达水平向量的对数似然； 例如 -0.465、0.111、-0.437。 计算的西格玛值为 0.4717，a = 0.8861。

    loglike <- function(x,theta,sigma,a) {
      sum(log(theta/(2*a)*(x<a)*(x>-a) + (1- 
      theta)*dnorm(x,mean=0.0,sd=sigma)))
    }

我想将 theta 的可能值向量输入到函数中，并用每个 theta 值绘制函数的结果。 所以我想到了一个 for 循环来旋转函数中 theta 的所有值，以从 loglike 函数的输出创建一个值列表。 然后根据对数结果创建一个 theta 图。

   theta <- seq(0,1,0.01)  

   result_ <- c()
   for (i in 1:length(theta)) {
      result <- loglike(x,i,sigma,a)
      result_[i] <- result[i]
   }

但这只是给了我警告信息：

   Warning messages:
   1: In log(theta/(2 * a) * (x < a) * (x > -a) + (1 - theta) * 
   dnorm(x,  :
     NaNs produced

我知道该函数有效，因为当我一次使用一个函数时，它会返回值，例如

    loglike(x,0,sigma,a)
    > -333.454234

我对 R 作为一种语言不是很有信心，所以非常感谢任何帮助。

Answer 1

感谢您添加 x、a 和 sigma 的值：这些使我能够找到您的错误。

在这个循环中：

  for (i in 1:length(theta)) {
      result <- loglike(x,i,sigma,a)
      result_[i] <- result[i]
   }

您将i作为 theta 的值输入到 loglike 函数中。 请记住， i是从 1 到 101 的循环的索引号，但是您需要循环每个循环的 theta 值，因此它应该是theta[i] 。

此外，因为您的loglike函数只返回一个数字，所以您将一个数字存储到result ，然后尝试将result[i]存储到result_[i] 。 当循环超过 i = 1 时，您正在尝试移动一个不存在的数字。

这是一个完整的工作示例：

loglike <- function(x, theta, sigma, a)
{
  sum(log((theta / (2*a)) * (x < a) * (x > -a) +
          (1 - theta) * dnorm(x, mean = 0.0, sd = sigma)))
}

x       <-  0.111
sigma   <-  0.4717
a       <-  0.8861
theta   <- seq(0, 1, 0.01)
result  <- numeric()

for (i in 1:length(theta)) result[i] <- loglike(x, theta[i], sigma, a)