model.frame.default（object、data、xlev = xlev）のエラー：可変長が異なる（ 'z'のために見つかりました）

Question

この質問はhereとhereの前に尋ねられました。データセットに欠落データがありません。model.frameからmodel.matrixを作成しようとしています。再現可能な例を次に示します。

dat = data.frame(x = rep(1:3, each = 5), 
    y = rnorm(15), 
    z = rep(c(1,1.1,1.3), each = 5)) 
# this works 
mt = model.matrix(model.frame(lm(y ~ x,dat))) 
# this produces the error 
mt = model.matrix(model.frame(lm(y ~ x + z,dat))) 

なぜこの問題を解決するのですか？あるいは、インターセプト列を使用して右辺値の行列を別の、効果的かつ堅牢な方法で取得するにはどうすればよいですか？

Answer 1

問題は、model.matrix()は、最初の引数として数式またはタームオブジェクトが必要であるということです。ソリューション：

model.matrix(y ~ x + z, model.frame(lm(y ~ x + z, dat)))