הידען - Hayadan

למ"ס: התמותה העודפת בגלל הקורונה מתחילה לבלוט בנתונים

אבי בליזובסקי ·

Summary: סך התמותה מסוף מרץ עד סוף אוקטובר גבוה ב-10.4% מהתמותה הצפויה ב-2020. כלומר, זהו עודף התמותה בתקופת הקורונה עד סוף חודש אוקטובר

סך התמותה מסוף מרץ עד סוף אוקטובר גבוה ב-10.4% מהתמותה הצפויה ב-2020. כלומר, זהו עודף התמותה בתקופת הקורונה עד סוף חודש אוקטובר

04 באפריל 2020. משאיות של הצבא האיטלקי מעבירות את גופות ההרוגים בברגמו (איטליה) עקב זיהום בנגיף קורונה העטרה למקדש השריפה בפרארה, איטליה. צילום: Shutterstock  filippo.rubin
04 באפריל 2020. משאיות של הצבא האיטלקי מעבירות את גופות ההרוגים בברגמו (איטליה) עקב זיהום בנגיף קורונה העטרה למקדש השריפה בפרארה, איטליה. צילום: Shutterstock filippo.rubin

הלשכה המרכזית לסטטיסטיקה הוציאה דו"ח על מספר המתים וסיבות המוות בישראל בתקופת הקורונה. מסתבר שבאוקטובר אחוז הנפטרים מקורונה ל-21% מכלל הנפטרים.

תקציר

לכל הגרפים והטבלאות באתר הלמ"ס

מתחילת 2020 עד סוף חודש אוקטובר נפטרו בישראל 40,315 תושבים.[1] בשנת 2019 נפטרו באותם חודשים 38,056 תושבים, הבדל של 2,259 נפטרים. בשלושת החודשים הראשונים של 2020, סך התמותה היה נמוך לעומת שלוש השנים הקודמות (2019-2017), כנראה בשל עונת שפעת מוקדמת שהשפיעה על התמותה בעיקר בסוף 2019 ורק מעט בתחילת 2020. בחודשים אפריל-מאי התמותה הייתה מעט גבוהה לעומת שנים קודמות. אך מחודש יולי, ובמיוחד מאוגוסט התמותה ב-2020 גבוהה יותר באופן משמעותי. בחודש ספטמבר, שיעור התמותה היה גבוה ב-19% משיעור התמותה בחודש זה ב-2019 (וכן מחודש ספטמבר בממוצע שלוש השנים האחרונות), ובחודש אוקטובר השיעור גבוה ב-24% לעומת 2019 (וב-21% לעומת ממוצע 2019-2017).

הגידול בתמותה הוא בהיקף קרוב מאוד להיקף התמותה הידועה מנגיף קורונה. מקרי התחלואה הראשונים בנגיף קורונה התגלו בישראל בסוף חודש פברואר 2020. מקרה התמותה הראשון היה ב-20 במרץ. מחודש מרץ עד סוף אוקטובר נפטרו בישראל לפי נתוני משרד הבריאות 2,537 תושבים מנגיף קורונה, 6.3% מסך הנפטרים בחודשים אלו. בחודשים הראשונים מאז כניסת  נגיף קורונה לישראל, מספר הנפטרים מהנגיף היה נמוך יחסית וכך גם האחוז מסך הנפטרים. ואולם מאמצע יולי אחוז זה נמצא במגמת עלייה, והוא הגיע ל-16%מסך הנפטרים בספטמבר ול-21% באוקטובר. עם זאת, בשבועות האחרונים של אוקטובר חלה ירידה באחוז הנפטרים מקורונה מסך הנפטרים - מ-24% בשבוע הראשון של אוקטובר ל-15% בשבוע האחרון של החודש.


<img src="https://www.hayadan.org.il/images/content3/2020/12/Depositphotos_375580302_l-2015-500x155.jpg" alt="קשישים בתקופת הקורונה. <a href="https://depositphotos.com/">המחשה: depositphotos.com
קשישים בתקופת הקורונה. המחשה: depositphotos.com

פטירות בקרב בני 70 ומעלה

פטירות לפי מין

פטירות לפי קבוצת אוכלוסייה

עודף תמותה

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

עודף תמותה ניתן לחישוב בצורות שונות, להלן יוצגו שתי שיטות:

עודף תמותה לפי שיטת P-score

כאשר בוחנים את ההבדל במספר הפטירות ב-2020 לעומת שנים קודמות בשיטת P-score, ניתן לראות את העודף בשבועות 16-13 (19/4/20-13/3/20) – עודף ממוצע של 11.5% ואת העלייה בעודף החל משבוע 31 (סוף יולי). החל משבוע 37 (שהתחיל ב-7/9/20) למעט השבוע האחרון של אוקטובר, העודף בכל השבועות גבוה מ-20%. העודף הממוצע בכל השבועות 44-13 הוא 13.9%.

