* tspmed_cert.do -- certification script for tspmed
* Expected values were computed independently with the Python reference
* implementation (econmediate 0.2.0) on the same data and formulas.
* Run this file from the folder that contains it; the ado-files are in the same folder.
version 14
clear all
set more off
adopath ++ "`c(pwd)'"
cscript "tspmed certification" adofile tspmed

use tspmed_cert_stationary, clear
tsset t

* ---------------------------------------------------------------- ols, prewhitened QS
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ols) simreps(199)
matrix B0 = (1.3878950770539593, 1.053648388986642, 0.558887134432134, 0.5950211185820942, 0.778431616721342, 0.30422944633848203, -0.512539779496867, 0.19269580151029939)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.009633852513791813, -0.005477416372386623, 0.002943780903178635, 0.0008922515577381271, -0.001339871556241911, -5.199283046433934e-05, 0.0021323094424559475, 0.0006188205051674029 \ -0.005477416372386623, 0.007450929888041327, -0.00040192651087973434, -0.003249663793101239, 0.0026877271314626767, -0.00020273678977558956, -0.0004780718011869901, 9.973421199679388e-05 \ 0.002943780903178635, -0.00040192651087973434, 0.004925285111007727, -0.0007970320161254224, 0.00011400576730924559, -0.00019871641894567043, 0.0013370030166880943, 0.0016487467204928839 \ 0.0008922515577381271, -0.003249663793101239, -0.0007970320161254224, 0.008401187290231987, -0.001327524712943518, -0.0001433895324625969, 0.00071175157221724, 0.0001670123247132736 \ -0.001339871556241911, 0.0026877271314626767, 0.00011400576730924559, -0.001327524712943518, 0.007545977214366561, -0.004027103052007016, -0.001018160723774128, -0.00048504915746149814 \ -5.199283046433934e-05, -0.00020273678977558956, -0.00019871641894567043, -0.0001433895324625969, -0.004027103052007016, 0.006457422837698472, 5.8612279703112614e-05, 9.482651481168241e-05 \ 0.0021323094424559475, -0.0004780718011869901, 0.0013370030166880943, 0.00071175157221724, -0.001018160723774128, 5.8612279703112614e-05, 0.006288895769171583, -0.002722770410370338 \ 0.0006188205051674029, 9.973421199679388e-05, 0.0016487467204928839, 0.0001670123247132736, -0.00048504915746149814, 9.482651481168241e-05, -0.002722770410370338, 0.009871482532013608)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.3878950770539593, 0.09815219057288313 \ 0.5337413302492273, 0.13555784663915196 \ 1.9216364073031866, 0.12658798568462262 \ 0.2777535480805533, 0.058721474121738594 \ 0.8201932188947091, 0.13153402559378033 \ -0.2864518886454818, 0.04992244618457268)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 70.8846734723574) < 1e-05
assert reldif(G[2,1], 89.66928760230732) < 1e-05
assert reldif(G[3,1], 0.9570085283869924) < 1e-05
assert reldif(G[4,1], 1.4234434578566657) < 1e-05
assert reldif(G[5,1], 6.655315030108993) < 1e-05
assert reldif(G[6,1], 16.690743430831734) < 1e-05
assert reldif(G[7,1], 0.3610727590960134) < 1e-05
assert reldif(G[8,1], -0.09164420372123201) < 1e-05
assert reldif(G[9,1], 2.967314459306559) < 1e-05
assert reldif(G[10,1], 0.2596514113111813) < 1e-05
assert reldif(G[11,1], 0.5620036545150647) < 1e-05
assert reldif(G[12,1], 0.07977838264674408) < 1e-05
assert reldif(G[13,1], 0.5394340629547885) < 1e-05
assert reldif(G[1,5], 37.64712935143924) < 1e-05
assert reldif(G[2,5], 36.91020965857951) < 1e-05
assert reldif(G[3,5], 1.262431708780018) < 1e-05
assert reldif(G[4,5], 3.6975934391142165) < 1e-05
assert reldif(G[5,5], 6.14914566306593) < 1e-05
assert reldif(G[6,5], 2.4210891037253157) < 1e-05
assert reldif(G[7,5], 0.9659318747343986) < 1e-05
assert reldif(G[8,5], 0.14516898724779798) < 1e-05
assert reldif(G[9,5], 2.9081668875606574) < 1e-05
assert reldif(G[10,5], 0.5827598022741863) < 1e-05
assert reldif(G[11,5], 0.428862457624154) < 1e-05
assert reldif(G[12,5], 0.044612006567490536) < 1e-05
assert reldif(G[13,5], 1.0156762027421355) < 1e-05
assert reldif(G[1,9], 64.36463104460047) < 1e-05
assert reldif(G[2,9], 75.23169779983512) < 1e-05
assert reldif(G[3,9], 1.0005616365733232) < 1e-05
assert reldif(G[4,9], 4.350548560631523) < 1e-05
assert reldif(G[5,9], 4.89153966769193) < 1e-05
assert reldif(G[6,9], 39.169590558122515) < 1e-05
assert reldif(G[7,9], 1.2218641426680361) < 1e-05
assert reldif(G[8,9], -0.009126522726126815) < 1e-05
assert reldif(G[9,9], 3.3420026082179044) < 1e-05
assert reldif(G[10,9], 0.3484300995865954) < 1e-05
assert reldif(G[11,9], 0.8037095953276027) < 1e-05
assert reldif(G[12,9], 0.08872313195156317) < 1e-05
assert reldif(G[13,9], 0.4864644081532931) < 1e-05

