2

我在一个学生研究职位，并且是 R 的新手。我问了一个类似的问题（发布在这里：MLE 问题）。我已经解决了最初的问题，但我遇到了更多关于此功能的问题。

我仍在使用此函数来尝试估计 theta[i]，在此处输入图像描述 其中每个其他变量都是当前已知的。

下面是我的代码：

  maxParam <- function(theta) {
    logl <- sum(for (i in 1:length(doses)) {
      sum(
        for (j in 1:LITTERS.M) {
          sum(
            for (k in 0:(litterResponses[i,j]-1)) {
              sum(log10(probabilityResponses[i] + k * theta[i]))
            }
            +
            for (k in 0:(litterSizes[i,j]-litterResponses[i,j]-1)) {
              sum(log10(1 - probabilityResponses[i] + k * theta[i]))
            }
            -
            for (k in 0:(litterSizes[i,j] - 1)) {
              sum(log10(1 + k * theta[i]))
            }
          )
        }
      )
    })

    return (-logl)
  }

  mle.fit <- mle(maxParam, start=list(theta=c(1,1,1,1,1,1)))
  print(mle.fit)

我被抛出的错误是：

错误：缺少参数“theta”，没有默认值

如果错误很愚蠢，我深表歉意，我对 R 知之甚少。

注意： 我使用 (1,1,1,1,1,1) 的向量作为 theta 的测试。这不是实际数据。剂量是对应于血清剂量水平的 6 向量。Litter Responses 是一个矩阵，描述了每窝每剂量对血清的反应。LitterSizes 是一个矩阵，描述了每窝每剂量的窝大小。LITTERS.M 是与血清接触的初始窝数。ProbabilityResponses 是一个向量，用于描述给定小鼠将受到血清影响的概率。

4

2 回答 2

3

该函数mle不接受初始起始值的向量。通过优化找到的每个参数都需要作为标量传递。将函数的声明更改为：

maxParam <- function(theta_1 = 1, theta_2 = 1, etc) {
    theta <- unlist(as.list(environment()))

    ... # rest of function follows
  }

whereetc表示根据需要将其替换theta_3 = 1为theta_4 = 1。mle然后可以使用以下方法调用该函数：

mle.fit <- mle(maxParam)
于 2016-09-15T01:35:34.330 回答
1

不久前我遇到了同样的问题。mle不接受向量参数，其中一样mle2。看看这个。

于 2016-10-24T03:46:03.217 回答