אומדי Realized Kernel: מדידת תנודתיות עם רעש מיקרו־מבני

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

השונות הממומשת הרגילה משתמשת רק במכפלת ההשהיה אפס: היא מסכמת את ריבועי התשואות. בדגימה צפופה מאוד, הרעש עלול להצטבר מהר יותר משהיא קולטת תנועה חדשה במחיר החבוי. Realized Kernel מוסיף לאיבר ההשהיה אפס אוטו־קווריאנסים משוקללים בהשהיות חיוביות ושליליות. תחת הנחות מתאימות ובהינתן רוחב פס גדל, האיברים בהשהיות שאינן אפס מסייעים לקזז את השפעת הרעש העיקרית תוך שמירת השינוי במחיר היעיל.

Barndorff-Nielsen, Hansen, Lunde ו־Shephard פיתחו Realized Kernels למדידת שינוי בדיעבד בנוכחות חיכוכים בשוק, במאמרם המקורי ב־Econometrica. השיטה היא חישוב משוקלל בקפידה של שונות ארוכת־טווח, ולא הוראה להשתמש בכל tick בלי לבדוק את הנתונים. מדד הבסיס הפשוט יותר מוסבר במדריך לחישוב תנודתיות ממומשת. כאן המוקד הוא תיקון הגרעין והבחירות הכרוכות בו.

נסמן ב־\(X_t\) את לוגריתם המחיר היעיל החבוי, וב־\(Y_{t_i}\) את לוגריתם המחיר המתועד בזמן \(t_i\):

\[ Y_{t_i}=X_{t_i}+\epsilon_{t_i}, \]

כאשר \(\epsilon_{t_i}\) מייצג רעש מיקרו־מבני. מודל בסיס נפוץ מניח שהמחיר החבוי הוא תהליך רציף בעל קצב שונות מקומי \(\sigma_t^2\), והיעד בפרק הזמן הקבוע \([0,T]\) הוא השונות המשולבת:

\[ IV_T=\int_0^T \sigma_t^2\,dt. \]

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

אם התהליך החבוי כולל קפיצות, השונות הריבועית כוללת גם את ריבועי הקפיצות. השאלה אם זהו היעד תלויה במודל ובשאלת המחקר. גזירה למחיר רציף מכוונת לשונות המשולבת; אין לתאר Realized Kernel כמפריד בין שונות רציפה לקפיצות אלא אם מוסיפים מבנה חסין לקפיצות. המדריך לשונות ריבועית מסביר את ההבדל בין שונות משולבת לשינוי הכולל במסלול.

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

עבור סדרה של \(n\) תשואות נצפות \(r_i=Y_{t_i}-Y_{t_{i-1}}\), נגדיר את סכום האוטו־קווריאנסים הלא־מנורמל בהשהיה \(h\):

\[ \Gamma_h=\sum_{i=h+1}^{n} r_i r_{i-h}, \qquad \Gamma_{-h}=\Gamma_h. \]

בהשהיה אפס, \(\Gamma_0=\sum_{i=1}^n r_i^2\), כלומר השונות הממומשת הרגילה. בהשהיה אחת, המכפלות מצמידות תשואות סמוכות. במודל של רעש מחיר חיבורי ובלתי תלוי, רכיב הרעש בתשואות שכנות מקבל סימנים מנוגדים: שגיאת תצפית חיובית מעלה תשואה אחת ומורידה את הבאה. כך עשוי להיווצר אוטו־קווריאנס שלילי בהשהיה הראשונה. דפוסים בהשהיות נוספות עשויים גם הם להראות כיצד הרעש משפיע על סדרת התשואות.

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

האוטו־קווריאנסים אינם מסמנים tick יחיד כ"רעש" או כ"אות"; הם צוברים מכפלות צולבות לאורך הפרק. \(\Gamma_1\) שלילי יכול להתאים לתנודות בין מחיר הקנייה והמכירה, אך אינו מוכיח לבדו מנגנון מסוים. אם לתשואות האמיתיות יש תלות סדרתית או שלרעש מבנה אחר, אותו אוטו־קווריאנס עשוי לשקף כמה השפעות. יש לפרש את האומד לפי המודל המוצהר ומשקלי הגרעין, ולא כאבחון ישיר של השוק מכל השהיה.