* ---------------------------------------------------------------- ols, Newey-West
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ols) vce(hac) simreps(199)
matrix B0 = (1.3878950770539593, 1.053648388986642, 0.558887134432134, 0.5950211185820942, 0.778431616721342, 0.30422944633848203, -0.512539779496867, 0.19269580151029939)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.008272497818928679, -0.00458090911487464, 0.0021326498853280877, 0.000723844767152199, -0.0015019037263447857, 0.0004878227081863506, 0.0009508355866111502, 0.0007841931057601831 \ -0.00458090911487464, 0.006232598148876879, 0.00018122615943771224, -0.0028183286010813355, 0.0022619179676434255, -0.0009375428762684524, 1.025743538769643e-06, 0.000207119738140007 \ 0.0021326498853280877, 0.00018122615943771224, 0.004050534423982788, -0.0010132511125228569, 0.0003454962794084097, -0.0005319533963808167, 0.0005442026600320679, 0.0008803149172127621 \ 0.000723844767152199, -0.0028183286010813355, -0.0010132511125228569, 0.0074995768010841565, -0.0009286582101716041, 6.722833163766662e-05, -4.4978761950474885e-05, 0.0010023219472366572 \ -0.0015019037263447857, 0.0022619179676434255, 0.0003454962794084097, -0.0009286582101716041, 0.005655182959787607, -0.003060172414293284, -0.0007959310692805986, -0.0001239535910594669 \ 0.0004878227081863506, -0.0009375428762684524, -0.0005319533963808167, 6.722833163766662e-05, -0.003060172414293284, 0.004827437001816519, 0.0002177136774303535, 0.00014403985829650252 \ 0.0009508355866111502, 1.025743538769643e-06, 0.0005442026600320679, -4.4978761950474885e-05, -0.0007959310692805986, 0.0002177136774303535, 0.00518091242301074, -0.0023205427935620642 \ 0.0007841931057601831, 0.000207119738140007, 0.0008803149172127621, 0.0010023219472366572, -0.0001239535910594669, 0.00014403985829650252, -0.0023205427935620642, 0.008349239589570482)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.3878950770539593, 0.09095327272466285 \ 0.5337413302492273, 0.12116778299506421 \ 1.9216364073031866, 0.10739649595605923 \ 0.2777535480805533, 0.05354625253347455 \ 0.8201932188947091, 0.11732578904554004 \ -0.2864518886454818, 0.04868850582797447)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 70.8846734723574) < 1e-05
assert reldif(G[2,1], 89.66928760230732) < 1e-05
assert reldif(G[3,1], 0.9570085283869924) < 1e-05
assert reldif(G[4,1], 1.4234434578566657) < 1e-05
assert reldif(G[5,1], 6.655315030108993) < 1e-05
assert reldif(G[6,1], 16.690743430831734) < 1e-05
assert reldif(G[7,1], 0.3610727590960134) < 1e-05
assert reldif(G[8,1], -0.09164420372123201) < 1e-05
assert reldif(G[9,1], 2.967314459306559) < 1e-05
assert reldif(G[10,1], 0.2596514113111813) < 1e-05
assert reldif(G[11,1], 0.5620036545150647) < 1e-05
assert reldif(G[12,1], 0.07977838264674408) < 1e-05
assert reldif(G[13,1], 0.5394340629547885) < 1e-05
assert reldif(G[1,5], 37.64712935143924) < 1e-05
assert reldif(G[2,5], 36.91020965857951) < 1e-05
assert reldif(G[3,5], 1.262431708780018) < 1e-05
assert reldif(G[4,5], 3.6975934391142165) < 1e-05
assert reldif(G[5,5], 6.14914566306593) < 1e-05
assert reldif(G[6,5], 2.4210891037253157) < 1e-05
assert reldif(G[7,5], 0.9659318747343986) < 1e-05
assert reldif(G[8,5], 0.14516898724779798) < 1e-05
assert reldif(G[9,5], 2.9081668875606574) < 1e-05
assert reldif(G[10,5], 0.5827598022741863) < 1e-05
assert reldif(G[11,5], 0.428862457624154) < 1e-05
assert reldif(G[12,5], 0.044612006567490536) < 1e-05
assert reldif(G[13,5], 1.0156762027421355) < 1e-05
assert reldif(G[1,9], 64.36463104460047) < 1e-05
assert reldif(G[2,9], 75.23169779983512) < 1e-05
assert reldif(G[3,9], 1.0005616365733232) < 1e-05
assert reldif(G[4,9], 4.350548560631523) < 1e-05
assert reldif(G[5,9], 4.89153966769193) < 1e-05
assert reldif(G[6,9], 39.169590558122515) < 1e-05
assert reldif(G[7,9], 1.2218641426680361) < 1e-05
assert reldif(G[8,9], -0.009126522726126815) < 1e-05
assert reldif(G[9,9], 3.3420026082179044) < 1e-05
assert reldif(G[10,9], 0.3484300995865954) < 1e-05
assert reldif(G[11,9], 0.8037095953276027) < 1e-05
assert reldif(G[12,9], 0.08872313195156317) < 1e-05
assert reldif(G[13,9], 0.4864644081532931) < 1e-05

* ---------------------------------------------------------------- ols, White, trend
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ols) vce(robust) trend simreps(199)
matrix B0 = (1.3898369897199585, 1.0519952292579675, 0.5615836215770832, 0.6001246520850474, 0.7761915194639748, 0.32000590060988876, -0.5078291661342602, 0.15952010966875957)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.006233636494825091, -0.0034350779272985976, 0.0017871967906987524, -0.0008003304397689328, -0.0004871947122568801, 0.00029613594504511323, 0.0006232749751585766, 0.0002914869478482602 \ -0.0034350779272985976, 0.004462287197837741, -9.931246947036585e-05, -0.0012619400625619914, 0.000919527844895379, -0.0005342949206741565, -3.023359300261539e-05, -3.8482967884109084e-05 \ 0.0017871967906987524, -9.931246947036585e-05, 0.0033103364061253834, -0.0005796816131338511, -5.306518839131758e-05, 2.321650860615622e-05, 0.0005989271961699502, 0.00014727204621673123 \ -0.0008003304397689328, -0.0012619400625619914, -0.0005796816131338511, 0.005527660904033309, -0.00041615863295089116, 0.00048304467052671365, 0.00010096175024826747, -0.0002447787328135214 \ -0.0004871947122568801, 0.000919527844895379, -5.306518839131758e-05, -0.00041615863295089116, 0.003815358731918326, -0.002187659738406594, -0.0001569175238837689, 7.118698323073e-06 \ 0.00029613594504511323, -0.0005342949206741565, 2.321650860615622e-05, 0.00048304467052671365, -0.002187659738406594, 0.004841824656155156, 7.11869832307299e-06, 0.0003892613195473059 \ 0.0006232749751585766, -3.023359300261539e-05, 0.0005989271961699502, 0.00010096175024826747, -0.0001569175238837689, 7.11869832307299e-06, 0.0038965223683872087, -0.0024722660480880343 \ 0.0002914869478482602, -3.8482967884109084e-05, 0.00014727204621673123, -0.0002447787328135214, 7.118698323073e-06, 0.0003892613195473059, -0.0024722660480880343, 0.006197681302442752)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.3898369897199585, 0.07895338178196099 \ 0.5313612332064462, 0.10038281184479943 \ 1.9211982229264049, 0.09400849107999108 \ 0.2765780370112288, 0.04430609784969075 \ 0.8165497754665944, 0.09171981119218846 \ -0.2851885422601481, 0.04172481413574041)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 70.8749835979347) < 1e-05
assert reldif(G[2,1], 89.43624804559359) < 1e-05
assert reldif(G[3,1], 0.9572384069372966) < 1e-05
assert reldif(G[4,1], 4.664488568534075) < 1e-05
assert reldif(G[5,1], 18.56154784843353) < 1e-05
assert reldif(G[6,1], 16.90797050047112) < 1e-05
assert reldif(G[7,1], 0.43700834373098363) < 1e-05
assert reldif(G[8,1], -0.10189134726333642) < 1e-05
assert reldif(G[9,1], 2.979374423411134) < 1e-05
assert reldif(G[10,1], 0.2718698583409457) < 1e-05
assert reldif(G[11,1], 0.6279652715783814) < 1e-05
assert reldif(G[12,1], 0.08384796563096741) < 1e-05
assert reldif(G[13,1], 0.508639757663494) < 1e-05
assert reldif(G[1,5], 37.354175923464986) < 1e-05
assert reldif(G[2,5], 36.81220002134435) < 1e-05
assert reldif(G[3,5], 1.266961498083285) < 1e-05
assert reldif(G[4,5], 3.939210255386627) < 1e-05
assert reldif(G[5,5], 7.214543626356818) < 1e-05
assert reldif(G[6,5], 2.0247511710803847) < 1e-05
assert reldif(G[7,5], 1.1698315018743315) < 1e-05
assert reldif(G[8,5], 0.15654971004669913) < 1e-05
assert reldif(G[9,5], 2.8805321075890533) < 1e-05
assert reldif(G[10,5], 0.504925185932061) < 1e-05
assert reldif(G[11,5], 0.7188209803700368) < 1e-05
assert reldif(G[12,5], 0.048491067590474946) < 1e-05
assert reldif(G[13,5], 1.4866877323754433) < 1e-05
assert reldif(G[1,9], 61.597440983721995) < 1e-05
assert reldif(G[2,9], 69.7110426950164) < 1e-05
assert reldif(G[3,9], 1.0226744524694416) < 1e-05
assert reldif(G[4,9], 6.007757375154271) < 1e-05
assert reldif(G[5,9], 13.483556038787953) < 1e-05
assert reldif(G[6,9], 35.203758898768115) < 1e-05
assert reldif(G[7,9], 0.8277194498349768) < 1e-05
assert reldif(G[8,9], 0.041093068103718974) < 1e-05
assert reldif(G[9,9], 3.269641440063374) < 1e-05
assert reldif(G[10,9], 0.19705524627438262) < 1e-05
assert reldif(G[11,9], 0.8707431968762392) < 1e-05
assert reldif(G[12,9], 0.0901888323671397) < 1e-05
assert reldif(G[13,9], 0.7656001057768498) < 1e-05

