他の 2 つの線形モデルの残差を使用して線形モデルを構築する問題に取り組んでいます。実際のデータセットを使用するよりも問題を簡単にここに置くことができるため、UN3 データセットを使用して問題を示しました。
ここに私のRコードがあります:
head(UN3)
m1.lgFert.purban <- lm(log(Fertility) ~ Purban, data=UN3)
m2.lgPPgdp.purban <- lm(log(PPgdp) ~ Purban, data=UN3)
m3 <- lm(residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban))
これが私が得ているエラーです:
> m3 <- lm(residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban))
Error in model.frame.default(formula = residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban), :
variable lengths differ (found for 'residuals(m2.lgPPgdp.purban)')
このエラーが実際に発生する理由がよくわかりません。ログ関連の問題であれば、最初の 2 つのモデルを構築しているときにエラーが発生するはずです。