עבור השהיה מרבית \(H\), מוסכמה נפוצה במדגם סופי היא

\[ \widehat{IV}_{RK}(H)=\Gamma_0+2\sum_{h=1}^{H} k\!\left(\frac{h}{H}\right)\Gamma_h. \]

הגורם 2 כולל את איברי ההשהיות השליליות הסימטריים. \(k(x)\) הוא משקל גרעין בתחום \(0\le x\le1\), עם \(k(0)=1\), והמשקלים דועכים לאפס סמוך לקצה רוחב הפס. מאמרים ותוכנות עשויים להשתמש במוסכמות מעט שונות לקצוות ולאינדוקס, למשל \(H+1\) במכנה. יש לציין את המוסכמה ולהשתמש בה בעקביות כדי לשחזר תוצאה.

גרעין Parzen הוא בחירה נפוצה שאינה שלילית:

\[ k(x)= \begin{cases} 1-6x^2+6x^3, & 0\le x\le \tfrac12,\\ 2(1-x)^3, & \tfrac12 < x\le1,\\ 0, & x>1. \end{cases} \]

הוא נותן משקל גדול להשהיות קצרות ומפחית בהדרגה את משקלן של הרחוקות. ההדרגתיות מועילה משום שרעש בתדירות גבוהה עשוי ליצור תלות בין תשואות, בעוד שאומדני השהיות רחוקות נעשים תנודתיים יותר. בניית Parzen הלא־שלילית נועדה להשאיר את אומד Realized Kernel החד־משתני לא־שלילי תחת המימוש המוגדר. לא לכל הגרעינים יש תכונה זו. גרעין flat-top עשוי להציע מאפייני הטיה נוחים, אך להפיק אומדן שלילי במדגם סופי.

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

<!-- learn:illustration --> <!-- Text-free concept: הדים סימטריים של תשואות דועכים בהדרגה כשהם מתרחקים ממרכז בהיר, וממחישים משקלי השהיה פוחתים בתוך רוחב פס סופי. המחשה רעיונית, לא נתוני שוק או פלט של אומד. -->

רוחב הפס \(H\) הוא ההשהיה הגדולה ביותר שנכללת בסכום. אם הוא קטן מדי, תלות סדרתית חשובה של רעש התצפית עלולה להישאר ללא טיפול. אם הוא גדול מדי, מצטרפות מכפלות השהיה תנודתיות יותר ועשויות להגדיל את שונות הדגימה. לכן רוחב הפס הוא בחירת כיול סטטיסטית, לא מספר התצפיות למחיקה ולא פרק זמן שמתאים לכל נכס.

בתיאוריה האסימפטוטית של Realized Kernel עם Parzen, רוחב הפס האופטימלי גדל בסדר גודל \(n^{3/5}\). ביטוי מעשי שמופיע בספרות הוא

\[ H^*=c^*\,\xi^{4/5}n^{3/5}, \qquad c^*_{\text{Parzen}}=3.5134, \]

כאשר \(n\) הוא מספר התשואות ברשת הדקה ו־\(\xi\) מסכם את רמת הרעש היחסית ואת הקוורטיות המשולבת של המחיר החבוי. ביטוי אחד לסקאלה הוא

\[ \xi^2=\frac{\omega^2}{\sqrt{T\int_0^T\sigma_u^4\,du}}, \]

כאשר \(\omega^2\) היא שונות רעש התצפית. בפועל הכמויות אינן ידועות ויש לאמוד אותן, לעיתים באמצעות אומדי־עזר לשונות הרעש ולקוורטיות המשולבת. הנוסחה מראה שמספר ההשהיות האופטימלי תלוי גם בגודל המדגם וגם ביחס רעש־לאות; היא אינה הגדרה קבועה שמתאימה לכולם.

