Spatial Autocorrelation
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Higher-Order Asymptotics
Higher-Order Asymptotics Todd Kuffner Washington University in St. Louis WHOA-PSI 2016 1 / 113 First- and Higher-Order Asymptotics Classical Asymptotics in Statistics: available sample size n ! 1 First-Order Asymptotic Theory: asymptotic statements that are correct to order O(n−1=2) Higher-Order Asymptotics: refinements to first-order results 1st order 2nd order 3rd order kth order error O(n−1=2) O(n−1) O(n−3=2) O(n−k=2) or or or or o(1) o(n−1=2) o(n−1) o(n−(k−1)=2) Why would anyone care? deeper understanding more accurate inference compare different approaches (which agree to first order) 2 / 113 Points of Emphasis Convergence pointwise or uniform? Error absolute or relative? Deviation region moderate or large? 3 / 113 Common Goals Refinements for better small-sample performance Example Edgeworth expansion (absolute error) Example Barndorff-Nielsen’s R∗ Accurate Approximation Example saddlepoint methods (relative error) Example Laplace approximation Comparative Asymptotics Example probability matching priors Example conditional vs. unconditional frequentist inference Example comparing analytic and bootstrap procedures Deeper Understanding Example sources of inaccuracy in first-order theory Example nuisance parameter effects 4 / 113 Is this relevant for high-dimensional statistical models? The Classical asymptotic regime is when the parameter dimension p is fixed and the available sample size n ! 1. What if p < n or p is close to n? 1. Find a meaningful non-asymptotic analysis of the statistical procedure which works for any n or p (concentration inequalities) 2. Allow both n ! 1 and p ! 1. 5 / 113 Some First-Order Theory Univariate (classical) CLT: Assume X1;X2;::: are i.i.d. -
Geomathematics in Hungarian Geology
JOURNAL OF HUNGARIAN GEOMATHEMATICS Volume 1 Geomathematics in Hungarian Geology George Bárdossy1 1Hungarian Academy of Sciences ([email protected]) Abstract The application of mathematical methods has a long tradition in Hungary. The main bases of geomathematics are the universities of the country, more closely the departments related to geology, such as general geology, stratigraphy, paleontology, structural geology, mineralogy, pet-rography, geochemistry, hydrogeology and applied geology. The Hungarian Geological Survey, the Geological Institute of Hungary and the Geochemical Research Laboratory of the Hungarian Academy of Sciences are institutions where geomathematical methods found broad applications. Finally, some mining and exploration companies, like the Hungarian Oil Company (MOL), the Bakony Bauxite Mining Company and others are regularly using geomathematical methods, mainly for the evaluation of exploration results, for deposit and reservoir modelling and for the estimation of resources. Keywords: geomathematics, applications. The application of mathematical methods has a long tradition in Hungary. The main bases of geomathematics are the universities of the country, more closely the departments related to geology, such as general geology, stratigraphy, paleontology, structural geology, mineralogy, pet-rography, geochemistry, hydrogeology and applied geology. The Hungarian Geological Survey, the Geological Institute of Hungary and the Geochemical Research Laboratory of the Hungarian Academy of Sciences are institutions where geomathematical methods found broad applications. Finally, some mining and exploration companies, like the Hungarian Oil Company (MOL), the Bakony Bauxite Mining Company and others are regularly using geomathematical methods, mainly for the evaluation of exploration results, for deposit and reservoir modelling and for the estimation of resources. In the following examples of recent geomathematical applications are listed. -
Applied Computing and Geosciences
APPLIED COMPUTING AND GEOSCIENCES AUTHOR INFORMATION PACK TABLE OF CONTENTS XXX . • Description p.1 • Editorial Board p.2 • Guide for Authors p.3 ISSN: 2590-1974 DESCRIPTION . Applied Computing & Geosciences is an online-only, open access journal focused on all aspects of computing in the geosciences. Like its companion title Computers & Geosciences, Applied Computing & Geosciences' mission is to advance and disseminate knowledge in all the related areas of at the interface between computer sciences and geosciences. Applied Computing & Geosciences publishes original articles, review articles and case-studies. Alongside welcoming direct submissions, the journal will benefit from an Article Transfer Service which will allow the author(s) to transfer their manuscript online from Computers & Geosciences thus saving authors time and effort spent on formatting and resubmitting. Applied Computing & Geosciences offers the community an innovative, efficient and flexible route for the publication of scientifically and ethically sound articles which address problems in the geosciences which includes: Algorithms; Artificial Intelligence; Computational Methods; Computer Graphics; Computer Visualization; Data Models; Data Processing; Database Retrieval; Distributed Systems; E-Geoscience; Geocomputation; Geographical Information Systems; Geoinformatics; Geomathematics; Image Analysis; Information Retrieval; Modelling; Near and Remote Sensing Data Analysis; Ontologies; Parallel Systems; Programming Languages; Remote Sensing; Simulation; Social Media; -
2006 Executive Committee CSPG 2005 Strategic Planning Session Geoscience Mixer 2005 William (Bill) Carruthers Gu
December 11/11/05 6:16 PM Page 1 Canadian Publication Mail Contract - 40070050 $3.00 VOLUME 32, ISSUE 11 DECEMBER 2005 ■ 2006 Executive Committee ■ CSPG 2005 Strategic Planning Session ■ Geoscience Mixer 2005 ■ William (Bill) Carruthers Gussow (1908-2005) ■ 2006 CSPG CSEG CWLS Joint Convention December 11/14/05 9:19 PM Page 2 December 11/11/05 6:16 PM Page 3 CSPG OFFICE #160, 540 - 5th Avenue SW Calgary,Alberta, Canada T2P 0M2 Tel:403-264-5610 Fax: 403-264-5898 Web: www.cspg.org Office hours: Monday to Friday, 8:30am to 4:00pm CONTENTS Business Manager:Tim Howard Email: [email protected] Office Manager: Deanna Watkins Email: [email protected] Communications Manager: Jaimè Croft Larsen Email: [email protected] Conventions Manager: Lori Humphrey-Clements ARTICLES Email: [email protected] Corporate Relations Manager: Kim MacLean 2006 EXECUTIVE COMMITTEE . 28 Email: [email protected] CSPG 2005 STRATEGIC PLANNING SESSION . 33 EDITORS/AUTHORS Please submit RESERVOIR articles to the CSPG GEOSCIENCE MIXER 2005 . 35 office. Submission deadline is the 23rd day of the month, two months prior to issue date. GOOGLING GEOMORPHOLOGY . 36 (e.g., January 23 for the March issue). To publish an article, the CSPG requires digital WILLIAM (BILL) CARRUTHERS GUSSOW (1908 - 2005) . 38 copies of the document. Text should be in Microsoft Word format and illustrations should 2006 CSPG CSEG CWLS JOINT CONVENTION . 43 be in TIFF format at 300 dpi. For additional information on manuscript preparation, refer to the Guidelines for Authors published in the CSPG Bulletin or contact the editor. -
Seismic Wavefield Imaging of Earth's Interior Across Scales
TECHNICAL REVIEWS Seismic wavefield imaging of Earth’s interior across scales Jeroen Tromp Abstract | Seismic full- waveform inversion (FWI) for imaging Earth’s interior was introduced in the late 1970s. Its ultimate goal is to use all of the information in a seismogram to understand the structure and dynamics of Earth, such as hydrocarbon reservoirs, the nature of hotspots and the forces behind plate motions and earthquakes. Thanks to developments in high- performance computing and advances in modern numerical methods in the past 10 years, 3D FWI has become feasible for a wide range of applications and is currently used across nine orders of magnitude in frequency and wavelength. A typical FWI workflow includes selecting seismic sources and a starting model, conducting forward simulations, calculating and evaluating the misfit, and optimizing the simulated model until the observed and modelled seismograms converge on a single model. This method has revealed Pleistocene ice scrapes beneath a gas cloud in the Valhall oil field, overthrusted Iberian crust in the western Pyrenees mountains, deep slabs in subduction zones throughout the world and the shape of the African superplume. The increased use of multi- parameter inversions, improved computational and algorithmic efficiency , and the inclusion of Bayesian statistics in the optimization process all stand to substantially improve FWI, overcoming current computational or data- quality constraints. In this Technical Review, FWI methods and applications in controlled- source and earthquake seismology are discussed, followed by a perspective on the future of FWI, which will ultimately result in increased insight into the physics and chemistry of Earth’s interior. -
The Method of Maximum Likelihood for Simple Linear Regression
08:48 Saturday 19th September, 2015 See updates and corrections at http://www.stat.cmu.edu/~cshalizi/mreg/ Lecture 6: The Method of Maximum Likelihood for Simple Linear Regression 36-401, Fall 2015, Section B 17 September 2015 1 Recapitulation We introduced the method of maximum likelihood for simple linear regression in the notes for two lectures ago. Let's review. We start with the statistical model, which is the Gaussian-noise simple linear regression model, defined as follows: 1. The distribution of X is arbitrary (and perhaps X is even non-random). 2. If X = x, then Y = β0 + β1x + , for some constants (\coefficients", \parameters") β0 and β1, and some random noise variable . 3. ∼ N(0; σ2), and is independent of X. 4. is independent across observations. A consequence of these assumptions is that the response variable Y is indepen- dent across observations, conditional on the predictor X, i.e., Y1 and Y2 are independent given X1 and X2 (Exercise 1). As you'll recall, this is a special case of the simple linear regression model: the first two assumptions are the same, but we are now assuming much more about the noise variable : it's not just mean zero with constant variance, but it has a particular distribution (Gaussian), and everything we said was uncorrelated before we now strengthen to independence1. Because of these stronger assumptions, the model tells us the conditional pdf 2 of Y for each x, p(yjX = x; β0; β1; σ ). (This notation separates the random variables from the parameters.) Given any data set (x1; y1); (x2; y2);::: (xn; yn), we can now write down the probability density, under the model, of seeing that data: n n (y −(β +β x ))2 Y 2 Y 1 − i 0 1 i p(yijxi; β0; β1; σ ) = p e 2σ2 2 i=1 i=1 2πσ 1See the notes for lecture 1 for a reminder, with an explicit example, of how uncorrelated random variables can nonetheless be strongly statistically dependent. -
Use of the Kurtosis Statistic in the Frequency Domain As an Aid In
lEEE JOURNALlEEE OF OCEANICENGINEERING, VOL. OE-9, NO. 2, APRIL 1984 85 Use of the Kurtosis Statistic in the FrequencyDomain as an Aid in Detecting Random Signals Absmact-Power spectral density estimation is often employed as a couldbe utilized in signal processing. The objective ofthis method for signal ,detection. For signals which occur randomly, a paper is to compare the PSD technique for signal processing frequency domain kurtosis estimate supplements the power spectral witha new methodwhich computes the frequency domain density estimate and, in some cases, can be.employed to detect their presence. This has been verified from experiments vith real data of kurtosis (FDK) [2] forthe real and imaginary parts of the randomly occurring signals. In order to better understand the detec- complex frequency components. Kurtosis is defined as a ratio tion of randomlyoccurring signals, sinusoidal and narrow-band of a fourth-order central moment to the square of a second- Gaussian signals are considered, which when modeled to represent a order central moment. fading or multipath environment, are received as nowGaussian in Using theNeyman-Pearson theory in thetime domain, terms of a frequency domain kurtosis estimate. Several fading and multipath propagation probability density distributions of practical Ferguson [3] , has shown that kurtosis is a locally optimum interestare considered, including Rayleigh and log-normal. The detectionstatistic under certain conditions. The reader is model is generalized to handle transient and frequency modulated referred to Ferguson'swork for the details; however, it can signals by taking into account the probability of the signal being in a be simply said thatit is concernedwith detecting outliers specific frequency range over the total data interval. -
Statistical Models in R Some Examples
Statistical Models Statistical Models in R Some Examples Steven Buechler Department of Mathematics 276B Hurley Hall; 1-6233 Fall, 2007 Statistical Models Outline Statistical Models Linear Models in R Statistical Models Regression Regression analysis is the appropriate statistical method when the response variable and all explanatory variables are continuous. Here, we only discuss linear regression, the simplest and most common form. Remember that a statistical model attempts to approximate the response variable Y as a mathematical function of the explanatory variables X1;:::; Xn. This mathematical function may involve parameters. Regression analysis attempts to use sample data find the parameters that produce the best model Statistical Models Linear Models The simplest such model is a linear model with a unique explanatory variable, which takes the following form. y^ = a + bx: Here, y is the response variable vector, x the explanatory variable, y^ is the vector of fitted values and a (intercept) and b (slope) are real numbers. Plotting y versus x, this model represents a line through the points. For a given index i,y ^i = a + bxi approximates yi . Regression amounts to finding a and b that gives the best fit. Statistical Models Linear Model with 1 Explanatory Variable ● 10 ● ● ● ● ● 5 y ● ● ● ● ● ● y ● ● y−hat ● ● ● ● 0 ● ● ● ● x=2 0 1 2 3 4 5 x Statistical Models Plotting Commands for the record The plot was generated with test data xR, yR with: > plot(xR, yR, xlab = "x", ylab = "y") > abline(v = 2, lty = 2) > abline(a = -2, b = 2, col = "blue") > points(c(2), yR[9], pch = 16, col = "red") > points(c(2), c(2), pch = 16, col = "red") > text(2.5, -4, "x=2", cex = 1.5) > text(1.8, 3.9, "y", cex = 1.5) > text(2.5, 1.9, "y-hat", cex = 1.5) Statistical Models Linear Regression = Minimize RSS Least Squares Fit In linear regression the best fit is found by minimizing n n X 2 X 2 RSS = (yi − y^i ) = (yi − (a + bxi )) : i=1 i=1 This is a Calculus I problem. -
A Statistical Test Suite for Random and Pseudorandom Number Generators for Cryptographic Applications
Special Publication 800-22 Revision 1a A Statistical Test Suite for Random and Pseudorandom Number Generators for Cryptographic Applications AndrewRukhin,JuanSoto,JamesNechvatal,Miles Smid,ElaineBarker,Stefan Leigh,MarkLevenson,Mark Vangel,DavidBanks,AlanHeckert,JamesDray,SanVo Revised:April2010 LawrenceE BasshamIII A Statistical Test Suite for Random and Pseudorandom Number Generators for NIST Special Publication 800-22 Revision 1a Cryptographic Applications 1 2 Andrew Rukhin , Juan Soto , James 2 2 Nechvatal , Miles Smid , Elaine 2 1 Barker , Stefan Leigh , Mark 1 1 Levenson , Mark Vangel , David 1 1 2 Banks , Alan Heckert , James Dray , 2 San Vo Revised: April 2010 2 Lawrence E Bassham III C O M P U T E R S E C U R I T Y 1 Statistical Engineering Division 2 Computer Security Division Information Technology Laboratory National Institute of Standards and Technology Gaithersburg, MD 20899-8930 Revised: April 2010 U.S. Department of Commerce Gary Locke, Secretary National Institute of Standards and Technology Patrick Gallagher, Director A STATISTICAL TEST SUITE FOR RANDOM AND PSEUDORANDOM NUMBER GENERATORS FOR CRYPTOGRAPHIC APPLICATIONS Reports on Computer Systems Technology The Information Technology Laboratory (ITL) at the National Institute of Standards and Technology (NIST) promotes the U.S. economy and public welfare by providing technical leadership for the nation’s measurement and standards infrastructure. ITL develops tests, test methods, reference data, proof of concept implementations, and technical analysis to advance the development and productive use of information technology. ITL’s responsibilities include the development of technical, physical, administrative, and management standards and guidelines for the cost-effective security and privacy of sensitive unclassified information in Federal computer systems. -
BULLETIN of DELTA STATE UNIVERSITY (USPS) 152-880) Cleveland, Mississippi 38733 Founded April 9, 1924 Opened September 25, 1925
2011-2012ug1 pages 1-65.qxp 6/3/2011 2:47 PM Page 1 BULLETIN OF DELTA STATE UNIVERSITY (USPS) 152-880) Cleveland, Mississippi 38733 Founded April 9, 1924 Opened September 25, 1925 UNDERGRADUATE ANNOUNCEMENTS EIGHTY-FOURTH SESSION 2011-2012 CATALOG This bulletin presents information which, at the time of preparation for printing, most ac- curately describes the course offerings, policies, procedures, regulations, and requirements of the University. However, it does not establish contractual relationships. The University reserves the right to alter or change any statement contained herein without prior notice. Volume 86 June 2011 Number 2 PUBLISHED ANNUALLY BY DELTA STATE UNIVERSITY CLEVELAND, MISSISSIPPI 38733 Entered as Second-Class Matter July 1926, at the Post Office at Cleveland, Mississippi, Under the Act of February 28, 1925. Postmaster: Send change of address to Delta State University, Cleveland, Mississippi 38733 2011-2012ug1 pages 1-65.qxp 6/3/2011 2:47 PM Page 2 TABLE OF CONTENTS Board of Trustees ..................................................................................................3 The University Calendar ......................................................................................4 General Information ............................................................................................6 Facilities ............................................................................................................11 Student Activities ................................................................................................18 -
Chapter 5 Statistical Models in Simulation
Chapter 5 Statistical Models in Simulation Banks, Carson, Nelson & Nicol Discrete-Event System Simulation Purpose & Overview The world the model-builder sees is probabilistic rather than deterministic. Some statistical model might well describe the variations. An appropriate model can be developed by sampling the phenomenon of interest: Select a known distribution through educated guesses Make estimate of the parameter(s) Test for goodness of fit In this chapter: Review several important probability distributions Present some typical application of these models ٢ ١ Review of Terminology and Concepts In this section, we will review the following concepts: Discrete random variables Continuous random variables Cumulative distribution function Expectation ٣ Discrete Random Variables [Probability Review] X is a discrete random variable if the number of possible values of X is finite, or countably infinite. Example: Consider jobs arriving at a job shop. Let X be the number of jobs arriving each week at a job shop. Rx = possible values of X (range space of X) = {0,1,2,…} p(xi) = probability the random variable is xi = P(X = xi) p(xi), i = 1,2, … must satisfy: 1. p(xi ) ≥ 0, for all i ∞ 2. p(x ) =1 ∑i=1 i The collection of pairs [xi, p(xi)], i = 1,2,…, is called the probability distribution of X, and p(xi) is called the probability mass function (pmf) of X. ٤ ٢ Continuous Random Variables [Probability Review] X is a continuous random variable if its range space Rx is an interval or a collection of intervals. The probability that X lies in the interval [a,b] is given by: b P(a ≤ X ≤ b) = f (x)dx ∫a f(x), denoted as the pdf of X, satisfies: 1. -
Ed Leith Cretaceous Menagerie Promotes Outreach at University of Manitoba
Canadian Publication Mail Contract - 40070050 $3.00 VOLUME 31, ISSUE 4 APRIL 2004 ■ Ed Leith Cretaceous Menagerie Promotes Outreach at University of Manitoba ■ Your Roadmap to Tradeshow Success ■ 90 Reasons to Exhibit AccuMapTM Comprehensive, trusted oil and gas information for the y. a w Western Canadian e Sedimentary Basin, Northern th f o Territories, and East Coast. p e From mineral rights to sur- t s face information, E&P data y r e to M&A, AccuMap offers v e unparalleled data and analy- u o sis tools. y h t i New ventures. Exploration. w e Economics. Exploitation. r ’ e Drilling. Production. w , Delivery. Acquisitions and s a Divestitures. Wherever your g & opportunities take you, IHS il o Energy supports you every f o step of the way. it u rs pu ur yo In AccuMap Hotkeys reference guide! Yours free in this issue of the Reservoir. CSPG OFFICE #160, 540 - 5th Avenue SW Calgary,Alberta, Canada T2P 0M2 Tel:403-264-5610 Fax: 403-264-5898 Email: [email protected] Web: www.cspg.org Office hours: Monday to Friday, 8:30am to 4:00pm CONTENTS Business Manager:Tim Howard Email: [email protected] Office Manager: Deanna Watkins Email: [email protected] Communications Manager: Jaimè Croft Email: [email protected] ARTICLES Conventions Manager: Lori Humphrey-Clements Email: [email protected] Corporate Relations Manager: Kim MacLean Email: [email protected] STAFF PROFILE: KIM MACLEAN, CORPORATE RELATIONS . .7 ED LEITH CRETACEOUS MENAGERIE PROMOTES OUTREACH EDITORS/AUTHORS AT UNIVERSITY OF MANITOBA . 22 Please submit RESERVOIR articles to the CSPG office.