Ideology and Social Networks in the U.S. Congress

Size: px
Start display at page:

Download "Ideology and Social Networks in the U.S. Congress"

Transcription

1 Ideology and Social Networks in the U.S. Congress James H. Fowler University of California, San Diego July 11, 2007

2 Improve connectedness scores (Fowler 2005) Cosponsorship is about 2 things: The idea The person behind the idea The goal is to create a Social Space Model that simultaneously estimates ideological and social ideal points. Place legislators and bills in policy space (Poole and Rosenthal 1997; Clinton, Jackman, and Rivers 2004) Place sponsors and cosponsors in social space (Hoff, Raftery, and Handcock 2002) Create alternative Social Utility Model

3 Utility from Policy Social Utility Other factors The Double-Center Operator n legislators choose whether to cosponsor m bills. Each bill j = 1...m is sponsored by legislator s(j) Each legislator i = 1...n chooses between the bill at ζ j and the status quo located at ψ j, both in R d { 1 if legislator i cosponosors bill j y ij = 0 otherwise (1)

4 Utility from Policy Social Utility Other factors The Double-Center Operator Legislator i receives utility for supporting policies close to her own ideal point x i in R d policy space: U bill status quo ij = x i ζ j 2 Uij = x i ψ j 2 U cosponsor ij = x i ζ j 2 + x i ψ j 2 (2)

5 Utility from Policy Social Utility Other factors The Double-Center Operator Legislator i also receives utility for supporting people close to her own social ideal point in R δ social space λ i = legislator i s ideal point in R δ social sender space (or hub space, Kleinberg 1999) αs(j) = the ideal point for bill j s sponsor in R δ social receiver space (or authority space) U social ij = λ i α s(j) 2 U ideological ij = x i x s(j) 2 U personal ij = λ i α s(j) 2 κ x i x s(j) 2 (3)

6 Utility from Policy Social Utility Other factors The Double-Center Operator Utility also affected by bill-specific factors φ j and legislator-specific factors ξ i : U ij = x i ζ j 2 + x i ψ j 2 λ i α s(j) 2 This simplifies to where κ x i x s(j) 2 + ξ i + φ j (4) U ij = 2x i β j + 2λ i α s(j) + η i + θ j (5) β j = ψ j ζ j κx s(j), η i = λ 2 i κx 2 i + ξ i θ j = ψ 2 j ζ 2 j x 2 s(j) κx2 s(j) + φ j

7 Utility from Policy Social Utility Other factors The Double-Center Operator Suppose observed cosponsorship means true utility of cosponsorship is high (1), while not observing one means that the true utility is low (0). True utility is observed cosponsorship decision minus an error term: Substituting, we get: U ij = y ij ν ij (6) y ij = 2x i β j + 2λ i α s(j) + η i + θ j + ν ij (7)

8 Utility from Policy Social Utility Other factors The Double-Center Operator Define the double-center operator D(.) for a matrix Z (Poole 2005; Clinton, Jackman, Rivers 2005): D(z ij ) = (z ij z i. z.j + z.. )/( 2) (8) Suppose ν ij is i.i.d. stable density and dimension-by-dimension means of x, β, λ, and α s(j) equal 0. Apply double-center operator in equation (8) to both sides of equation (7): D(y ij ) = x i β j λ i α s(j) + ɛ ij (9)

9 Utility from Policy Social Utility Other factors The Double-Center Operator By construction, social distance (λ i α s(j) ) is uncorrelated with ideological distance (x i β j ). Therefore, singular value decomposition (SVD) provides the best d dimensional approximation of x i and β j (Frobenius-Perron): D(Y) = XΣB (10)

10 Utility from Policy Social Utility Other factors The Double-Center Operator Suppose each legislator k sponsors n k bills. We can define an n by n matrix such that: µ ik = 1 n k s(j)=k D(y ij ) = x i βk λ i α k + υ ik (11) where β k = 1 n k s(j)=k β j and υ ik is a stable density (but not i.i.d. across sponsors). Rearranging yields another SVD formulation: XΣB M = ΛTA (12)

11 Utility from Policy Social Utility Other factors The Double-Center Operator Standard errors via Metropolis-Hastings! If errors ɛ ij in equation (9) are normally distributed, then: L(x i, β j, λ i, α s(j) y ij ) = ( ) 2 D(yij ) x i β j + λ i α s(j) i,j,(i j) (13)

12 Utility Add Error and Double Center Estimate dyad-specific mean utility transfers γ ik from legislator k to legislator i. These utilities can be thought of as social influence vote-trading specific instances of monetary transfers from one legislator to another (such as PAC contributions to a friend s campaign) To control for legislative activity, assume γ i. = γ.s(j) = γ.. = 0.

13 Utility Add Error and Double Center Personal utility from equation (3) becomes: U personal ij = γ is(j) κ x i x s(j) 2 (14) The full utility equation from (4) becomes: which simplifies to: U ij = x i ζ j 2 + x i ψ j 2 + γ is(j) κ x i x s(j) 2 + ξ i + φ j (15) U ij = 2x i β j + γ is(j) + η i + θ j (16)

14 Utility Add Error and Double Center Establishing a relationship between cosponsorship and the utility it provides yields an analogue to equation (7): y ij = 2x i β j + γ is(j) + η i + θ j + ν ij (17) Double-centering both sides yields: D(y ij ) = x i β j γ is(j) 2 + ɛ ij (18)

15 Utility Add Error and Double Center The SVD provides the best d dimensional approximation of x i and β j : D(Y) = XΣB (19) Suppose each legislator k sponsors n k bills. We can define an n by n matrix such that: µ ik = 1 n k s(j)=k D(y ij ) = x i βk γ ik 2 Rearranging yields the social utility matrix: (20) Γ = 2(XΣB M) (21)

16 Utility Add Error and Double Center Standard errors via Metropolis-Hastings! L(x i, β j y ij ) = i,j (D(y ij ) x i β j ) 2 (22)

17 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Generate true data with independent random draw from a uniform distribution for each legislator (x) bill (ζ) status quo (ψ) and either: legislator s sender space (λ) and receiver space (α) for Social Space Model legislator dyad s specific social utility (γ) for Social Utility Model

18 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Parameterize social weight: U total ij = ωu social ij + (1 ω)u ideological ij (23) Normally-disributed error term Choose a constant threshold U to generate specific rate of cosponsorship such that { 1 if Uij > U y ij = 0 otherwise (24) Option to change the utilities from quadratic u 2 to gaussian σφ(u).

19 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Baseline configuration: 100 legislators, each sponsors on average 10 bills (for 1000 total bills) rate of cosponsorship is set to 0.1 (Fowler 2005) the effect of social utility is approximately equal to ideological utility (ω = 0.5) no error in the utility utilities are quadratic Use Spearman rank correlation (doesn t matter).

20 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Legislator Ideal Point Bill Parameter Ideological Utility Social Utility ρ ρ ρ NOM ρ ρ ρ NOM ρ ρ ρ NOM ρ λ α Baseline , , , , , , , Bills per Legislator 0.966, , , , , , , Bill per Legislator 0.935, , , , , , , Legislators (1 Bill Each) 0.977, , , , , , , Social Weight 0.992, , , , , , , Social Weight 0.997, , , , , , ,0.022

21 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Legislator Ideal Point Bill Parameter Ideological Utility Social Utility ρ ρ ρ NOM ρ ρ ρ NOM ρ ρ ρ NOM ρ λ α Baseline , , , , , , , Rate of Cosponsorship 0.974, , , , , , , Rate of Cosponsorship 0.965, , , , , , , Standard Deviation Error 0.972, , , , , , , Standard Deviation Error 0.968, , , , , , ,0.843 Gaussian Utility 0.871, , , , , , ,0.829

22 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Legislator Ideal Point Bill Parameter Ideological Utility Social Utility ρ ρ ρ NOM ρ ρ ρ NOM ρ ρ ρ NOM ρ γ Baseline , , , , , , , Bills per Legislator 0.987, , , , , , , Bill per Legislator 0.947, , , , , , , Legislators (1 Bill Each) 0.991, , , , , , , Social Weight 0.997, , , , , , , Social Weight 0.997, , , , , , ,0.020

23 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Legislator Ideal Point Bill Parameter Ideological Utility Social Utility ρ ρ ρ NOM ρ ρ ρ NOM ρ ρ ρ NOM ρ γ Baseline , , , , , , , Rate of Cosponsorship 0.990, , , , , , , Rate of Cosponsorship 0.993, , , , , , , Standard Deviation Error 0.991, , , , , , , Standard Deviation Error 0.990, , , , , , ,0.721 Gaussian Utility 0.910, , , , , , ,0.732

24 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Social Space Model Social Utility Model Congress Years ρ x ρ β ρ λ α ρ x ρ β ρ γ 93rd , , , , , , th , , , , , , th , , , , , , th , , , , , , th , , , , , , th , , , , , , th , , , , , , th , , , , , ,0.933

25 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Social Space Model Social Utility Model Congress Years ρ x ρ β ρ λ α ρ x ρ β ρ γ 101st , , , , , , nd , , , , , , rd , , , , , , th , , , , , , th , , , , , , th , , , , , , th , , , , , , th , , , , , ,0.925

26 The Process Data Generation Baseline Configuration and Evlauation Social Space Model Results Social Utility Model Results U.S. Senate Results Social Space Model Social Utility Model Real Data, 99th Senate Singular Value Singular Value Singular Value Eigenvalue rank Eigenvalue rank Eigenvalue rank is Superior for Cosponsorship Outlier values in left panel indicate two structured dimensions in data generated by the Social Space Model Center panel shows lack of structure in data generated by the Social Utility Model Right panel shows a lack of structure in data taken from the 99th Senate

27 We can recover social utilities without ideology! SVD procedure (surprisingly) works better than W-NOMINATE Social Utility Model works best on real data Can generalize to simultaneous modeling of any bipartite network with a unipartite projection.

SPATIAL VOTING (MULTIPLE DIMENSIONS)

SPATIAL VOTING (MULTIPLE DIMENSIONS) SPATIAL VOTING (MULTIPLE DIMENSIONS) 1 Assumptions Alternatives are points in an n-dimensional space. Examples for 2D: Social Issues and Economic Issues Domestic Spending and Foreign Spending Single-peaked

More information

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Friday April 1 ± ǁ

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Friday April 1 ± ǁ . α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω Friday April 1 ± ǁ 1 Chapter 5. Photons: Covariant Theory 5.1. The classical fields 5.2. Covariant

More information

Statistics 202: Data Mining. c Jonathan Taylor. Week 2 Based in part on slides from textbook, slides of Susan Holmes. October 3, / 1

Statistics 202: Data Mining. c Jonathan Taylor. Week 2 Based in part on slides from textbook, slides of Susan Holmes. October 3, / 1 Week 2 Based in part on slides from textbook, slides of Susan Holmes October 3, 2012 1 / 1 Part I Other datatypes, preprocessing 2 / 1 Other datatypes Document data You might start with a collection of

More information

Part I. Other datatypes, preprocessing. Other datatypes. Other datatypes. Week 2 Based in part on slides from textbook, slides of Susan Holmes

Part I. Other datatypes, preprocessing. Other datatypes. Other datatypes. Week 2 Based in part on slides from textbook, slides of Susan Holmes Week 2 Based in part on slides from textbook, slides of Susan Holmes Part I Other datatypes, preprocessing October 3, 2012 1 / 1 2 / 1 Other datatypes Other datatypes Document data You might start with

More information

Hierarchical Modeling for Univariate Spatial Data

Hierarchical Modeling for Univariate Spatial Data Hierarchical Modeling for Univariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Spatial Domain 2 Geography 890 Spatial Domain This

More information

CSSS/STAT/SOC 321 Case-Based Social Statistics I. Levels of Measurement

CSSS/STAT/SOC 321 Case-Based Social Statistics I. Levels of Measurement CSSS/STAT/SOC 321 Case-Based Social Statistics I Levels of Measurement Christopher Adolph Department of Political Science and Center for Statistics and the Social Sciences University of Washington, Seattle

More information

Estimating the Distribution of Voter Preferences Using Partially Aggregated Voting Data

Estimating the Distribution of Voter Preferences Using Partially Aggregated Voting Data Estimating the Distribution of Voter Preferences Using Partially Aggregated Voting Data James M. Snyder, Jr. Department of Political Science and Department of Economics Massachusetts Institute of Technology

More information

Unit 5: Centrality. ICPSR University of Michigan, Ann Arbor Summer 2015 Instructor: Ann McCranie

Unit 5: Centrality. ICPSR University of Michigan, Ann Arbor Summer 2015 Instructor: Ann McCranie Unit 5: Centrality ICPSR University of Michigan, Ann Arbor Summer 2015 Instructor: Ann McCranie What does centrality tell us? We often want to know who the most important actors in a network are. Centrality

More information

Stat 542: Item Response Theory Modeling Using The Extended Rank Likelihood

Stat 542: Item Response Theory Modeling Using The Extended Rank Likelihood Stat 542: Item Response Theory Modeling Using The Extended Rank Likelihood Jonathan Gruhl March 18, 2010 1 Introduction Researchers commonly apply item response theory (IRT) models to binary and ordinal

More information

Hierarchical Modelling for Univariate Spatial Data

Hierarchical Modelling for Univariate Spatial Data Hierarchical Modelling for Univariate Spatial Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department

More information

sparse and low-rank tensor recovery Cubic-Sketching

sparse and low-rank tensor recovery Cubic-Sketching Sparse and Low-Ran Tensor Recovery via Cubic-Setching Guang Cheng Department of Statistics Purdue University www.science.purdue.edu/bigdata CCAM@Purdue Math Oct. 27, 2017 Joint wor with Botao Hao and Anru

More information

Department of Statistics. Bayesian Modeling for a Generalized Social Relations Model. Tyler McCormick. Introduction.

Department of Statistics. Bayesian Modeling for a Generalized Social Relations Model. Tyler McCormick. Introduction. A University of Connecticut and Columbia University A models for dyadic data are extensions of the (). y i,j = a i + b j + γ i,j (1) Here, y i,j is a measure of the tie from actor i to actor j. The random

More information

CSE 1400 Applied Discrete Mathematics Definitions

CSE 1400 Applied Discrete Mathematics Definitions CSE 1400 Applied Discrete Mathematics Definitions Department of Computer Sciences College of Engineering Florida Tech Fall 2011 Arithmetic 1 Alphabets, Strings, Languages, & Words 2 Number Systems 3 Machine

More information

Hierarchical Modelling for Univariate Spatial Data

Hierarchical Modelling for Univariate Spatial Data Spatial omain Hierarchical Modelling for Univariate Spatial ata Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A.

More information

Further Maths A2 (M2FP2D1) Assignment ψ (psi) A Due w/b 19 th March 18

Further Maths A2 (M2FP2D1) Assignment ψ (psi) A Due w/b 19 th March 18 α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω The mathematician s patterns, like the painter s or the poet s, must be beautiful: the ideas, like the colours or the words, must fit together in a harmonious

More information

Recommender systems, matrix factorization, variable selection and social graph data

Recommender systems, matrix factorization, variable selection and social graph data Recommender systems, matrix factorization, variable selection and social graph data Julien Delporte & Stéphane Canu stephane.canu@litislab.eu StatLearn, april 205, Grenoble Road map Model selection for

More information

and in each case give the range of values of x for which the expansion is valid.

and in each case give the range of values of x for which the expansion is valid. α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)

More information

On an Additive Semigraphoid Model for Statistical Networks With Application to Nov Pathway 25, 2016 Analysis -1 Bing / 38Li,

On an Additive Semigraphoid Model for Statistical Networks With Application to Nov Pathway 25, 2016 Analysis -1 Bing / 38Li, On an Additive Semigraphoid Model for Statistical Networks With Application to Pathway Analysis - Bing Li, Hyunho Chun & Hongyu Zhao Kim Youngrae SNU Stat. Multivariate Lab Nov 25, 2016 On an Additive

More information

CONVEX OPTIMIZATION LEARNING OF FAITHFUL EUCLIDEAN DISTANCE REPRESENTATIONS IN NONLINEAR DIMENSIONALITY REDUCTION

CONVEX OPTIMIZATION LEARNING OF FAITHFUL EUCLIDEAN DISTANCE REPRESENTATIONS IN NONLINEAR DIMENSIONALITY REDUCTION CONVEX OPTIMIZATION LEARNING OF FAITHFUL EUCLIDEAN DISTANCE REPRESENTATIONS IN NONLINEAR DIMENSIONALITY REDUCTION CHAO DING AND HOU-DUO QI Abstract. Classical multidimensional scaling only works well when

More information

Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.

Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we

More information

Fixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility

Fixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility American Economic Review: Papers & Proceedings 2016, 106(5): 400 404 http://dx.doi.org/10.1257/aer.p20161082 Fixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility By Gary Chamberlain*

More information

Concentration Ellipsoids

Concentration Ellipsoids Concentration Ellipsoids ECE275A Lecture Supplement Fall 2008 Kenneth Kreutz Delgado Electrical and Computer Engineering Jacobs School of Engineering University of California, San Diego VERSION LSECE275CE

More information

Dimension reduction, PCA & eigenanalysis Based in part on slides from textbook, slides of Susan Holmes. October 3, Statistics 202: Data Mining

Dimension reduction, PCA & eigenanalysis Based in part on slides from textbook, slides of Susan Holmes. October 3, Statistics 202: Data Mining Dimension reduction, PCA & eigenanalysis Based in part on slides from textbook, slides of Susan Holmes October 3, 2012 1 / 1 Combinations of features Given a data matrix X n p with p fairly large, it can

More information

Bayesian non-parametric model to longitudinally predict churn

Bayesian non-parametric model to longitudinally predict churn Bayesian non-parametric model to longitudinally predict churn Bruno Scarpa Università di Padova Conference of European Statistics Stakeholders Methodologists, Producers and Users of European Statistics

More information

Advanced Quantitative Research Methodology Lecture Notes: January Ecological 28, 2012 Inference1 / 38

Advanced Quantitative Research Methodology Lecture Notes: January Ecological 28, 2012 Inference1 / 38 Advanced Quantitative Research Methodology Lecture Notes: Ecological Inference 1 Gary King http://gking.harvard.edu January 28, 2012 1 c Copyright 2008 Gary King, All Rights Reserved. Gary King http://gking.harvard.edu

More information

Modeling conditional distributions with mixture models: Theory and Inference

Modeling conditional distributions with mixture models: Theory and Inference Modeling conditional distributions with mixture models: Theory and Inference John Geweke University of Iowa, USA Journal of Applied Econometrics Invited Lecture Università di Venezia Italia June 2, 2005

More information

Linear Algebra and Dirac Notation, Pt. 3

Linear Algebra and Dirac Notation, Pt. 3 Linear Algebra and Dirac Notation, Pt. 3 PHYS 500 - Southern Illinois University February 1, 2017 PHYS 500 - Southern Illinois University Linear Algebra and Dirac Notation, Pt. 3 February 1, 2017 1 / 16

More information

Two-sample hypothesis testing for random dot product graphs

Two-sample hypothesis testing for random dot product graphs Two-sample hypothesis testing for random dot product graphs Minh Tang Department of Applied Mathematics and Statistics Johns Hopkins University JSM 2014 Joint work with Avanti Athreya, Vince Lyzinski,

More information

THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR

THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR WEN LI AND MICHAEL K. NG Abstract. In this paper, we study the perturbation bound for the spectral radius of an m th - order n-dimensional

More information

FINM 331: MULTIVARIATE DATA ANALYSIS FALL 2017 PROBLEM SET 3

FINM 331: MULTIVARIATE DATA ANALYSIS FALL 2017 PROBLEM SET 3 FINM 331: MULTIVARIATE DATA ANALYSIS FALL 2017 PROBLEM SET 3 The required files for all problems can be found in: http://www.stat.uchicago.edu/~lekheng/courses/331/hw3/ The file name indicates which problem

More information

Fast Estimation of Ideal Points with Massive Data

Fast Estimation of Ideal Points with Massive Data Fast Estimation of Ideal Points with Massive Data Kosuke Imai James Lo Jonathan Olmsted First Draft: December, 014 This Draft: February 5, 015 Abstract Estimating ideological preferences across time and

More information

Complex Networks CSYS/MATH 303, Spring, Prof. Peter Dodds

Complex Networks CSYS/MATH 303, Spring, Prof. Peter Dodds Complex Networks CSYS/MATH 303, Spring, 2011 Prof. Peter Dodds Department of Mathematics & Statistics Center for Complex Systems Vermont Advanced Computing Center University of Vermont Licensed under the

More information

Chapter 6. Intuitionistic Fuzzy PROMETHEE Technique. AIDS stands for acquired immunodeficiency syndrome. AIDS is the final stage. 6.

Chapter 6. Intuitionistic Fuzzy PROMETHEE Technique. AIDS stands for acquired immunodeficiency syndrome. AIDS is the final stage. 6. Chapter 6 Intuitionistic Fuzzy PROMETHEE Technique 6.1 Introduction AIDS stands for acquired immunodeficiency syndrome. AIDS is the final stage of HIV infection, and not everyone who has HIV advances to

More information

A Mixture Modeling Approach for Empirical Testing of Competing Theories

A Mixture Modeling Approach for Empirical Testing of Competing Theories A Mixture Modeling Approach for Empirical Testing of Competing Theories Kosuke Imai Dustin Tingley Princeton University January 8, 2010 Imai and Tingley (Princeton) Mixture Modeling Gakusyuin 2010 1 /

More information

Trade and Inequality: From Theory to Estimation

Trade and Inequality: From Theory to Estimation Trade and Inequality: From Theory to Estimation Elhanan Helpman Oleg Itskhoki Marc Muendler Stephen Redding Harvard Princeton UC San Diego Princeton MEF Italia Dipartimento del Tesoro September 2014 1

More information

Spatial Process Estimates as Smoothers: A Review

Spatial Process Estimates as Smoothers: A Review Spatial Process Estimates as Smoothers: A Review Soutir Bandyopadhyay 1 Basic Model The observational model considered here has the form Y i = f(x i ) + ɛ i, for 1 i n. (1.1) where Y i is the observed

More information

Influence Maximization in Social Networks: An Ising-model-based Approach

Influence Maximization in Social Networks: An Ising-model-based Approach Influence Maximization in Social etworks: An Ising-model-based Approach Shihuan Liu, Lei Ying, and Srinivas Shakkottai Department of Electrical and Computer Engineering, Iowa State University Email: {liush08,

More information

last name ID 1 c/cmaker/cbreaker 2012 exam version a 6 pages 3 hours 40 marks no electronic devices SHOW ALL WORK

last name ID 1 c/cmaker/cbreaker 2012 exam version a 6 pages 3 hours 40 marks no electronic devices SHOW ALL WORK last name ID 1 c/cmaker/cbreaker 2012 exam version a 6 pages 3 hours 40 marks no electronic devices SHOW ALL WORK 8 a b c d e f g h i j k l m n o p q r s t u v w x y z 1 b c d e f g h i j k l m n o p q

More information

Distributed Detection via Gaussian Running Consensus: Large Deviations Asymptotic Analysis

Distributed Detection via Gaussian Running Consensus: Large Deviations Asymptotic Analysis Distributed Detection via Gaussian Running Consensus: Large Deviations Asymptotic Analysis Dragana Bajović, Du san Jakovetić, João Xavier, Bruno Sinopoli and José M. F. Moura Abstract We study, by large

More information

Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US

Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Gerdie Everaert 1, Lorenzo Pozzi 2, and Ruben Schoonackers 3 1 Ghent University & SHERPPA 2 Erasmus

More information

Agent-Based Methods for Dynamic Social Networks. Duke University

Agent-Based Methods for Dynamic Social Networks. Duke University Agent-Based Methods for Dynamic Social Networks Eric Vance Institute of Statistics & Decision Sciences Duke University STA 395 Talk October 24, 2005 Outline Introduction Social Network Models Agent-based

More information

Writing Game Theory in L A TEX

Writing Game Theory in L A TEX Writing Game Theory in L A TEX Thiago Silva First Version: November 22, 2015 This Version: November 13, 2017 List of Figures and Tables 1 2x2 Matrix: Prisoner s ilemma Normal-Form Game............. 3 2

More information

Sticky Leverage. João Gomes, Urban Jermann & Lukas Schmid Wharton School and UCLA/Duke. September 28, 2013

Sticky Leverage. João Gomes, Urban Jermann & Lukas Schmid Wharton School and UCLA/Duke. September 28, 2013 Sticky Leverage João Gomes, Urban Jermann & Lukas Schmid Wharton School and UCLA/Duke September 28, 213 Introduction Models of monetary non-neutrality have traditionally emphasized the importance of sticky

More information

I. Relationship with previous work

I. Relationship with previous work x x i t j J t = {0, 1,...J t } j t (p jt, x jt, ξ jt ) p jt R + x jt R k k ξ jt R ξ t T j = 0 t (z i, ζ i, G i ), ζ i z i R m G i G i (p j, x j ) i j U(z i, ζ i, x j, p j, ξ j ; G i ) = u(ζ i, x j,

More information

Mathematics Review Exercises. (answers at end)

Mathematics Review Exercises. (answers at end) Brock University Physics 1P21/1P91 Mathematics Review Exercises (answers at end) Work each exercise without using a calculator. 1. Express each number in scientific notation. (a) 437.1 (b) 563, 000 (c)

More information

Y1 Double Maths Assignment λ (lambda) Exam Paper to do Core 1 Solomon C on the VLE. Drill

Y1 Double Maths Assignment λ (lambda) Exam Paper to do Core 1 Solomon C on the VLE. Drill α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Nature is an infinite sphere of which the centre is everywhere and the circumference nowhere Blaise Pascal Y Double Maths Assignment λ (lambda) Tracking

More information

Estimation of Treatment Effects under Essential Heterogeneity

Estimation of Treatment Effects under Essential Heterogeneity Estimation of Treatment Effects under Essential Heterogeneity James Heckman University of Chicago and American Bar Foundation Sergio Urzua University of Chicago Edward Vytlacil Columbia University March

More information

Statistical Inference on Large Contingency Tables: Convergence, Testability, Stability. COMPSTAT 2010 Paris, August 23, 2010

Statistical Inference on Large Contingency Tables: Convergence, Testability, Stability. COMPSTAT 2010 Paris, August 23, 2010 Statistical Inference on Large Contingency Tables: Convergence, Testability, Stability Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu COMPSTAT

More information

Day two: Linear Algebra

Day two: Linear Algebra Day two: Linear Algebra Methods camp instructors September 7th, 2016 1 / 94 Outline Vectors Basic notation Multiplying by a scalar Addition and subtraction (sidenote on conformability) Two forms of vector

More information

Random Effects Models for Network Data

Random Effects Models for Network Data Random Effects Models for Network Data Peter D. Hoff 1 Working Paper no. 28 Center for Statistics and the Social Sciences University of Washington Seattle, WA 98195-4320 January 14, 2003 1 Department of

More information

Exogeneity tests and weak identification

Exogeneity tests and weak identification Cireq, Cirano, Départ. Sc. Economiques Université de Montréal Jean-Marie Dufour Cireq, Cirano, William Dow Professor of Economics Department of Economics Mcgill University June 20, 2008 Main Contributions

More information

Jed Chou. April 13, 2015

Jed Chou. April 13, 2015 of of CS598 AGB April 13, 2015 Overview of 1 2 3 4 5 Competing Approaches of Two competing approaches to species tree inference: Summary methods: estimate a tree on each gene alignment then combine gene

More information

ELEC633: Graphical Models

ELEC633: Graphical Models ELEC633: Graphical Models Tahira isa Saleem Scribe from 7 October 2008 References: Casella and George Exploring the Gibbs sampler (1992) Chib and Greenberg Understanding the Metropolis-Hastings algorithm

More information

Bindel, Fall 2016 Matrix Computations (CS 6210) Notes for

Bindel, Fall 2016 Matrix Computations (CS 6210) Notes for 1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.

More information

Granger Causality and Dynamic Structural Systems 1

Granger Causality and Dynamic Structural Systems 1 Granger Causality and Dynamic Structural Systems 1 Halbert White and Xun Lu Department of Economics, University of California, San Diego December 10, 2009 1 forthcoming, Journal of Financial Econometrics

More information

Contents. basic algebra. Learning outcomes. Time allocation. 1. Mathematical notation and symbols. 2. Indices. 3. Simplification and factorisation

Contents. basic algebra. Learning outcomes. Time allocation. 1. Mathematical notation and symbols. 2. Indices. 3. Simplification and factorisation basic algebra Contents. Mathematical notation and symbols 2. Indices 3. Simplification and factorisation 4. Arithmetic of algebraic fractions 5. Formulae and transposition Learning outcomes In this workbook

More information

The Return of the Gibson Paradox

The Return of the Gibson Paradox The Return of the Gibson Paradox Timothy Cogley, Thomas J. Sargent, and Paolo Surico Conference to honor Warren Weber February 212 Introduction Previous work suggested a decline in inflation persistence

More information

On the complexity of stoquastic Hamiltonians

On the complexity of stoquastic Hamiltonians On the complexity of stoquastic Hamiltonians Ian Kivlichan December 11, 2015 Abstract Stoquastic Hamiltonians, those for which all off-diagonal matrix elements in the standard basis are real and non-positive,

More information

How to Use the Internet for Election Surveys

How to Use the Internet for Election Surveys How to Use the Internet for Election Surveys Simon Jackman and Douglas Rivers Stanford University and Polimetrix, Inc. May 9, 2008 Theory and Practice Practice Theory Works Doesn t work Works Great! Black

More information

ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS

ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS ECO 513 Fall 2009 C. Sims HIDDEN MARKOV CHAIN MODELS 1. THE CLASS OF MODELS y t {y s, s < t} p(y t θ t, {y s, s < t}) θ t = θ(s t ) P[S t = i S t 1 = j] = h ij. 2. WHAT S HANDY ABOUT IT Evaluating the

More information

Analytically tractable processes on networks

Analytically tractable processes on networks University of California San Diego CERTH, 25 May 2011 Outline Motivation 1 Motivation Networks Random walk and Consensus Epidemic models Spreading processes on networks 2 Networks Motivation Networks Random

More information

Modeling homophily and stochastic equivalence in symmetric relational data

Modeling homophily and stochastic equivalence in symmetric relational data Modeling homophily and stochastic equivalence in symmetric relational data Peter D. Hoff Departments of Statistics and Biostatistics University of Washington Seattle, WA 98195-4322. hoff@stat.washington.edu

More information

High dimensional Ising model selection

High dimensional Ising model selection High dimensional Ising model selection Pradeep Ravikumar UT Austin (based on work with John Lafferty, Martin Wainwright) Sparse Ising model US Senate 109th Congress Banerjee et al, 2008 Estimate a sparse

More information

Chapter 9: Hypothesis Testing Sections

Chapter 9: Hypothesis Testing Sections 1 / 22 : Hypothesis Testing Sections Skip: 9.2 Testing Simple Hypotheses Skip: 9.3 Uniformly Most Powerful Tests Skip: 9.4 Two-Sided Alternatives 9.5 The t Test 9.6 Comparing the Means of Two Normal Distributions

More information

Discrete Variables and Gradient Estimators

Discrete Variables and Gradient Estimators iscrete Variables and Gradient Estimators This assignment is designed to get you comfortable deriving gradient estimators, and optimizing distributions over discrete random variables. For most questions,

More information

Parameterizing Exponential Family Models for Random Graphs: Current Methods and New Directions

Parameterizing Exponential Family Models for Random Graphs: Current Methods and New Directions Carter T. Butts p. 1/2 Parameterizing Exponential Family Models for Random Graphs: Current Methods and New Directions Carter T. Butts Department of Sociology and Institute for Mathematical Behavioral Sciences

More information

Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017

Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Put your solution to each problem on a separate sheet of paper. Problem 1. (5106) Let X 1, X 2,, X n be a sequence of i.i.d. observations from a

More information

Distributed Detection over Random Networks: Large Deviations Performance Analysis

Distributed Detection over Random Networks: Large Deviations Performance Analysis Distributed Detection over Random Networks: Large Deviations Performance Analysis Dragana Bajović, Du san Jakovetić, João Xavier, Bruno Sinopoli and José M. F. Moura arxiv:02.4668v [cs.it] 2 Dec 200 Abstract

More information

ESTIMATION: ONE DIMENSION (PART 1)

ESTIMATION: ONE DIMENSION (PART 1) ESTIMATION: ONE DIMENSION (PART 1) 1 Goals of this Unit Examine how we can use real world data to estimate the location of ideal points (and cut points). Interest group scores. Optimal classification.

More information

Bayesian Melding. Assessing Uncertainty in UrbanSim. University of Washington

Bayesian Melding. Assessing Uncertainty in UrbanSim. University of Washington Bayesian Melding Assessing Uncertainty in UrbanSim Hana Ševčíková University of Washington hana@stat.washington.edu Joint work with Paul Waddell and Adrian Raftery University of Washington UrbanSim Workshop,

More information

The Singular Value Decomposition

The Singular Value Decomposition The Singular Value Decomposition Philippe B. Laval KSU Fall 2015 Philippe B. Laval (KSU) SVD Fall 2015 1 / 13 Review of Key Concepts We review some key definitions and results about matrices that will

More information

R k. t + 1. n E t+1 = ( 1 χ E) W E t+1. c E t+1 = χ E Wt+1 E. Γ E t+1. ) R E t+1q t K t. W E t+1 = ( 1 Γ E t+1. Π t+1 = P t+1 /P t

R k. t + 1. n E t+1 = ( 1 χ E) W E t+1. c E t+1 = χ E Wt+1 E. Γ E t+1. ) R E t+1q t K t. W E t+1 = ( 1 Γ E t+1. Π t+1 = P t+1 /P t R k E 1 χ E Wt E n E t+1 t t + 1 n E t+1 = ( 1 χ E) W E t+1 c E t+1 = χ E Wt+1 E t + 1 q t K t Rt+1 E 1 Γ E t+1 Π t+1 = P t+1 /P t W E t+1 = ( 1 Γ E t+1 ) R E t+1q t K t Π t+1 Γ E t+1 K t q t q t K t j

More information

Detection And Estimation Of Block Structure In Spatial Weight Matrix

Detection And Estimation Of Block Structure In Spatial Weight Matrix Detection And Estimation Of Block Structure In Spatial Weight Matrix Clifford Lam and Pedro CL Souza 2 Department of Statistics, London School of Economics and Political Science 2 Department of Economics,

More information

Reflections in Hilbert Space III: Eigen-decomposition of Szegedy s operator

Reflections in Hilbert Space III: Eigen-decomposition of Szegedy s operator Reflections in Hilbert Space III: Eigen-decomposition of Szegedy s operator James Daniel Whitfield March 30, 01 By three methods we may learn wisdom: First, by reflection, which is the noblest; second,

More information

ECO Class 6 Nonparametric Econometrics

ECO Class 6 Nonparametric Econometrics ECO 523 - Class 6 Nonparametric Econometrics Carolina Caetano Contents 1 Nonparametric instrumental variable regression 1 2 Nonparametric Estimation of Average Treatment Effects 3 2.1 Asymptotic results................................

More information

Using SVD to Recommend Movies

Using SVD to Recommend Movies Michael Percy University of California, Santa Cruz Last update: December 12, 2009 Last update: December 12, 2009 1 / Outline 1 Introduction 2 Singular Value Decomposition 3 Experiments 4 Conclusion Last

More information

Dynamic Factor Models and Factor Augmented Vector Autoregressions. Lawrence J. Christiano

Dynamic Factor Models and Factor Augmented Vector Autoregressions. Lawrence J. Christiano Dynamic Factor Models and Factor Augmented Vector Autoregressions Lawrence J Christiano Dynamic Factor Models and Factor Augmented Vector Autoregressions Problem: the time series dimension of data is relatively

More information

Improving Judicial Ideal Point Estimates with a More Realistic Model of Opinion Content

Improving Judicial Ideal Point Estimates with a More Realistic Model of Opinion Content Improving Judicial Ideal Point Estimates with a More Realistic Model of Opinion Content Kevin M. Quinn Department of Government Harvard University kevin quinn@harvard.edu Jong Hee Park Department of Political

More information

Polar Decomposition of a Matrix

Polar Decomposition of a Matrix Polar Decomposition of a Matrix Garrett Buffington May 4, 2014 Table of Contents 1 The Polar Decomposition What is it? Square Root Matrix The Theorem Table of Contents 1 The Polar Decomposition What is

More information

COMS 4721: Machine Learning for Data Science Lecture 18, 4/4/2017

COMS 4721: Machine Learning for Data Science Lecture 18, 4/4/2017 COMS 4721: Machine Learning for Data Science Lecture 18, 4/4/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University TOPIC MODELING MODELS FOR TEXT DATA

More information

Learning and Global Dynamics

Learning and Global Dynamics Learning and Global Dynamics James Bullard 10 February 2007 Learning and global dynamics The paper for this lecture is Liquidity Traps, Learning and Stagnation, by George Evans, Eran Guse, and Seppo Honkapohja.

More information

CSI 445/660 Part 3 (Networks and their Surrounding Contexts)

CSI 445/660 Part 3 (Networks and their Surrounding Contexts) CSI 445/660 Part 3 (Networks and their Surrounding Contexts) Ref: Chapter 4 of [Easley & Kleinberg]. 3 1 / 33 External Factors ffecting Network Evolution Homophily: basic principle: We tend to be similar

More information

DSGE Methods. Estimation of DSGE models: Maximum Likelihood & Bayesian. Willi Mutschler, M.Sc.

DSGE Methods. Estimation of DSGE models: Maximum Likelihood & Bayesian. Willi Mutschler, M.Sc. DSGE Methods Estimation of DSGE models: Maximum Likelihood & Bayesian Willi Mutschler, M.Sc. Institute of Econometrics and Economic Statistics University of Münster willi.mutschler@uni-muenster.de Summer

More information

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Wednesday March 30 ± ǁ

. α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω. Wednesday March 30 ± ǁ . α β γ δ ε ζ η θ ι κ λ μ Aμ ν(x) ξ ο π ρ ς σ τ υ φ χ ψ ω. Α Β Γ Δ Ε Ζ Η Θ Ι Κ Λ Μ Ν Ξ Ο Π Ρ Σ Τ Υ Φ Χ Ψ Ω Wednesday March 30 ± ǁ 1 Chapter 5. Photons: Covariant Theory 5.1. The classical fields 5.2. Covariant

More information

Problem Set Solutions for the video course A Mathematics Course for Political and Social Research

Problem Set Solutions for the video course A Mathematics Course for Political and Social Research Problem Set Solutions for the video course A Mathematics Course for Political and Social Research David A Siegel July 0, 05 Associate Professor, Department of Political Science, Duke University, Durham,

More information

Multi-stage convex relaxation approach for low-rank structured PSD matrix recovery

Multi-stage convex relaxation approach for low-rank structured PSD matrix recovery Multi-stage convex relaxation approach for low-rank structured PSD matrix recovery Department of Mathematics & Risk Management Institute National University of Singapore (Based on a joint work with Shujun

More information

The University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013.

The University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013. The University of Texas at Austin Department of Electrical and Computer Engineering EE381V: Large Scale Learning Spring 2013 Assignment Two Caramanis/Sanghavi Due: Tuesday, Feb. 19, 2013. Computational

More information

CCP Estimation. Robert A. Miller. March Dynamic Discrete Choice. Miller (Dynamic Discrete Choice) cemmap 6 March / 27

CCP Estimation. Robert A. Miller. March Dynamic Discrete Choice. Miller (Dynamic Discrete Choice) cemmap 6 March / 27 CCP Estimation Robert A. Miller Dynamic Discrete Choice March 2018 Miller Dynamic Discrete Choice) cemmap 6 March 2018 1 / 27 Criteria for Evaluating Estimators General principles to apply when assessing

More information

Empirical Macroeconomics

Empirical Macroeconomics Empirical Macroeconomics Francesco Franco Nova SBE April 18, 2018 Francesco Franco Empirical Macroeconomics 1/23 Invertible Moving Average A difference equation interpretation Consider an invertible MA1)

More information

Robust Testing and Variable Selection for High-Dimensional Time Series

Robust Testing and Variable Selection for High-Dimensional Time Series Robust Testing and Variable Selection for High-Dimensional Time Series Ruey S. Tsay Booth School of Business, University of Chicago May, 2017 Ruey S. Tsay HTS 1 / 36 Outline 1 Focus on high-dimensional

More information

Community detection in time-dependent, multiscale, and multiplex networks

Community detection in time-dependent, multiscale, and multiplex networks Community detection in time-dependent, multiscale, and multiplex networks P. J. Mucha, T. Richardson, K. Macon, M. A. Porter, J.-P. Onnela Jukka-Pekka Onnela Harvard University MERSIH, November 13, 2009

More information

Linear Algebra Review

Linear Algebra Review January 29, 2013 Table of contents Metrics Metric Given a space X, then d : X X R + 0 and z in X if: d(x, y) = 0 is equivalent to x = y d(x, y) = d(y, x) d(x, y) d(x, z) + d(z, y) is a metric is for all

More information

4sec 2xtan 2x 1ii C3 Differentiation trig

4sec 2xtan 2x 1ii C3 Differentiation trig A Assignment beta Cover Sheet Name: Question Done Backpack Topic Comment Drill Consolidation i C3 Differentiation trig 4sec xtan x ii C3 Differentiation trig 6cot 3xcosec 3x iii C3 Differentiation trig

More information

MATH 581D FINAL EXAM Autumn December 12, 2016

MATH 581D FINAL EXAM Autumn December 12, 2016 MATH 58D FINAL EXAM Autumn 206 December 2, 206 NAME: SIGNATURE: Instructions: there are 6 problems on the final. Aim for solving 4 problems, but do as much as you can. Partial credit will be given on all

More information

Bargaining, Information Networks and Interstate

Bargaining, Information Networks and Interstate Bargaining, Information Networks and Interstate Conflict Erik Gartzke Oliver Westerwinter UC, San Diego Department of Political Sciene egartzke@ucsd.edu European University Institute Department of Political

More information

(Part 1) High-dimensional statistics May / 41

(Part 1) High-dimensional statistics May / 41 Theory for the Lasso Recall the linear model Y i = p j=1 β j X (j) i + ɛ i, i = 1,..., n, or, in matrix notation, Y = Xβ + ɛ, To simplify, we assume that the design X is fixed, and that ɛ is N (0, σ 2

More information

Polynomial Chaos and Karhunen-Loeve Expansion

Polynomial Chaos and Karhunen-Loeve Expansion Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the

More information

. D CR Nomenclature D 1

. D CR Nomenclature D 1 . D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the

More information

Network Inference From Grouped Observations. Using Hub Models. Supplementary Material

Network Inference From Grouped Observations. Using Hub Models. Supplementary Material Statistica Sinica: Supplement 1 Network Inference From Grouped Observations Using Hub Models YUNPENG ZHAO 1 AND CHARLES WEKO 2 George Mason University 1 and United States Army 2 Supplementary Material

More information

ERRATA. ± Lines are calculated before ( ) or after (+) the Anchor. If the Anchor is a page, t and b indicate, respectively, top and bottom.

ERRATA. ± Lines are calculated before ( ) or after (+) the Anchor. If the Anchor is a page, t and b indicate, respectively, top and bottom. ERRATA 22/02/2010 Fundamentals of Neutrino Physics and Astrophysics C. Giunti and C.. Kim Oxford University Press publication date: 15 March 2007; 728 pages ± Lines are calculated before or after + the

More information