※上記の広告は60日以上更新のないWIKIに表示されています。更新することで広告が下部へ移動します。

変数変換の妥当性


### 元データ y を z = f(y) と変換する ###
### y, z それぞれにモデルを当てはめて AIC を比較 ###
fm.y <- lm(y ~ x)
AIC.y <- AIC(fm.y)
#
fm.z <- lm(z ~ x)
### 対数変換 z = f(y) = log(y) の場合の例 ###
### df/dy = 1/y を計算し、-2∑log|df/df| を算出する ###
AIC.z <- AIC(fm.z) - 2 * sum(log(abs(1/y)))
### AIC.y, AIC.z の小さいほうを採用する ###