עם זאת, נמצאו הבדלים לפי גילים. בגיל 19-0 התמותה בכל שנת 2020 נמוכה יותר מהתמותה בשנים קודמות בצורה משמעותית. בגיל 59-20 יש תנודתיות, אך סה"כ התמותה דומה לשנים הקודמות. בגיל 70 ומעלה, בדומה למגמה בכלל האוכלוסייה, קיים עודף ממוצע של 13% בשבועות 17-13. החל משבוע 31 חלה עלייה בעודף, ומשבוע 33 העודף גבוה מ-30%. בשבועיים האחרונים של אוקטובר העודף נמוך יותר.

בהשוואה בין ישראל למדינות רבות שחלקן מוצגות בתרשים 5, בולט כי עודף התמותה בישראל בחודשים הראשונים של מגפת הקורונה היה נמוך, אך משבוע 31 (סוף יולי) ובמיוחד משבוע 37 (תחילת ספטמבר) ועד שבוע 42 (אמצע אוקטובר) עודף התמותה בישראל גבוה יותר ממרבית המדינות. בסוף אוקטובר קיימת ירידה בולטת בעודף התמותה בישראל, לעומת כמה מדינות שבהן ישנה עלייה.

עודף תמותה לפי המודל של צוות "אקלים וקורונה"

כאמור, מדד ה-P-score המבוסס רק על מספר הפטירות ומשווה אותן למספר הפטירות בשנים קודמות, אינו מביא בחשבון שינויים בגודל האוכלוסייה ובהרכבה, ואינו מתחשב במגמה הקיימת בתמותה. ישראל מאופיינת בגידול אוכלוסייה גבוה יחסית למדינות אחרות, באוכלוסייה מזדקנת (כלומר חלקן של קבוצות הגיל המבוגרות עולה) ובמגמה כללית של ירידה בתמותה לאורך השנים. כדי לבחון באופן מדויק יותר את עודף התמותה בישראל, פותח מודל לחישוב מספר הנפטרים שהיה צפוי בשנת 2020 ללא מגפה בהתחשב בגורמים אלו, וכן בגורמי אקלים והימצאות שפעת. מספר הפטירות הנצפה הושווה אליו. המודל פותח על ידי קבוצת חוקרים בצוות "אקלים וקורונה", אחת הקבוצות שהוקמו במסגרת - Academia IL Collective Impact: Covid19.[3]

המודל[4] מתבסס על הנתונים האלה לשנים 2019-2000:

  1. פטירות יומיות לפי שש קבוצות גיל.
  2.  אומדני אוכלוסייה ממוצעת בכל שנה.
  3. נתוני אקלים יומיים מתחנת בית דגן – טמפרטורה - מקסימום ומינימום, קרינה - ממוצע ומקסימום יומי, ערכים מחושבים של לחות מוחלטת ושל נקודת הטל.[5]
  4. שיעור הפניות הגולמי בשבוע למרפאות "מכבי שירותי בריאות" בשל תחלואה דמוית שפעת כאינדיקציה לתחלואה בשפעת.

תוצאות המודל מצביעות על כך שמתחילת 2020 עד סוף אוקטובר נצפה עודף תמותה של כמעט 2,300 פטירות (6% לעומת הצפוי). ואולם אם בוחנים את העודף רק לאחר מקרה התמותה הראשון מנגיף קורונה (משבוע 13) עד סוף אוקטובר, העודף המתקבל הוא של מעט יותר מ-2,700 פטירות – 10.4% יותר מהצפוי. יצוין כי כל עודף התמותה ב-2020 הוא לאחר שבוע 13 שכן ב-12 השבועות הראשונים של השנה היה חוסר בתמותה של כ-440 פטירות (3.6%) לעומת הצפוי. משבוע 13 ועד סוף יולי עודף התמותה היה של כ-370 פטירות, כך שמרבית העודף הוא החל בחודשים אוגוסט – ספטמבר.

לוח א - עודף פטירות בכלל התקופה, שבועות 44-1, לפי גיל

(שבוע 1 מתחיל ב-30.12.19, שבוע 44 מתחיל ב-26.10.20)

קבוצת גיל מספר פטירות צפוי (רווח ניבוי 95%) מספר פטירות נצפה עודף פטירות
כלל האוכלוסייה 38,444 (38,060-38,829) 40,741 2,297*
19-0 807 (751-864) 638 169-*
59-20 4,003 (3,878-4,127) 4,137 134*
69-60 4,912 (4,774-5,050) 5,093 181*
79-70 7,717 (7,544-7,890) 8,485 768*
89-80 12,781 (12,559-13,003) 13,511 730*
90+ 8,224 (8,046-8,402) 8,877 653*

* מובהק סטטיסטית ברמת מובהקות 0.05>p

לוח ב- עודף פטירות לפני תחילת התמותה מנגיף קורונה, שבועות 12-1, לפי גיל

(שבוע 1 מתחיל ב-30.12.19,  שבוע 12 מתחיל ב-16.3.20)