התכנון האסימפטוטי דורש שרוחב הפס יגדל אך יישאר קטן ביחס לגודל המדגם, לרוב בניסוח \(H\to\infty\) ו־\(H/n\to0\). כלל \(n^{3/5}\) מאזן את ההטיה המובילה והשונות תחת המודל ותנאי הגרעין במאמר. במקטע נתונים סופי, אומדן plug-in עלול להיות לא יציב אם אומד־העזר לרעש או לקוורטיות גרוע. דווחו את \(H\), מוסכמת האינדוקס, ערכי העזר וניתוח הרגישות לחלופות סבירות.

נבחן סדרת תשואות היפותטית קטנה ומכוונת, ביחידות נקודות בסיס:

\[ (r_1,r_2,r_3,r_4)=(1,-1,1,-1)\ \mathrm{bp}. \]

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

\[ \Gamma_0=1^2+(-1)^2+1^2+(-1)^2=4\ \mathrm{bp}^2. \]

המכפלות בהשהיה הראשונה הן \((-1),(-1),(-1)\), ולכן \(\Gamma_1=-3\ \mathrm{bp}^2\). בהשהיה השנייה \(r_3r_1=1\) ו־\(r_4r_2=1\), ולכן \(\Gamma_2=2\ \mathrm{bp}^2\).

נקבע \(H=2\) רק כדי להמחיש את הנוסחה. בגרעין Parzen, \(k(1/2)=1-6(1/4)+6(1/8)=1/4\), ואילו \(k(1)=0\). לכן

\[ \widehat{IV}_{RK}=4+2\left(\tfrac14\right)(-3)+2(0)(2) =2.5\ \mathrm{bp}^2. \]

השונות הממומשת הרגילה היא \(4\ \mathrm{bp}^2\); האיבר השלילי המשוקלל בהשהיה הראשונה מקטין כאן את אומדן הגרעין. סכום ההשהיה השנייה אינו תורם, כי משקל Parzen שלו אפס בגבול. התוצאה אינה שלילית, אך החשבון לבדו אינו מוכיח שהיא שווה לשינוי החבוי במחיר היעיל. לא נצפה המסלול החבוי ולא חושב רווח אי־ודאות, ולכן אין כאן תחזית או מסקנת מסחר. השורש הריבועי, כ־1.58 נקודות בסיס, הוא רק סקאלת סטיית התקן של הפרק ההיפותטי הזה.

הדוגמה מראה גם מדוע צריך לציין את רוחב הפס ואת כלל האינדוקס. שינוי ב־\(H\), במוסכמת הקצוות או בתשואות משנה את המשקלים או המכפלות הנכללות. אפשר לשחזר את התוצאה כי כל תשואה, סכום השהיה, משקל ויחידה מוצגים.

סכומי השהיות לא־מנורמלים עלולים לתת משקל לא־מידתי לתצפיות בקצוות, במיוחד אם למחיר הרשום הראשון או האחרון יש שגיאת רעש גדולה. מאמרים מעשיים על Realized Kernels דנים במיצוע מקומי, המכונה לעיתים jittering, להפחתת השפעות קצה. במימוש נפוץ מחליפים מחיר קצה בממוצע מקומי של תצפיות סמוכות לפני חישוב התשואות. יש לתעד את אורך החלון והאם ממוצעים קצה אחד או את שניהם.

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

המחקר המעשי משנת 2009 על Realized Kernels דן בנתוני עסקאות וציטוטים, השפעות קצה, מגמות מקומיות ובגרעין Parzen הלא־שלילי. הלקח המעשי הוא ש"עמיד לרעש" אינו "אדיש לנתונים". רכיב מחיר המשתנה בהדרגה לאורך חלון קצר, שינוי בחיכוכי השוק או אופן בניית הציטוטים עשויים לערער את פרשנות הבסיס. Jittering עשוי להפחית בעיית קצה אחת, אך אינו מתקן סדרה שאינה מיושרת או תצפית שגויה באמצע המדגם.