* ---------------------------------------------------------------- ardl, AIC lags
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ardl) simreps(199)
matrix B0 = (1.2879531937652264, 1.0250670201483552, 0.33606920269491236, 0.49544580811097205, -0.4749358940083863, 0.7758067250034982, 0.22592783701154234, -0.7288902520818865, -0.39140824300767935, 0.20737531065367967, -0.5006676377902359)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.025981086309116552, -0.015089192001409592, 0.012125906752674064, 0.0012868900458161588, -0.0015714008687892792, -6.722648079395303e-06, -0.002772577718121093, -0.0013393495306074379, 0.005414715826097153, 0.002520011437378357, 0.0025387253199563786 \ -0.015089192001409592, 0.022045949457238936, -0.003635948613125179, -0.008165245015870713, 0.0006283739643143251, 0.0026574303608744056, 0.000840442950489057, 0.0014843916291574744, -0.0026625099427956085, -0.0029034959409011546, -0.0010714336283286275 \ 0.012125906752674064, -0.003635948613125179, 0.02011569195455084, -0.00012408484404609663, -0.0027838547299120064, -0.0005045974093144282, -0.00043998015616466803, -0.0014366306185250573, 0.0015255877022497294, 0.005368220373214688, 0.003993921314805924 \ 0.0012868900458161588, -0.008165245015870713, -0.00012408484404609663, 0.01885024988566632, -0.0013997048812938527, -0.002387819903392944, 0.00048156796312399784, -0.00047762124802059355, 0.0010751390012947303, 0.0017214421263382102, -0.0002590618817568624 \ -0.0015714008687892792, 0.0006283739643143251, -0.0027838547299120064, -0.0013997048812938527, 0.0027555916823735673, 0.000463819690441287, 2.6195127839492636e-05, 0.0007027495677778535, 0.00041448986002219016, -0.0005219836640540053, -0.0002208045792015098 \ -6.722648079395303e-06, 0.0026574303608744056, -0.0005045974093144282, -0.002387819903392944, 0.000463819690441287, 0.012949597793740808, -0.007291001884016568, 0.0017109794581140767, 4.275191970331604e-05, -0.0009457560886482344, 0.0005856232525989224 \ -0.002772577718121093, 0.000840442950489057, -0.00043998015616466803, 0.00048156796312399784, 2.6195127839492636e-05, -0.007291001884016568, 0.013426265956061539, 0.000538975369422625, -8.757878767763643e-05, -0.001828482053948056, -0.0005695951294186919 \ -0.0013393495306074379, 0.0014843916291574744, -0.0014366306185250573, -0.00047762124802059355, 0.0007027495677778535, 0.0017109794581140767, 0.000538975369422625, 0.004546143249495062, -0.00027806499386646267, 0.001625756967507485, -0.00016194030841597837 \ 0.005414715826097153, -0.0026625099427956085, 0.0015255877022497294, 0.0010751390012947303, 0.00041448986002219016, 4.275191970331604e-05, -8.757878767763643e-05, -0.00027806499386646267, 0.012864256504821881, -0.009665956376272385, -0.00031646437480463234 \ 0.002520011437378357, -0.0029034959409011546, 0.005368220373214688, 0.0017214421263382102, -0.0005219836640540053, -0.0009457560886482344, -0.001828482053948056, 0.001625756967507485, -0.009665956376272385, 0.0319319260320705, 0.0051191676342032385 \ 0.0025387253199563786, -0.0010714336283286275, 0.003993921314805924, -0.0002590618817568624, -0.0002208045792015098, 0.0005856232525989224, -0.0005695951294186919, -0.00016194030841597837, -0.00031646437480463234, 0.0051191676342032385, 0.00469274004841836)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.2879531937652264, 0.1611864954439136 \ 0.663713631654583, 0.19102358056811908 \ 1.9516668254198093, 0.18218616249804395 \ 0.34007527463700993, 0.08176820248852287 \ 0.7952538878103903, 0.17635913238384168 \ -0.1315402561558073, 0.06429075610584908)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 3.7734003605572597) < 1e-05
assert reldif(G[2,1], 2.8001265571843064) < 1e-05
assert reldif(G[3,1], 1.9950714831582643) < 1e-05
assert reldif(G[4,1], 5.158577749291415) < 1e-05
assert reldif(G[5,1], 55.16278395023806) < 1e-05
assert reldif(G[6,1], 0.058678040194921066) < 1e-05
assert reldif(G[7,1], 0.509092203479423) < 1e-05
assert reldif(G[8,1], -0.05384286415512015) < 1e-05
assert reldif(G[9,1], 2.8053972434949324) < 1e-05
assert reldif(G[10,1], 0.05561612908849032) < 1e-05
assert reldif(G[11,1], 0.9012259595488495) < 1e-05
assert reldif(G[12,1], 0.04502413967525365) < 1e-05
assert reldif(G[13,1], 0.705220547044973) < 1e-05
assert reldif(G[15,1], 15.555064972604356) < 1e-05
assert reldif(G[16,1], 15.555064972604356) < 1e-05
assert reldif(G[17,1], -8.527895028812312) < 1e-05
assert reldif(G[18,1], -8.527895028812312) < 1e-05
assert reldif(G[1,5], 1.7925600025306963) < 1e-05
assert reldif(G[2,5], 0.7047581545787641) < 1e-05
assert reldif(G[3,5], 1.9561269408439987) < 1e-05
assert reldif(G[4,5], 8.66335477957138) < 1e-05
assert reldif(G[5,5], 21.825474610372716) < 1e-05
assert reldif(G[6,5], 2.77289508921628) < 1e-05
assert reldif(G[7,5], 1.1295247886791735) < 1e-05
assert reldif(G[8,5], 0.15459720503253843) < 1e-05
assert reldif(G[9,5], 3.118954668909463) < 1e-05
assert reldif(G[10,5], 0.11000383217694583) < 1e-05
assert reldif(G[11,5], 0.29611957475956113) < 1e-05
assert reldif(G[12,5], 0.04510837536550649) < 1e-05
assert reldif(G[13,5], 1.8316008411359341) < 1e-05
assert reldif(G[15,5], 36.09403491362372) < 1e-05
assert reldif(G[16,5], 36.09403491362372) < 1e-05
assert reldif(G[17,5], -10.131262209147394) < 1e-05
assert reldif(G[18,5], -10.131262209147394) < 1e-05
assert reldif(G[1,9], 3.578804703371358) < 1e-05
assert reldif(G[2,9], 3.257352537557663) < 1e-05
assert reldif(G[3,9], 1.9459482765001577) < 1e-05
assert reldif(G[4,9], 3.2447675972669865) < 1e-05
assert reldif(G[5,9], 16.304965660728694) < 1e-05
assert reldif(G[6,9], 0.5671862417796697) < 1e-05
assert reldif(G[7,9], 0.07789983952096675) < 1e-05
assert reldif(G[8,9], -0.03815262842002954) < 1e-05
assert reldif(G[9,9], 2.9589399702185477) < 1e-05
assert reldif(G[10,9], 0.6221952683767443) < 1e-05
assert reldif(G[11,9], 0.6035063292284003) < 1e-05
assert reldif(G[12,9], 0.04056671707293663) < 1e-05
assert reldif(G[13,9], 1.1541405318747435) < 1e-05
assert reldif(G[15,9], 27.92620890121037) < 1e-05
assert reldif(G[16,9], 27.92620890121037) < 1e-05
assert reldif(G[17,9], -9.012149881678354) < 1e-05
assert reldif(G[18,9], -9.012149881678354) < 1e-05

