برمجة R / التقييم والملكية والنتائج القابلة للإعادة / الدرس 13

الصيغ والتباينات والنماذج المعممة

افحص الترميز ومقياس الاستجابة قبل تفسير معامل أو تنبؤ.

وقت القراءة التقريبي: 4 دقيقةتاريخ مراجعة المحتوى 2026-10-28

في هذه الصفحة

أهداف الدرس

توسيع حدود الصيغة وتحديد أعمدة الترميز.
تمييز عائلة GLM عن دالة الربط.
فحص أعداد الاستجابة الثنائية ومقياس التنبؤ.

الصيغة لغة نمذجة

تضم a*b في الصيغة a وb وتفاعلهما a:b، ولا تحسب الضرب العنصري ببساطة. تطلب I(a*b) الحساب المعتاد داخل حد. يغير حذف الثابت ترميز الأعمدة وتفسيرها. ويختلف تحويل الاستجابة عن اختيار ربط أو إضافة إزاحة ثابتة. ابن مصفوفة نموذج صغيرة وافحصها قبل تفسير تركيب معقد؛ فهي تربط الرموز بالمتنبئ العددي وتبين دخول التأثيرات والتفاعلات والمراجع المقصودة.

تحدد التباينات المقارنات الفئوية

يعلن المثال مستويين وترميز معالجة صريحاً. يحتوي صف التحكم ثابتاً واحداً ومؤشراً صفراً، ويحتوي صف المعالجة واحدين. يشير ثابت النموذج الخطي بهذا الترميز إلى متوسط التحكم الشرطي ويمثل المؤشر فرقاً على مقياس المتنبئ الخطي. قد تمثل ترميزات أخرى مقارنات مختلفة دون تغيير فضاء الدوال في نموذج مكافئ كامل الرتبة. يغير الترميز معنى المعامل؛ فاحفظه عند مشاركة الملاءمة وبناء مصفوفة بيانات جديدة.

r برنامج توضيحي بلغة R الأساسية؛ لم تتوفر بيئة R أثناء المراجعة

d <- data.frame(group = factor(c("control", "treated", "control"),
                               levels = c("control", "treated")))
X <- model.matrix(~ group, data = d,
                  contrasts.arg = list(group = contr.treatment(c("control", "treated"), base = 1)))
cat("columns=", ncol(X), "\n", sep = "")
cat("control row=", paste(X[1, ], collapse = ","), "\n", sep = "")
cat("treated row=", paste(X[2, ], collapse = ","), "\n", sep = "")

الناتج المتوقع

columns=2
control row=1,0
treated row=1,1

يعطي الترميز الصريح ثابتاً ومؤشر معالجة، ويعلن المرجع بدلاً من وراثته من الخيارات العامة.

يفصل النموذج المعمم العائلة والربط

يجمع GLM متنبئاً خطياً وعائلة توزيع وربطاً يربط المتوسط الشرطي بالمتنبئ. يعبر الثنائي ذو الربط اللوغاريتمي للأرجحية عنها خطياً؛ ويعيد predict من نوع response احتمالاً بينما يبقى المقياس الآخر لوغاريتم أرجحية. يستخدم Poisson غالباً ربطاً لوغاريتمياً للعد، لكن التعرض والتشتت والاعتماد أسئلة مستقلة. لا يثبت اسم العائلة الافتراضات. حدد الوحدة والمتوسط والتباين والتعرض قبل الاختيار من مجال الخرج وحده.

تحفظ الاستجابة المجمعة مقامها

يقدم المثال نجاحاً وفشلاً في عمودين. يجمع نموذج الثابت ثلاثة نجاحات في ثماني تجارب، فيعطي 0.375 ولوغاريتم أرجحية log(3/5). يهمل متوسط نسبتين متساوي الوزن المقامات عند اختلاف عدد التجارب. افحص عدم سلبية الأعداد ومجموع التجارب وترتيب العمودين والفقد. قد يصعب الفصل أو ندرة الفئات ملاءمة لوجستية ويمنع التعرف؛ فتحتاج التحذيرات والمعاملات المتطرفة تشخيصاً. المثال واجهة صغيرة ولا يثبت ملاءمة نموذج واقعي معقد.

r برنامج توضيحي بلغة R الأساسية؛ لم تتوفر بيئة R أثناء المراجعة

d <- data.frame(success = c(2, 1), failure = c(2, 3))
fit <- glm(cbind(success, failure) ~ 1, data = d,
           family = binomial(link = "logit"), na.action = na.fail)
cat(sprintf("logit=%.6f\n", unname(coef(fit)[1])))
p <- predict(fit, type = "response")
cat("probabilities=", paste(sprintf("%.3f", p), collapse = ","), "\n", sep = "")

الناتج المتوقع

logit=-0.510826
probabilities=0.375,0.375

يجمع نموذج الثابت ثلاثة نجاحات من ثماني تجارب؛ ويختلف احتمال الاستجابة عن مقياس المتنبئ الخطي.

يحفظ التنبؤ العقد الكامل

احفظ المستويات والترميز والإزاحات والتحويلات والوحدات مع النموذج. الإزاحة مساهمة معلومة بمعامل ثابت وليست متنبئاً عاديّاً يقدر أثره. احفظ معلمات التحويل المتعلمة من التدريب ولا تعِد تعلمها من الجديد. افحص التحذيرات والتقارب والرتبة والبواقي قبل تفسير المعامل. اعرض عدم اليقين والحدود وفق التصميم. النتائج مشتقة تحليلياً؛ ولم ينفذ محلل صيغة أو باني ترميز أو محسن GLM أو تنبؤ R أثناء المراجعة.

الأخطاء الشائعة وإصلاحها

فُسر معامل لوجستي كفرق احتمال

لماذا يحدث؟: المعامل على مقياس لوغاريتم الأرجحية.
كيف تصلحه؟: حول التنبؤ وفق الربط وحدد المقارنة.
كيف تتجنب تكراره؟: سم مقياس كل خرج.

تغير النموذج بعد ترتيب مرجع العامل

لماذا يحدث؟: تغير معنى أعمدة الترميز.
كيف تصلحه؟: احفظ المستويات والتباينات أو أعد تعريفها صراحة.
كيف تتجنب تكراره؟: افحص مصفوفة صغيرة قبل مقارنة المعاملات.

تساوى وزن نسب ذات أعداد تجارب مختلفة عرضاً

لماذا يحدث؟: أهمل مقام الاستجابة الثنائية.
كيف تصلحه؟: قدم النجاح والفشل أو تمثيلاً موزوناً مكافئاً موثقاً.
كيف تتجنب تكراره؟: افحص أعداد النتيجة وإجمالي التجارب.

اختبر فهمك

1. وسع y ~ ab وقارنه بـ y ~ I(ab).

اعرض الحل المشروح

يضم الأول a وb والتفاعل مع الثابت الافتراضي. يضم الثاني حاصل الضرب الحسابي كحد واحد مع ثابت. تختلف الأعمدة عموماً للمتغيرات العددية؛ فافحص model.matrix ولا تعتمد على التشابه البصري.

2. نجاح واحد من تجربتين وتسعة من عشر؛ ما الاحتمال المجمع لنموذج ثابت؟