The sensitivity Package
|
|
- Christopher Robbins
- 5 years ago
- Views:
Transcription
1 The sensitivity Package January 5, 2007 Version Date Title Sensitivity Analysis Author Gilles Pujol (this package was originally developped at Commissariat a l Energie Atomique CEA, Service d Etudes et de Simulation du Comportement des Combustibles CEA/DEN/CAD/DEC/SESC, France ; the author would like to thank : Patrick Obry, Bertrand Iooss, Claire Cannamela and Frederic Michel). Maintainer Bertrand Iooss <bertrand.iooss@cea.fr> Depends R, boot This package allows to perform sensitivity analyses within the R environment. Implemented methods are : linear sensitivity analysis (SRC, PCC, rank analysis), the screening method of Morris, the Sobol global sensitivity indices (two methods of estimation), and the FAST method. License CeCILL version 2 R topics documented: fast morris INTRODUCTION sobol srcpcc tell testmodels Index 12 1
2 2 fast fast Fourier Amplitude Sensitivity Test Usage fast is the implementation of the Fourier Amplitude Sensitivity Test. fast(method = "saltelli99", model = NULL, factors, n, M = 4, omega = NULL, q = NULL, q.arg = NULL,...) Arguments method model factors n M omega q Details q.arg the method: "saltelli99" only the model the number of factors, or their names the sample size the interference parameter the set of frequencies the names of the quantile functions for the factors distributions the quantile parameters... any other arguments for model which are passed unchanged each time it is called The method "saltelli99" is the so-called extended fast method wich provides estimations of both first order and total indices at a low computational cost. model is a function or a predictor (a class with a predict method) computing the response y based on the sample given by x. If no model is specified, the indices will be computed when one gives the response. factors could either be a single number or a vector of character strings. n is the length of the discretization of the s-space (for computing Fourier coefficients) and M is the number of harmonics to sum (for computing partial variances). If the set of frequencies omega is not given, the function use the set recommended by the corresponding method. For the method "saltelli99", the first frequency is the greater, associated with the input variable to assess, and the other frequencies are associated with the complementary set. If q and q.args are not given, the factors will be considerd uniform on [0,1]. q is a list of character strings giving the names of the quantile functions (one for each factor), such as qunif, qnorm... q could also be a single character string (the same for all). q.arg is a list of lists, each list being additional parameters for the corresponding quantile function. For example, the parameters of the quantile qunif could be (min=1, max=2) giving an uniform distribution on [1,2]. If q is a single character string, then q.arg must be a single list.
3 morris 3 Value fast returns an object of class "fast". An object of class "fast" is a list containing the following components: x y S St the factor sample the response the estimations of the first-order indices the estimations of the total indices (method "saltelli99") Computational cost For the method "saltelli99", the number of model evaluations is p n where p is the number of factors. References Saltelli, A., Tarantola, S. and Chan, K., 1999, A quantitative, model independent method for global sensitivity analysis of model output. Technometrics, 41, Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis. Wiley. Cukier, R. I., Levine, H. B. and Schuler, K. E., 1978, Nonlinear sensitivity analysis of multiparameter model systems. J. Comput. Phys., 26, Examples # Test case : the non-monotonic Ishigami function sa <- fast(model = ishigami.fun, factors = 3, n = 1000, q = "qunif", q.arg = list(min = -pi, max = pi)) print(sa) plot(sa) morris The Morris OAT Screening Method morris is the implementation of the Morris OAT Screening method. This function generates the Morris design of experiments and computes the measures of sensitivity µ and σ. Usage morris(model = NULL, factors, levels, R, jump = NULL, min = 0, max = 1, scale = TRUE, optim = NULL,...)
4 4 morris Arguments model factors levels R jump min max Details Value scale optim the model the number of factors, or their names the number of levels of the design grid the number of repetitions of the design, i.e. the number of elementary effect computed per factor the grid jump coefficient the minimum values for the factors the maximum values for the factors logical. If TRUE, the input and output data are scaled optimization of the design for better coverage of the space (cf Campolongo 2005), not documented yet (for informations feel free to ask the maintainer)... any other arguments for model which are passed unchanged each time it is called model is a function or a predictor (a class with a predict method) computing the response y based on the sample given by x. If no model is specified, the indices will be computed when one gives the response. factors could either be a single number or a vector of character strings. The number of levels is not necessary the same for each space coordinate. It is the case when levels is a single integer. min and max are boundaries of the region of experimentation. They can be single values (the same for each factor) or vectors. jump is such that: i = jump i max i min i levels i 1 If jump is given as NULL and the number of levels is even (for each component), then jump has the value recommended by Morris: jump = levels/2. If jump is a single value, then it is taken the same for each coordinate. morris returns an object of class "morris". An object of class "morris" is a list containing the following components: x y ee mu sigma the design of experiments (input sample) the response the matrix of the elementary effects the estimations of the µ index the estimations of the σ index
5 INTRODUCTION 5 Computational cost The number of model evaluations is (p + 1) R where p is the number of factors. References Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis. Wiley. Morris, M. D., 1991, Factorial sampling plans for preliminary computational experiments. Technometrics, 33, Examples # Test case : the non-monotonic function of Morris sa <- morris(model = morris.fun, factors = 20, levels = 4, R = 4) print(sa) plot(sa) INTRODUCTION Package sensitivity : Sensitivity Analysis The sensitivity package implements sensitivity analysis methods: linear and monotonic sensitivity analysis (SRC, PCC, SRRC, PRCC), the screening method of Morris, and non-linear global sensitivity analysis (the Sobol indices, the FAST method). The functions of this package generate the design of experiments (depending on the method of analysis) and compute the sensitivity indices based on the model inputs and outputs. All sensitivity indices can be estimated with the bootstrap technique which allows to estimate the bias, and basic bootstrap confidence intervals. Text and graphical outputs display the results of the analysis. Details The approach applied when performing a sensitivity analysis (SA) is as follows: step 1 The model is defined: it is a function that returns the (real) ouput values (called responses), corresponding to a sample of (real) input parameters (called factors). step 2 A sensitivity analysis method is chosen. Parameters of this method must be in accordance with objectives and technical constraints (like computational time). step 3 A design of experiments (DOE) corresponding to the SA method is generated. step 4 The model is evaluated on the DOE values. step 5 The sensitivity indices are computed, based on input and output values. step 6 Post-treatments...
6 6 sobol The sensitivity package allows to follow this methodology: (step 1) The model can be internal or external to R. If internal, it can be a function that takes an unique matrix or data.frame parameter and returns a numeric vector. It can also be a predictor, i.e. an object wich can be called with the predict method. One should note that all the responses must be computed by a single call to the model function (then, the model can be vectorized). If the model is external it does not have to be interfaced with R: the user won t have to give a model to the function. Then, it will stop just after generating the DOE. The responses have to be computed by the user, whithin R or not. Calculations will start again when the user gives the corresponding responses (via the tell function). The four next steps depend upon the type of the model: For internal models: (step 2-5) sa <- method(model, parameters...) For external models: (step 2-3) sa <- method(model = NULL, parameters...) (step 4) external to R (or not), and the result is loaded by the user in the y variable (step 5) tell(sa, y) method should be the name of a SA function, such as srcpcc, morris, sobol, or fast. These function create the object sa of class "srcpcc", "morris", "sobol", or "fast". For further information on these function, see the corresponding documentation. Finally, for displaying the results of the analysis: (step 6) print(sa); plot(sa) References Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis. Wiley, See Also srcpcc morris sobol fast tell testmodels sobol Sobol Non-linear Sensitivity Analysis sobol is the implementation of the Monte Carlo estimation of the Sobol indices. Usage sobol(method = "sobol93", model = NULL, x1, x2, max.order = 1, nboot = 0, conf = 0.95,...)
7 sobol 7 Arguments method model x1 x2 Details Value max.order nboot conf the method: "sobol93" or "saltelli02" the model the first random sample the second random sample the maximum order of indices to compute (method "sobol93") the number of bootstrap replicates the confidence level for bootstrap confidence intervals... any other arguments for model which are passed unchanged each time it is called Two methods. The method "sobol93" computes all the Sobol indices (coming from the HDMR- ANOVA decomposition) from order 1 to order given by the argument max.order. The method "saltelli02" computes both first order and total indices at a reduced computational cost. model is a function or a predictor (a class with a predict method) computing the response y based on the sample given by x. If no model is specified, the indices will be computed when one gives the response. The Monte Carlo estimation requires two independent random samples x1 and x2. They must have the same dimensions. sobol returns an object of class "sobol". An object of class "sobol" is a list containing the following components: x y S St the factor sample the response used the estimations of the Sobol sensitivity indices the estimations of the total indices (method "saltelli02") Computational cost For the method "sobol93", the number of model evaluations is n (N + 1) where n is the size of the samples x1 and x2, and N is the number of indices to estimate. For the method "saltelli02", the number of model evaluations is n (p + 2) where p is the number of factors (for the estimation of 2p indices). References Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis. Wiley. Sobol, I. M., 1993, Sensitivity analysis for non-linear mathematical model. Math. Modelling Comput. Exp., 1, Saltelli, A., 2002, Making best use of model evaluations to compute sensitivity indices. Computer Physics Communication, 145,
8 8 srcpcc Examples # Test case : the non-monotonic Sobol g-function # The method of sobol requires 2 samples # There are 8 factors, all following the uniform distribution # on [0,1] n < x <- data.frame(matrix(nr = 2 * n, nc = 8)) for (i in 1:8) x[, i] <- runif(2 * n) # sensitivity analysis sa <- sobol(model = sobol.fun, x1 = x[1:n,], x2 = x[(n+1):(2*n),], max.order = 2, nboot = 10 print(sa) #plot(sa) srcpcc Linear Sensitivity Analysis srcpcc computes the standardized regression coefficients (SRC) and the partial correlation coefficients (PCC). Analysis can be done on the ranks; then the indices are the standardized rank regression coefficients (SRRC) and the partial rank correlation coefficients (PRCC). Usage srcpcc(model = NULL, x, pcc = TRUE, rank = FALSE, nboot = 0, conf = 0.95,...) Arguments model x pcc Details rank nboot conf the model the input sample logical. If TRUE, the P(R)CCs are computed logical. If TRUE, the analysis is done on the ranks the number of bootstrap replicates the confidence level for bootstrap confidence intervals... any other arguments for model which are passed unchanged each time it is called model is a function or a predictor (a class with a predict method) computing the response y based on the sample given by x. If no model is specified, the indices will be computed when one gives the response.
9 tell 9 Value srcpcc returns an object of class "srcpcc". An object of class "srcpcc" is a list containing the following components: y src pcc the response the estimations of the SRC indices (or SRRC if rank analysis is requested) if requested, the estimations of the PCC indices (or PRCC if rank analysis is requested) Computational cost The number of model evaluations is n where n is the size of the sample x. References Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis, Wiley. Examples # linear model : Y = X1 + X2 + X3 model1 <- function(x) x[, 1] + x[, 2] + x[, 3] # a 100-sample with X1 ~ U(0.5, 1.5) # X2 ~ U(1.5, 4.5) # X3 ~ U(4.5, 13.5) n <- 100 x <- data.frame(x1 = runif(n, 0.5, 1.5), X2 = runif(n, 1.5, 4.5), X3 = runif(n, 4.5, 13.5)) # sensitivity analysis sa <- srcpcc(model = model1, x = x, nboot = 100) print(sa) par(mfrow = c(1,2)) plot(sa, ask = FALSE) tell Computation Of Sensitivity Indices When The Model Is External tell is used to tell a sensitivity analysis object the results of the simulations. It is used when the model is not given when parametring the sensitivity analysis (whith the model=null argument). For example, it is the case when the model is external to R.
10 10 testmodels Usage tell(sa, y = NULL) Arguments sa y the sensitivity analysis object the response Details Value sa is an object returned by a sensitivity analysis function, such as srcpcc, morris, sobol,... y should be a numeric vector. tell doesn t return anything. It does the sensitivity analysis, and stores all the results in the variable sa. Examples # Example of the FAST method # (one should note the call with model = NULL) sa <- fast(model = NULL, factors = 8, n = 1000) # at this stage, only the design of experiment (sa$x) was generated # the response is computed "manually": y <- sobol.fun(sa$x) # at this place could be a # call to an external code # then, the sensitivity analysis: tell(sa, y) print(sa) # Remark : because the model is a simple R function, # this example is equivalent to : ## Not run: sa <- fast(model = sobol.fun, factors = 8, n = 1000) testmodels Test Models For Sensitivity Analysis These functions are standard testcase for sensitivity analysis benchmarks. There are: the g-function of Sobol, the function of Ishigami and the function of Morris.
11 testmodels 11 Usage sobol.fun(x) ishigami.fun(x) morris.fun(x) Arguments x the matrix or data.frame containing the input values. Value All these functions return a numeric vector containig the values of the function. These functions are vectorized. References Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis. Wiley,
12 Index Topic misc fast, 1 INTRODUCTION, 5 morris, 3 sobol, 6 srcpcc, 8 tell, 9 testmodels, 10 fast, 1, 6 INTRODUCTION, 5 ishigami.fun (testmodels), 10 morris, 3, 6, 10 morris.fun (testmodels), 10 plot.fast.saltelli99 (fast), 1 plot.morris (morris), 3 plot.sobol.saltelli02 (sobol), 6 plot.sobol.sobol93 (sobol), 6 plot.srcpcc (srcpcc), 8 print.fast.saltelli99 (fast), 1 print.morris (morris), 3 print.sobol.saltelli02 (sobol), 6 print.sobol.sobol93 (sobol), 6 print.srcpcc (srcpcc), 8 sensitivity (INTRODUCTION), 5 sensitivity-package (INTRODUCTION), 5 sobol, 6, 6, 10 sobol.fun (testmodels), 10 srcpcc, 6, 8, 10 tell, 6, 9 tell.fast.saltelli99 (fast), 1 tell.morris (morris), 3 tell.sobol.saltelli02 (sobol), 6 tell.sobol.sobol93 (sobol), 6 tell.srcpcc (srcpcc), 8 testmodels, 6, 10 12
Package sensitivity. R topics documented: February 15, Version Date Title Sensitivity Analysis
Version 1.6-1 Date 2012-12-28 Title Sensitivit Analsis Package sensitivit Februar 15, 2013 Author Gilles Pujol, Bertrand Iooss, Alexandre Janon Maintainer Bertrand Iooss Depends R (>=
More informationPackage pse. June 11, 2017
Type Package Package pse June 11, 2017 Title Parameter Space Exploration with Latin Hypercubes Version 0.4.7 Date 2017-06-11 Author Andre Chalom, Paulo Inacio Knegt Lopez de Prado Maintainer Andre Chalom
More informationSimulation: Solving Dynamic Models ABE 5646 Week 12, Spring 2009
Simulation: Solving Dynamic Models ABE 5646 Week 12, Spring 2009 Week Description Reading Material 12 Mar 23- Mar 27 Uncertainty and Sensitivity Analysis Two forms of crop models Random sampling for stochastic
More informationPackage CompModSA. February 14, 2012
Package CompModSA February 14, 2012 Type Package Title Sensitivity Analysis for Complex Computer Models Version 1.3.1 Date 2011-11-05 Author Curtis Storlie , with contributions from
More informationDesign of Experiments
Seite 1 von 1 Design of Experiments Module Overview In this module, you learn how to create design matrices, screen factors, and perform regression analysis and Monte Carlo simulation using Mathcad. Objectives
More informationUniversity of Leeds. School of Process, Environmental and Materials Engineering Energy and Resources Research Institute GUI-HDMR
University of Leeds School of Process, Environmental and Materials Engineering Energy and Resources Research Institute User Documentation (Draft) GUI-HDMR Version 1.1 Developer: Tilo Ziehn Contact: Professor
More informationPackage sensitivity. September 23, 2017
Version 1.15.0 Package sensitivit September 23, 2017 Title Global Sensitivit Analsis of Model Outputs Author Gilles Pujol, Bertrand Iooss, Alexandre Janon with contributions from Khalid Boumhaout, Sebastien
More informationBeta-Regression with SPSS Michael Smithson School of Psychology, The Australian National University
9/1/2005 Beta-Regression with SPSS 1 Beta-Regression with SPSS Michael Smithson School of Psychology, The Australian National University (email: Michael.Smithson@anu.edu.au) SPSS Nonlinear Regression syntax
More informationSTATISTICS (STAT) Statistics (STAT) 1
Statistics (STAT) 1 STATISTICS (STAT) STAT 2013 Elementary Statistics (A) Prerequisites: MATH 1483 or MATH 1513, each with a grade of "C" or better; or an acceptable placement score (see placement.okstate.edu).
More informationIntroduction to Statistical Analyses in SAS
Introduction to Statistical Analyses in SAS Programming Workshop Presented by the Applied Statistics Lab Sarah Janse April 5, 2017 1 Introduction Today we will go over some basic statistical analyses in
More informationResources for statistical assistance. Quantitative covariates and regression analysis. Methods for predicting continuous outcomes.
Resources for statistical assistance Quantitative covariates and regression analysis Carolyn Taylor Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC January 24, 2017 Department
More informationPackage citools. October 20, 2018
Type Package Package citools October 20, 2018 Title Confidence or Prediction Intervals, Quantiles, and Probabilities for Statistical Models Version 0.5.0 Maintainer John Haman Functions
More informationWe deliver Global Engineering Solutions. Efficiently. This page contains no technical data Subject to the EAR or the ITAR
Numerical Computation, Statistical analysis and Visualization Using MATLAB and Tools Authors: Jamuna Konda, Jyothi Bonthu, Harpitha Joginipally Infotech Enterprises Ltd, Hyderabad, India August 8, 2013
More informationGAUSSIAN PROCESS RESPONSE SURFACE MODELING AND GLOBAL SENSITIVITY ANALYSIS USING NESSUS
UNCECOMP 2017 2 nd ECCOMAS Thematic Conference on International Conference on Uncertainty Quantification in Computational Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou (eds.) Rhodes
More informationLecture 25: Review I
Lecture 25: Review I Reading: Up to chapter 5 in ISLR. STATS 202: Data mining and analysis Jonathan Taylor 1 / 18 Unsupervised learning In unsupervised learning, all the variables are on equal standing,
More informationVisualization tools for uncertainty and sensitivity analyses on thermal-hydraulic transients
Joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo 013 (SNA + MC 013) La Cité des Sciences et de l Industrie, Paris, France, October 7-31, 013 Visualization tools
More informationPackage acebayes. R topics documented: November 21, Type Package
Type Package Package acebayes November 21, 2018 Title Optimal Bayesian Experimental Design using the ACE Algorithm Version 1.5.2 Date 2018-11-21 Author Antony M. Overstall, David C. Woods & Maria Adamou
More information2014 Stat-Ease, Inc. All Rights Reserved.
What s New in Design-Expert version 9 Factorial split plots (Two-Level, Multilevel, Optimal) Definitive Screening and Single Factor designs Journal Feature Design layout Graph Columns Design Evaluation
More informationGetting started with simulating data in R: some helpful functions and how to use them Ariel Muldoon August 28, 2018
Getting started with simulating data in R: some helpful functions and how to use them Ariel Muldoon August 28, 2018 Contents Overview 2 Generating random numbers 2 rnorm() to generate random numbers from
More informationBluman & Mayer, Elementary Statistics, A Step by Step Approach, Canadian Edition
Bluman & Mayer, Elementary Statistics, A Step by Step Approach, Canadian Edition Online Learning Centre Technology Step-by-Step - Minitab Minitab is a statistical software application originally created
More informationPackage uclaboot. June 18, 2003
Package uclaboot June 18, 2003 Version 0.1-3 Date 2003/6/18 Depends R (>= 1.7.0), boot, modreg Title Simple Bootstrap Routines for UCLA Statistics Author Maintainer
More informationThe Bolstad Package. July 9, 2007
The Bolstad Package July 9, 2007 Version 0.2-12 Date 2007-09-07 Title Bolstad functions Author James Curran Maintainer James M. Curran A set of
More informationThe simpleboot Package
The simpleboot Package April 1, 2005 Version 1.1-1 Date 2005-03-31 LazyLoad yes Depends R (>= 2.0.0), boot Title Simple Bootstrap Routines Author Maintainer Simple bootstrap
More informationInteractive Building Design Space Exploration Using Regionalized Sensitivity Analysis Østergård, Torben; Jensen, Rasmus Lund; Maagaard, Steffen
Aalborg Universitet Interactive Building Design Space Exploration Using Regionalized Sensitivity Analysis Østergård, Torben; Jensen, Rasmus Lund; Maagaard, Steffen Published in: DOI (link to publication
More informationSTATS PAD USER MANUAL
STATS PAD USER MANUAL For Version 2.0 Manual Version 2.0 1 Table of Contents Basic Navigation! 3 Settings! 7 Entering Data! 7 Sharing Data! 8 Managing Files! 10 Running Tests! 11 Interpreting Output! 11
More informationBootstrapped and Means Trimmed One Way ANOVA and Multiple Comparisons in R
Bootstrapped and Means Trimmed One Way ANOVA and Multiple Comparisons in R Another way to do a bootstrapped one-way ANOVA is to use Rand Wilcox s R libraries. Wilcox (2012) states that for one-way ANOVAs,
More informationTECHNICAL REPORT Use of Sobol's Quasirandom Sequence Generator for Integration of Modified Uncertainty Importance Measure
Journal of NUCLEAR SCIENCE and TECHNOLOGY, 32 [11], pp.1164-1173 (November 1995). TECHNICAL REPORT Use of Sobol's Quasirandom Sequence Generator for Integration of Modified Uncertainty Importance Measure
More informationFathom Dynamic Data TM Version 2 Specifications
Data Sources Fathom Dynamic Data TM Version 2 Specifications Use data from one of the many sample documents that come with Fathom. Enter your own data by typing into a case table. Paste data from other
More informationPackage deming. June 19, 2018
Package deming June 19, 2018 Title Deming, Thiel-Sen and Passing-Bablock Regression Maintainer Terry Therneau Generalized Deming regression, Theil-Sen regression and Passing-Bablock
More informationJMP Book Descriptions
JMP Book Descriptions The collection of JMP documentation is available in the JMP Help > Books menu. This document describes each title to help you decide which book to explore. Each book title is linked
More informationUsing PC SAS/ASSIST* for Statistical Analyses
Using PC SAS/ASSIST* for Statistical Analyses Margaret A. Nemeth, Monsanto Company lptroductjon SAS/ASSIST, a user friendly, menu driven applications system, is available on several platforms. This paper
More informationPredict Outcomes and Reveal Relationships in Categorical Data
PASW Categories 18 Specifications Predict Outcomes and Reveal Relationships in Categorical Data Unleash the full potential of your data through predictive analysis, statistical learning, perceptual mapping,
More informationPackage lcc. November 16, 2018
Type Package Title Longitudinal Concordance Correlation Version 1.0 Author Thiago de Paula Oliveira [aut, cre], Rafael de Andrade Moral [aut], John Hinde [aut], Silvio Sandoval Zocchi [aut, ctb], Clarice
More informationSensitivity analysis of the design portal frames of steel industrial buildings
Sensitivity analysis of the design portal frames of steel industrial buildings Trong Ha Nguyen 1,* and Xuan Hung Dang 2 1 Faculty of Civil Engineering, Vinh University, 182 LeDuan Vinh City, Vietnam 2
More informationBootstrapping Methods
Bootstrapping Methods example of a Monte Carlo method these are one Monte Carlo statistical method some Bayesian statistical methods are Monte Carlo we can also simulate models using Monte Carlo methods
More informationData analysis using Microsoft Excel
Introduction to Statistics Statistics may be defined as the science of collection, organization presentation analysis and interpretation of numerical data from the logical analysis. 1.Collection of Data
More informationPaper PO-06. Gone are the days when social and behavioral science researchers should simply report obtained test statistics (e.g.
Paper PO-06 CI_MEDIATE: A SAS Macro for Computing Point and Interval Estimates of Effect Sizes Associated with Mediation Analysis Thanh V. Pham, University of South Florida, Tampa, FL Eun Kyeng Baek, University
More informationThe relaimpo Package
The relaimpo Package October 1, 2007 Title Relative importance of regressors in linear models Version 1.2-2 Date 2007-09-30 Author Ulrike Groemping Description relaimpo provides several metrics for assessing
More informationMacros and ODS. SAS Programming November 6, / 89
Macros and ODS The first part of these slides overlaps with last week a fair bit, but it doesn t hurt to review as this code might be a little harder to follow. SAS Programming November 6, 2014 1 / 89
More informationCDAA No. 4 - Part Two - Multiple Regression - Initial Data Screening
CDAA No. 4 - Part Two - Multiple Regression - Initial Data Screening Variables Entered/Removed b Variables Entered GPA in other high school, test, Math test, GPA, High school math GPA a Variables Removed
More informationSampling and Monte-Carlo Integration
Sampling and Monte-Carlo Integration Sampling and Monte-Carlo Integration Last Time Pixels are samples Sampling theorem Convolution & multiplication Aliasing: spectrum replication Ideal filter And its
More informationStatistical Matching using Fractional Imputation
Statistical Matching using Fractional Imputation Jae-Kwang Kim 1 Iowa State University 1 Joint work with Emily Berg and Taesung Park 1 Introduction 2 Classical Approaches 3 Proposed method 4 Application:
More informationarxiv: v2 [stat.ap] 9 Jun 2008
Global sensitivity analysis of computer models with functional inputs arxiv:0802.1009v2 [stat.ap] 9 Jun 2008 Bertrand IOOSS and Mathieu RIBATET Submitted to: Reliability Engineering and System Safety for
More informationData Analysis and Solver Plugins for KSpread USER S MANUAL. Tomasz Maliszewski
Data Analysis and Solver Plugins for KSpread USER S MANUAL Tomasz Maliszewski tmaliszewski@wp.pl Table of Content CHAPTER 1: INTRODUCTION... 3 1.1. ABOUT DATA ANALYSIS PLUGIN... 3 1.3. ABOUT SOLVER PLUGIN...
More informationPackage fbroc. August 29, 2016
Type Package Package fbroc August 29, 2016 Title Fast Algorithms to Bootstrap Receiver Operating Characteristics Curves Version 0.4.0 Date 2016-06-21 Implements a very fast C++ algorithm to quickly bootstrap
More informationMaths PoS: Year 7 HT1. Students will colour code as they work through the scheme of work. Students will learn about Number and Shape
Maths PoS: Year 7 HT1 Students will learn about Number and Shape Number: Use positive and negative numbers in context and position them on a number line. Recall quickly multiplication facts up to 10 10
More informationExtending the capabilities of agriculture simulators using R, an introduction to the apsimr package by Bryan Stanfill
1 Extending the capabilities of agriculture simulators using R, an introduction to the apsimr package by Bryan Stanfill Introduction The Agricultural Production System simulator (APSIM) is a widely used,
More informationSTATISTICS FOR PSYCHOLOGISTS
STATISTICS FOR PSYCHOLOGISTS SECTION: JAMOVI CHAPTER: USING THE SOFTWARE Section Abstract: This section provides step-by-step instructions on how to obtain basic statistical output using JAMOVI, both visually
More informationStat 528 (Autumn 2008) Density Curves and the Normal Distribution. Measures of center and spread. Features of the normal distribution
Stat 528 (Autumn 2008) Density Curves and the Normal Distribution Reading: Section 1.3 Density curves An example: GRE scores Measures of center and spread The normal distribution Features of the normal
More informationUncertainty Quantification and Sensitivity Analysis of Reservoir Forecasts with Machine Learning
CS 9 Proect Final Report - ihoon Park Uncertainty Quantification and Sensitivity Analysis of Reservoir Forecasts with Machine Learning ihoon Park (hpark3@stanford.edu) Introduction Successful development
More informationThe grplasso Package
The grplasso Package June 27, 2007 Type Package Title Fitting user specified models with Group Lasso penalty Version 0.2-1 Date 2007-06-27 Author Lukas Meier Maintainer Lukas Meier
More informationSimulation Supported POD Methodology and Validation for Automated Eddy Current Procedures
4th International Symposium on NDT in Aerospace 2012 - Th.1.A.1 Simulation Supported POD Methodology and Validation for Automated Eddy Current Procedures Anders ROSELL, Gert PERSSON Volvo Aero Corporation,
More informationSimulation studies. Patrick Breheny. September 8. Monte Carlo simulation Example: Ridge vs. Lasso vs. Subset
Simulation studies Patrick Breheny September 8 Patrick Breheny BST 764: Applied Statistical Modeling 1/17 Introduction In statistics, we are often interested in properties of various estimation and model
More informationRobust Linear Regression (Passing- Bablok Median-Slope)
Chapter 314 Robust Linear Regression (Passing- Bablok Median-Slope) Introduction This procedure performs robust linear regression estimation using the Passing-Bablok (1988) median-slope algorithm. Their
More informationRecent advances in Metamodel of Optimal Prognosis. Lectures. Thomas Most & Johannes Will
Lectures Recent advances in Metamodel of Optimal Prognosis Thomas Most & Johannes Will presented at the Weimar Optimization and Stochastic Days 2010 Source: www.dynardo.de/en/library Recent advances in
More informationLecture: Simulation. of Manufacturing Systems. Sivakumar AI. Simulation. SMA6304 M2 ---Factory Planning and scheduling. Simulation - A Predictive Tool
SMA6304 M2 ---Factory Planning and scheduling Lecture Discrete Event of Manufacturing Systems Simulation Sivakumar AI Lecture: 12 copyright 2002 Sivakumar 1 Simulation Simulation - A Predictive Tool Next
More informationMinitab 17 commands Prepared by Jeffrey S. Simonoff
Minitab 17 commands Prepared by Jeffrey S. Simonoff Data entry and manipulation To enter data by hand, click on the Worksheet window, and enter the values in as you would in any spreadsheet. To then save
More informationAn Introduction to the Bootstrap
An Introduction to the Bootstrap Bradley Efron Department of Statistics Stanford University and Robert J. Tibshirani Department of Preventative Medicine and Biostatistics and Department of Statistics,
More informationMultiple runs of the same parameter combination with R package pse
Multiple runs of the same parameter combination with R package pse Chalom, A. Prado, P.I. Version 0.3.1, November 23, 2013 This document presents an extension to the practical tutorial found in the vignette
More informationStatistics (STAT) Statistics (STAT) 1. Prerequisites: grade in C- or higher in STAT 1200 or STAT 1300 or STAT 1400
Statistics (STAT) 1 Statistics (STAT) STAT 1200: Introductory Statistical Reasoning Statistical concepts for critically evaluation quantitative information. Descriptive statistics, probability, estimation,
More informationOverview. Background. Locating quantitative trait loci (QTL)
Overview Implementation of robust methods for locating quantitative trait loci in R Introduction to QTL mapping Andreas Baierl and Andreas Futschik Institute of Statistics and Decision Support Systems
More informationR Programming: Worksheet 6
R Programming: Worksheet 6 Today we ll study a few useful functions we haven t come across yet: all(), any(), `%in%`, match(), pmax(), pmin(), unique() We ll also apply our knowledge to the bootstrap.
More informationGeneral Factorial Models
In Chapter 8 in Oehlert STAT:5201 Week 9 - Lecture 2 1 / 34 It is possible to have many factors in a factorial experiment. In DDD we saw an example of a 3-factor study with ball size, height, and surface
More informationLecture 7: Linear Regression (continued)
Lecture 7: Linear Regression (continued) Reading: Chapter 3 STATS 2: Data mining and analysis Jonathan Taylor, 10/8 Slide credits: Sergio Bacallado 1 / 14 Potential issues in linear regression 1. Interactions
More informationComputer Experiments. Designs
Computer Experiments Designs Differences between physical and computer Recall experiments 1. The code is deterministic. There is no random error (measurement error). As a result, no replication is needed.
More informationThe glmmml Package. August 20, 2006
The glmmml Package August 20, 2006 Version 0.65-1 Date 2006/08/20 Title Generalized linear models with clustering A Maximum Likelihood and bootstrap approach to mixed models. License GPL version 2 or newer.
More informationExponential Random Graph Models for Social Networks
Exponential Random Graph Models for Social Networks ERGM Introduction Martina Morris Departments of Sociology, Statistics University of Washington Departments of Sociology, Statistics, and EECS, and Institute
More informationSOM+EOF for Finding Missing Values
SOM+EOF for Finding Missing Values Antti Sorjamaa 1, Paul Merlin 2, Bertrand Maillet 2 and Amaury Lendasse 1 1- Helsinki University of Technology - CIS P.O. Box 5400, 02015 HUT - Finland 2- Variances and
More informationR Programming Basics - Useful Builtin Functions for Statistics
R Programming Basics - Useful Builtin Functions for Statistics Vectorized Arithmetic - most arthimetic operations in R work on vectors. Here are a few commonly used summary statistics. testvect = c(1,3,5,2,9,10,7,8,6)
More informationPackage CALIBERrfimpute
Type Package Package CALIBERrfimpute June 11, 2018 Title Multiple Imputation Using MICE and Random Forest Version 1.0-1 Date 2018-06-05 Functions to impute using Random Forest under Full Conditional Specifications
More informationIn this computer exercise we will work with the analysis of variance in R. We ll take a look at the following topics:
UPPSALA UNIVERSITY Department of Mathematics Måns Thulin, thulin@math.uu.se Analysis of regression and variance Fall 2011 COMPUTER EXERCISE 2: One-way ANOVA In this computer exercise we will work with
More informationWeek 4: Simple Linear Regression III
Week 4: Simple Linear Regression III Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Goodness of
More informationRolling Markov Chain Monte Carlo
Rolling Markov Chain Monte Carlo Din-Houn Lau Imperial College London Joint work with Axel Gandy 4 th July 2013 Predict final ranks of the each team. Updates quick update of predictions. Accuracy control
More informationSome methods for the quantification of prediction uncertainties for digital soil mapping: Universal kriging prediction variance.
Some methods for the quantification of prediction uncertainties for digital soil mapping: Universal kriging prediction variance. Soil Security Laboratory 2018 1 Universal kriging prediction variance In
More informationApplied Regression Modeling: A Business Approach
i Applied Regression Modeling: A Business Approach Computer software help: SAS SAS (originally Statistical Analysis Software ) is a commercial statistical software package based on a powerful programming
More informationThe partial Package. R topics documented: October 16, Version 0.1. Date Title partial package. Author Andrea Lehnert-Batar
The partial Package October 16, 2006 Version 0.1 Date 2006-09-21 Title partial package Author Andrea Lehnert-Batar Maintainer Andrea Lehnert-Batar Depends R (>= 2.0.1),e1071
More informationLAB #2: SAMPLING, SAMPLING DISTRIBUTIONS, AND THE CLT
NAVAL POSTGRADUATE SCHOOL LAB #2: SAMPLING, SAMPLING DISTRIBUTIONS, AND THE CLT Statistics (OA3102) Lab #2: Sampling, Sampling Distributions, and the Central Limit Theorem Goal: Use R to demonstrate sampling
More informationEvaluation Metrics. (Classifiers) CS229 Section Anand Avati
Evaluation Metrics (Classifiers) CS Section Anand Avati Topics Why? Binary classifiers Metrics Rank view Thresholding Confusion Matrix Point metrics: Accuracy, Precision, Recall / Sensitivity, Specificity,
More informationDesign for sensitivity analysis, in Chapman and Hall Handbook of Design of Experiments
Design for sensitivity analysis, in Chapman and Hall Handbook of Design of Experiments Andrea Saltelli, William Becker Joint Research Centre, Institute for the Protection and Security of the Citizen, Ispra
More informationPackage simtool. August 29, 2016
Type Package Package simtool August 29, 2016 Title Conduct Simulation Studies with a Minimal Amount of Source Code. Version 1.0.3 Date 2014-09-07 Author Marsel Scheer Maintainer
More informationAn introduction to design of computer experiments
An introduction to design of computer experiments Derek Bingham Statistics and Actuarial Science Simon Fraser University Department of Statistics and Actuarial Science Outline Designs for estimating a
More informationThe ctest Package. January 3, 2000
R objects documented: The ctest Package January 3, 2000 bartlett.test....................................... 1 binom.test........................................ 2 cor.test.........................................
More informationExpectation-Maximization Methods in Population Analysis. Robert J. Bauer, Ph.D. ICON plc.
Expectation-Maximization Methods in Population Analysis Robert J. Bauer, Ph.D. ICON plc. 1 Objective The objective of this tutorial is to briefly describe the statistical basis of Expectation-Maximization
More informationTopics in Machine Learning
Topics in Machine Learning Gilad Lerman School of Mathematics University of Minnesota Text/slides stolen from G. James, D. Witten, T. Hastie, R. Tibshirani and A. Ng Machine Learning - Motivation Arthur
More informationThe lhs Package. October 22, 2006
The lhs Package October 22, 2006 Type Package Title Latin Hypercube Samples Version 0.3 Date 2006-10-21 Author Maintainer Depends R (>= 2.0.1) This package
More informationPackage intccr. September 12, 2017
Type Package Package intccr September 12, 2017 Title Semiparametric Competing Risks Regression under Interval Censoring Version 0.2.0 Author Giorgos Bakoyannis , Jun Park
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 13 Random Numbers and Stochastic Simulation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright
More informationThe theory of the linear model 41. Theorem 2.5. Under the strong assumptions A3 and A5 and the hypothesis that
The theory of the linear model 41 Theorem 2.5. Under the strong assumptions A3 and A5 and the hypothesis that E(Y X) =X 0 b 0 0 the F-test statistic follows an F-distribution with (p p 0, n p) degrees
More informationPackage actyr. August 5, 2015
Type Package Package actyr August 5, 2015 Title Replication Package for Abbring, Campbell, Tilly, Yang (2014): Very Simple Markov Perfect Industry Dynamics Version 1.0 Date 2015-08-01 Author Jan Tilly
More informationPsychology 282 Lecture #21 Outline Categorical IVs in MLR: Effects Coding and Contrast Coding
Psychology 282 Lecture #21 Outline Categorical IVs in MLR: Effects Coding and Contrast Coding In the previous lecture we learned how to incorporate a categorical research factor into a MLR model by using
More informationSensitivity and Uncertainty Analysis of Expected Loss Costs in Hurricane Models
Florida International University FIU Digital Commons FIU Electronic Theses and Dissertations University Graduate School 6-18-2012 Sensitivity and Uncertainty Analysis of Expected Loss Costs in Hurricane
More informationChapter 3. Bootstrap. 3.1 Introduction. 3.2 The general idea
Chapter 3 Bootstrap 3.1 Introduction The estimation of parameters in probability distributions is a basic problem in statistics that one tends to encounter already during the very first course on the subject.
More informationMS in Applied Statistics: Study Guide for the Data Science concentration Comprehensive Examination. 1. MAT 456 Applied Regression Analysis
MS in Applied Statistics: Study Guide for the Data Science concentration Comprehensive Examination. The Part II comprehensive examination is a three-hour closed-book exam that is offered on the second
More informationSIMULATION METHOD TO ESTIMATE THE UNCERTAINTIES OF ISO SPECIFICATIONS
SIMULATION METHOD TO ESTIMATE THE UNCERTAINTIES OF ISO SPECIFICATIONS Jean-Marc Linares, Jean Michel Sprauel To cite this version: Jean-Marc Linares, Jean Michel Sprauel. SIMULATION METHOD TO ESTIMATE
More informationPackage texteffect. November 27, 2017
Version 0.1 Date 2017-11-27 Package texteffect November 27, 2017 Title Discovering Latent Treatments in Text Corpora and Estimating Their Causal Effects Author Christian Fong
More informationPackage glmmml. R topics documented: March 25, Encoding UTF-8 Version Date Title Generalized Linear Models with Clustering
Encoding UTF-8 Version 1.0.3 Date 2018-03-25 Title Generalized Linear Models with Clustering Package glmmml March 25, 2018 Binomial and Poisson regression for clustered data, fixed and random effects with
More informationCOPYRIGHTED MATERIAL CONTENTS
PREFACE ACKNOWLEDGMENTS LIST OF TABLES xi xv xvii 1 INTRODUCTION 1 1.1 Historical Background 1 1.2 Definition and Relationship to the Delta Method and Other Resampling Methods 3 1.2.1 Jackknife 6 1.2.2
More informationCalibration and emulation of TIE-GCM
Calibration and emulation of TIE-GCM Serge Guillas School of Mathematics Georgia Institute of Technology Jonathan Rougier University of Bristol Big Thanks to Crystal Linkletter (SFU-SAMSI summer school)
More informationRolling Markov Chain Monte Carlo
Rolling Markov Chain Monte Carlo Din-Houn Lau Imperial College London Joint work with Axel Gandy 4 th September 2013 RSS Conference 2013: Newcastle Output predicted final ranks of the each team. Updates
More informationSTAT 311 (3 CREDITS) VARIANCE AND REGRESSION ANALYSIS ELECTIVE: ALL STUDENTS. CONTENT Introduction to Computer application of variance and regression
STAT 311 (3 CREDITS) VARIANCE AND REGRESSION ANALYSIS ELECTIVE: ALL STUDENTS. CONTENT Introduction to Computer application of variance and regression analysis. Analysis of Variance: one way classification,
More information