* ---------------------------------------------------------------- ardl, fixed lags
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ardl) ardllags(2 1) simreps(199)
matrix B0 = (1.2854990754530156, 1.1059735985636536, 0.33198122001469976, 0.5352984121711878, -0.4576007346285045, 0.8138118869101406, 0.2461644908336899, -0.6768470520275972, -0.39319542650811706, 0.21115765768469433, -0.5074960122018759)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.029065253480003085, -0.01654487590429487, 0.012815877463079675, 0.0008435252266059631, -0.001381288427770928, 0.0005303915379658364, -0.0027629934218154177, -0.0012067015516693917, 0.005891104082651586, 0.0016044975048399184, 0.0025171502868530683 \ -0.01654487590429487, 0.02367398200581214, -0.0036719707939250794, -0.008997276758641768, 0.0006537194747686537, 0.0020503120856518713, 0.0011960522716760557, 0.0009257384555584635, -0.0025044535377720086, -0.0025658868710613093, -0.0005114085573773473 \ 0.012815877463079675, -0.0036719707939250794, 0.020621260224830115, -0.0005594627809434555, -0.002582927045034004, -0.0003067720719210327, -0.0007871719997772959, -0.0014726939452240227, 0.0019880543518856156, 0.004079461748064224, 0.0034745738073149917 \ 0.0008435252266059631, -0.008997276758641768, -0.0005594627809434555, 0.020881740977275715, -0.001056734636260515, -0.002049712568502935, -0.00010862134311573378, -8.33368763115048e-05, 0.0012800615082958701, 0.0020880279278140106, -0.0006743196000748975 \ -0.001381288427770928, 0.0006537194747686537, -0.002582927045034004, -0.001056734636260515, 0.002610208221118058, 0.00036511147135023455, 1.640965222143641e-05, 0.0003838386163003435, 0.0003794150190582925, -0.0002889674425109932, 0.00011193063106185831 \ 0.0005303915379658364, 0.0020503120856518713, -0.0003067720719210327, -0.002049712568502935, 0.00036511147135023455, 0.013943240370835617, -0.007925777712800081, 0.0009845957558654935, -0.0001126777353746347, -0.002036004303389154, 0.0006517578271159436 \ -0.0027629934218154177, 0.0011960522716760557, -0.0007871719997772959, -0.00010862134311573378, 1.640965222143641e-05, -0.007925777712800081, 0.01494851529809184, 0.0006275788406722299, 1.3379751977783294e-05, -0.0023296536135593185, -0.0005186349694970768 \ -0.0012067015516693917, 0.0009257384555584635, -0.0014726939452240227, -8.33368763115048e-05, 0.0003838386163003435, 0.0009845957558654935, 0.0006275788406722299, 0.003621848577825141, -0.00035985783667250796, 0.001021876592890064, -0.0001712394078050629 \ 0.005891104082651586, -0.0025044535377720086, 0.0019880543518856156, 0.0012800615082958701, 0.0003794150190582925, -0.0001126777353746347, 1.3379751977783294e-05, -0.00035985783667250796, 0.012475744877619186, -0.008967379234364421, -7.934269435537733e-05 \ 0.0016044975048399184, -0.0025658868710613093, 0.004079461748064224, 0.0020880279278140106, -0.0002889674425109932, -0.002036004303389154, -0.0023296536135593185, 0.001021876592890064, -0.008967379234364421, 0.029658949574679107, 0.003615047127771096 \ 0.0025171502868530683, -0.0005114085573773473, 0.0034745738073149917, -0.0006743196000748975, 0.00011193063106185831, 0.0006517578271159436, -0.0005186349694970768, -0.0001712394078050629, -7.934269435537733e-05, 0.003615047127771096, 0.005691597181085434)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.2854990754530156, 0.17048534681327285 \ 0.7695209637235204, 0.2040826125577915 \ 2.055020039176536, 0.196958152856008 \ 0.374459104560301, 0.0812084622022938 \ 0.9000544611198853, 0.1908530347817185 \ -0.1305334973963649, 0.06359289978065456)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 3.7862719411561248) < 1e-05
assert reldif(G[2,1], 3.3213551731820177) < 1e-05
assert reldif(G[3,1], 1.9952000881479204) < 1e-05
assert reldif(G[4,1], 4.959288427863079) < 1e-05
assert reldif(G[5,1], 42.29722192868966) < 1e-05
assert reldif(G[6,1], 0.19234880507910224) < 1e-05
assert reldif(G[7,1], 0.7814267535631316) < 1e-05
assert reldif(G[8,1], -0.04558736787224777) < 1e-05
assert reldif(G[9,1], 2.7405600041336577) < 1e-05
assert reldif(G[10,1], 0.16868409029191186) < 1e-05
assert reldif(G[11,1], 0.7092179351385055) < 1e-05
assert reldif(G[12,1], 0.043447143893205475) < 1e-05
assert reldif(G[13,1], 0.7155041248156309) < 1e-05
assert reldif(G[15,1], 14.5580256425416) < 1e-05
assert reldif(G[16,1], 14.5580256425416) < 1e-05
assert reldif(G[17,1], -8.22724145665211) < 1e-05
assert reldif(G[18,1], -8.22724145665211) < 1e-05
assert reldif(G[1,5], 4.5879765059641135) < 1e-05
assert reldif(G[2,5], 4.00022369019463) < 1e-05
assert reldif(G[3,5], 1.9940148278305105) < 1e-05
assert reldif(G[4,5], 8.946876541336689) < 1e-05
assert reldif(G[5,5], 18.283366410747863) < 1e-05
assert reldif(G[6,5], 3.52254338820382) < 1e-05
assert reldif(G[7,5], 1.5468040937383718) < 1e-05
assert reldif(G[8,5], 0.1730264650921167) < 1e-05
assert reldif(G[9,5], 3.173026264653031) < 1e-05
assert reldif(G[10,5], 0.0745434737910146) < 1e-05
assert reldif(G[11,5], 0.3277178726273779) < 1e-05
assert reldif(G[12,5], 0.04419599815047226) < 1e-05
assert reldif(G[13,5], 2.0134238135155034) < 1e-05
assert reldif(G[15,5], 36.061460271904004) < 1e-05
assert reldif(G[16,5], 36.061460271904004) < 1e-05
assert reldif(G[17,5], -10.231731877748825) < 1e-05
assert reldif(G[18,5], -10.231731877748825) < 1e-05
assert reldif(G[1,9], 5.22823335496496) < 1e-05
assert reldif(G[2,9], 2.981761069322219) < 1e-05
assert reldif(G[3,9], 1.9401691274535602) < 1e-05
assert reldif(G[4,9], 3.4578319721553132) < 1e-05
assert reldif(G[5,9], 18.012153510837784) < 1e-05
assert reldif(G[6,9], 0.47218444402001847) < 1e-05
assert reldif(G[7,9], 0.1115147759291213) < 1e-05
assert reldif(G[8,9], -0.045186465782407224) < 1e-05
assert reldif(G[9,9], 2.948770264985546) < 1e-05
assert reldif(G[10,9], 0.6858892677776337) < 1e-05
assert reldif(G[11,9], 0.7382449158323631) < 1e-05
assert reldif(G[12,9], 0.04440871307865041) < 1e-05
assert reldif(G[13,9], 1.139106041197817) < 1e-05
assert reldif(G[15,9], 21.47031437738245) < 1e-05
assert reldif(G[16,9], 21.47031437738245) < 1e-05
assert reldif(G[17,9], -7.9138655994588385) < 1e-05
assert reldif(G[18,9], -7.9138655994588385) < 1e-05

