Introduction

This is a protype forecast report. The goal is to assist users in evaluating current forecasts in the context of past forecasts.

The report includes the following key components:

  1. Mean Area, Production, and Yield over the years 2007 to 2017:This provides context for interpreting the grain data

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

  3. Yield Forecast, (based on May 2021 Precip, NDVI, Et0 ) expressed as percent of mean yield over the period 2007 - 2017.

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

Summary Figures

Mean area, production, and yields for the years 2007 to 2017.



Out of Sample Forecast Error (MAPE)

Mean Absolute Percent Error (MAPE) calculated based on out of sample seasonfal forecasts conducted over 2001-2014. Lower scores indicate greater accuracy. Forecasts are based on model type MODEL1



Yield Forecast for May 2021

Forecast values expressed as percent of Mean Yields over the Years 2007 - 2017.

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 lables for reference.

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)
Awdal 0.77 0.23 1.32
Hiiraan 0.80 0.54 1.05
Middle Shabelle 0.87 0.54 1.20
Lower Shabelle 0.83 0.69 0.97
Gedo 0.99 0.84 1.14
Middle Juba 0.80 0.43 1.17
Lower Juba 0.67 0.45 0.90

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↩︎

LS0tCnRpdGxlOiAiRm9yZWNhc3QgUmVwb3J0IHdpdGggQW5hbG9nIFllYXJzLVNvbWFsaWEiCm91dHB1dDoKICBodG1sX25vdGVib29rOgogICAgdG9jOiB5ZXMKZmlnX3dpZHRoOiA3CmZpZ19oZWlnaHQ6IDYKZmlnX2NhcHRpb246IHRydWUKLS0tCgpgYGB7cixldmFsPVRSVUUsZWNobz1GQUxTRSx3YXJuaW5nPUZBTFNFLG1lc3NhZ2U9RkFMU0V9CiMtLS0tLS0tLS0tLS0tLS0tLS0tQmFzZSBTZXR1cC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCnJtKGxpc3Q9bHMoKSkKCiMtLS1TZXQgUHJvamVjdCBEaXJlY3RvcmllcwpkaXJCYXNlPC0nL1ZvbHVtZXMvR29vZ2xlRHJpdmUvTXkgRHJpdmUvJwoKZGlyQmFzZTI8LScvVm9sdW1lcy9Hb29nbGVEcml2ZS9TaGFyZWQgZHJpdmVzL0NIQyBUZWFtIERyaXZlIC8nCgojLVByb2plY3QgRGlyZWN0b3JpZXMKZGlyUHJvajwtcGFzdGUwKGRpckJhc2UyLCdwcm9qZWN0X21hY2hpbmVfbGVhcm5pbmdfZm9yZWNhc3RpbmcvJykgI3Byb2plY3QgZGlyZWN0b3J5CgpkaXJWaWV3ZXI8LXBhc3RlMChkaXJQcm9qLCd2aWV3ZXIvJykKZGlyVmlld2VyT3V0U3RhdGljPC1wYXN0ZTAoZGlyVmlld2VyLCd2aWV3ZXJfc3RhdGljX3NoYXBlcy8nKQpkaXJWaWV3ZXJEeW5hbWljPC1wYXN0ZTAoZGlyVmlld2VyLCd2aWV3ZXJfZHluYW1pY19zaGFwZXMvJykKCmRpclJlcG9ydDwtcGFzdGUwKGRpclByb2osJ2ZvcmVjYXN0X3JlcG9ydGluZy8nKQpkaXJSZXBvcnRSZGF0YTwtcGFzdGUwKGRpclJlcG9ydCwnZm9yZWNhc3RfcmVwb3J0aW5nX1JkYXRhLycpCgpsaWJyYXJ5KHN0cmluZ3IpCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShkcGx5cikKbGlicmFyeShyYXN0ZXIpCmxpYnJhcnkocmdkYWwpCmxpYnJhcnkobWdjdikKbGlicmFyeSh0aWR5cikKbGlicmFyeShsdWJyaWRhdGUpCmxpYnJhcnkoc2YpCmxpYnJhcnkocm1hcHNoYXBlcikKbGlicmFyeSh2aXJpZGlzKQpsaWJyYXJ5KHNjYWxlcykKbGlicmFyeShwbG90bHkpCmxpYnJhcnkoZm9yY2F0cykKbGlicmFyeShrbml0cikKbGlicmFyeShrYWJsZUV4dHJhKQpsaWJyYXJ5KHNoaW55KQojPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PQoKI1BhcmFtZXRlcnMKQ1VSUkVOVF9ZRUFSPC0yMDIxCk1PTlRIPC01CkRFS0FEPC0zCk1PREVMPC0nRVQnCkNPVU5UUlk8LSdTb21hbGlhJwojQ1JPUF9BUkVBPC0yICAjUGVyY2VudApBTkFMT0dfWUVBUlM8LWMoMTk5OSwyMDAwLDIwMDEsMjAwOCwyMDA5LDIwMTEsMjAxMiwyMDE3KQoKbW9udGhfbmFtZTwtbW9udGgubmFtZVtNT05USF0gI21vbnRoIHRoZSBwcm9kdWN0IGlzIGJhc2VkIG9uCgojLS1Mb2FkIEV4aXN0aW5nIFBsb3RzCnNldHdkKGRpclJlcG9ydFJkYXRhKQpsb2FkKGZpbGU9cGFzdGUwKCcyMF9mb3JlY2FzdF9yZXBvcnRpbmdfbWFpbl9wbG90cycsbW9udGhfbmFtZSxDT1VOVFJZLCcuUmRhdGEnKSkKbG9hZChmaWxlPXBhc3RlMCgnMDFfZm9yZWNhc3RfcmVwb3J0X2Fnc3RhdG1hcHNfJyxDT1VOVFJZLCcuUmRhdGEnKSkKYGBgCiMgSW50cm9kdWN0aW9uClRoaXMgaXMgYSBwcm90eXBlIGZvcmVjYXN0IHJlcG9ydC4gKipUaGUgZ29hbCBpcyB0byBhc3Npc3QgdXNlcnMgaW4gZXZhbHVhdGluZyBjdXJyZW50IGZvcmVjYXN0cyBpbiB0aGUgY29udGV4dCBvZiBwYXN0IGZvcmVjYXN0cy4gKioKClRoZSByZXBvcnQgaW5jbHVkZXMgdGhlIGZvbGxvd2luZyBrZXkgY29tcG9uZW50czoKCjEuICoqKk1lYW4gQXJlYSwgUHJvZHVjdGlvbiwgYW5kIFlpZWxkIG92ZXIgdGhlIHllYXJzIGByIGxpc192YXJzX3JlcG9ydCRtaW5fYWdgIHRvIGByIGxpc192YXJzX3JlcG9ydCRtYXhfYWdgOioqKlRoaXMgcHJvdmlkZXMgY29udGV4dCBmb3IgaW50ZXJwcmV0aW5nIHRoZSBncmFpbiBkYXRhCgoyLiAqKipIaXN0b3JpY2FsIE91dCBvZiBTYW1wbGUgRm9yZWNhc3QgRXJyb3IgQXZlcmFnZWQgT3ZlciBhIDEwIFllYXIrIFBlcmlvZCoqKi4gV2Ugc2hvdyB0aGUgKioqTSoqKmVhbiAqKipBKioqYnNvbHV0ZSAqKipQKioqZXJjZW50ICoqKkUqKipycm9yICgqKipNQVBFKioqKS4KCjMuICoqKllpZWxkIEZvcmVjYXN0KioqLCAoYmFzZWQgb24gYHIgbGlzX3ZhcnNfcmVwb3J0JG1vbnRoX25hbWVgIGByIGxpc192YXJzX3JlcG9ydCRtYXhfZXZhcl95ZWFyYCBgciBsaXNfdmFyc19yZXBvcnQkdmFyX25hbWVgICkgZXhwcmVzc2VkIGFzICoqcGVyY2VudCBvZiBtZWFuIHlpZWxkKiogb3ZlciB0aGUgcGVyaW9kIGByIGxpc192YXJzX3JlcG9ydCRtaW5fYWdgIC0gYHIgbGlzX3ZhcnNfcmVwb3J0JG1heF9hZ2AuCgo0LiAqKipGb3JlY2FzdHMgYW5kIEZvcmVjYXN0IEVycm9yIGluIEFuYWxvZyBZZWFycyoqKiBXZSBzaG93IHNwZWNpZmljIGZvcmVjYXN0cyBhbmQgZm9yZWNhc3QgZXJyb3IgZm9yIHllYXJzIHdlcmUgY2xpbWF0b2xvZ2ljYWxseSBzaW1pbGFyIHRvIHRob3NlIGluIHRoaXMgcmVwb3J0LgoKCgojIFN1bW1hcnkgRmlndXJlcwoKIyMjIyBNZWFuIGFyZWEsIHByb2R1Y3Rpb24sIGFuZCB5aWVsZHMgZm9yIHRoZSB5ZWFycyBgciBsaXNfdmFyc19yZXBvcnQkbWluX2FnYCB0byBgciBsaXNfdmFyc19yZXBvcnQkbWF4X2FnYC4KCioqKgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRX0KcDFhbGwKYGBgCioqKgoKCgoKCiMjIyMgT3V0IG9mIFNhbXBsZSBGb3JlY2FzdCBFcnJvciAoTUFQRSkKTWVhbiBBYnNvbHV0ZSBQZXJjZW50IEVycm9yIChNQVBFKSBjYWxjdWxhdGVkIGJhc2VkIG9uIG91dCBvZiBzYW1wbGUgc2Vhc29uZmFsIGZvcmVjYXN0cyBjb25kdWN0ZWQgb3ZlciAyMDAxLTIwMTQuICoqTG93ZXIgc2NvcmVzIGluZGljYXRlIGdyZWF0ZXIgYWNjdXJhY3kqKi4gRm9yZWNhc3RzIGFyZSBiYXNlZCBvbiBtb2RlbCB0eXBlIGBNT0RFTGBeW0V4dHJhL0V4dGVuZGVkIFRyZWVzLiBBIHR5cGUgb2YgUmFuZG9tIEZvcmVzdCBNb2RlbF0KCioqKgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRX0KcDIKYGBgCioqKgoKCiMgWWllbGQgRm9yZWNhc3QgZm9yIGByIGxpc192YXJzX3JlcG9ydCRtb250aF9uYW1lYCBgciBsaXNfdmFyc19yZXBvcnQkbWF4X2V2YXJfeWVhcmAgCkZvcmVjYXN0IHZhbHVlcyBleHByZXNzZWQgYXMgcGVyY2VudCBvZiBNZWFuIFlpZWxkcyBvdmVyIHRoZSBZZWFycyBgciBsaXNfdmFyc19yZXBvcnQkbWluX2FnYCAtIGByIGxpc192YXJzX3JlcG9ydCRtYXhfYWdgLiAKClRoZSBmaWd1cmUgc2hvd3MgKioqcHJlZGljdGVkIHBlcmNlbnQgb2YgbWVhbioqKiAoY2VudGVyKSBhcyB3ZWxsIGFzIGxvd2VyIChsZWZ0KSBhbmQgaGlnaGVyIChyaWdodCkgKioqcHJlZGljdGVkIHBlcmNlbnQgb2YgbWVhbioqKiBpbnRlcnZhbHMuCgoqKioKYGBge3IsZWNobz1GQUxTRSx3YXJuaW5nPUZBTFNFLG1lc3NhZ2U9RkFMU0V9CiNwMzwtZ2dwbG90bHkocD1wMyx0b29sdGlwPWMoJ2Rpc3RyaWN0JywncGVyX29mX21lYW4nKSkKcDEKYGBgCioqKgoKIyMjIyBTdGF0aWMgYW5kIER5bmFtaWMgVmVyc2lvbiBvZiBNYWluIEZvcmVjYXN0ClJvbGwgb3ZlciB0aGUgcG9seWdvbiBib3JkZXJzIHRvIGdldCB0aGUgZGlzdHJpY3QgbmFtZSBhbmQgJSBvZiBtZWFuIGZvcmVjYXN0IHZhbHVlLgoKYGBge3IsZWNobz1GQUxTRSx3YXJuaW5nPUZBTFNFLG1lc3NhZ2U9RkFMU0UsZmlnLmFsaWduPSdjZW50ZXInLGZpZy53aWR0aD04LGZpZy5oZWlnaHQ9OH0KI3AxTHMKcDFMPC1nZ3Bsb3RseShwMUwsdG9vbHRpcD1jKCdhZG1pbjEnLCd2YWx1ZScpKSAlPiUgbGF5b3V0KGxlZ2VuZCA9IGxpc3Qob3JpZW50YXRpb24gPSAiaCIsIHggPSAwLjQsIHkgPSAtMC4yKSkKcDFMCmBgYAojIyMjIFN0YXRpYyBWZXJzaW9uClRoaXMgbWFwIHNob3dzIHRoZSBtYWluICUgb2YgbWVhbiBmb3JlY2FzdCB2YWx1ZSBhbG9uZyB3aXRoIGRpc3RyaWN0IGxhYmxlcyBmb3IgcmVmZXJlbmNlLgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRSxmaWcuYWxpZ249J2NlbnRlcicsZmlnLndpZHRoPTgsZmlnLmhlaWdodD04fQpwMUxzCgpgYGAKCiMjIyMgVGhpcyB0YWJsZSBzaG93cyB0aGUgZm9yZWNhc3QgcGVyY2VudGFnZSBvZiBtZWFuIHZhbHVlcyBpbiB0aGUgYWJvdmUgdGFibGUgYWxvbmcgd2l0aCB0aGUgbWVhbiB5aWVsZCB2YWx1ZXMgZnJvbSB0aGUgZmlyc3QgZmlndXJlLiAKKioqCmBgYHtyLGVjaG89RkFMU0Usd2FybmluZz1GQUxTRSxtZXNzYWdlPUZBTFNFfQp0MTwta25pdHI6OmthYmxlKGR0YWIsIGNhcHRpb24gPSAnVGFibGUgb2YgTWVhbiBZaWVsZHMgYW5kIFByZWRpY3RlZCBQZXJjZW50IG9mIE1lYW4gVmFsdWVzJykKdDE8LWthYmxlX3N0eWxpbmcodDEsYm9vdHN0cmFwX29wdGlvbnMgPSBjKCJzdHJpcGVkIiwgImhvdmVyIiwiY29uZGVuc2VkIikpCnNjcm9sbF9ib3godDEsIGhlaWdodCA9ICczMDBweCcsIHdpZHRoID0gJzEwMCUnLAogIGJveF9jc3MgPSAiYm9yZGVyOiAxcHggc29saWQgI2RkZDsgcGFkZGluZzogMXB4OyAiLCBleHRyYV9jc3MgPSBOVUxMLAogIGZpeGVkX3RoZWFkID0gVFJVRSkKCgpgYGAKKioqCgojIEFuYWxvZyBZZWFyIEZvcmVjYXN0cwpZaWVsZCBmb3JlY2FzdHMgaW4gYW5hbG9nIHllYXJzLiA8LS1ERVNDUklQVElPTiBPRiBBTkFMT0cgWUVBUiBQUk9DRVNTLS0+LiAgCioqKgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRSxmaWcud2lkdGg9Ny41LGZpZy5oZWlnaHQ9Ny41fQpwMwpgYGAKKioqCgojIEFuYWxvZyBZZWFyIEZvcmVjYXN0cyBFcnJvcnMKRm9yZWNhc3QgZXJyb3JzIGluIGFuYWxvZyB5ZWFycy4gSWYgb2JzZXJ2ZWQgZGF0YSBpcyBub3QgYXZhaWxhYmxlIGluIGEgZ2l2ZW4geWVhciB3ZSBjYW5ub3QgY2FsY3VsYXRlIGZvcmVjYXN0IGVycm9ycy4gVmFsdWVzIGFyZSBleHByZXNzZWQgYSBwZXJjZW50YWdlICBvZiBvYnNlcnZlZCB5aWVsZHMgaW4gYSBnaXZlbiB5ZWFyIF8odClfOgoKKioqCiQkXGZyYWN7KG9ic2VydmVkX3sodCl9LWZvcmVjYXN0X3sodCl9KX17b2JzZXJ2ZWRfeyh0KX19JCQKKioqCgoqKipQb3NpdGl2ZSAoKykqKiogdmFsdWVzIGluZGljYXRlIGFuIF91bmRlciBwcmVkaWN0aW9uXy4gKioqTmVnYXRpdmUgKC0pKioqIHZhbHVlcyBpbmRpY2F0ZSBhbiBfb3ZlciBwcmVkaWN0aW9uXy4KCioqKgpgYGB7cixlY2hvPUZBTFNFLHdhcm5pbmc9RkFMU0UsbWVzc2FnZT1GQUxTRSxmaWcud2lkdGg9Ny41LGZpZy5oZWlnaHQ9Ny41fQpwNApgYGA=