Temporally-Reweighted Chinese Restaurant Process Mixtures for Clustering, Imputing, and Forecasting Multivariate Time Series

Total Page:16

File Type:pdf, Size:1020Kb

Temporally-Reweighted Chinese Restaurant Process Mixtures for Clustering, Imputing, and Forecasting Multivariate Time Series Temporally-Reweighted Chinese Restaurant Process Mixtures for Clustering, Imputing, and Forecasting Multivariate Time Series Feras A. Saad Vikash K. Mansinghka Probabilistic Computing Project Probabilistic Computing Project Massachusetts Institute of Technology Massachusetts Institute of Technology Abstract there are tens or hundreds of time series. One chal- lenge in these settings is that the data may reflect un- derlying processes with widely varying, non-stationary This article proposes a Bayesian nonparamet- dynamics [13]. Another challenge is that standard ric method for forecasting, imputation, and parametric approaches such as state-space models and clustering in sparsely observed, multivariate vector autoregression often become statistically and time series data. The method is appropri- numerically unstable in high dimensions [20]. Models ate for jointly modeling hundreds of time se- from these families further require users to perform ries with widely varying, non-stationary dy- significant custom modeling on a per-dataset basis, or namics. Given a collection of N time se- to search over a large set of possible parameter set- ries, the Bayesian model first partitions them tings and model configurations. In econometrics and into independent clusters using a Chinese finance, there is an increasing need for multivariate restaurant process prior. Within a clus- methods that exploit sparsity, are computationally ef- ter, all time series are modeled jointly us- ficient, and can accurately model hundreds of time se- ing a novel \temporally-reweighted" exten- ries (see introduction of [15], and references therein). sion of the Chinese restaurant process mix- ture. Markov chain Monte Carlo techniques This paper presents a nonparametric Bayesian method are used to obtain samples from the poste- for multivariate time series that aims to address some rior distribution, which are then used to form of the above challenges. The model is based on two ex- predictive inferences. We apply the tech- tensions to Dirichlet process mixtures. First, we intro- nique to challenging forecasting and imputa- duce a recurrent version of the Chinese restaurant pro- tion tasks using seasonal flu data from the US cess mixture to capture temporal dependences. Sec- Center for Disease Control and Prevention, ond, we add a hierarchical prior to discover groups of demonstrating superior forecasting accuracy time series whose underlying dynamics are modeled and competitive imputation accuracy as com- jointly. Unlike autoregressive models, our approach is pared to multiple widely used baselines. We designed to interpolate in regimes where it has seen further show that the model discovers inter- similar history before, and reverts to a broad prior in pretable clusters in datasets with hundreds of previously unseen regimes. This approach does not time series, using macroeconomic data from sacrifice predictive accuracy, when there is sufficient the Gapminder Foundation. signal to make a forecast or impute missing data. We apply the method to forecasting flu rates in 10 US regions using flu, weather, and Twitter data from the 1 Introduction US Center for Disease Control and Prevention. Quan- titative results show that the method outperforms sev- Multivariate time series data is ubiquitous, arising in eral Bayesian and non-Bayesian baselines, including domains such as macroeconomics, neuroscience, and Facebook Prophet, multi-output Gaussian processes, public health. Unfortunately, forecasting, imputation, seasonal ARIMA, and the HDP-HMM. We also show and clustering problems can be difficult to solve when competitive imputation accuracy with widely used sta- st tistical techniques. Finally, we apply the method to Proceedings of the 21 International Conference on Ar- clustering hundreds of macroeconomic time series from tificial Intelligence and Statistics (AISTATS) 2018, Lan- zarote, Spain. PMLR: Volume 84. Copyright 2018 by the Gapminder, detecting meaningful clusters of countries author(s). whose data exhibit coherent temporal patterns. Temporally-Reweighted Chinese Restaurant Process Mixtures for Multivariate Time Series 2 Related Work 3 Temporally-Reweighted Chinese Restaurant Process Mixture Model The temporally-reweighted Chinese restaurant process (TRCRP) mixture we introduce in Section3 can be We first outline the notations and basic setup as- directly seen as a time series extension to a family of sumed throughout this paper. Let xn : n = 1;:::;N nonparametric Bayesian regression models for cross- denote a collection of N discrete-timef series, whereg th n sectional data [18, 35, 28, 22, 23]. These methods the first T variables of the n time series is x1:T = n n n operate on an exchangeable data sequence xi with (x1 ; x2 ; : : : ; xT ). Slice notation is used to index sub- f g n n n exogenous covariates yi ; the prior CRP cluster prob- sequences of variables, so that x = (x ; : : : ; x ) f g t1:t2 t1 t2 ability p(zi = k) for each observation xi is reweighted for t1 < t2. Superscript n will be often be omit- based on yi. Our method extends this idea to a time ted when discussing a single time series. The re- series xt ; the prior CRP cluster probability p(zt = k) mainder of this section develops a generative pro- f g for xt is now reweighted based on the p previous val- cess for the joint distribution of all random variables n ues xt 1:t p. Moreover, the hierarchical extension in xt : t = 1; : : : ; N; n = 1;:::;N in the N time series, − − Section 3.4 coincides with CrossCat [21], when all tem- whichf we proceed to describe ing stages. poral dependencies are removed (by setting p = 0). Temporal extensions to the Dirichlet process have been 3.1 Background: CRP representation of previously used in the context of dynamic clustering Dirichlet process mixture models [38,1]. The latter work derives a recurrent CRP as the limit of a finite dynamic mixture model. Unlike Our approach is based on a temporal extension of the the method in this paper, those models are used for standard Dirichlet process mixture (DPM), which we clustering batched data and dynamic topic modeling review briefly. First consider the standard DPM in the [7], rather than data analysis tasks such as forecasting non-temporal setting [11], with concentration α and or imputation in real-valued, multivariate time series. base measure πΘ. The joint distribution of a sequence of m exchangeable random variables (x ; : : : ; x ) is: For multivariate time series, recent nonparametric 1 m Bayesian methods include using the dependent Dirich- P DP (α; πΘ); θ∗ P P; xj θ∗ F ( θ∗): let process for dynamic density estimation [31]; hier- ∼ j j ∼ j j ∼ · j j archical DP priors over the state in hidden Markov models [HDP-HMM; 12, 19]; Pitman-Yor mixtures of The DPM can be represented in terms of the Chi- non-linear state-space models for clustering [26]; and nese restaurant process [2]. As P is almost-surely dis- crete, the m draws θ∗ P contain repeated val- DP mixtures [8] and Polya trees [27] for modeling noise j ∼ distributions. As nonparametric Bayesian extensions ues, thereby inducing a clustering among data xj. Let of state-space models, all of these approaches specify λF be the hyperparameters of πΘ, θk be the unique f g priors that fall under distinct model classes to the one values among the θj∗ , and zj denote the cluster as- developed in this paper. They typically encode para- signment of xj which satisfies θj∗ = θzj . Define njk metric assumptions (such as linear autoregression and to be the number of observations xi with zi = k for hidden-state transition matrices), or integrate explicit i < j. Using the conditional distribution of zj given previous cluster assignments z1:j 1, the joint distribu- specifications of underlying temporal dynamics such − as seasonality, trends, and time-varying functionals. tion of exchangeable data sequence (x1; x2;::: ) in the Our method instead builds purely empirical models CRP mixture model can be described sequentially: and uses simple infinite mixtures to detect patterns iid in the data, without relying on dataset-specific cus- θk πΘ( λF ) tomizations. As a multivariate interpolator, the TR- f g ∼ · j CRP mixture is best applied to time series where there Pr [zj = k z1:j 1; α](j = 1; 2;::: ) is no structural theory of the temporal dynamics, and ( j − njk if 1 k max (z1:j 1) where there is sufficient statistical signal in the history ≤ ≤ − / α if k = max (z1:j 1) + 1 of the time series to inform probable future values. − (1) To the best of our knowledge, this paper presents the first multivariate, nonparametric Bayesian model that xj zj; θk F ( θzj ) j f g ∼ ·| provides strong baseline results without specifying cus- tom dynamics on a problem-specific basis; and that The CRP mixture model (1), and algorithms for pos- has been benchmarked against multiple Bayesian and terior inference, have been studied extensively for non- non-Bayesian techniques to cluster, impute, and fore- parametric modeling in a variety of statistical applica- cast sparsely observed real-world time series data. tions (for a survey see [37], and references therein). Saad and Mansinghka α; λG values are a-priori more likely to have similar distribu- tions for generating xt1 and xt2 , because G increases the probability that zt = zt = k , so that xt and f 1 2 g 1 λF xt2 are both drawn from F ( θk). z1 z2 z3 z4 ·| ··· Figure1 shows a graphical model for the TRCRP mix- ture (2) with window size p = 1. The model pro- ceeds as follows: first assume the initial p observations x x x x x θk 0 1 2 3 4 (x p+1; : : : ; x0) are fixed or have a known joint dis- ··· − k=1; 2;::: tribution. At step t, the generative process samples a cluster assignment zt, whose probability of joining Figure 1: Graphical model for the TRCRP mixture cluster k is a product of (i) the CRP probability for in a single time series x = (x ; x ;::: ) with lagged zt = k given all previous cluster assignments z1:t 1, 1 2 f g − and (ii) the \cohesion" term G(xt p:t 1; λG;Dtk).
Recommended publications
  • A Class of Measure-Valued Markov Chains and Bayesian Nonparametrics
    Bernoulli 18(3), 2012, 1002–1030 DOI: 10.3150/11-BEJ356 A class of measure-valued Markov chains and Bayesian nonparametrics STEFANO FAVARO1, ALESSANDRA GUGLIELMI2 and STEPHEN G. WALKER3 1Universit`adi Torino and Collegio Carlo Alberto, Dipartimento di Statistica e Matematica Ap- plicata, Corso Unione Sovietica 218/bis, 10134 Torino, Italy. E-mail: [email protected] 2Politecnico di Milano, Dipartimento di Matematica, P.zza Leonardo da Vinci 32, 20133 Milano, Italy. E-mail: [email protected] 3University of Kent, Institute of Mathematics, Statistics and Actuarial Science, Canterbury CT27NZ, UK. E-mail: [email protected] Measure-valued Markov chains have raised interest in Bayesian nonparametrics since the seminal paper by (Math. Proc. Cambridge Philos. Soc. 105 (1989) 579–585) where a Markov chain having the law of the Dirichlet process as unique invariant measure has been introduced. In the present paper, we propose and investigate a new class of measure-valued Markov chains defined via exchangeable sequences of random variables. Asymptotic properties for this new class are derived and applications related to Bayesian nonparametric mixture modeling, and to a generalization of the Markov chain proposed by (Math. Proc. Cambridge Philos. Soc. 105 (1989) 579–585), are discussed. These results and their applications highlight once again the interplay between Bayesian nonparametrics and the theory of measure-valued Markov chains. Keywords: Bayesian nonparametrics; Dirichlet process; exchangeable sequences; linear functionals
    [Show full text]
  • Examining Employment Relations in the Ethnic Chinese Restaurant Sector Within the UK Context
    Examining employment relations in the ethnic Chinese restaurant sector within the UK context By: Xisi Li A thesis submitted in partial fulfilment of the requirements for the degree of Doctor of Philosophy The University of Sheffield Faculty of Social Science Management School December 2017 Declaration No portion of the work referred to in the thesis has been submitted in support of an application for another degree or qualification of this or any other university or other institution of learning. Abstract Studies of employment relations in ethnic minority small firms have long been focused on the South Asian and Black communities. While the richness of these accounts has contributed much to our understanding of employment relations in small firms both relating to members of minority communities and more widely, there remains scope for engaging with a greater diversity of minority ethnic communities in the UK context. Specifically, there has not been any extensive research focusing on the ethnic Chinese community. The PhD thesis aims to examine employment relations in the ethnic Chinese restaurant sector within the UK context to address the current research gap. The research is located within a rich ethnographic tradition. The fieldwork for the current study consisted of participant observation of restaurant work over a period of seven months spent in two ethnic Chinese restaurants in Sheffield. The researcher worked as a full-time front area waiter and a full-time kitchen assistant. The field work enabled the researcher to develop a nuanced understanding of workplace behaviours. By focusing on four different aspects – the product market, the labour market, multi-cultural workforces and the informality-ethnicity interaction, this research thoroughly demonstrates how shop-floor behaviours and employment relations in the two case study firms are influenced by a variety of factors.
    [Show full text]
  • Markov Process Duality
    Markov Process Duality Jan M. Swart Luminy, October 23 and 24, 2013 Jan M. Swart Markov Process Duality Markov Chains S finite set. RS space of functions f : S R. ! S For probability kernel P = (P(x; y))x;y2S and f R define left and right multiplication as 2 X X Pf (x) := P(x; y)f (y) and fP(x) := f (y)P(y; x): y y (I do not distinguish row and column vectors.) Def Chain X = (Xk )k≥0 of S-valued r.v.'s is Markov chain with transition kernel P and state space S if S E f (Xk+1) (X0;:::; Xk ) = Pf (Xk ) a:s: (f R ) 2 P (X0;:::; Xk ) = (x0;:::; xk ) , = P[X0 = x0]P(x0; x1) P(xk−1; xk ): ··· µ µ µ Write P ; E for process with initial law µ = P [X0 ]. x δx x 2 · P := P with δx (y) := 1fx=yg. E similar. Jan M. Swart Markov Process Duality Markov Chains Set k µ k x µk := µP (x) = P [Xk = x] and fk := P f (x) = E [f (Xk )]: Then the forward and backward equations read µk+1 = µk P and fk+1 = Pfk : In particular π invariant law iff πP = π. Function h harmonic iff Ph = h h(Xk ) martingale. , Jan M. Swart Markov Process Duality Random mapping representation (Zk )k≥1 i.i.d. with common law ν, take values in (E; ). φ : S E S measurable E × ! P(x; y) = P[φ(x; Z1) = y]: Random mapping representation (E; ; ν; φ) always exists, highly non-unique.
    [Show full text]
  • 1 Introduction Branching Mechanism in a Superprocess from a Branching
    数理解析研究所講究録 1157 巻 2000 年 1-16 1 An example of random snakes by Le Gall and its applications 渡辺信三 Shinzo Watanabe, Kyoto University 1 Introduction The notion of random snakes has been introduced by Le Gall ([Le 1], [Le 2]) to construct a class of measure-valued branching processes, called superprocesses or continuous state branching processes ([Da], [Dy]). A main idea is to produce the branching mechanism in a superprocess from a branching tree embedded in excur- sions at each different level of a Brownian sample path. There is no clear notion of particles in a superprocess; it is something like a cloud or mist. Nevertheless, a random snake could provide us with a clear picture of historical or genealogical developments of”particles” in a superprocess. ” : In this note, we give a sample pathwise construction of a random snake in the case when the underlying Markov process is a Markov chain on a tree. A simplest case has been discussed in [War 1] and [Wat 2]. The construction can be reduced to this case locally and we need to consider a recurrence family of stochastic differential equations for reflecting Brownian motions with sticky boundaries. A special case has been already discussed by J. Warren [War 2] with an application to a coalescing stochastic flow of piece-wise linear transformations in connection with a non-white or non-Gaussian predictable noise in the sense of B. Tsirelson. 2 Brownian snakes Throughout this section, let $\xi=\{\xi(t), P_{x}\}$ be a Hunt Markov process on a locally compact separable metric space $S$ endowed with a metric $d_{S}(\cdot, *)$ .
    [Show full text]
  • The Globalization of Chinese Food ANTHROPOLOGY of ASIA SERIES Series Editor: Grant Evans, University Ofhong Kong
    The Globalization of Chinese Food ANTHROPOLOGY OF ASIA SERIES Series Editor: Grant Evans, University ofHong Kong Asia today is one ofthe most dynamic regions ofthe world. The previously predominant image of 'timeless peasants' has given way to the image of fast-paced business people, mass consumerism and high-rise urban conglomerations. Yet much discourse remains entrenched in the polarities of 'East vs. West', 'Tradition vs. Change'. This series hopes to provide a forum for anthropological studies which break with such polarities. It will publish titles dealing with cosmopolitanism, cultural identity, representa­ tions, arts and performance. The complexities of urban Asia, its elites, its political rituals, and its families will also be explored. Dangerous Blood, Refined Souls Death Rituals among the Chinese in Singapore Tong Chee Kiong Folk Art Potters ofJapan Beyond an Anthropology of Aesthetics Brian Moeran Hong Kong The Anthropology of a Chinese Metropolis Edited by Grant Evans and Maria Tam Anthropology and Colonialism in Asia and Oceania Jan van Bremen and Akitoshi Shimizu Japanese Bosses, Chinese Workers Power and Control in a Hong Kong Megastore WOng Heung wah The Legend ofthe Golden Boat Regulation, Trade and Traders in the Borderlands of Laos, Thailand, China and Burma Andrew walker Cultural Crisis and Social Memory Politics of the Past in the Thai World Edited by Shigeharu Tanabe and Charles R Keyes The Globalization of Chinese Food Edited by David Y. H. Wu and Sidney C. H. Cheung The Globalization of Chinese Food Edited by David Y. H. Wu and Sidney C. H. Cheung UNIVERSITY OF HAWAI'I PRESS HONOLULU Editorial Matter © 2002 David Y.
    [Show full text]
  • Immigration and Restaurants in Chicago During the Era of Chinese Exclusion, 1893-1933
    University of South Carolina Scholar Commons Theses and Dissertations Summer 2019 Exclusive Dining: Immigration and Restaurants in Chicago during the Era of Chinese Exclusion, 1893-1933 Samuel C. King Follow this and additional works at: https://scholarcommons.sc.edu/etd Recommended Citation King, S. C.(2019). Exclusive Dining: Immigration and Restaurants in Chicago during the Era of Chinese Exclusion, 1893-1933. (Doctoral dissertation). Retrieved from https://scholarcommons.sc.edu/etd/5418 This Open Access Dissertation is brought to you by Scholar Commons. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of Scholar Commons. For more information, please contact [email protected]. Exclusive Dining: Immigration and Restaurants in Chicago during the Era of Chinese Exclusion, 1893-1933 by Samuel C. King Bachelor of Arts New York University, 2012 Submitted in Partial Fulfillment of the Requirements For the Degree of Doctor of Philosophy in History College of Arts and Sciences University of South Carolina 2019 Accepted by: Lauren Sklaroff, Major Professor Mark Smith, Committee Member David S. Shields, Committee Member Erica J. Peters, Committee Member Yulian Wu, Committee Member Cheryl L. Addy, Vice Provost and Dean of the Graduate School Abstract The central aim of this project is to describe and explicate the process by which the status of Chinese restaurants in the United States underwent a dramatic and complete reversal in American consumer culture between the 1890s and the 1930s. In pursuit of this aim, this research demonstrates the connection that historically existed between restaurants, race, immigration, and foreign affairs during the Chinese Exclusion era.
    [Show full text]
  • Recruitment Costs for Chinese Workers in EU Member States Case
    An Initiative funded by the Partnership Instrument (PI) of the European Union EU-China Dialogue on Migration and Mobility Support Project 支持中欧人员往来和移民领域对话项目 Recruitment costs for Chinese workers in EU Member States Case study: Chinese chefs in Germany Recruitment costs for Chinese workers in EU Member States Case study: Chinese chefs in Germany Sabrina Kouba Copyright © 2018 International Labour Organization and International Organization for Migration First published 2018 Publications of the International Labour Office (ILO) enjoy copyright under Protocol 2 of the Universal Copyright Convention. Nevertheless, short excerpts from them may be reproduced without authorization, on condition that the source is indicated. For rights of adaptation, reproduction or translation, application should be made to ILO Publications (Rights and Licensing), International Labour Office, CH-1211 Geneva 22, Switzerland, or by email: [email protected]. The International Labour Office welcomes such applications. Libraries, institutions and other users registered with a reproduction rights organization may make copies in accordance with the licences issued to them for this purpose. Visit https://www.ifrro.org/ to find the reproduction rights organization in your country. ILO ISBN: 978-92-2-131678-7 (print); 978-92-2-131679-4 (web pdf) IOM ISBN: 978-92-9068-764-1 The designations employed in this publication, which are in conformity with United Nations practice, and the presentation of material therein do not imply the expression of any opinion whatsoever on the part of the ILO or the International Organization for Migration (IOM) concerning the legal status of any country, area, territory or city, or of its authorities, or concerning the delimitation of its frontiers or boundaries.
    [Show full text]
  • 1 Markov Chain Notation for a Continuous State Space
    The Metropolis-Hastings Algorithm 1 Markov Chain Notation for a Continuous State Space A sequence of random variables X0;X1;X2;:::, is a Markov chain on a continuous state space if... ... where it goes depends on where is is but not where is was. I'd really just prefer that you have the “flavor" here. The previous Markov property equation that we had P (Xn+1 = jjXn = i; Xn−1 = in−1;:::;X0 = i0) = P (Xn+1 = jjXn = i) is uninteresting now since both of those probabilities are zero when the state space is con- tinuous. It would be better to say that, for any A ⊆ S, where S is the state space, we have P (Xn+1 2 AjXn = i; Xn−1 = in−1;:::;X0 = i0) = P (Xn+1 2 AjXn = i): Note that it is okay to have equalities on the right side of the conditional line. Once these continuous random variables have been observed, they are fixed and nailed down to discrete values. 1.1 Transition Densities The continuous state analog of the one-step transition probability pij is the one-step tran- sition density. We will denote this as p(x; y): This is not the probability that the chain makes a move from state x to state y. Instead, it is a probability density function in y which describes a curve under which area represents probability. x can be thought of as a parameter of this density. For example, given a Markov chain is currently in state x, the next value y might be drawn from a normal distribution centered at x.
    [Show full text]
  • Simulation of Markov Chains
    Copyright c 2007 by Karl Sigman 1 Simulating Markov chains Many stochastic processes used for the modeling of financial assets and other systems in engi- neering are Markovian, and this makes it relatively easy to simulate from them. Here we present a brief introduction to the simulation of Markov chains. Our emphasis is on discrete-state chains both in discrete and continuous time, but some examples with a general state space will be discussed too. 1.1 Definition of a Markov chain We shall assume that the state space S of our Markov chain is S = ZZ= f:::; −2; −1; 0; 1; 2;:::g, the integers, or a proper subset of the integers. Typical examples are S = IN = f0; 1; 2 :::g, the non-negative integers, or S = f0; 1; 2 : : : ; ag, or S = {−b; : : : ; 0; 1; 2 : : : ; ag for some integers a; b > 0, in which case the state space is finite. Definition 1.1 A stochastic process fXn : n ≥ 0g is called a Markov chain if for all times n ≥ 0 and all states i0; : : : ; i; j 2 S, P (Xn+1 = jjXn = i; Xn−1 = in−1;:::;X0 = i0) = P (Xn+1 = jjXn = i) (1) = Pij: Pij denotes the probability that the chain, whenever in state i, moves next (one unit of time later) into state j, and is referred to as a one-step transition probability. The square matrix P = (Pij); i; j 2 S; is called the one-step transition matrix, and since when leaving state i the chain must move to one of the states j 2 S, each row sums to one (e.g., forms a probability distribution): For each i X Pij = 1: j2S We are assuming that the transition probabilities do not depend on the time n, and so, in particular, using n = 0 in (1) yields Pij = P (X1 = jjX0 = i): (Formally we are considering only time homogenous MC's meaning that their transition prob- abilities are time-homogenous (time stationary).) The defining property (1) can be described in words as the future is independent of the past given the present state.
    [Show full text]
  • Chapter 8: Markov Chains
    149 Chapter 8: Markov Chains 8.1 Introduction So far, we have examined several stochastic processes using transition diagrams and First-Step Analysis. The processes can be written as {X0, X1, X2,...}, where Xt is the state at time t. A.A.Markov On the transition diagram, Xt corresponds to 1856-1922 which box we are in at step t. In the Gambler’s Ruin (Section 2.7), Xt is the amount of money the gambler possesses after toss t. In the model for gene spread (Section 3.7), Xt is the number of animals possessing the harmful allele A in generation t. The processes that we have looked at via the transition diagram have a crucial property in common: Xt+1 depends only on Xt. It does not depend upon X0, X1,...,Xt−1. Processes like this are called Markov Chains. Example: Random Walk (see Chapter 4) none of these steps matter for time t+1 ? time t+1 time t ? In a Markov chain, the future depends only upon the present: NOT upon the past. 150 .. ...... .. .. ...... .. ....... ....... ....... ....... ... .... .... .... .. .. .. .. ... ... ... ... ....... ..... ..... ....... ........4................. ......7.................. The text-book image ...... ... ... ...... .... ... ... .... ... ... ... ... ... ... ... ... ... ... ... .... ... ... ... ... of a Markov chain has ... ... ... ... 1 ... .... ... ... ... ... ... ... ... .1.. 1... 1.. 3... ... ... ... ... ... ... ... ... 2 ... ... ... a flea hopping about at ... ..... ..... ... ...... ... ....... ...... ...... ....... ... ...... ......... ......... .. ........... ........ ......... ........... ........
    [Show full text]
  • Chop Suey As Imagined Authentic Chinese Food: the Culinary Identity of Chinese Restaurants in the United States
    UC Santa Barbara Journal of Transnational American Studies Title Chop Suey as Imagined Authentic Chinese Food: The Culinary Identity of Chinese Restaurants in the United States Permalink https://escholarship.org/uc/item/2bc4k55r Journal Journal of Transnational American Studies, 1(1) Author Liu, Haiming Publication Date 2009-02-16 DOI 10.5070/T811006946 Peer reviewed eScholarship.org Powered by the California Digital Library University of California Chop Suey as Imagined Authentic Chinese Food: The Culinary Identity of Chinese Restaurants in the United States HAIMING LIU Introduction In the small hours of one morning in 1917, John Doe, a white laborer, strolled into the Dragon Chop Suey House at 630 West Sixth Street, Los Angeles, and ordered chicken chop suey. The steaming bowl was set before Mr. Doe by a grinning Japanese. “I won’t eat it,” barked Mr. Doe, “There’s no poultry in it.” The flying squad was called in and was happily annoyed at this midnight incident. The officers offered to act as a jury and demanded sample bowls of chop suey. The Japanese owner declined and Mr. Doe was free to go.1 The laborer demanded real meat, the officers wanted free meals, and the owner of this Chinese restaurant was actually Japanese, but everyone was thoroughly familiar with the concept of chop suey. As this story shows, by 1917 chop suey was a well‐known restaurant meal in America. Food is a cultural tradition. The popularity of Chinese restaurants reflects how an Asian cuisine was transplanted and developed in American society. Chinese migration was a transnational flow of people, social networks, and cultural values.
    [Show full text]
  • Markov Random Fields and Stochastic Image Models
    Markov Random Fields and Stochastic Image Models Charles A. Bouman School of Electrical and Computer Engineering Purdue University Phone: (317) 494-0340 Fax: (317) 494-3358 email [email protected] Available from: http://dynamo.ecn.purdue.edu/»bouman/ Tutorial Presented at: 1995 IEEE International Conference on Image Processing 23-26 October 1995 Washington, D.C. Special thanks to: Ken Sauer Suhail Saquib Department of Electrical School of Electrical and Computer Engineering Engineering University of Notre Dame Purdue University 1 Overview of Topics 1. Introduction (b) Non-Gaussian MRF's 2. The Bayesian Approach i. Quadratic functions ii. Non-Convex functions 3. Discrete Models iii. Continuous MAP estimation (a) Markov Chains iv. Convex functions (b) Markov Random Fields (MRF) (c) Parameter Estimation (c) Simulation i. Estimation of σ (d) Parameter estimation ii. Estimation of T and p parameters 4. Application of MRF's to Segmentation 6. Application to Tomography (a) The Model (a) Tomographic system and data models (b) Bayesian Estimation (b) MAP Optimization (c) MAP Optimization (c) Parameter estimation (d) Parameter Estimation 7. Multiscale Stochastic Models (e) Other Approaches (a) Continuous models 5. Continuous Models (b) Discrete models (a) Gaussian Random Process Models 8. High Level Image Models i. Autoregressive (AR) models ii. Simultaneous AR (SAR) models iii. Gaussian MRF's iv. Generalization to 2-D 2 References in Statistical Image Modeling 1. Overview references [100, 89, 50, 54, 162, 4, 44] 4. Simulation and Stochastic Optimization Methods [118, 80, 129, 100, 68, 141, 61, 76, 62, 63] 2. Type of Random Field Model 5. Computational Methods used with MRF Models (a) Discrete Models i.
    [Show full text]