* ---------------------------------------------------------------- dml, blocked cross-fitting
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(dml) simreps(199)
matrix B0 = (1.4003813235899187, 1.041287760782942, 0.5659391086585296, 0.771762727296859, -0.4993046254191984)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.009636754765625477, -0.005552584131608523, 0.002838512478875303, -0.0015184713596711511, 0.0021464132745194113 \ -0.005552584131608523, 0.007691403358326579, -0.00043249160533678516, 0.002669779957817504, -0.00064569023467824 \ 0.002838512478875303, -0.00043249160533678516, 0.004688672546114945, -0.0001384817624901642, 0.0013852288496505392 \ -0.0015184713596711511, 0.002669779957817504, -0.0001384817624901642, 0.007625436131705088, -0.0013063271881473585 \ 0.0021464132745194113, -0.00064569023467824, 0.0013852288496505392, -0.0013063271881473585, 0.006352315327633602)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.4003813235899187, 0.09816697390330219 \ 0.5210510675038604, 0.13391974039888743 \ 1.9214323910937792, 0.12423074907547527 \ 0.2711784551561823, 0.05849877486145741 \ 0.8036270821626826, 0.13092081732438435 \ -0.28257601465882215, 0.04919969124230578)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 70.789090726391) < 1e-05
assert reldif(G[2,1], 89.78566883930161) < 1e-05
assert reldif(G[3,1], 0.9511512067815807) < 1e-05
assert reldif(G[4,1], 1.5081950548198886) < 1e-05
assert reldif(G[5,1], 5.870800219034411) < 1e-05
assert reldif(G[6,1], 15.95709878827789) < 1e-05
assert reldif(G[7,1], 0.6386602422234027) < 1e-05
assert reldif(G[8,1], -0.09908062250269171) < 1e-05
assert reldif(G[9,1], 2.8515294499526487) < 1e-05
assert reldif(G[10,1], 2.6008653899407697) < 1e-05
assert reldif(G[11,1], 0.6051716494745218) < 1e-05
assert reldif(G[12,1], 0.07293587868028761) < 1e-05
assert reldif(G[13,1], 0.5474488578497423) < 1e-05
assert reldif(G[1,5], 39.97292050903345) < 1e-05
assert reldif(G[2,5], 38.06451976580917) < 1e-05
assert reldif(G[3,5], 1.230580883063193) < 1e-05
assert reldif(G[4,5], 0.8693282254618184) < 1e-05
assert reldif(G[5,5], 1.7139079574323424) < 1e-05
assert reldif(G[6,5], 2.240910674647541) < 1e-05
assert reldif(G[7,5], 1.391971644797827) < 1e-05
assert reldif(G[8,5], 0.18009295231733524) < 1e-05
assert reldif(G[9,5], 2.9375869244220776) < 1e-05
assert reldif(G[10,5], 1.322418379249308) < 1e-05
assert reldif(G[11,5], 0.7184187361187678) < 1e-05
assert reldif(G[12,5], 0.045060068022520716) < 1e-05
assert reldif(G[13,5], 0.3644005675441257) < 1e-05
assert reldif(G[1,9], 67.36809670498) < 1e-05
assert reldif(G[2,9], 82.78474686623159) < 1e-05
assert reldif(G[3,9], 0.9719155560899738) < 1e-05
assert reldif(G[4,9], 3.785781626612361) < 1e-05
assert reldif(G[5,9], 3.830945432598182) < 1e-05
assert reldif(G[6,9], 41.908388032736994) < 1e-05
assert reldif(G[7,9], 1.1574173114039394) < 1e-05
assert reldif(G[8,9], -0.020274879895032294) < 1e-05
assert reldif(G[9,9], 3.33085915262502) < 1e-05
assert reldif(G[10,9], 0.5340326540473579) < 1e-05
assert reldif(G[11,9], 0.8647642413662024) < 1e-05
assert reldif(G[12,9], 0.09774857213153776) < 1e-05
assert reldif(G[13,9], 0.40549308284611574) < 1e-05

use tspmed_cert_coint, clear
tsset t

