Displaying 2 results from an estimated 2 matches for "pred4".
Did you mean:
preds
2006 May 27
1
Recommended package nlme: bug in predict.lme when an independent variable is a polynomial (PR#8905)
...ly(Orthodont$age, 3)
> mm2 <- model.matrix(~ poly(age, 3, coefs = attr(p, "coefs")) + Sex, data =
Newdata)
>
> data.frame(pred1 = predict(fm, level = 0, newdata = Newdata),
+ pred2 = mm1 %*% fixef(fm),
+ pred3 = head(predict(fm, level = 0)),
+ pred4 = mm2 %*% fixef(fm))
pred1 pred2 pred3 pred4
1 18.61469 18.61469 23.13079 23.13079
2 23.23968 23.23968 24.11227 24.11227
3 29.90620 29.90620 25.59375 25.59375
4 36.19756 36.19756 27.03819 27.03819
5 18.61469 18.61469 23.13079 23.13079
6 23.23968 23.23968 24.11227 24.11227
Best regard...
2006 May 30
0
(PR#8905) Recommended package nlme: bug in predict.lme when an independent variable is a polynomial
...)) + Sex,=
data =3D
> >> Newdata)
> >>>
> >>> data.frame(pred1 =3D predict(fm, level =3D 0, newdata =3D Newdata),
> >> + pred2 =3D mm1 %*% fixef(fm),
> >> + pred3 =3D head(predict(fm, level =3D 0)),
> >> + pred4 =3D mm2 %*% fixef(fm))
> >> pred1 pred2 pred3 pred4
> >> 1 18.61469 18.61469 23.13079 23.13079
> >> 2 23.23968 23.23968 24.11227 24.11227
> >> 3 29.90620 29.90620 25.59375 25.59375
> >> 4 36.19756 36.19756 27.03819 27.03819
> >> 5 1...