Introduction

This is the Somalia December Sorghum forecast report. The goal is to assist users in evaluating current forecasts in the context of past forecasts and to translate the forecast component into key assumptions about food security.

The report includes the following key components:

  1. A verbal summary of assumptions based on the statistics in this forecast

  2. Mean Area, Production, and Yield over the years 2012 to 2022:This provides context for interpreting the grain data

  3. Historical Out of Sample Forecast Error Averaged Over a 10 Year+ Period. We show the Mean Absolute Percent Error (MAPE).

  4. Yield Forecast, (based on December 2021 Precip, NDVI, Et0 ) expressed as percent of mean yield over the period 2012 - 2022.

  5. Forecasts and Forecast Error in Analog Years We show specific forecasts and forecast error for years were climatologically similar to those in this report.

Assumption Statements [NOT YET UPDATED FOR THIS OND REPORT]

Average Forecast Error For This Point in the Season

At this point in the season historical forecast error, is on average, below 50% for 46 out of the 46 admin units used in the forecast.



Summary Figures

Mean area, production, and yields for the years 2012 to 2022.



Out of Sample Forecast Error (MAPE)

Mean Absolute Percent Error (MAPE) calculated based on historical out of sample seasonal forecasts. Lower scores indicate greater accuracy. Forecasts are based on model type MODEL1



Yield Forecast for December 2021

Forecast values expressed as percent of Mean Yields over the Years 2012 - 2022.

The figure shows predicted percent of mean (center) as well as lower (left) and higher (right) predicted percent of mean intervals.



Static and Dynamic Version of Main Forecast

Roll over the polygon borders to get the district name and % of mean forecast value.

Static Version

This map shows the main % of mean forecast value along with district labels for reference.

Discrete Map

This map shows the main % of mean forecast binned into discrete values. Averages are based on the most recent 10 year period of observed yields: 2012 - 2022

This table shows the forecast percentage of mean values in the above table along with the mean yield values from the first figure.


Table of Mean Yields and Predicted Percent of Mean Values
District % of mean % of mean (low) % of mean (high)
Hargeysa 1.29 0.88 1.70
Gabiley 1.29 0.95 1.63
Marka 0.79 0.50 1.08
Afgooye 1.12 0.79 1.46
Qoryooley 1.08 0.89 1.27
Wanla Weyn 1.15 0.94 1.37
Baydhabo 0.89 0.65 1.12
Buur Hakaba 0.95 0.58 1.32
Diinsoor 0.82 0.53 1.11
Qansax Dheere 0.66 0.32 1.00
Tiyeglow 1.22 0.83 1.61
Waajid 0.79 0.33 1.25
Garbahaarey 1.30 1.22 1.38
Baardheere 0.79 0.44 1.14
Bu'aale 0.35 0.00 0.70
Saakow 0.61 0.22 1.01
Afmadow 0.71 0.54 0.88

Analog Year Forecasts

Yield forecasts in analog years. <–DESCRIPTION OF ANALOG YEAR PROCESS–>.
***


Analog Year Forecasts Errors

Forecast errors in analog years. If observed data is not available in a given year we cannot calculate forecast errors. Values are expressed a percentage of observed yields in a given year (t):


\[\frac{(observed_{(t)}-forecast_{(t)})}{observed_{(t)}}\] ***

Positive (+) values indicate an under prediction. Negative (-) values indicate an over prediction.



  1. Extra/Extended Trees. A type of Random Forest Model↩︎