* ---------------------------------------------------------------- fmols
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(fmols) simreps(199)
matrix B0 = (1.3746346085899688, 0.9948323642790176, 0.39792046357191535, 0.5386609917025, 0.8117259105812914, 0.29414843914431416, -0.4901467021411046, 0.2915842715560044)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.008326946780593819, -0.006574253423015292, 0.005590169602287607, 0.00021550034576224257, 4.796353428721227e-05, -1.4097526286590162e-05, 4.674322512333849e-05, -1.3738850872629598e-05 \ -0.006574253423015438, 0.007500454956938377, -0.001383677513909903, -0.0018518480924990693, 3.870370050691924e-19, -2.447677319634029e-19, 3.7718984073772395e-19, -2.3854024454458736e-19 \ 0.005590169602287695, -0.0013836775139099155, 0.009327617112961723, -0.002300226670469198, -3.135099857484578e-19, 1.7837027311137435e-19, -3.0553352791939443e-19, 1.7383209880718188e-19 \ 0.00021550034576229624, -0.0018518480924991055, -0.0023002266704691917, 0.001475677430277659, -1.4097526286590617e-05, 2.6411306420519714e-05, -1.3738850872630042e-05, 2.573933844039045e-05 \ 4.796353428720795e-05, 6.108049796855515e-18, -5.52705763328556e-19, -1.4097526286592198e-05, 0.0008161742910432642, -0.00023989138193031894, 0.00030213189563514013, -8.880313771767491e-05 \ -1.4097526286588896e-05, -1.938020821373602e-18, 2.0895442752638548e-19, 2.641130642052018e-05, -0.00023989138193031902, 0.0004494295429567796, -8.880313771767492e-05, 0.0001663701016536535 \ 4.6743225123334284e-05, 5.9526461292303084e-18, -5.386435821749142e-19, -1.3738850872631585e-05, 0.00030213189563514013, -8.880313771767491e-05, 0.0005236826171398618, -0.00015392171512531542 \ -1.3738850872628367e-05, -1.8887128501565825e-18, 2.0363811782294415e-19, 2.573933844039091e-05, -8.880313771767492e-05, 0.0001663701016536535, -0.00015392171512531544, 0.00028836786683728347)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.3746346085899688, 0.09125210562255416 \ 0.6124918038358906, 0.09702507125494422 \ 1.9871264124258594, 0.041478339811682155 \ 0.30822991431539987, 0.04697139732823676 \ 0.8075312067701245, 0.07582745300760112 \ -0.19503940293423386, 0.048206040627205377)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 37.5232820218507) < 1e-05
assert reldif(G[2,1], 41.75403536450219) < 1e-05
assert reldif(G[3,1], 1.2330761316193242) < 1e-05
assert reldif(G[4,1], 12.749625057370253) < 1e-05
assert reldif(G[5,1], 31.9980994779772) < 1e-05
assert reldif(G[6,1], 4.075697817375617) < 1e-05
assert reldif(G[7,1], 0.5540831205289457) < 1e-05
assert reldif(G[8,1], 0.05690244554959566) < 1e-05
assert reldif(G[9,1], 2.798868081248662) < 1e-05
assert reldif(G[10,1], 1.0171705179744213) < 1e-05
assert reldif(G[11,1], 0.7926039336998487) < 1e-05
assert reldif(G[12,1], 0.07362752829893054) < 1e-05
assert reldif(G[13,1], 0.8239904683060887) < 1e-05
assert reldif(G[20,1], -10.500363839319661) < 1e-05
assert reldif(G[1,5], 72.26288511213954) < 1e-05
assert reldif(G[2,5], 95.50577411783689) < 1e-05
assert reldif(G[3,5], 0.9326046239457821) < 1e-05
assert reldif(G[4,5], 0.11902497715600235) < 1e-05
assert reldif(G[5,5], 9.489170250093364) < 1e-05
assert reldif(G[6,5], 8.302222994146113) < 1e-05
assert reldif(G[7,5], 2.2821891552701437) < 1e-05
assert reldif(G[8,5], 0.17223039363812706) < 1e-05
assert reldif(G[9,5], 3.3183030082974487) < 1e-05
assert reldif(G[10,5], 3.199038010544043) < 1e-05
assert reldif(G[11,5], 1.6361513271884185) < 1e-05
assert reldif(G[12,5], 0.06938942176957641) < 1e-05
assert reldif(G[13,5], 1.0962817834064884) < 1e-05
assert reldif(G[20,5], -8.696537883060257) < 1e-05
assert reldif(G[1,9], 64.72818611847167) < 1e-05
assert reldif(G[2,9], 109.04669586952798) < 1e-05
assert reldif(G[3,9], 1.0041830884847927) < 1e-05
assert reldif(G[4,9], 4.9652299525044885) < 1e-05
assert reldif(G[5,9], 6.806764554977233) < 1e-05
assert reldif(G[6,9], 15.37705409740983) < 1e-05
assert reldif(G[7,9], 1.105411920730461) < 1e-05
assert reldif(G[8,9], -0.13198430000097805) < 1e-05
assert reldif(G[9,9], 2.807993996155514) < 1e-05
assert reldif(G[10,9], 3.1234336930751265) < 1e-05
assert reldif(G[11,9], 0.718910665586823) < 1e-05
assert reldif(G[12,9], 0.08741535619735266) < 1e-05
assert reldif(G[13,9], 4.875543837814318) < 1e-05
assert reldif(G[20,9], -6.913339070191363) < 1e-05

* ---------------------------------------------------------------- fmols, no prewhitening
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(fmols) noprewhite simreps(199)
matrix B0 = (1.390936287275193, 0.9969411976394761, 0.4317173777567468, 0.5285236844548172, 0.7935888714186282, 0.3001214971279613, -0.49116808423808067, 0.29150499892698617)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.005641650571639874, -0.004454170485212906, 0.0037874366635570916, 0.0001460052142631138, 1.7512104196899953e-05, -5.147188440513313e-06, 2.7242824120962954e-05, -8.007258740910224e-06 \ -0.004454170485213005, 0.005081688055697278, -0.0009374654651939973, -0.0012546591355652383, 1.4131219605797178e-19, -8.936785184583035e-20, 2.1983327988853865e-19, -1.390257071640321e-19 \ 0.003787436663557152, -0.0009374654651940058, 0.006319622042021225, -0.0015584433829452674, -1.1446653418656386e-19, 6.512528433895443e-20, -1.780706432258883e-19, 1.0131259197212362e-19 \ 0.00014600521426315016, -0.0012546591355652628, -0.001558443382945263, 0.0009997970009228661, -5.147188440513479e-06, 9.643108183871e-06, -8.007258740910484e-06, 1.500136689908007e-05 \ 1.7512104196898374e-05, 2.2301276599398898e-18, -2.0179999371350253e-19, -5.1471884405140575e-06, 0.0007148271840991548, -0.00021010326215459635, 0.0002563177578964837, -7.533736584188608e-05 \ -5.147188440512852e-06, -7.075963659480511e-19, 7.629195307697004e-20, 9.64310818387117e-06, -0.0002101032621545964, 0.0003936223649388818, -7.53373658418861e-05, 0.00014114236878972795 \ 2.72428241209605e-05, 3.4693132774867872e-18, -3.139315341283534e-19, -8.007258740911381e-06, 0.0002563177578964837, -7.533736584188608e-05, 0.0005982236126727896, -0.00017583093552724533 \ -8.007258740909508e-06, -1.1007771041910405e-18, 1.1868409622006359e-19, 1.5001366899080337e-05, -7.53373658418861e-05, 0.00014114236878972795, -0.00017583093552724538, 0.0003294141555057062)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.390936287275193, 0.07511092178582421 \ 0.5791156425403783, 0.08061241644746817 \ 1.9700519298155712, 0.0375280622700725 \ 0.293959582372325, 0.03916831301259492 \ 0.7911614399054474, 0.06253657409222299 \ -0.2120457973650691, 0.040448497405212475)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 37.67009236082639) < 1e-05
assert reldif(G[2,1], 42.3009917551009) < 1e-05
assert reldif(G[3,1], 1.2287894976371234) < 1e-05
assert reldif(G[4,1], 13.018493328877275) < 1e-05
assert reldif(G[5,1], 31.75739221309582) < 1e-05
assert reldif(G[6,1], 4.27222076397106) < 1e-05
assert reldif(G[7,1], 0.42767980381675286) < 1e-05
assert reldif(G[8,1], 0.031832100932266516) < 1e-05
assert reldif(G[9,1], 2.8072073155785673) < 1e-05
assert reldif(G[10,1], 1.0171705179744213) < 1e-05
assert reldif(G[11,1], 0.7926039336998487) < 1e-05
assert reldif(G[12,1], 0.07362752829893054) < 1e-05
assert reldif(G[13,1], 0.8239904683060887) < 1e-05
assert reldif(G[20,1], -10.477529932743828) < 1e-05
assert reldif(G[1,5], 70.88331362780292) < 1e-05
assert reldif(G[2,5], 93.39827876019253) < 1e-05
assert reldif(G[3,5], 0.9375747254041547) < 1e-05
assert reldif(G[4,5], 0.06836227863745026) < 1e-05
assert reldif(G[5,5], 10.58990825094895) < 1e-05
assert reldif(G[6,5], 7.820722155334909) < 1e-05
assert reldif(G[7,5], 2.3152882425267065) < 1e-05
assert reldif(G[8,5], 0.16836465097622533) < 1e-05
assert reldif(G[9,5], 3.331321127819084) < 1e-05
assert reldif(G[10,5], 3.199038010544043) < 1e-05
assert reldif(G[11,5], 1.6361513271884185) < 1e-05
assert reldif(G[12,5], 0.06938942176957641) < 1e-05
assert reldif(G[13,5], 1.0962817834064884) < 1e-05
assert reldif(G[20,5], -8.694415139149958) < 1e-05
assert reldif(G[1,9], 64.83785279155907) < 1e-05
assert reldif(G[2,9], 109.25553105208206) < 1e-05
assert reldif(G[3,9], 1.0036253577593464) < 1e-05
assert reldif(G[4,9], 4.922126702679187) < 1e-05
assert reldif(G[5,9], 6.752656965372833) < 1e-05
assert reldif(G[6,9], 15.408166783707966) < 1e-05
assert reldif(G[7,9], 1.1192914709573363) < 1e-05
assert reldif(G[8,9], -0.13412732313702988) < 1e-05
assert reldif(G[9,9], 2.8104664672798627) < 1e-05
assert reldif(G[10,9], 3.1234336930751265) < 1e-05
assert reldif(G[11,9], 0.718910665586823) < 1e-05
assert reldif(G[12,9], 0.08741535619735266) < 1e-05
assert reldif(G[13,9], 4.875543837814318) < 1e-05
assert reldif(G[20,9], -6.906062104899487) < 1e-05

