Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' — when running a nested function

huangapple
go评论236阅读模式
英文:

Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' -- when running a nested function

问题

问题：运行get_pi_ij()时出现错误：

Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any'

Called from: as.vector(data)

这个函数的第一步是将生成的alphas和beta_prelims转换为与c匹配的矩阵，以便它们可以一起计算。这是出现问题的地方，我无法找出问题所在。如果我在前面的函数中使用<<-将alphas和betas保存到全局环境中，并将它们替换为在有问题的函数中的那些，那么它就可以正常工作。因此，我认为问题可能与我在矩阵创建中调用函数的方式有关。

  1. get_pi_ij <- function() {
  2.   alphas <-
  3.     matrix(get_alpha(),
  4.            nrow = length(get_alpha()),
  5.            ncol = length(get_alpha()),
  6.            byrow = FALSE)
  7.   betas <-
  8.     matrix(get_beta_prelim(),
  9.            nrow = length(get_beta_prelim()),
  10.            ncol = length(get_beta_prelim()),
  11.            byrow = TRUE) 
  12.   pi_ij <- exp(alphas + betas + gamma * c)
  13.   return(pi_ij)
  14. }
  15. get_pi_ij()

我添加了完整的代码，因为它不太长，前面的部分只是定义。这样更容易测试。

一切都正常工作，如预期那样

  1. size <- 18 
  2. gamma <- -0.07 
  4. c <- structure(c(0, 4, 8, 12, 16, 4, 4, 8, 12, 16, 8, 8, 8, 12, 16, 
  5.             12, 12, 16, 16, 4, 0, 4, 8, 12, 8, 4, 4, 8, 12, 12, 8, 8, 12, 
  6.             16, 12, 12, 16, 16, 8, 4, 0, 4, 8, 12, 8, 4, 8, 12, 16, 12, 8, 
  7.             12, 16, 16, 12, 16, 16, 12, 8, 4, 0, 4, 12, 8, 4, 4, 8, 16, 12, 
  8.             8, 8, 12, 16, 12, 12, 16, 16, 12, 8, 4, 0, 16, 12, 8, 4, 4, 16, 
  9.             12, 8, 8, 8, 16, 12, 12, 16, 4, 8, 12, 12, 16, 0, 4, 8, 12, 16, 
  10.             4, 4, 8, 12, 16, 8, 8, 12, 12, 4, 4, 8, 8, 12, 4, 0, 4, 8, 12, 
  11.             8, 4, 4, 8, 12, 8, 8, 12, 12, 8, 4, 4, 4, 8, 8, 4, 0, 4, 8, 12, 
  12.             8, 4, 8, 12, 12, 8, 12, 12, 12, 8, 8, 4, 4, 12, 8, 4, 0, 4, 12, 
  13.             8, 4, 4, 8, 12, 8, 8, 12, 16, 12, 12, 8, 4, 16, 12, 8, 4, 0, 
  14.             16, 12, 8, 4, 4, 12, 8, 8, 12, 8, 12, 16, 16, 16, 4, 8, 12, 12, 
  15.             16, 0, 4, 8, 12, 16, 4, 8, 12, 8, 8, 8, 12, 12, 12, 4, 4, 8, 
  16.             8, 12, 4, 0, 4, 8, 12, 4, 4, 8, 8, 8, 8, 8, 8, 8, 8, 4, 4, 4, 
  17.             8, 8, 4, 0, 4, 8, 8, 4, 8, 8, 12, 12, 12, 8, 8, 12, 8, 8, 4, 
  18.             4, 12, 8, 4, 0, 4, 8, 4, 4, 8, 16, 16, 16, 12, 8, 16, 12, 12, 
  19.             8, 4, 16, 12, 8, 4, 0, 12, 8, 4, 8, 12, 12, 16, 16, 16, 8, 8, 
  20.             12, 12, 12, 4, 4, 8, 8, 12, 0, 4, 8, 4, 12, 12, 12, 12, 12, 8, 
  21.             8, 8, 8, 8, 8, 4, 4, 4, 8, 4, 0, 4, 4, 16, 16, 16, 12, 12, 12, 
  22.             12, 12, 8, 8, 12, 8, 8, 4, 4, 8, 4, 0, 4, 16, 16, 16, 16, 16, 
  23.             12, 12, 12, 12, 12, 8, 8, 8, 8, 8, 4, 4, 4, 0), 
  24.           .Dim = c(19L, 19L), 
  25.           .Dimnames = list(NULL, c("X1", "X
  27. <details>
  28. <summary>英文:</summary>
  30. Problem: running `get_pi_ij()` gives the error: 
  32. `Error in as.vector(x, mode) : cannot coerce type &#39;closure&#39; to vector of type &#39;any&#39;`
  34. `Called from: as.vector(data)`
  36. The first thing this function does is to make the resulting `alphas` and `beta_prelims` into matrixes that match `c` so that they can be calculated together. This is where something goes wrong, and I have not been able to figure out what. If I use `&lt;&lt;-` to save the `alphas` and `betas` to the global environment in the prior functions for alphas and betas and replace that with those in the faulty function, it works. So I assume it has to do with how I call the functions inside the matrix creation. 
  38.     get_pi_ij &lt;- function() {
  39.       alphas &lt;-
  40.         matrix(get_alpha(),
  41.                nrow = length(get_alpha()),
  42.                ncol = length(get_alpha()),
  43.                byrow = FALSE)
  44.       betas &lt;-
  45.         matrix(get_beta_prelim,
  46.                nrow = length(get_beta_prelim()),
  47.                ncol = length(get_beta_prelim()),
  48.                byrow = TRUE) 
  49.       pi_ij &lt;- exp(alphas + betas + gamma * c)
  50.       return(pi_ij)
  51.     }
  52.     get_pi_ij()
  54. I added the full code cause it&#39;s not too long and the first parts are just definitions. Makes it easier to test it. 
  56. **Everything up to the final function works as it is supposed to**
  58.     size &lt;- 18 
  59.     gamma &lt;- -0.07 
  61.     c &lt;- structure(c(0, 4, 8, 12, 16, 4, 4, 8, 12, 16, 8, 8, 8, 12, 16, 
  62.                 12, 12, 16, 16, 4, 0, 4, 8, 12, 8, 4, 4, 8, 12, 12, 8, 8, 12, 
  63.                 16, 12, 12, 16, 16, 8, 4, 0, 4, 8, 12, 8, 4, 8, 12, 16, 12, 8, 
  64.                 12, 16, 16, 12, 16, 16, 12, 8, 4, 0, 4, 12, 8, 4, 4, 8, 16, 12, 
  65.                 8, 8, 12, 16, 12, 12, 16, 16, 12, 8, 4, 0, 16, 12, 8, 4, 4, 16, 
  66.                 12, 8, 8, 8, 16, 12, 12, 16, 4, 8, 12, 12, 16, 0, 4, 8, 12, 16, 
  67.                 4, 4, 8, 12, 16, 8, 8, 12, 12, 4, 4, 8, 8, 12, 4, 0, 4, 8, 12, 
  68.                 8, 4, 4, 8, 12, 8, 8, 12, 12, 8, 4, 4, 4, 8, 8, 4, 0, 4, 8, 12, 
  69.                 8, 4, 8, 12, 12, 8, 12, 12, 12, 8, 8, 4, 4, 12, 8, 4, 0, 4, 12, 
  70.                 8, 4, 4, 8, 12, 8, 8, 12, 16, 12, 12, 8, 4, 16, 12, 8, 4, 0, 
  71.                 16, 12, 8, 4, 4, 12, 8, 8, 12, 8, 12, 16, 16, 16, 4, 8, 12, 12, 
  72.                 16, 0, 4, 8, 12, 16, 4, 8, 12, 8, 8, 8, 12, 12, 12, 4, 4, 8, 
  73.                 8, 12, 4, 0, 4, 8, 12, 4, 4, 8, 8, 8, 8, 8, 8, 8, 8, 4, 4, 4, 
  74.                 8, 8, 4, 0, 4, 8, 8, 4, 8, 8, 12, 12, 12, 8, 8, 12, 8, 8, 4, 
  75.                 4, 12, 8, 4, 0, 4, 8, 4, 4, 8, 16, 16, 16, 12, 8, 16, 12, 12, 
  76.                 8, 4, 16, 12, 8, 4, 0, 12, 8, 4, 8, 12, 12, 16, 16, 16, 8, 8, 
  77.                 12, 12, 12, 4, 4, 8, 8, 12, 0, 4, 8, 4, 12, 12, 12, 12, 12, 8, 
  78.                 8, 8, 8, 8, 8, 4, 4, 4, 8, 4, 0, 4, 4, 16, 16, 16, 12, 12, 12, 
  79.                 12, 12, 8, 8, 12, 8, 8, 4, 4, 8, 4, 0, 4, 16, 16, 16, 16, 16, 
  80.                 12, 12, 12, 12, 12, 8, 8, 8, 8, 8, 4, 4, 4, 0), 
  81.               .Dim = c(19L, 19L), 
  82.               .Dimnames = list(NULL, c(&quot;X1&quot;, &quot;X2&quot;, &quot;X3&quot;, &quot;X4&quot;, &quot;X5&quot;, 
  83.                                        &quot;X6&quot;, &quot;X7&quot;, &quot;X8&quot;, &quot;X9&quot;, &quot;X10&quot;, 
  84.                                        &quot;X11&quot;, &quot;X12&quot;, &quot;X13&quot;, &quot;X14&quot;, &quot;X15&quot;, 
  85.                                        &quot;X16&quot;, &quot;X17&quot;, &quot;X18&quot;, &quot;X19&quot;)))
  87.     set.seed(12345) 
  88.     h_share &lt;- diff(c(0, sort(runif(size)), 1)) 
  89.     e_share &lt;- diff(c(0, sort(runif(size)), 1)) 
  91.     alpha &lt;- numeric(size + 1)
  92.     beta &lt;- numeric(size + 1)
  93.     
  94.     get_beta_prelim &lt;- function() {
  95.       a_matrix &lt;- exp(alpha + t(c) * gamma)
  96.       beta_prelim &lt;- log(e_share) - log(colSums(a_matrix))
  97.       return(beta_prelim)
  98.     }
  99.     
  100.     get_beta &lt;- function() {
  101.       beta &lt;- get_beta_prelim() - get_beta_prelim()[[1]]
  102.       return(beta)
  103.     }
  104.     
  105.     get_alpha_prelim &lt;- function() {
  106.       b_matrix &lt;- t(exp(beta + t(c) * gamma))
  107.       alpha_prelim &lt;- log(h_share) - log(rowSums(b_matrix))
  108.       return(alpha_prelim)
  109.     }
  110.     
  111.     get_alpha &lt;- function() {
  112.       alpha &lt;- get_alpha_prelim() - get_alpha_prelim()[[1]]
  113.       return(alpha)
  114.     }
  116.     get_pi_ij &lt;- function() {
  117.       alphas &lt;-
  118.         matrix(get_alpha(),
  119.                nrow = length(get_alpha()),
  120.                ncol = length(get_alpha()),
  121.                byrow = FALSE)
  122.       betas &lt;-
  123.         matrix(get_beta_prelim,
  124.                nrow = length(get_beta_prelim()),
  125.                ncol = length(get_beta_prelim()),
  126.                byrow = TRUE) 
  127.       pi_ij &lt;- exp(alphas + betas + gamma * c)
  128.       return(pi_ij)
  129.     }
  130.     get_pi_ij()
  132. </details>
  135. # 答案1
  136. **得分**: 0
  138. ```r
  139. get_pi_ij <- function() {
  140.   alphas <-
  141.     matrix(get_alpha(),
  142.            nrow = length(get_alpha()),
  143.            ncol = length(get_alpha()),
  144.            byrow = FALSE)
  145.   betas <-
  146.     matrix(get_beta_prelim(), # your were missing "()"
  147.            nrow = length(get_beta_prelim()),
  148.            ncol = length(get_beta_prelim()),
  149.            byrow = TRUE) 
  150.   pi_ij <- exp(alphas + betas + gamma * c)
  151.   return(pi_ij)
  152. }
英文: 
  1. get_pi_ij &lt;- function() {
  2.   alphas &lt;-
  3.     matrix(get_alpha(),
  4.            nrow = length(get_alpha()),
  5.            ncol = length(get_alpha()),
  6.            byrow = FALSE)
  7.   betas &lt;-
  8.     matrix(get_beta_prelim(), # your were missing &quot;()&quot;
  9.            nrow = length(get_beta_prelim()),
  10.            ncol = length(get_beta_prelim()),
  11.            byrow = TRUE) 
  12.   pi_ij &lt;- exp(alphas + betas + gamma * c)
  13.   return(pi_ij)
  14. }

huangapple
  • 本文由 发表于 2023年6月1日 18:49:04
  • 转载请务必保留本文链接：https://go.coder-hub.com/76381116.html
  • function
  • matrix
  • nested-function
  • r
  • vector
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定