LS0tCnRpdGxlOiAiRm9yZWNhc3QgUmVwb3J0IHdpdGggQW5hbG9nIFllYXJzLVNvbWFsaWEgT05EIgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazoKICAgIHRvYzogeWVzCmZpZ193aWR0aDogNwpmaWdfaGVpZ2h0OiA2CmZpZ19jYXB0aW9uOiB0cnVlCi0tLQoKYGBge3IsZXZhbD1UUlVFLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFLHJlc3VsdHM9J2hpZGUnfQojLS0tLS0tLS0tLS0tLS0tLS0tLUJhc2UgU2V0dXAtLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQpybShsaXN0PWxzKCkpCgojLS0tU2V0IFByb2plY3QgRGlyZWN0b3JpZXMKZGlyQmFzZTwtJy9Wb2x1bWVzL0dvb2dsZURyaXZlL015IERyaXZlLycKCmRpckJhc2UyPC0nL1ZvbHVtZXMvR29vZ2xlRHJpdmUvU2hhcmVkIGRyaXZlcy9DSEMgVGVhbSBEcml2ZSAvJwoKIy1Qcm9qZWN0IERpcmVjdG9yaWVzCmRpclByb2o8LXBhc3RlMChkaXJCYXNlMiwncHJvamVjdF9tYWNoaW5lX2xlYXJuaW5nX2ZvcmVjYXN0aW5nLycpICNwcm9qZWN0IGRpcmVjdG9yeQoKZGlyVmlld2VyPC1wYXN0ZTAoZGlyUHJvaiwndmlld2VyLycpCmRpclZpZXdlck91dFN0YXRpYzwtcGFzdGUwKGRpclZpZXdlciwndmlld2VyX3N0YXRpY19zaGFwZXMvJykKZGlyVmlld2VyRHluYW1pYzwtcGFzdGUwKGRpclZpZXdlciwndmlld2VyX2R5bmFtaWNfc2hhcGVzLycpCgpkaXJSZXBvcnQ8LXBhc3RlMChkaXJQcm9qLCdmb3JlY2FzdF9yZXBvcnRpbmcvJykKZGlyUmVwb3J0UmRhdGE8LXBhc3RlMChkaXJSZXBvcnQsJ2ZvcmVjYXN0X3JlcG9ydGluZ19SZGF0YS8nKQoKbGlicmFyeShzdHJpbmdyKQpsaWJyYXJ5KGdncGxvdDIpCmxpYnJhcnkoZHBseXIpCmxpYnJhcnkocmFzdGVyKQpsaWJyYXJ5KHJnZGFsKQpsaWJyYXJ5KG1nY3YpCmxpYnJhcnkodGlkeXIpCmxpYnJhcnkobHVicmlkYXRlKQpsaWJyYXJ5KHNmKQpsaWJyYXJ5KHJtYXBzaGFwZXIpCmxpYnJhcnkodmlyaWRpcykKbGlicmFyeShzY2FsZXMpCmxpYnJhcnkocGxvdGx5KQpsaWJyYXJ5KGZvcmNhdHMpCmxpYnJhcnkoa25pdHIpCmxpYnJhcnkoa2FibGVFeHRyYSkKbGlicmFyeShzaGlueSkKIz09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0KCiNQYXJhbWV0ZXJzCkNVUlJFTlRfWUVBUjwtMjAyMQpNT05USDwtMTIKREVLQUQ8LTMKTU9ERUw8LSdFVCcKQ09VTlRSWTwtJ1NvbWFsaWEnClBST0RVQ1Q8LSdTb3JnaHVtJwojQ1JPUF9BUkVBPC0yICAjUGVyY2VudApBTkFMT0dfWUVBUlM8LWMoMTk5OSwyMDAwLDIwMDEsMjAwOCwyMDA5LDIwMTEsMjAxMiwyMDE3KQoKbW9udGhfbmFtZTwtbW9udGgubmFtZVtNT05USF0gI21vbnRoIHRoZSBwcm9kdWN0IGlzIGJhc2VkIG9uCgojLS1Mb2FkIEV4aXN0aW5nIFBsb3RzCnNldHdkKGRpclJlcG9ydFJkYXRhKQpsb2FkKGZpbGU9cGFzdGUwKCcyMF9mb3JlY2FzdF9yZXBvcnRpbmdfbWFpbl9wbG90cycsbW9udGhfbmFtZSxDT1VOVFJZLCdfJyxQUk9EVUNULCcuUmRhdGEnKSkKbG9hZChmaWxlPXBhc3RlMCgnMDFfZm9yZWNhc3RfcmVwb3J0X2Fnc3RhdG1hcHNfJyxDT1VOVFJZLCdfJyxQUk9EVUNULCcuUmRhdGEnKSkKYGBgCiMgSW50cm9kdWN0aW9uClRoaXMgaXMgdGhlIGByIENPVU5UUllgIGByIG1vbnRoX25hbWVgIGByIFBST0RVQ1RgIGZvcmVjYXN0IHJlcG9ydC4gKipUaGUgZ29hbCBpcyB0byBhc3Npc3QgdXNlcnMgaW4gZXZhbHVhdGluZyBjdXJyZW50IGZvcmVjYXN0cyBpbiB0aGUgY29udGV4dCBvZiBwYXN0IGZvcmVjYXN0cyoqIGFuZCAqKnRvIHRyYW5zbGF0ZSB0aGUgZm9yZWNhc3QgY29tcG9uZW50IGludG8ga2V5IGFzc3VtcHRpb25zIGFib3V0IGZvb2Qgc2VjdXJpdHkuKioKClRoZSByZXBvcnQgaW5jbHVkZXMgdGhlIGZvbGxvd2luZyBrZXkgY29tcG9uZW50czoKCjEuICoqKkEgdmVyYmFsIHN1bW1hcnkgb2YgYXNzdW1wdGlvbnMgYmFzZWQgb24gdGhlIHN0YXRpc3RpY3MgaW4gdGhpcyBmb3JlY2FzdCoqKiAKCjIuICoqKk1lYW4gQXJlYSwgUHJvZHVjdGlvbiwgYW5kIFlpZWxkIG92ZXIgdGhlIHllYXJzIGByIGxpc192YXJzX3JlcG9ydCRtaW5fYWdgIHRvIGByIGxpc192YXJzX3JlcG9ydCRtYXhfYWdgOioqKlRoaXMgcHJvdmlkZXMgY29udGV4dCBmb3IgaW50ZXJwcmV0aW5nIHRoZSBncmFpbiBkYXRhCgozLiAqKipIaXN0b3JpY2FsIE91dCBvZiBTYW1wbGUgRm9yZWNhc3QgRXJyb3IgQXZlcmFnZWQgT3ZlciBhIDEwIFllYXIrIFBlcmlvZCoqKi4gV2Ugc2hvdyB0aGUgKioqTSoqKmVhbiAqKipBKioqYnNvbHV0ZSAqKipQKioqZXJjZW50ICoqKkUqKipycm9yICgqKipNQVBFKioqKS4KCjQuICoqKllpZWxkIEZvcmVjYXN0KioqLCAoYmFzZWQgb24gYHIgbGlzX3ZhcnNfcmVwb3J0JG1vbnRoX25hbWVgIGByIGxpc192YXJzX3JlcG9ydCRtYXhfZXZhcl95ZWFyYCBgciBsaXNfdmFyc19yZXBvcnQkdmFyX25hbWVgICkgZXhwcmVzc2VkIGFzICoqcGVyY2VudCBvZiBtZWFuIHlpZWxkKiogb3ZlciB0aGUgcGVyaW9kIGByIGxpc192YXJzX3JlcG9ydCRtaW5fYWdgIC0gYHIgbGlzX3ZhcnNfcmVwb3J0JG1heF9hZ2AuCgo1LiAqKipGb3JlY2FzdHMgYW5kIEZvcmVjYXN0IEVycm9yIGluIEFuYWxvZyBZZWFycyoqKiBXZSBzaG93IHNwZWNpZmljIGZvcmVjYXN0cyBhbmQgZm9yZWNhc3QgZXJyb3IgZm9yIHllYXJzIHdlcmUgY2xpbWF0b2xvZ2ljYWxseSBzaW1pbGFyIHRvIHRob3NlIGluIHRoaXMgcmVwb3J0LgoKIyBBc3N1bXB0aW9uIFN0YXRlbWVudHMgW05PVCBZRVQgVVBEQVRFRCBGT1IgVEhJUyBPTkQgUkVQT1JUXQoKIyMjIyBBdmVyYWdlIEZvcmVjYXN0IEVycm9yIEZvciBUaGlzIFBvaW50IGluIHRoZSBTZWFzb24KQXQgdGhpcyBwb2ludCBpbiB0aGUgc2Vhc29uICoqKmhpc3RvcmljYWwgZm9yZWNhc3QgZXJyb3IsIGlzIG9uIGF2ZXJhZ2UsIGJlbG93IDUwJSBmb3IgNDYgb3V0IG9mIHRoZSA0NiBhZG1pbiB1bml0cyoqKiB1c2VkIGluIHRoZSBmb3JlY2FzdC4KCioqKgoKCgoqKioKCiMgU3VtbWFyeSBGaWd1cmVzCgojIyMjIE1lYW4gYXJlYSwgcHJvZHVjdGlvbiwgYW5kIHlpZWxkcyBmb3IgdGhlIHllYXJzIGByIGxpc192YXJzX3JlcG9ydCRtaW5fYWdgIHRvIGByIGxpc192YXJzX3JlcG9ydCRtYXhfYWdgLgoKKioqCmBgYHtyLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFfQpwMWFsbApgYGAKKioqCgoKCgoKIyMjIyBPdXQgb2YgU2FtcGxlIEZvcmVjYXN0IEVycm9yIChNQVBFKQpNZWFuIEFic29sdXRlIFBlcmNlbnQgRXJyb3IgKE1BUEUpIGNhbGN1bGF0ZWQgYmFzZWQgb24gaGlzdG9yaWNhbCBvdXQgb2Ygc2FtcGxlIHNlYXNvbmFsIGZvcmVjYXN0cy4gKipMb3dlciBzY29yZXMgaW5kaWNhdGUgZ3JlYXRlciBhY2N1cmFjeSoqLiBGb3JlY2FzdHMgYXJlIGJhc2VkIG9uIG1vZGVsIHR5cGUgYE1PREVMYF5bRXh0cmEvRXh0ZW5kZWQgVHJlZXMuIEEgdHlwZSBvZiBSYW5kb20gRm9yZXN0IE1vZGVsXQoKKioqCmBgYHtyLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFfQpwMgpgYGAKKioqCgoKIyBZaWVsZCBGb3JlY2FzdCBmb3IgYHIgbGlzX3ZhcnNfcmVwb3J0JG1vbnRoX25hbWVgIGByIGxpc192YXJzX3JlcG9ydCRtYXhfZXZhcl95ZWFyYCAKRm9yZWNhc3QgdmFsdWVzIGV4cHJlc3NlZCBhcyBwZXJjZW50IG9mIE1lYW4gWWllbGRzIG92ZXIgdGhlIFllYXJzIGByIGxpc192YXJzX3JlcG9ydCRtaW5fYWdgIC0gYHIgbGlzX3ZhcnNfcmVwb3J0JG1heF9hZ2AuIAoKVGhlIGZpZ3VyZSBzaG93cyAqKipwcmVkaWN0ZWQgcGVyY2VudCBvZiBtZWFuKioqIChjZW50ZXIpIGFzIHdlbGwgYXMgbG93ZXIgKGxlZnQpIGFuZCBoaWdoZXIgKHJpZ2h0KSAqKipwcmVkaWN0ZWQgcGVyY2VudCBvZiBtZWFuKioqIGludGVydmFscy4KCioqKgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRSxmaWcud2lkdGg9OCxmaWcuaGVpZ2h0PTh9CiNwMzwtZ2dwbG90bHkocD1wMyx0b29sdGlwPWMoJ2Rpc3RyaWN0JywncGVyX29mX21lYW4nKSkKcDEKYGBgCioqKgoKIyMjIyBTdGF0aWMgYW5kIER5bmFtaWMgVmVyc2lvbiBvZiBNYWluIEZvcmVjYXN0ClJvbGwgb3ZlciB0aGUgcG9seWdvbiBib3JkZXJzIHRvIGdldCB0aGUgZGlzdHJpY3QgbmFtZSBhbmQgJSBvZiBtZWFuIGZvcmVjYXN0IHZhbHVlLgoKYGBge3IsZWNobz1GQUxTRSx3YXJuaW5nPUZBTFNFLG1lc3NhZ2U9RkFMU0UsZmlnLmFsaWduPSdjZW50ZXInLGZpZy53aWR0aD04LGZpZy5oZWlnaHQ9OH0KI3AxTHMKcDFMPC1nZ3Bsb3RseShwMUwsdG9vbHRpcD1jKCdhZG1pbjInLCd2YWx1ZScpKSAlPiUgbGF5b3V0KGxlZ2VuZCA9IGxpc3Qob3JpZW50YXRpb24gPSAiaCIsIHggPSAwLjQsIHkgPSAtMC4yKSkKcDFMCmBgYAojIyMjIFN0YXRpYyBWZXJzaW9uClRoaXMgbWFwIHNob3dzIHRoZSBtYWluICUgb2YgbWVhbiBmb3JlY2FzdCB2YWx1ZSBhbG9uZyB3aXRoIGRpc3RyaWN0IGxhYmVscyBmb3IgcmVmZXJlbmNlLgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRSxmaWcuYWxpZ249J2NlbnRlcicsZmlnLndpZHRoPTgsZmlnLmhlaWdodD04fQpwMUxzCgpgYGAKIyMjIyBEaXNjcmV0ZSBNYXAKVGhpcyBtYXAgc2hvd3MgdGhlIG1haW4gJSBvZiBtZWFuIGZvcmVjYXN0IGJpbm5lZCBpbnRvIGRpc2NyZXRlIHZhbHVlcy4gKipBdmVyYWdlcyoqIGFyZSBiYXNlZCBvbiB0aGUgbW9zdCByZWNlbnQgMTAgeWVhciBwZXJpb2Qgb2Ygb2JzZXJ2ZWQgeWllbGRzOiBgciBsaXNfdmFyc19yZXBvcnQkbWluX2FnYCAtIGByIGxpc192YXJzX3JlcG9ydCRtYXhfYWdgCmBgYHtyLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFLGZpZy5hbGlnbj0nY2VudGVyJyxmaWcud2lkdGg9OCxmaWcuaGVpZ2h0PTh9CnAxTGMKCmBgYAoKIyMjIyBUaGlzIHRhYmxlIHNob3dzIHRoZSBmb3JlY2FzdCBwZXJjZW50YWdlIG9mIG1lYW4gdmFsdWVzIGluIHRoZSBhYm92ZSB0YWJsZSBhbG9uZyB3aXRoIHRoZSBtZWFuIHlpZWxkIHZhbHVlcyBmcm9tIHRoZSBmaXJzdCBmaWd1cmUuIAoqKioKYGBge3IsZWNobz1GQUxTRSx3YXJuaW5nPUZBTFNFLG1lc3NhZ2U9RkFMU0V9CnQxPC1rbml0cjo6a2FibGUoZHRhYiwgY2FwdGlvbiA9ICdUYWJsZSBvZiBNZWFuIFlpZWxkcyBhbmQgUHJlZGljdGVkIFBlcmNlbnQgb2YgTWVhbiBWYWx1ZXMnKQp0MTwta2FibGVfc3R5bGluZyh0MSxib290c3RyYXBfb3B0aW9ucyA9IGMoInN0cmlwZWQiLCAiaG92ZXIiLCJjb25kZW5zZWQiKSkKc2Nyb2xsX2JveCh0MSwgaGVpZ2h0ID0gJzMwMHB4Jywgd2lkdGggPSAnMTAwJScsCiAgYm94X2NzcyA9ICJib3JkZXI6IDFweCBzb2xpZCAjZGRkOyBwYWRkaW5nOiAxcHg7ICIsIGV4dHJhX2NzcyA9IE5VTEwsCiAgZml4ZWRfdGhlYWQgPSBUUlVFKQoKCmBgYAoqKioKCiMgQW5hbG9nIFllYXIgRm9yZWNhc3RzCllpZWxkIGZvcmVjYXN0cyBpbiBhbmFsb2cgeWVhcnMuIDwtLURFU0NSSVBUSU9OIE9GIEFOQUxPRyBZRUFSIFBST0NFU1MtLT4uICAKKioqCmBgYHtyLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFLGZpZy53aWR0aD03LjUsZmlnLmhlaWdodD03LjV9CnAzCmBgYAoqKioKCiMgQW5hbG9nIFllYXIgRm9yZWNhc3RzIEVycm9ycwpGb3JlY2FzdCBlcnJvcnMgaW4gYW5hbG9nIHllYXJzLiBJZiBvYnNlcnZlZCBkYXRhIGlzIG5vdCBhdmFpbGFibGUgaW4gYSBnaXZlbiB5ZWFyIHdlIGNhbm5vdCBjYWxjdWxhdGUgZm9yZWNhc3QgZXJyb3JzLiBWYWx1ZXMgYXJlIGV4cHJlc3NlZCBhIHBlcmNlbnRhZ2UgIG9mIG9ic2VydmVkIHlpZWxkcyBpbiBhIGdpdmVuIHllYXIgXyh0KV86CgoqKioKJCRcZnJhY3sob2JzZXJ2ZWRfeyh0KX0tZm9yZWNhc3Rfeyh0KX0pfXtvYnNlcnZlZF97KHQpfX0kJAoqKioKCioqKlBvc2l0aXZlICgrKSoqKiB2YWx1ZXMgaW5kaWNhdGUgYW4gX3VuZGVyIHByZWRpY3Rpb25fLiAqKipOZWdhdGl2ZSAoLSkqKiogdmFsdWVzIGluZGljYXRlIGFuIF9vdmVyIHByZWRpY3Rpb25fLgoKKioqCmBgYHtyLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFLGZpZy53aWR0aD03LjUsZmlnLmhlaWdodD03LjV9CnA0CmBgYA==