קבוצת גיל מספר פטירות צפוי (רווח ניבוי 95%) מספר פטירות נצפה עודף פטירות
כלל האוכלוסייה 12,219 (12,002-12,436) 11,779 440-*
19-0 241 (210 – 272) 204 -37*
59-20 1,158 1091 – 1225)) 1,153 5-
69-60 1,533 (1,456-1,610) 1,459  74-
79-70 2,464 (2366-2,561) 2,323 141-*
89-80 4,162 4,035 – 4,289)) 3,968 194-*
90+ 2,661 2,559 – 2,762)) 2,672 11

* מובהק סטטיסטית ברמת מובהקות 0.05>p

לוח ג - עודף פטירות בתקופת הקורונה, שבועות 44-13 לפי גיל

(שבוע 13 מתחיל ב-23.3.20, שבוע 44 מתחיל ב-26.10.20)

קבוצת מספר פטירות צפוי מספר פטירות מספר פטירות עודף פטירות
גיל (רווח ניבוי 95%) נצפה מנגיף קורונה מספר אחוז
כלל האוכלוסייה 26,225 (25,908-26,543) 28,962 2,537 2,737* 10.4
19-0 566 (519-614) 434 3 132-* -23.3
59-20 2,845 (2,740-2,950) 2,984 169 139* 4.9
69-60 3,379 (3,265-3,494) 3,634 309 255* 7.5
79-70 5,253 (5,111-5,396) 6,162 620 909* 17.3
89-80 8,619 (8,436-8,801) 9,543 882 924* 10.7
90+ 5,564 (5,417-5,710) 6,205 552 641* 11.5

* מובהק סטטיסטית ברמת מובהקות 0.05>p

עודף התמותה שונה בקבוצות הגיל השונות והוא קיים בעיקר בקרב בני 70 ומעלה. בתקופת הקורונה (שבועות 44-13), קבוצת הגיל שבה העודף הוא הגדול ביותר היא 79-70 – עודף של 17.3%. בקרב בני 89-80 ו-90 ומעלה העודף גבוה מ-10%. בקרב צעירים עד גיל 20 היה חוסר תמותה גם בתקופת הקורונה (23%).

הפטירות מנגיף קורונה מסבירות כמעט את כל עודף התמותה, בכלל הגילים - 94% מהעודף מוסבר על ידי נגיף קורונה. אך גם כאן נמצא הבדל בין קבוצות הגיל:


[1] מקור נתוני הפטירות – רשות האוכלוסין וההגירה, מרשם האוכלוסין; נתוני אוקטובר אינם סופיים ועשויים לעלות עוד.

[2] Excess mortality during the Coronavirus pandemic (COVID-19)

[3] תשתית אקדמית ייחודית המשותפת לכלל אוניברסיטאות המחקר בישראל שהוקמה לאחר פרוץ מגפת הקורונה ומטרתה להעניק לקובעי מדיניות תשתית מחקרית רלוונטית ואמינה לקבלת החלטות מושכלות

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

[4] כדי לחקור את דפוסי התמותה היומית (עקומת תמותה יומית בסיסית) והקשר עם הגורמים המסבירים נבנה מודל של רגרסיה פואסונית מסוג GAM Generalized Additive Model )).  לכל אחת משש  קבוצות הגיל הותאם מודל רגרסיה בנפרד, כאשר המשתנה התלוי היה שיעור התמותה. כדי לבחון מגמות לאורך זמן בשיעורי התמותה הוכנסה השנה הכרונולוגית כמשתנה מסביר. תמותה מתאפיינת על ידי מחזוריות שנתית הדומה לזו של מזג האוויר, עם תמותה גבוהה בחורף ונמוכה בקיץ. מחזוריות זו קשורה לשינויי מזג האוויר ולתופעות התחלואה המתלוות אליהם, אך היא אינה חופפת להם בדיוק ולכן חושב לה מודל בנפרד. מעבר למחזוריות, ייתכנו גם קשרים ממוקדים עם מזג אוויר לא שגרתי, למשל גל חום קיצוני, במיוחד בעונות המעבר, הקשור לתוספת תמותה. לכן, המודל כולל גם אפשרות להשפעות מסוג זה. המודל משקף גם תנודות שבועיות בשיעורי תמותה, עם ירידה ברורה בשבת וירידה מתונה יותר ביום שישי. כדי  לנבא את מספר הפטירות הצפוי בשנת 2020, נעשה שימוש במודלים הבסיסיים שנבנו. הניבוי חושב לכל יום ולכל קבוצת גיל. הניבויים היומיים סוכמו כדי לקבל ניבויים שבועיים וכן חושבו רווחי ניבוי ברמת סמך של 95%. בכל שבוע חושב הפער בין הערך הצפוי לערך הנצפה וחושב סיכום תקופתי. ניבויים שבועיים לכלל האוכלוסייה חושבו מסכום הניבויים של כל קבוצות הגיל.

[5] הטמפרטורה שבה מתקרר האוויר בלחץ קבוע ובתכולת אדי מים קבועה עד הגיעו לרוויה, כלומר למצב שבו אדי המים חייבים לעבור עיבוי ולהפוך מגז לנוזל.