* ---------------------------------------------------------------- fmols, fixed bandwidth
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(fmols) fmbw(6) simreps(199)
matrix B0 = (1.391011785722009, 1.001533562935819, 0.437532382839839, 0.5256792303827444, 0.7948915194596091, 0.29874158119039773, -0.49067701999309216, 0.29253616369669316)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.006704712411293185, -0.005293474268744392, 0.004501107128610561, 0.00017351712229473896, 4.72090428648856e-05, -1.387576484183164e-05, 4.8819409380625395e-05, -1.4349086598121378e-05 \ -0.005293474268744509, 0.006039235600415836, -0.0011141130170736412, -0.0014910758068664744, 3.8094871101856887e-19, -2.409174078166692e-19, 3.939434046452605e-19, -2.491354377334622e-19 \ 0.004501107128610632, -0.0011141130170736514, 0.007510434721501286, -0.0018521024227301014, -3.085783100790934e-19, 1.7556441563943072e-19, -3.191043480556981e-19, 1.815531635390721e-19 \ 0.00017351712229478217, -0.0014910758068665036, -0.0018521024227300962, 0.0011881897461991807, -1.3875764841832087e-05, 2.5995842788765694e-05, -1.4349086598121843e-05, 2.6882597364478607e-05 \ 4.7209042864881364e-05, 6.0119669862920696e-18, -5.440114132624251e-19, -1.387576484183365e-05, 0.0008155866809495812, -0.00023971867053894488, 0.000315456225221768, -9.271944808533601e-05 \ -1.3875764841830401e-05, -1.9075347425692157e-18, 2.0566746534628887e-19, 2.5995842788766172e-05, -0.00023971867053894496, 0.00044910597317672295, -9.271944808533604e-05, 0.00017370719548525306 \ 4.8819409380621004e-05, 6.217043593249905e-18, -5.625684038501267e-19, -1.4349086598123454e-05, 0.000315456225221768, -9.2719448085336e-05, 0.0006416487709138136, -0.00018859453435080817 \ -1.4349086598120093e-05, -1.9726034220134374e-18, 2.1268307039718228e-19, 2.6882597364479084e-05, -9.271944808533602e-05, 0.00017370719548525303, -0.00018859453435080822, 0.00035332638753171434)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.391011785722009, 0.08188230829442747 \ 0.5814234498695203, 0.08781017080977675 \ 1.9724352355915293, 0.041472475940511594 \ 0.29477441863643883, 0.04258863246986955 \ 0.7961105356318493, 0.0680734631553433 \ -0.21468708576232892, 0.043944005348269186)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 37.83109938085457) < 1e-05
assert reldif(G[2,1], 42.46781421133927) < 1e-05
assert reldif(G[3,1], 1.2280500185544125) < 1e-05
assert reldif(G[4,1], 13.06134222398524) < 1e-05
assert reldif(G[5,1], 31.57633104824339) < 1e-05
assert reldif(G[6,1], 4.292380447583924) < 1e-05
assert reldif(G[7,1], 0.4242286657470786) < 1e-05
assert reldif(G[8,1], 0.02333835586904282) < 1e-05
assert reldif(G[9,1], 2.803249410336728) < 1e-05
assert reldif(G[10,1], 1.0171705179744213) < 1e-05
assert reldif(G[11,1], 0.7926039336998487) < 1e-05
assert reldif(G[12,1], 0.07362752829893054) < 1e-05
assert reldif(G[13,1], 0.8239904683060887) < 1e-05
assert reldif(G[20,1], -10.471582807737637) < 1e-05
assert reldif(G[1,5], 70.90025182485942) < 1e-05
assert reldif(G[2,5], 93.45269681433845) < 1e-05
assert reldif(G[3,5], 0.9376449967914604) < 1e-05
assert reldif(G[4,5], 0.06868160288405012) < 1e-05
assert reldif(G[5,5], 10.343545104562086) < 1e-05
assert reldif(G[6,5], 7.812585632472714) < 1e-05
assert reldif(G[7,5], 2.352007829628086) < 1e-05
assert reldif(G[8,5], 0.16886614985267076) < 1e-05
assert reldif(G[9,5], 3.33561353306953) < 1e-05
assert reldif(G[10,5], 3.199038010544043) < 1e-05
assert reldif(G[11,5], 1.6361513271884185) < 1e-05
assert reldif(G[12,5], 0.06938942176957641) < 1e-05
assert reldif(G[13,5], 1.0962817834064884) < 1e-05
assert reldif(G[20,5], -8.695972221409844) < 1e-05
assert reldif(G[1,9], 64.76530551436832) < 1e-05
assert reldif(G[2,9], 109.05779651598085) < 1e-05
assert reldif(G[3,9], 1.0040722667597848) < 1e-05
assert reldif(G[4,9], 4.920897863275962) < 1e-05
assert reldif(G[5,9], 6.725961959493234) < 1e-05
assert reldif(G[6,9], 15.329259287453198) < 1e-05
assert reldif(G[7,9], 1.116296725636902) < 1e-05
assert reldif(G[8,9], -0.13391211932515187) < 1e-05
assert reldif(G[9,9], 2.8106193190607494) < 1e-05
assert reldif(G[10,9], 3.1234336930751265) < 1e-05
assert reldif(G[11,9], 0.718910665586823) < 1e-05
assert reldif(G[12,9], 0.08741535619735266) < 1e-05
assert reldif(G[13,9], 4.875543837814318) < 1e-05
assert reldif(G[20,9], -6.912643007447881) < 1e-05

