CAL 1976 1 1
ALLOCATE 30 2015:1
OPEN DATA
DATA(FORMAT=XLS,ORG=OBS,VERBOSE,missing=-9999) /

set LnCumMW = Log(CumMW)
Set LnPVCost   = log(PVCost)
set LnBatCost = log(BatCost)
set LnBatCap = log(BatCap)
set LnCarbonTx = log(carbontx)
set LnGDP = log(GDP)
set LnAusCoal = log(AusCoal)


Statistics(noprint) LnCumMW
set LnCumMWa = (LnCumMW-%mean)/(Sqrt(%variance))
Statistics(noprint) LnPVCost
set LnPVCosta = (LnPVCost-%mean)/(Sqrt(%variance))
Statistics(noprint) LnBatCost
set LnBatCosta = (LnBatCost-%mean)/(Sqrt(%variance))
Statistics(noprint) LnBatCap
set LnBatCapa = (LnBatCap-%mean)/(Sqrt(%variance))
Statistics(noprint) LnGDP
set LnGDPa = (LnGDP-%mean)/(Sqrt(%variance))
Statistics(noprint) Lncarbontx
set Lncarbontxa = (Lncarbontx-%mean)/(Sqrt(%variance))
Statistics(noprint) LnAusCoal
set LnAusCoala = (LnAusCoal-%mean)/(Sqrt(%variance))


linreg LnPVCost / resid
   # constant LnCumMW

linreg LnBatCost / resid
   # constant LnBatCap


source cats.src
@cats(lags=2, dettrend=cimean)
#  LnCumMWa LnPVCosta LnAusCoala LnGDPa
@cats(lags=2, exog, dettrend=cimean)
#  LnCumMWa LnPVCosta
#  LnAusCoala LnGDPa