From 13a30a0a4e95f1a71577441c20298cd87256334f Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@jrc.ec.europa.eu>
Date: Wed, 17 Jan 2018 19:33:39 +0100
Subject: [PATCH] add a new case for ratios between endogenous variables, like
 trade balance to GDP.

---
 matlab/utilities/dataset/quarterly2annual.m | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/matlab/utilities/dataset/quarterly2annual.m b/matlab/utilities/dataset/quarterly2annual.m
index 8663a09ff..13c6bafd9 100644
--- a/matlab/utilities/dataset/quarterly2annual.m
+++ b/matlab/utilities/dataset/quarterly2annual.m
@@ -14,6 +14,7 @@ function [ya, yass, gya, gyass] = quarterly2annual(y,yss,GYTREND0,type,islog,aux
 %          5 annual price as quantity weighted average
 %          6 annual quantity from average price
 %          7 annual nominal from Q real and deflator
+%          8 annual ratio [e.g. trade balance to GDP]
 % islog    0 level (default)
 %          1 log-level
 %          2 growth rate Q frequency
@@ -120,6 +121,19 @@ switch type
     yn = (y+yss).*(yaux+yauxss) - yss.*yauxss;
     [ya, yass] = quarterly2annual(yn,yss.*yauxss,GYTREND0+GYTREND0aux,typeaux,0,0);
     GYTREND0=GYTREND0+GYTREND0aux;
+
+    case 8
+        % numerator
+        yn = y;
+        [yna, ynass] = quarterly2annual(yn,yss,GYTREND0,typeaux(1),0,0);
+        % denominator
+        yd = yaux;
+        [yda, ydass] = quarterly2annual(yd,yauxss,GYTREND0aux,typeaux(2),0,0);
+        % ratio
+        yass = ynass/ydass;
+        ya = (yna+ynass)./(yda+ydass)-yass;    
+        GYTREND0 = GYTREND0 - GYTREND0aux;
+
   otherwise
     error('Wrong type input')
 end