* ---------------------------------------------------------------- dols
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(dols) leads(2) dlags(2) simreps(199)
matrix B0 = (1.451561845936293, 0.9244171331015689, 0.4494374898861969, 0.5459382604963713, 0.7934647195620541, 0.3022355949567148, -0.4854476985800698, 0.290642173518142)
matrix B1 = e(b)
assert mreldif(B1, B0) < 1e-06
matrix V0 = (0.021582423294596496, -0.018472966092061484, 0.013332613131573818, 0.0011552490735992718, -0.0005582164335664181, 0.0013141787779733095, -0.000331753713404216, 0.0001547641519420869 \ -0.018472966092061484, 0.020756388471934027, -0.004563487156259019, -0.004576597089777619, 0.00029470160534945966, -0.0015072932697836573, -8.812099881747281e-05, 8.989960983287837e-05 \ 0.013332613131573818, -0.004563487156259019, 0.02049625626698131, -0.0046829669708823365, -0.0013256455353373177, 0.00035023432839530373, -0.0015917938940022349, 0.00047542505827002056 \ 0.0011552490735992718, -0.004576597089777619, -0.0046829669708823365, 0.0029697124527805433, 0.0002577413408918154, 0.00027768112744800446, 0.00046265853712033867, -0.00023999060819034868 \ -0.0005582164335664181, 0.00029470160534945966, -0.0013256455353373177, 0.0002577413408918154, 0.0010309744032074542, -0.00020717226960279865, 0.0006404957489427492, -5.898963133686027e-05 \ 0.0013141787779733095, -0.0015072932697836573, 0.00035023432839530373, 0.00027768112744800446, -0.00020717226960279865, 0.0006896423593483339, -9.511754472608463e-07, 3.3219337585858046e-05 \ -0.000331753713404216, -8.812099881747281e-05, -0.0015917938940022349, 0.00046265853712033867, 0.0006404957489427492, -9.511754472608463e-07, 0.0009916200157070731, -0.00016541575806474816 \ 0.0001547641519420869, 8.989960983287837e-05, 0.00047542505827002056, -0.00023999060819034868, -5.898963133686027e-05, 3.3219337585858046e-05, -0.00016541575806474816, 0.00020042291736596725)
matrix V1 = e(V)
assert mreldif(V1, V0) < 1e-06
matrix E0 = (1.451561845936293, 0.1469095752207231 \ 0.5153139861539368, 0.1590000853426109 \ 1.9668758320902295, 0.05721301607262699 \ 0.26199619607217633, 0.07823581422057131 \ 0.7334923812747944, 0.11992195753945485 \ -0.21817839512085768, 0.07566391658732925)
matrix E1 = e(effects)
matrix E1 = E1[1..., 1..2]
assert mreldif(E1, E0) < 1e-06
matrix G = e(diag)
assert reldif(G[1,1], 58.02133609728411) < 1e-05
assert reldif(G[2,1], 57.80908960785067) < 1e-05
assert reldif(G[3,1], 1.0699661179813544) < 1e-05
assert reldif(G[4,1], 31.013514521431805) < 1e-05
assert missing(G[5,1])
assert reldif(G[6,1], 22.613942627183185) < 1e-05
assert reldif(G[7,1], 0.06225066260069513) < 1e-05
assert reldif(G[8,1], -0.014425629942947759) < 1e-05
assert reldif(G[9,1], 3.0725646532031483) < 1e-05
assert reldif(G[10,1], 2.834070138118511) < 1e-05
assert reldif(G[11,1], 0.8950503957105963) < 1e-05
assert reldif(G[12,1], 0.0639215823084695) < 1e-05
assert reldif(G[13,1], 0.9371059738910824) < 1e-05
assert reldif(G[20,1], -9.052552883957944) < 1e-05
assert reldif(G[1,5], 71.4546258306175) < 1e-05
assert reldif(G[2,5], 95.06703223137161) < 1e-05
assert reldif(G[3,5], 0.9289638624223177) < 1e-05
assert reldif(G[4,5], 7.639085767739157) < 1e-05
assert reldif(G[5,5], 87.30896317249284) < 1e-05
assert reldif(G[6,5], 6.57990862242578) < 1e-05
assert reldif(G[7,5], 2.103973587456984) < 1e-05
assert reldif(G[8,5], 0.16028469562906394) < 1e-05
assert reldif(G[9,5], 3.321463529014133) < 1e-05
assert reldif(G[10,5], 4.1683687516348575) < 1e-05
assert reldif(G[11,5], 2.1129075336931087) < 1e-05
assert reldif(G[12,5], 0.06034395464727493) < 1e-05
assert reldif(G[13,5], 0.9664140086975266) < 1e-05
assert reldif(G[20,5], -8.596199238680668) < 1e-05
assert reldif(G[1,9], 61.791107140693505) < 1e-05
assert reldif(G[2,9], 98.02196318388967) < 1e-05
assert reldif(G[3,9], 1.0204637440971138) < 1e-05
assert reldif(G[4,9], 19.44898003640116) < 1e-05
assert reldif(G[5,9], 88.04194985324081) < 1e-05
assert reldif(G[6,9], 18.249572867975218) < 1e-05
assert reldif(G[7,9], 0.886798007082826) < 1e-05
assert reldif(G[8,9], -0.1296634116521152) < 1e-05
assert reldif(G[9,9], 2.8599299904618065) < 1e-05
assert reldif(G[10,9], 3.196698249311144) < 1e-05
assert reldif(G[11,9], 0.47529202478593374) < 1e-05
assert reldif(G[12,9], 0.06671470165902282) < 1e-05
assert reldif(G[13,9], 2.268271963935086) < 1e-05
assert reldif(G[20,9], -6.911232966522234) < 1e-05

* ---------------------------------------------------------------- simulated p-values are valid probabilities
use tspmed_cert_stationary, clear
tsset t
set seed 3
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ardl) simreps(199)
matrix G = e(diag)
forvalues r = 15/18 {
    assert G[`r',4] >= 0 & G[`r',4] <= 1
}
assert inlist(G[19,1], 0, 0.5, 1)
assert G[11,4] >= 0 & G[11,4] <= 1 & G[12,4] >= 0 & G[12,4] <= 1

* ---------------------------------------------------------------- bootstraps run
use tspmed_cert_coint, clear
tsset t
set seed 4
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(fmols) bootstrap(50) nodots nodiag
assert "`e(bmethod)'" == "sieve" & el(e(effects),3,5) < . & el(e(effects),3,6) < .
use tspmed_cert_stationary, clear
tsset t
tspmed y, treatment(d) mediators(m1 m2) controls(x) estimator(ols) bootstrap(50) nodots nodiag
assert "`e(bmethod)'" == "block" & el(e(effects),3,2) < .

* ---------------------------------------------------------------- error checks
rcof "tspmed y, treatment(d) mediators(m1) estimator(bogus)" == 198
rcof "tspmed y, treatment(d) mediators(m1) estimator(dml)" == 198
rcof "tspmed y, treatment(d) mediators(m1) estimator(ols) bmethod(sieve)" == 198
rcof "tspmed y, treatment(d) mediators(m1) vce(bogus)" == 198
rcof "tspmed y, treatment(d) mediators(m1) simreps(10)" == 198
di as res "tspmed certification completed"