שיטת הדלתה

שיטת הדלתה בסטטיסטיקה היא תוצאה המאפשרת את אמידת השונות של פונקציה של אמד לפרמטר, כאשר התפלגותו האסימפטוטית של האמד היא נורמלית וסטיית התקן של האמד ידועה או ניתנת לאמידה.

היסטוריה

סקירה של ההיסטוריה של שיטת הדלתה ניתנה על ידי והר הוף ומכתב תשובה מאת פורטנוי. הרעיון שליו מבוססת שיטת הדלתה היה ידוע כבר במאה התשע עשרה. ב-1838 הוצג הרעיון הבסיס על ידי בסל, וב-1838 הוא הוצג בספרו של איירי. השימוש הסטטיסטי הראשון בשיטה נעשה ככל הנראה ב-1928 על ידי קלי. הניסוח הפורמלי הוצג ב-1935 על ידי דוב.

שיטת הדלתה

משפט

תהי $X_{n}$ סדרה של משתנים מקריים בעלת התפלגות אסימפטוטית נורמלית, ובאופן פורמלי ${{\sqrt {n}}[X_{n}-\theta ]\,{\xrightarrow {D}}\,{\mathcal {N}}(0,\sigma ^{2})},$, כאשר $\theta$ ו-$\sigma ^{2}$ קבועים ממשיים, ו-${\xrightarrow {D}}$ מציין התכנסות בהתפלגות.

כן תהא $g$ פונקציה הגזירה בנקודה $\theta$ כך ש-$g'(\theta )$ רציפה ו-$g'(\theta )\neq 0$.

אזי: ${{\sqrt {n}}[g(X_{n})-g(\theta )]\,{\xrightarrow {D}}\,{\mathcal {N}}(0,\sigma ^{2}\cdot [g'(\theta )]^{2})}$.

הוכחה:

על פי משפט הערך הממוצע של לגראנז' (משפט ערך הביניים) קיים ${\tilde {\theta }}$ כך ש-$g(X_{n})=g(\theta )+g'({\tilde {\theta }})(X_{n}-\theta )$, כך ש-${\tilde {\theta }}$ נמצא בין $X_{n}$ ובין $\theta$. (זהו למעשה פיתוח טיילור של $g$ סביב $\theta$.)

נסדר מחדש את האיברים ונכפיל ב-${\sqrt {n}}$, ונקבל כי ${\sqrt {n}}[g(X_{n})-g(\theta )]=g'\left({\tilde {\theta }}\right){\sqrt {n}}[X_{n}-\theta ]$.

מכיוון ש-$X_{n}\,{\xrightarrow {P}}\,\theta$ (כאשר ${\xrightarrow {P}}$ מציין התכנסות בהסתברות) מקבלים כי ${\tilde {\theta }}{\xrightarrow {P}}\,\theta$, ומכיוון ש-$g'(\theta )$ רציפה אנו מקבלים כי גם $g'({\tilde {\theta }}){\xrightarrow {p}}g'(\theta )$.

כזכור, על פי תנאי המשפט ${{\sqrt {n}}[X_{n}-\theta ]\,{\xrightarrow {D}}\,{\mathcal {N}}(0,\sigma ^{2})},$, ולכן על פי משפט סלוצקי מקבלים באופן מיידי כי ${{\sqrt {n}}[g(X_{n})-g(\theta )]\,{\xrightarrow {D}}\,{\mathcal {N}}(0,\sigma ^{2}\cdot [g'(\theta )]^{2})}$.

הערה: ניתן להוכיח כי השגיאה בקירוב שואפת בהסתברות לאפס. כן קיימת גרסה רב ממדית.

דוגמאות

ריבוע התוחלת

על פי משפט הגבול המרכזי, הממוצע של $n$ משתנים מקריים בלתי תלויים ושווי התפלגות בעלי תוחלת $\mu$ ושונות סופית וחיובית $\sigma ^{2}$ מתפלג אסימפטוטית נורמלית, כלומר ${\sqrt {n}}[{\bar {X}}_{n}-\mu ]{\xrightarrow {D}}N(0,\sigma ^{2})$. תהי $g(x)=x^{2}$, ולכן $g'(x)=2x$. מכאן נקבל כי ל-${\bar {X}}_{n}^{2}$ יש התפלגות אסימפטוטית נורמלית עם תוחלת $\mu ^{2}$ ושונות ${\frac {4\mu ^{2}\sigma ^{2}}{n}}$.

רגרסיה לוגיסטית

נתבונן במודל which(minu==0)$\log {\frac {\pi (x)}{1-\pi (x)}}=\alpha +\beta \ x$, ונסמן ב-${\hat {\beta }}_{n}$ את אמדן הנראות המקסימלית ל-$\beta$. מכיוון שזהו אמד נראות מקסימלית ידוע כי התפלגותו היא אסימפטוטית נורמלית עם תוחלת $\beta$ ושגיאת תקן $s>0$, כאשר $s$ נאמדת על ידי שימוש באינפורמציה של פישר. חוקרים מתעניינים בדרך כלל בערך $OR=e^{\beta }$ שמפורש כיחס הסיכויים של $Y$ בהינתן $X$. מכיוון ש-$g(x)=e^{x}$ היא פונקציה רציפה שנגזרתה רציפה, נקבל כי לאמדן יחס הסיכויים $e^{\hat {\beta }}$ יש התפלגות אסימפטוטית נורמלית עם תוחלת $e^{\beta }$, ולאמוד את שגיאת תקן שלו על ידי $se^{\hat {\beta }}$.

מכאן נוכל לקבל כי רווח סמך ברמת סמך $100(1-\alpha )\%$ ליחס הסיכויים הוא $e^{\hat {\beta }}\ \pm \ Z_{\frac {\alpha }{2}}se^{\hat {\beta }}$. כן נוכל לבדוק את ההשערה כי יחס הסיכויים שווה ל-1 על ידי הסטטיסטי $Z={\frac {e^{\hat {\beta }}-1}{se^{\hat {\beta }}}}$.

ראו גם

• רגרסיה לוגיסטית
• משפט הגבול המרכזי
• התכנסות בהתפלגות
• התכנסות בהסתברות
• משפט סלוצקי

לקריאה נוספת

• Lehmann, E. L.; Casella, G. (2006). Theory of point estimation, 2nd edition. Springer Science & Business Media. ISBN 0-387-98502-6.

קישורים חיצוניים

• Standard error of a nonlinear function of estimated regression coefficients - R alr package

• How can I estimate the standard error of transformed regression parameters in R using the delta method? - UCLA: Statistical Consulting Group.