r - R 要替换的项目数不是替换长度的倍数

Question

我写了一个函数，它必须创建一个 rxk 矩阵并用代码中给出的公式填充它。

mCreate<-function(r,k)
{
 result<-array(0,dim=c(r,k))
  for(R in 0:(r-1)) 
  {
    for(K in 0:(k-1))
     {
      result[R][K] = 5*R^K-2*K^R
      }
   }
   return(result)
}

当我尝试对此进行测试时，我总是得到“要替换的 R 个项目不是替换长度的倍数”。我不明白为什么，我在谷歌上搜索过，我总是得到信息，因为我要去出境。

Answer 1

我会用outer这个代替......

mCreate <- function( r , c ){
  r <- seq_len(r)
  c <- seq_len(c)
  m <- outer( r , c , FUN = function(x,y) 5*x^y-2*y^x )
  return( m )
}

mCreate( 4 , 4 )
#     [,1] [,2] [,3] [,4]
#[1,]    3    1   -1   -3
#[2,]    8   12   22   48
#[3,]   13   29   81  277
#[4,]   18   48  158  768

#  Sanity check:
mCreate( 1 , 1 )
#     [,1]
#[1,]    3

#  Safety check:
mCreate( 0 ,0 )
<0 x 0 matrix>