以下のコードでボックスコックス変換を実行しようとしている:
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
しかし、このエラーが出続けます:
Error in lm.fit(x,y,offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases Calls:
よろしくお願いします!
エラーメッセージ
lm.fit(x,y,offset = offset, singular.ok = singular.ok, ...) : 0 (非NA)の場合
というエラーメッセージは、lm(y ~ x)
コマンドによって生成されます;
以下はその例です:
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
エラーはデータ中のNA'または不正な変換によって引き起こされる可能性があります。
#From the mtcars dataset
mpg.reg3 <- lm(mpg ~ cylinders + displacement + horsepower + weight + acceleration + year + origin, data=Auto, na.action=na.exclude)
na.action=引数に注目してください。これを
na.excludeに設定すると、
lm関数はデータ中のNA'を無視することができます。もう一つのオプションは
na.omit`で、これは少し違った動作をします。
もう1つの問題は、データの変換が間違っている可能性がある。