Independent Component Models for Replicated Point Processes

Total Page:16

File Type:pdf, Size:1020Kb

Independent Component Models for Replicated Point Processes Independent component models for replicated point processes Daniel Gervini September 25, 2018 Abstract We propose a semiparametric independent-component model for the intensity func- tions of a point process. When independent replications of the process are available, we show that the estimators are consistent and asymptotically normal. We study the finite-sample behavior of the estimators by simulation, and as an example of application we analyze the spatial distribution of street robberies in the city of Chicago. Key words: Cox process; latent-variable model; Poisson process; spline smoothing. 1 Introduction Point processes in time and space have a broad range of applications, in diverse areas such as neuroscience, ecology, finance, astronomy, seismology, and many others. Exam- ples are given in classic textbooks like Cox and Isham (1980), Diggle (2013), Møller and Waagepetersen (2004), Streit (2010), and Snyder and Miller (1991), and in the papers cited below. However, the point-process literature has mostly focused on the single-realization case, such as the distribution of trees in a single forest (Jalilian et al., 2013) or the dis- tribution of cells in a single tissue sample (Diggle et al., 2006). But situations where several replications of a process are available are increasingly common. Among the few papers proposing statistical methods for replicated point processes we can cite Diggle et arXiv:1506.00137v1 [stat.ME] 30 May 2015 al. (1991), Baddeley et al. (1993), Diggle et al. (2000), Bell and Grunwald (2004), Landau et al. (2004), Wager et al. (2004), and Pawlas (2011). However, these papers mostly propose estimators for summary statistics of the processes (the so-called F , G and K statistics) rather than for the intensity functions that characterize the processes. When several replications of a process are available, it is possible to estimate the inten- sity functions directly thanks to the possibility of “borrowing strength” across replications. This is the basic idea behind many functional data methods that are applied to stochastic processes which, individually, are only sparsely sampled (James et al., 2000; Yao et al., 2005; Gervini, 2009). Functional Data Analysis has mostly focused on continuous-time processes; little work has been done on discrete point processes. We can mention Bouzas 1 et al. (2006, 2007) and Fern´andez-Alcal´aet al. (2012), which have rather limited scopes since they only estimate the mean of Cox processes, and Wu et al. (2013), who estimate the mean and the principal components of independent and identically distributed realizations of a Cox process, but the method of Wu et al. (2013) is based on kernel estimators of the covariance function that are not easily extended beyond i.i.d. replications and time- dependent processes. The semiparametric method we propose in this paper, on the other hand, can be applied just as easily to temporal or spatial processes, and can be extended to non-i.i.d. situations like ANOVA models or more complex data structures like marked point processes and multivariate processes, even though we will not go as far in this paper. As an example of application we will analyze the spatial distribution of street robberies in Chicago in the year 2014, where each day is seen as a replication of the process. The method we propose in this paper fits a non-negative independent-component model for the intensity functions, and provides estimators of the intensity functions for the individual replications as by-products. We establish the consistency and asymptotic normality of the component estimators in Section 3, study their finite-sample behavior by simulation in Section 4, and analyze the Chicago crime data in Section 5. 2 Independent-component model for intensity process A point process X is a random countable set in a space , where is usually R for S S temporal processes and R2 or R3 for spatial processes (Møller and Waagepetersen, 2004, ch. 2; Streit, 2010, ch. 2). Locally finite processes are those for which #(X B) < ∩ ∞ with probability one for any bounded B . For such processes we can define the count ⊆ S function N(B) = #(X B). Given a locally integrable function λ : [0, ), i.e. a ∩ S → ∞ function such that λ(t)dt < for any bounded B , the process X is a Poisson B ∞ ⊆ S process with intensityR function λ(t) if (i) N(B) follows a Poisson distribution with rate λ(t)dt for any bounded B , and (ii) conditionally on N(B) = m, the m points in B ⊆ S X B are independent and identically distributed with density λ˜(t) = λ(t)/ λ for any R ∩ B bounded B . ⊆ S R Let X = X B for a given B. Then the density function of X at x = t ,...,t B ∩ B B { 1 m} is f(x ) = f(m)f(t ,...,t m) (1) B 1 m| λ(t)dt m m = exp λ(t)dt { B } λ˜(t ). − m! × j B R j=1 Z Y Since the realizations of XB are sets, not vectors, what we mean by ‘density’ is the following: 2 if is the family of locally finite subsets of , i.e. N S = A : #(A B) < for all bounded B , N { ⊆ S ∩ ∞ ⊆ S} then for any F ⊆N ∞ P (X F ) = I( t ,...,t F )f( t ,...,t )dt dt B ∈ · · · { 1 m}∈ { 1 m} 1 · · · m m=0 ZB ZB X∞ exp λ(t)dt m = − B I( t ,...,t F ) λ(t ) dt dt , m! · · · { 1 m}∈ { j } 1 · · · m m=0 R B B j=1 X Z Z Y and, more generally, for any function h : [0, ) N → ∞ ∞ E h(X ) = h( t ,...,t )f( t ,...,t )dt dt . (2) { B } · · · { 1 m} { 1 m} 1 · · · m m=0 B B X Z Z A function h on is, essentially, a function well defined on m for any integer m which is N S invariant under permutation of the coordinates; for example, h( t ,...,t )= m t /m. { 1 m} j=1 j Single realizations of point processes are often modeled as Poisson processesP with fixed λs. But for replicated point processes a single intensity function λ rarely provides an ade- quate fit for all replications. It is more reasonable to assume that the λs are subject-specific and treat them as latent random effects. These processes are called doubly stochastic or Cox processes (Møller and Waagepetersen, 2004, ch. 5; Streit, 2010, ch. 8). We can think of a doubly stochastic process as a pair (X, Λ) where X Λ = λ is a Poisson process | with intensity λ, and Λ is a random function that takes values on the space of non- F negative locally integrable functions on . This provides an appropriate framework for the S type of applications we have in mind: the n replications can be seen as i.i.d. realizations (X1, Λ1),..., (Xn, Λn) of a doubly stochastic process (X, Λ), where X is observable but Λ is not. In this paper we will assume that all X s are observed on a common region B of ; i S the method could be extended to Xis observed on non-conformal regions Bi at the expense of a higher computational complexity. For many applications we can think of the intensity process Λ as a finite combination of independent components, p Λ(t)= Ukφk(t), (3) Xk=1 where φ ,...,φ are functions in and U ,...,U are independent nonnegative random 1 p F 1 p variables. This is the functional equivalent of the multivariate Independent Component decomposition (Hyv¨arinen et al., 2001). For identifiability we assume that B φk(t)dt = 1 R 3 p for all k, i.e. that the φks are density functions. Then B Λ(t)dt = k=1 Uk and, if we p define S = k=1 Uk and Wk = Uk/S, we have R P P p ˜ ˜ Λ(t)= SΛ(t), with Λ(t)= Wkφk(t). (4) Xk=1 The intensity process Λ(t), then, can be seen as the product of an intensity factor S ˜ and a scatter factor Λ(t) which is a convex combination of the φks. Conditionally on Λ(t) = λ(t) the count function N(B) has a Poisson distribution with rate S = s, and conditionally on N(B) = m the m points in XB are independent identically distributed realizations with density λ˜(t) which is determined by the Wks. In general S and W are not independent, but if the Uks are assumed to follow independent Gamma distributions, say U Gamma(α , β) with a common β, then S and W are independent with S k ∼ k ∼ Gamma( α , β) and W Dirichlet(α ,...,α ) (Bilodeau and Brenner, 1999, ch. 3). k ∼ 1 p t The marginalP density of XB, with a slight abuse of notation in writing xB = (m, ), can be expressed as f(xB) = f(m, t, s, w) ds dw (5) Z Z = f(t m, w)f(m s)f(s, w) ds dw | | Z Z with m p f(t m, w)= w φ (t ). | k k j jY=1 Xk=1 ˜ Since λ(t) is a mixture of the φks, for computational convenience we introduce indicator ′ vectors y ,..., y such that, for each j, y = 1 for some k and y ′ = 0 for all k = k, 1 m jk jk 6 with P (Y = 1 m, w)= w . In words, y indicates which φ the point t is coming from. jk | k j k j Conditionally on m and w, the yjs are independent Multinomial(1, w). Then, collecting all yjs into a single y for notational simplicity, (5) can also be written as f(xB) = f(m, t, y, s, w) ds dw y Z Z X = f(t m, y)f(y m, w)f(m s)f(s, w) ds dw y | | | Z Z X with m p f (t m, y)= φ (t )yjk | k j jY=1 kY=1 4 and m p y f(y m, w)= w jk .
Recommended publications
  • Associations Between Air Temperature and Daily Suicide Counts in Astana, Kazakhstan
    Medicina (Kaunas) 2013;49(8):379-85 379 Associations Between Air Temperature and Daily Suicide Counts in Astana, Kazakhstan Andrej M. Grjibovski1, 2, Gulmira Kozhakhmetova3, Aliya Kosbayeva4, Bettina Menne5 1Department of International Public Health, Norwegian Institute of Public Health, Norway, 2International School of Public Health, Northern State Medical University, Russia, 3Medical University of Astana, Kazakhstan, 4WHO Office in Kazakhstan, Kazakhstan, 5WHO European Centre for Environment and Health, Bonn Office, Germany Key Words: suicides; temperature; apparent temperature; Central Asia; Kazakhstan. Summary. Background and Objective. Seasonal variations in suicide mortality and its associa- tion with ambient air temperature have been observed in many countries. However, the evidence from Central Asia is scarce. The aim of the study was to assess the relationship between 4 indicators of air temperature and daily suicide counts in Astana, Kazakhstan. Material and Methods. The daily counts of suicides (ICD-10 codes, X60–X84) for the popu- lation of Astana in 2005–2010 were collected using death certificates and medical records at the Municipal Bureau of Forensic Medicine. Associations between the number of cases and mean, maxi- mum, mean apparent, and maximum apparent temperatures were studied using negative binomial regression models controlling for the effects of month, year, weekends, holidays, wind velocity, baro- metric pressure, and relative humidity. Results. Altogether, there were 685 suicides in Astana in 2005–2010. A clear seasonal pattern with the peak in summer was observed. In crude analyses, significant associations between suicide counts and all 4 temperatures were found. After adjustment for other variables, only apparent tem- peratures remained significantly associated with the outcome.
    [Show full text]
  • Mapping Slavery in the Nineteenth Century
    Mapping Slavery in the Nineteenth Century Commercial lithographer Henry S. Graham printed this choropleth map showing the distribution of the slave population in September 1861. The map shows in graphic terms the density of the slave population in the Southern states, based on figures from the 1860 census. Although the development of this map was a collaborative government effort, cartographers working for Edwin Hergesheimer, U.S. Coast Survey Drafting Division, created it. The development of this map was revolutionary for its time for several reasons. First, it was among the first of its kind, initiating a trend of statistical cartography in the United States that allowed the thematic mapping of larger social, political, and cultural trends. Second, this map represented an early use of statistical information from the census. Third, new techniques in shading developed by Hergesheimer were a path-breaking application of these new techniques to human geography. Finally, its makers went as far to use “moral statistics” in order to affect political change. This map was more than a collaborative effort initiated by the federal government; it is a reflection of the beliefs of the cartographers and administrators who made it. The United States Coast Survey took the statistics from the 8th Census, supervised by Joseph Camp Griffith Kennedy, who was pro-Union and created controversy regarding the information to be gathered about slaves. Kennedy wanted to include _________________________ NOAA Office of Coast Survey www.nauticalcharts.noaa.gov/history/CivilWar slaves by name, but Congress refused. Alexander D. Bache, the Coast Survey Superintendent, was also pro-Union, and allowed his staff to undertake such a map.
    [Show full text]
  • The Scientific Rationality of Early Statistics, 1833–1877
    The Scientific Rationality of Early Statistics, 1833–1877 Yasuhiro Okazawa St Catharine’s College This dissertation is submitted for the degree of Doctor of Philosophy. November 2018 Declaration Declaration This dissertation is the result of my own work and includes nothing which is the outcome of work done in collaboration except as declared in the Preface and specified in the text. It is not substantially the same as any that I have submitted, or, is being concurrently submitted for a degree or diploma or other qualification at the University of Cambridge or any other University or similar institution except as declared in the Preface and specified in the text. I further state that no substantial part of my dissertation has already been submitted, or, is being concurrently submitted for any such degree, diploma or other qualification at the University of Cambridge or any other University or similar institution except as declared in the Preface and specified in the text It does not exceed the prescribed word limit of 80,000 words for the Degree Committee of the Faculty of History. Yasuhiro Okazawa 13 November 2018 i Thesis Summary The Scientific Rationality of Early Statistics, 1833–1877 Yasuhiro Okazawa Summary This thesis examines the activities of the Statistical Society of London (SSL) and its contribution to early statistics—conceived as the science of humans in society—in Britain. The SSL as a collective entity played a crucial role in the formation of early statistics, as statisticians envisaged early statistics as a collaborative scientific project and prompted large-scale observation, which required cooperation among numerous statistical observers.
    [Show full text]
  • An Examination of Cities' Resilence to Violent Crime: a Contextual Analysis of Victim-Offender Convergence
    The University of Southern Mississippi The Aquila Digital Community Dissertations Summer 8-2013 An Examination of Cities' Resilence to Violent Crime: A Contextual Analysis of Victim-Offender Convergence Vanessa Hatch Woodward University of Southern Mississippi Follow this and additional works at: https://aquila.usm.edu/dissertations Part of the Criminology and Criminal Justice Commons, Legal Writing and Research Commons, and the Other Legal Studies Commons Recommended Citation Woodward, Vanessa Hatch, "An Examination of Cities' Resilence to Violent Crime: A Contextual Analysis of Victim-Offender Convergence" (2013). Dissertations. 210. https://aquila.usm.edu/dissertations/210 This Dissertation is brought to you for free and open access by The Aquila Digital Community. It has been accepted for inclusion in Dissertations by an authorized administrator of The Aquila Digital Community. For more information, please contact [email protected]. The University of Southern Mississippi AN EXAMINATION OF CITIES’ RESILIENCE TO VIOLENT CRIME: A CONTEXTUAL ANALYSIS OF VICTIM-OFFENDER CONVERGENCE by Vanessa Hatch Woodward Abstract of a Dissertation Submitted to the Graduate School of The University of Southern Mississippi in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy August 2013 ABSTRACT AN EXAMINATION OF CITIES’ RESILIENCE TO VIOLENT CRIME: A CONTEXTUAL ANALYSIS OF VICTIM-OFFENDER CONVERGENCE by Vanessa Hatch Woodward August 2013 In his 2012 Presidential Address to the American Society of Criminology, Robert Sampson purported that causality can only be reached when social science researchers accept that individual actions are dependent on social context. He referred to this as contextual causality and argued that future research needed to focus on how to measure and/or reoperationalize community measures of crime.
    [Show full text]
  • Spatial Multivariate Analysis of Guerry's Data in R
    Spatial multivariate analysis of Guerry's data in St´ephane Dray Universit´ede Lyon ; universit´eLyon 1 ; CNRS UMR 5558 Laboratoire de Biom´etrie et Biologie Evolutive 43 boulevard du 11 novembre 1918 Villeurbanne F-69622, France [email protected] http://pbil.univ-lyon1.fr/members/dray December 1, 2009 Abstract This vignette indicates how to perform the analyses described in Dray and Jombart (submitted) of data in the Guerry package, derived from Andr´e-Michel Guerry's (1833) Essai sur la Statistique Morale de la France, using . It demonstrates some classical methods for analysis of multi- variate spatial data that focus either on the multivariate aspect or on the spatial one, as well as some more modern methods that attempt to integrate geographical and multivariate aspects simultaneously. Contents 1 Introduction 3 1.1 Installation and loading of required packages . .3 1.2 Preliminary steps . .3 2 Standard approaches 3 2.1 Multivariate analysis . .4 2.2 Spatial autocorrelation . .7 2.2.1 The spatial weighting matrix . .7 2.2.2 Moran's Coefficient . .8 2.2.3 Moran scatterplot . 10 2.3 Toward an integration of multivariate and geographical aspects . 11 3 Spatial multivariate analysis 12 3.1 Spatial partition . 12 3.2 Spatial explanatory variables . 16 1 3.2.1 Trend surface of geographic coordinates . 16 3.2.2 Moran's eigenvector maps . 18 3.3 Spatial graph and weighting matrix . 20 4 Conclusions 24 2 1 Introduction A recent study (Friendly, 2007) revived Andr´e-Michel Guerry's (1833) Essai sur la Statistique Morale de la France.
    [Show full text]
  • Independent Component Models for Replicated Point Processes
    Independent component models for replicated point processes Daniel Gervini Department of Mathematical Sciences University of Wisconsin–Milwaukee August 9, 2016 Abstract We propose a semiparametric independent-component model for the intensity func- tions of replicated point processes. We show that the maximum likelihood estimators of the model parameters are consistent and asymptotically normal when the number of replications goes to in…nity. The …nite-sample behavior of the estimators is stud- ied by simulation. As an example of application, we analyze the temporal variation in the spatial distribution of street robberies in Chicago. Key words: Doubly-stochastic process; functional data analysis; latent-variable model; Poisson process; spline smoothing. 1 Introduction Point processes in time and space have a broad range of applications, in diverse areas such as neuroscience, ecology, …nance, astronomy, seismology, and many others. Examples are given in classic textbooks like Cox and Isham (1980), Diggle (2013), Møller and Waagepetersen (2004), Streit (2010), and Snyder and Miller (1991), and in the papers cited below. However, the point-process literature has mostly focused on single-realization cases, such as the distribution of trees in a single forest (Jalilian et al., 2013) or the distribution of cells in a single tissue sample (Diggle et al., 2006). Situations where several replications of a process are available are increasingly common. Among the few papers proposing statistical methods for replicated point processes we can cite Diggle et al. (1991), Baddeley et al. (1993), Diggle et al. (2000), Bell and Grunwald (2004), Landau et al. (2004), Wager et al. (2004), and Pawlas (2011). However, these papers propose estimators for summary statistics of the processes (the so-called F , G and K statistics) rather than the intensity functions, which would be more informative.
    [Show full text]
  • Statistical Analysis of Spatial Crime Data
    CHAPTER 33 Statistical Analysis of Spatial Crime Data WIM BERNASCO AND HENK ELFFERS INTRODUCTION Social scientists become increasingly aware of the relevance of space and place (Goodchild et al. 2000), and criminology is no exception. As a matter of fact, the geography of crime has been a focal concern of criminologists from the very start of the discipline. In nine- teenth century Europe, the “moral statistics” of pioneers, Guerry and Quetelet, empirically demonstrated that crime varied across geographical regions. They not only produced maps that visualized these differences, but also studied statistical relations between crime, poverty, and education. During the early twentieth century, researchers associated with the University of Chicago studied how crime and other social problems varied across urban communities, again mapping the geographical patterns and using community characteristics to explain these distributions. The work on juvenile delinquency has become a classic example (Shaw and McKay 1942). Interest in the geography of crime decreased somewhat between 1950 and 1980, to some extend possibly because the ecological fallacy (Robinson 1950) hampered the interpretation of aggregated data. Later, the link between communities and crime was revitalized (Bursik and Grasmick 1993; Sampson et al. 1997). In recent years, a new concern has emerged with micro units of place such as addresses or street segments (Eck and Weisburd 1995; St. Jean 2007). From the very start, geographical criminology has been an area of research where methodological and statistical innovations were either developed or adopted early. In their times, Guerry and Quetelet were pioneers and innovators, and their work is said to have been the launching pad for much of modern social science (Beirne 1987; Friendly 2007).
    [Show full text]
  • A Local Indicator of Multivariate Spatial Association
    A Local Indicator of Multivariate Spatial Association: Extending Geary’s c.∗ Luc Anselin Center for Spatial Data Science University of Chicago [email protected] November 9, 2017 ∗This research was funded in part by Award 1R01HS021752-01A1 from the Agency for Healthcare Research and Quality (AHRQ), “Advancing spatial evaluation methods to improve healthcare efficiency and quality.” Comments by Sergio Rey, Julia Koschinsky and the referees on an earlier draft are greatly appreciated. Abstract This paper extends the application of the Local Geary c statistic to a multivariate context. The statistic is conceptualized as a weighted distance in multivariate attribute space between an observation and its geographical neighbors. Inference is based on a conditional permutation approach. The interpretation of significant univariate Local Geary statistics is clarified and the differences with a multivariate case outlined. An empirical illustration uses Guerry’s classic data on moral statistics in 1830s France. Key words: spatial clusters, LISA, multivariate spatial association, Local Geary c, spatial data science 1 Introduction An important and growing component of geographical analysis is a focus on the local, reflected in new methods to deal with local spatial autocorrelation and local spatial heterogeneity (for general overviews, see, e.g., Unwin 1996, Fotheringham 1997, Unwin and Unwin 1998, Fother- ingham and Brunsdon 1999, Boots 2002, Boots and Okabe 2007, Lloyd 2010). Specifically, considerable interest has been devoted to local indicators of spatial association (LISA) since the original LISA framework was outlined in Anselin (1995, 1996), building upon the initial work by Getis and Ord (1992, 1996), and Ord and Getis (1995). Since the local statistics form the basis for a hypothesis test for (local) spatial randomness, they are strictly speaking outside the scope of exploratory data analysis (EDA) as outlined by, among others, Tukey (1977) and Good (1983), and its spatial counterparts (ESDA).
    [Show full text]
  • On the Nineteenth-Century Origins of Significance Testing and P-Hacking
    On the nineteenth-century origins of significance testing and p-hacking Glenn Shafer, Rutgers University The Game-Theoretic Probability and Finance Project Working Paper #55 First posted July 18, 2019. Last revised June 11, 2020. Project web site: http://www.probabilityandfinance.com Abstract Although the names significance test, p-value, and confidence interval came into use only in the 20th century, the methods they name were already used and abused in the 19th century. Knowledge of this earlier history can help us evaluate some of the ideas for improving statistical testing and estimation currently being discussed. This article recounts first the development of statistical testing and estima- tion after Laplace's discovery of the central limit theorem and then the sub- sequent transmission of these ideas into the English-language culture of math- ematical statistics in the early 20th century. I argue that the earlier history casts doubt on the efficacy of many of the competing proposals for improving on significance tests and p-values and for forestalling abuses. Rather than fur- ther complicate the way we now teach statistics, we should leave aside most of the 20th-century embellishments and emphasize exploratory data analysis and the idea of testing probabilities by betting against them. 1 Introduction 1 2 Laplace's theorem 3 2.1 Laplace's discovery . .3 2.2 Direct and inverse probability . .4 2.3 Laplacean and Gaussian least squares . .5 2.4 Seeing p-hacking in France . .7 2.5 The disappearance of Laplace's theorem in France . .8 3 Practical certainty 8 3.1 La limite de l'´ecart ..........................9 3.2 Tables of the normal distribution .
    [Show full text]
  • Report No. 30 Crime Statistics: User Perspectives Report No
    Crime Statistics: User Perspectives Statistics Commission Report No. 30 Crime Statistics: User Perspectives Report No. 30 Statistics Commission Report No. 30 Crime Statistics: User Perspectives September 2006 Report by the Statistics Commission Incorporating: Review of Crime Statistics by Matrix Research and Consultancy and Professor Tim Hope Statistics Commission Artillery House 11-19 Artillery Row London SW1P 1RT 020 7273 8008 www.statscom.org.uk © Crown Copyright 2006 Contents Page Part 1: Crime Statistics – User Perspectives 3 A report by the Statistics Commission Section 1: Recommendations 5 Section 2: Background to the report 7 Section 3: Rationale for recommendations 11 Part 2: Review of Crime Statistics 20 Report to the Commission by Matrix and Professor Tim Hope Executive Summary 25 Introduction 35 Approach 37 Context 39 Groups with an interest in crime statistics 49 Purposes for which crime statistics are used 55 Fitness for Purpose 63 Potential solutions 83 Findings and recommendations 95 Appendices 99 iii iv Preface By the Chairman of the Statistics Commission To produce a report on crime statistics is to travel a crowded road, one that may be in danger of intellectual grid-lock as ever more views from experts, commentators and the media run up against one another. Thus our decision to produce this report at this time may require some explanation, and perhaps justification. The Statistics Commission exists to help ensure official statistics are trustworthy and responsive to public needs. Crime statistics are the quintessential ‘official figures’, a measure both of society and of government, telling us something about the social hazards we face and something about the success of government and public services in containing those hazards.
    [Show full text]
  • Downloaded from Manchesterhive.Com at 09/26/2021 07:14:08AM Via Free Access
    2 All the world’s a stage: Paris 1855 n 1855 Parisians believed that their city was the centre of the world. On 15 IMay of that year emperor Napoleon III opened the second World’s Fair, which would attract over five million visitors. To Napoleon, this exposition was the international affirmation of his reign. Charles Louis Napoleon Bonaparte, the third son of Louis Napoleon Bonaparte, king of Holland, was elected presi- dent of the new French Republic in 1848. In 1852, he abandoned his republican ideals and had himself crowned emperor Napoleon III of France. He ordered the construction of the Palais de l’Industrie, a magnificent structure of glass and iron between the Champs-Élysées and the river Seine. The design was remi- niscent of the Crystal Palace in London, which had been built for the Great Exhibition of 1851. The colossal Palais provided accommodation for the 21,779 industrial exhibitors. At 250 metres long, 108 metres wide and 35 metres high, it was one of the largest, if not the most elegant, modern structures in Paris according to the Baedeker of 1878. The exposition’s 2,175 fine art entries were housed in a separate building, the Palais des Beaux-Arts. The construction of the palaces was an integral part of the grandiose urban renewal project that the emperor asked prefect Georges Haussmann to carry out in 1853. The visitors to the World’s Fair witnessed the beginning of a demolition and construction craze that would grip the city for years to come. Napoleon also employed less peaceful means to raise the prestige of his empire.
    [Show full text]
  • Lecture 14: Classification, Statistical Sins
    Lecture 14: Classification, Statistical Sins 6.0002 LECTURE 14 1 Announcements §Reading ◦ Chapter 21 §Course evaluations ◦ Online evaluation now through noon on Friday, December 16 §Will be making study code for final exam available later today 6.0002 LECTURE 14 2 Compare to KNN Results (from Monday) Average of 10 80/20 splits using KNN (k=3) Average of 10 80/20 splits LR Accuracy = 0.744 Accuracy = 0.804 Sensitivity = 0.629 Sensitivity = 0.719 Specificity = 0.829 Specificity = 0.859 Pos. Pred. Val. = 0.728 Pos. Pred. Val. = 0.767 Average of LOO testing using KNN (k=3) Average of LOO testing using LR Accuracy = 0.769 Accuracy = 0.786 Sensitivity = 0.663 Sensitivity = 0.705 Specificity = 0.842 Specificity = 0.842 Pos. Pred. Val. = 0.743 Pos. Pred. Val. = 0.754 Performance not much difference Logistic regression slightly better Logistic regression provides ins ight about variables 6.0002 LECTURE 14 3 Looking at F eature Weights model.classes_ = ['Died' 'Survived'] For label Survived Be wary of reading too C1 = 1.66761946545 much into the weights C2 = 0.460354552452 Features are often C3 = -0.50338282535 correlated age = -0.0314481062387 male gender = -2.39514860929 L1 regression tends to drive one variable to zero L2 (default) regression spreads weights across variables 6.0002 LECTURE 14 4 Correlated Features, a n Example §c1 + c2 + c3 = 1 ◦ I.e., values are not independent st ◦ Is being in 1 class good, or being in the other classes bad? §Suppose we eliminate c1? 6.0002 LECTURE 14 5 Comparative Results Original Features Modified Features Average of 20 80/20 split s LR Average of 20 80 /20 split s LR Accuracy = 0.778 Accuracy = 0.779 Sensitivity = 0.687 Sensitivity = 0.674 Specificity = 0.842 Specificity = 0.853 Pos.
    [Show full text]