1. ホーム
  2. r

[解決済み] lm.fit(x,y,offset = offset, singular.ok,...) 0 非NAケースでboxcox式で計算するとエラーになる。

2022-03-12 21:38:39

質問

次のコードでボックスコックス変換を実行しようとしています。

urban1 <- subset(ski,urban <= 4,na.rm=TRUE)
ski$gender <- as.numeric((as.character(ski$gender)),na.rm=TRUE)
urban1 <- as.numeric((as.character(urban1)))
x <- (ski$gender*urban1)
y <- ski$EPSI.
bc <- boxcox(y ~ x) 
(trans <- bc$x[which.max(bc$y)]) 
model3 <- lm(y ~ x) 
model3new <- lm(y^trans ~ x)
ski$EPSI. <- ski$EPSI. + 1

しかし、このエラーが出続けています。

lm.fit(x,y,offset = offset, singular.ok = singular.ok, ...) でエラー : 0 (non-NA) cases Calls: ... eval -> eval -> boxcar -> boxcar.formula -> lm -> lm.fit 実行停止中。

ありがとうございました。

どのように解決するのですか?

エラーメッセージ

<ブロッククオート

lm.fit(x,y,offset = offset, singular.ok = singular.ok, ...) : 0 (非NA)の場合

が生成されます。 lm(y ~ x) コマンドを使用する場合、変数 x または y (またはその両方)はNAのみです。
以下はその例である。

n <- 10
x <- rnorm(n,1)
y <- rep(NA,n)
lm(y ~ x)

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
  0 (non-NA) cases

あなたのコードで、私はテストすることをお勧めします(ちょうどあなたの lm コマンドを使用して、変数の1つにすべてのNAが使用されているかどうかを確認します。

all(is.na(x))
all(is.na(y))
all(is.na(y^trans))

私の例では

all(is.na(y))
[1] TRUE