从 function 中的 if 语句中抑制 NULL output

Question

我正在尝试使用 replicate() 获取值矩阵，但是当 if() 语句返回 NULL 时，它会创建一个列表而不是矩阵。 我在这里阅读了 web 页和问题，但似乎无法找到有用的东西。 我尝试了 invisible 和 sink 的变体，但仍然无法获得不返回 NULL 值的 output。

这是一个为第 5 个条目提供 NULL 值的示例。 如何让 if() 语句不返回任何内容，包括 NULL？

set.seed(10)
reps <- 10

f_myfun2 <- function(i){
  x1 <- rep(1:4, each=5)
  x2 <- rep(1:5, times=4)
  n <- length(x1)
  y <- 0.20 + 0.30*x1 + 0.7*x2 + 0.50*rnorm(n)
  cis <- confint(lm( y ~ x1*x2 ))
  int_lower <- cis[4,1]
  int_upper <- cis[4,2]
  if(int_lower > 0 | int_upper < 0){
    # I don't want to return anything, including a NULL value
    # tried various things including invisible, sink, etc
  }
  else{
    cis2 <- confint(lm( y ~ x1 + x2))
    c(cis2[2,], cis2[3,])
  }
}


sims <- replicate(reps, f_myfun2(1))
sims  # [[5]] is NULL rather than just missing
str(sims) # now it's a list rather than a matrix without the NULL
is(sims)

Answer 1

如果您的最终目标是将它们全部合并到一个数据帧/矩阵中，则您不必担心这些NULL值在合并时会自动删除。

set.seed(10)
sims <- replicate(reps, f_myfun2(1))
result <- do.call(rbind, sims)
result
#            2.5 %    97.5 %     2.5 %    97.5 %
# [1,]  0.26428935 0.5953288 0.5766971 0.8384067
# [2,]  0.20815417 0.5763875 0.5750148 0.8661288
# [3,]  0.16616864 0.5437533 0.5694641 0.8679710
# [4,]  0.05366132 0.5055326 0.5076891 0.8649247
# [5,]  0.26292580 0.6246576 0.6441824 0.9301565
# [6,]  0.21173249 0.5972766 0.5297051 0.8345045
# [7,] -0.01499442 0.5077043 0.5399975 0.9532272
# [8,]  0.23639931 0.5871463 0.5270991 0.8043890
# [9,] -0.09686737 0.2529322 0.4378804 0.7144212