النماذج الخطية ومصفوفات التصميم

d <- data.frame(x = 0:3, y = c(1, 3, 5, 7))
fit <- lm(y ~ x, data = d, na.action = na.fail)
cat(sprintf("intercept=%.1f; slope=%.1f\n", coef(fit)[1], coef(fit)[2]))
predicted <- predict(fit, newdata = data.frame(x = c(4, 5)))
cat("predictions=", paste(round(predicted, 6), collapse = ","), "\n", sep = "")

intercept=1.0; slope=2.0
predictions=9,11

تقع النقاط الاصطناعية على y=1+2x تماماً؛ وهذا يفحص الواجهة ولا يثبت افتراضات الاستدلال أو الجودة الخارجية.

d <- data.frame(x = 0:3, y = c(1, 3, 5, 7))
X <- model.matrix(~ x + I(2*x), data = d)
cat("columns=", ncol(X), "; rank=", qr(X)$rank, "\n", sep = "")
cat("duplicated slope information=", isTRUE(all.equal(X[, 3], 2*X[, 2])), "\n", sep = "")

columns=3; rank=2
duplicated slope information=TRUE

لا تحدد الأعمدة المتناسبة مع الثابت ثلاثة معاملات مستقلة؛ ويكفي فحص الرتبة دون عرض ملاءمة غير مستقرة.