התיאוריה הקלאסית מניחה הנחות על תהליך המחיר החבוי, הדגימה, הרעש, סדירות הגרעין ורוחב הפס. מבנים מסוימים של Realized Kernel עמידים לצורות רחבות יותר של רעש בעל תלות סדרתית ולזמני תצפית אנדוגניים מהדוגמה הפשוטה של רעש בלתי תלוי. עם זאת, העמידות תלויה במשפט המתאים: צריך לבחור גרעין ורוחב פס העומדים בתנאים הרלוונטיים. המחקר של Aït-Sahalia, Mykland ו־Zhang על רעש מיקרו־מבני תלוי מסביר מדוע יש להתייחס לתלות הרעש במפורש, ולא להניח שהנוסחה הבסיסית מכסה אוטומטית כל הרחבה.

בחירת הגרעין חשובה. משתמשים לעיתים ב־Parzen לאומד חד־משתני לא־שלילי. משקלי flat-top עשויים לכוון לתכונות הטיה אחרות ולהניב אומד שלילי במדגם סופי. אל תחליפו לגרעין HAC שרירותי רק משום שהוא מוכר: ספרות Realized Kernel מסבירה שכמה אפשרויות מוכרות, ובהן Bartlett במבנה שנבחן, אינן מספקות את אותה תוצאת עקיבות. אם האומד שלילי, ציינו את הגרעין ואת כלל המדגם הסופי; אל תחתכו אותו לאפס בשקט ואל תכנו אותו שונות פיזיקלית שלילית.

האומדן יורש גם את מגבלות ה־proxy למחיר ואת הגדרת היעד. קפיצות עשויות להיות חלק מהשונות הריבועית; חריגים יכולים לשלוט במכפלות; תצפיות לא־סדירות או אסינכרוניות יכולות לשנות את הגודל הנאמד; והתנודתיות או הרעש עשויים להשתנות במהלך הפרק. Realized Kernel חד־משתני אינו פותר בעצמו דגימה רב־משתנית שאינה מסונכרנת, אינו מסיר קפיצות ואינו מפריד תנודות bid–ask מכל שאר מקורות התלות. אומד עמיד הוא עמיד במסגרת מודל מוצהרת, לא בפני כל פגם נתונים אפשרי.

דיווח שניתן לשחזור צריך לציין את סדרת המחירים ואת תכנון הדגימה, חלון התצפית, יחידות התשואה, פונקציית הגרעין, ההשהיה המרבית \(H\), מוסכמת הקצוות, כלל המיצוע בקצוות וכל ניקוי. אם נעשה שימוש ברוחב פס plug-in, דווחו את אומדני הרעש והקוורטיות או פרטים מספיקים לשחזורם. ציינו אם התוצאה היא אומדן לשונות משולבת, לשונות ריבועית הכוללת קפיצות או למדד תנודתיות שעבר טרנספורמציה.

האומדן מתאר שינוי בפרק הזמן שנדגם. הוא אינו אומדן לתנודתיות עתידית אלא אם מודל חיזוי נפרד ממפה מדידות היסטוריות לאופק עתידי ונבדק מחוץ למדגם. הוא גם אינו אומדן למרווח מצוטט או לעלות ביצוע: מדריך מרווח bid–ask של Roll משתמש באוטו־קווריאנס תשואה כדי לאמוד גודל אחר. המדריך לתנודתיות ממומשת דו־שלבית מציג גישה אחרת העמידה לרעש, עם מבנה תיקון שונה.

השאלה המרכזית אינה רק אם להשתמש ב־Realized Kernel. השאלה היא אם נתוני המחיר, תכנית הדגימה, היעד, הגרעין ורוחב הפס מתאימים מספיק לשאלת המחקר כדי לתמוך בפרשנות המדווחת. טבלת רגישות ומגבלות מודל מפורשות מקלות על בדיקת השיקול. גם אומדן זהיר נשאר מדידה היסטורית עם אי־ודאות דגימה; הוא אינו מוכיח כלל מסחר רווחי ואינו מבטיח שאפקטים מיקרו־מבניים נעלמו.