r - タイプ「NULL」の非(リストまたはベクトル)に適用されるis.na()はどういう意味ですか?

Question

NAのないdata.frameからフォワード手順でCoxモデルを選択したいです。サンプルデータは次のとおりです。

test <- data.frame(
  x_1   = runif(100,0,1),
  x_2   = runif(100,0,5),
  x_3   = runif(100,10,20),
  time  = runif(100,50,200),
  event = c(rep(0,70),rep(1,30))
)

この表には何の意味もありませんが、とにかくモデルを構築しようとすると:

modeltest <- coxph(Surv(time, event) ~1, test)
modeltest.forward <- step(
  modeltest, 
  data      = test, 
  direction = "forward", 
  scope     = list(lower = ~ 1, upper = ~ x_1 + x_2 + x_3)
)

転送は最初のステップで終了し、次のように述べています。

in is.na(fit$coefficients) : is.na() はタイプ 'NULL' の非 (リストまたはベクトル) に適用されます

(3回)

上位モデルを変更しようとしましたが、試してみupper = ~ 1ましたが、警告が残ります。わかりません: NA がなく、ベクトルはすべて数値です (確認しました)。人々が同じ問題を抱えているかどうかを検索しましたが、見つけることができたのは、ベクターの名前またはクラスによる問題だけでした。

コードの何が問題になっていますか?