• No results found

Corrigendum to “Bayesian reduced rank regression in econometrics” [J. Econometrics 75 (1996) 121–146]

N/A
N/A
Protected

Academic year: 2021

Share "Corrigendum to “Bayesian reduced rank regression in econometrics” [J. Econometrics 75 (1996) 121–146]"

Copied!
4
0
0

Loading.... (view fulltext now)

Full text

(1)

http://www.diva-portal.org

Postprint

This is the accepted version of a paper published in Journal of Econometrics. This paper has been peer-reviewed but does not include the final publisher proof-corrections or journal pagination.

Citation for the original published paper (version of record): Karlsson, S. (2017)

Corrigendum to “Bayesian reduced rank regression in econometrics” [J. Econometrics 75 (1996) 121–146].

Journal of Econometrics, 201(1): 170-171

https://doi.org/10.1016/j.jeconom.2012.10.005

Access to the published version may require subscription. N.B. When citing this work, cite the original published paper.

Permanent link to this version:

(2)

Corrigendum to ”Bayesian reduced rank regression in

econometrics” [J. Econometrics, 75 (1996) 121-146]

Sune Karlsson

¨

Orebro University

October 10, 2012

Department of Statistics, ¨Orebro University Business School, 701 82 ¨Orebro, Sweden. E-mail:

Sune.Karlsson@oru.se

(3)

Geweke (1996) studied the reduced rank regression model Y = XΘ + ZA + E where

Y is an n × L matrix of dependent variables, X and Z contains p and k explanatory

variables and Θ and A are parameter matrices where Θ is assumed to have reduced

rank q < min (L, p) . The rows of E are assumed to be independent normal with mean

zero and variance matrix Σ. Under the reduced rank assumption Θ can be factored into

Θ = ΨΦ with Ψ a p × q matrix and Φ a q × L matrix, both of rank q. To identify the

model Geweke considers two normalizations, Φ = (Iq, Φ∗) with Ψ unrestricted

(normal-ization 1) and Ψ0 = (Iq, Ψ∗0) with Φ unrestricted (normalization 2). Geweke provides

full conditional posterior distributions for an informative prior where Σ is distributed

as inverse Wishart with parameter matrix S and v degrees of freedom, iW (S, v) , and

the elements of A, (Ψ, Φ∗) or (Φ∗, Ψ) are independent normal with mean zero and

vari-ance 1/τ2. Geweke also considers a flat, improper, prior on A, (Ψ, Φ

) or (Φ∗, Ψ) which

corresponds to setting τ = 0 in the independent normals.

The full conditional posteriors for the improper prior are correctly stated in Geweke

(1996) while the full conditional posteriors are incorrect in several cases for the proper

informative prior. The incorrect expressions for the parameters of the posterior

distribu-tions have, unfortunately, been picked up in the literature and used to construct Gibbs

samplers with incorrect stationary distributions. Geweke (2004) developed a method for

checking the correctness of posterior simulators and detected problems with the Gibbs

sampler coded up for the 1996 paper but failed to connect this with the incorrect

expres-sions for the full conditional posteriors.

The mean of the full conditional posterior for Ψ in normalization 1 given in equation

(11) of Geweke is incorrect. Let ψ = vec (Ψ) , the correct full conditional posterior is

normal, ψ|Y, A, Φ, Σ ∼ N ψ, Vψ , for

Vψ = τ2I + ΦΣ−1Φ0⊗ X0X −1 , ψ = Vψvec X0Y∗Σ−1Φ0  (1) where Y∗ = Y − ZA. 2

(4)

The mean and variance of the full conditional posterior for Φ∗ in normalization 1 given

in equation (13) of Geweke are incorrect. Let φ∗ = vec (Φ∗) , the correct full conditional

posterior is normal, φ∗|Y, A, Ψ, Σ ∼ Nφ∗, Vφ∗

 , for Vφ∗ =τ2I + Σ22⊗ Ψ0X0XΨ −1 , φ∗ = Vφ∗vecΨ0X0 Y∗ 1Σ 12− XΨΣ12+ Y∗ 2Σ 22 (2)

where (Y∗1, Y2∗) = Y∗ partitions Y∗ into n × q and n × (L − q) matrices and Σ12 and Σ22

are the upper right q × (L − q) and lower right (L − q) × (L − q) submatrices of Σ−1.

The mean and variance of the full conditional posterior for Ψ∗in normalization 2 given

in equation (15) of Geweke are incorrect. Let ψ∗ = vec (Ψ∗) , the correct full conditional

posterior is normal, ψ∗|Y, A, Φ, Σ ∼ Nψ∗, Vψ∗

 for

Vψ∗ = τ2I + ΦΣ−1Φ0 ⊗ X02X2

−1

, ψ∗ = Vψ∗vec X02(Y∗− X1Φ) Σ−1Φ0 (3)

where X1 contains the first q columns of X and X2 the remaining columns

Additional details and derivations are provided in Karlsson (2012). The correct

ex-pressions for the improper prior in Geweke (1996) can be obtained as special cases by

setting τ = 0 in (1), (2) and (3).

References

Geweke, J. (1996), ‘Bayesian reduced rank regression in econometrics’, Journal of

Econo-metrics 75, 121–146.

Geweke, J. (2004), ‘Getting it right: Joint distribution tests of posterior simulators’,

Journal of the American Statistical Association 99, 799–804.

Karlsson, S. (2012), Conditional posteriors for the reduced rank regression model, Working

Papers 2012:11, ¨Orebro University Business School.

References

Related documents

När testpersonen börjar bli trött så kommer det märkas att denne inte alltid orkar ända fram till linjen innan nästa beepljud drar igång.. Tp måste fortfarande nudda linjen

Faktorer såsom kommunikation mellan enheter och individens förmåga att förhålla sig till i dialog benämns som gynnande för att skapa en högre kompetens inom verksamheter

Hvor omfattende muggsoppangrepet er, vil også avhenge av andre faktorer som treslag, om trelasten er infi sert av muggsopp eller sporer før impregnering, mye muggsporer i

Flera arbetsterapeuter menar därför att det är viktigt att inte vara för många personer samtidigt när de ska möta denna klientgrupp.. För att undvika detta belyser ett

konkurrensmässiga, för förr eller senare när det här arbetet går framåt och vi kommer att se att vi får företag som till och med kan slå sig för bröstet och säga” hos

In order to evaluate the LC-PCR assay with clinical material, specimens from men and women attending the STD Clinic were analysed at the local laboratory in Örebro, Sweden using

Johan Holmgren and Pernilla Ivehammar, Making headway towards a better understanding of service frequency valuations: a study of how the relative valuation of train service frequency

The answer we have come up with in this chapter is that by regarding a basic supply of local bus services as merit goods, and providing a level of service above the basic minimum