UP School of Statistics Student Council Education and Research
|
|
- Beverly Chase
- 6 years ago
- Views:
Transcription
1 w UP School of Statistics Student Council Education and Research erho.weebly.com 0 erhomyhero@gmail.com f /erhoismyhero S133_HOA_001 Statistics 133 Bayesian Statistical Inference Use of R for Bayesian Analysis 1. Let X 1, X 2,, X n ~Be ( p) and suppose p ~Beta(1, 1). Find the MSE of the following estimators for the target parameter p : a. ^p MLE b. ^p B Note: Some lines of code may appear as two lines in this handout but they must be coded in R as one line #Create a sequence containing values of p for plotting p<-seq(from=0,to=1,by=0.01) #Set n to a value of 10 n<-10 #From the formula of the MSE of the MLE derived above mse.mle<-p*(1-p)/n #From the formula of the MSE of the Bayes estimator derived above mse.bayes<-n*p*(1-p)/(n+2)^2+((1-2*p)/(n+2))^2 Plotting Codes are displayed chronologically a. Method 1 plot(x=p,y=mse.mle,type='l',col='1',lwd="2",xlab="p",ylab="ms E",las=1) i. x = specify values on x-axis ii. y = specify values on y-axis iii. type = specify plot type (type?plot in R to see different plot types). Specified above is l (small letter L) indicating line type plot iv. lwd = specifiy line thickness v. col = specifies color of plotting vi. xlab = label of x-axis 1
2 vii. ylab = label of y-axis viii. las = 1 makes the y-values in the y-axis appear horizontally lines(x=p,y=mse.bayes,col='2',lwd="2",lty=2) lines function allows one to superimpose a plot to a pre-existing plot. lty = specifies the line type. lty = 2 means the dashed line is used type?par and look for lty to explore different line types Resulting plot b. Method 2 plot(x=p,y=mse.mle,type="l",axes="f",xlab="p",ylab="mse",lwd= "2") Axes = F means that the plot will appear without x and y axes as we will customize the tick marks of the axes to be shown below lines(x=p,y=mse.bayes,col='2',lwd="2",lty=2) This is the same lines function as the first method. This will superimpose the Bayes MSE graph to the MLE MSE graph box() axis(side=1,at=seq(from=0,to=1,by=0.1)) axis(side=2,at=seq(from=0,to=0.05,by=0.005),las=1) i. box() function wraps the plot around a box since the axes = F option in plot removes the box. ii. axis function creates customizable axes 2
3 iii. side=1 creates the x-axis or places tick marks on the bottom of the box and side=2 creates the y- axis or places tick marks on the left of the box. (side = 3 and side = 4 puts tick marks on the top and right, respectively) iv. at option defines the sequence of numbers we wish to put in the axes Resulting plot We now have more values in the x-axis compared to the previous plot. c. Adding legends (Optional) legend(x=0.35,y=0.01,legend=c("mle","bayes"),col=c(1,2),lty=c (1,2),lwd="1") i. place the legend function after the plotting functions such as plot, lines, box and axes ii. x = and y = specifies the coordinates of the top left of the legend box iii. legend = specifies the legend labels iv. col = specifies the color of the legend lines v. lty = specifies the line type of the legend lines vi. lwd = specifies the thickness of the legend lines vii. legend label MSE will have color 1 and lty 1 while label Bayes will have color 2 and lty 2 Using plotting method 2, we add labels p<-seq(from=0,to=1,by=0.01) n<-10 mse.mle<-p*(1-p)/n mse.bayes<-n*p*(1-p)/(n+2)^2+((1-2*p)/(n+2))^2 plot(x=p,y=mse.mle,type="l",axes="f",xlab="p",ylab="mse",lwd= "2") lines(x=p,y=mse.bayes,col='2',lwd="2",lty=2) box() axis(side=1,at=seq(from=0,to=1,by=0.1)) 3
4 axis(side=2,at=seq(from=0,to=0.05,by=0.005),las=1) legend(x=0.35,y=0.01,legend=c("mle","bayes"),col=c(1,2),lty=c (1,2),lwd="1") Resulting plot 2. Let X 1, X 2,, X n ~ Po(θ ) and suppose θ ~Gamma(r, λ ). Find the MSE of the following estimators for the target parameter θ : a. ^θ MLE b. ^θ B Suppose we use Ga(2, 1) as prior Codes theta<-seq(from=0,to=5,by=0.5) n<-6 lambda<-1 r<-2 mse.mle<-theta/n mse.bayes<-n*theta/(lambda+n)^2+((r-(lambda*theta))/ (lambda+n))^2 plot(x=theta,y=mse.bayes,type="l",lwd=2,col=1,xlab='theta',ylab= 'mse',las=1) lines(x=theta,y=mse.mle,lwd=2,col=1,lty=2) legend(x=0.1,y=0.7,legend=c("bayes","frequentist"),lty=c(1,2),lw d=2) 4
5 Resulting plot Bayesian Confidence Interval Estimation Note: Check if the following packages are installed: TeachingDemos, coda, VGAM Use of the Special Parametric Family of Distributions dbeta(x,shape1,shape2) returns the density or point in the graph provided the quantile and parameters x specifies quantile where we want to find the density shape1 specifies alpha parameter in Beta distribution shape2 specifies beta parameter in Beta distribution pbeta(q,shape1,shape2,lower.tail = TRUE) returns the probability or area to the left of the specified quantile q specifies quantile qbeta(p,shape1,shape2,lower.tail = TRUE) the inverse of the pbeta function. Returns the quantile for the specified left-tail probability p specifies the lower tail probability rbeta(n,shape1,shape2) Creates a random sample of size n for a Beta distribution with paramters shape1 and shape2 dnorm, pnorm, qnorm, rnorm are used for the normal distribution dgamma, pgamma, qgamma, rgamma are used for the gamma distribution 5
6 Graphical Presentation of the Special Parametric Distributions Consider a normal distribution with mean 0 and variance 1. dnorm(x = 1, mean = 0, sd = 1) Result: dnorm returns this value ( ) Quantile (x) = 1 pnorm(q = 1, mean = 0, sd =1, lower.tail = TRUE) Result: pnorm returns this area ( ) Quantile (q) = 1 6
7 qnorm(p = 0.95, mean = 0, sd = 1, lower.tail = TRUE) Result: (approx ) p = 0.95 qnorm returns this value ( ) Equal Tail Credible Intervals Suppose that the probability of contracting a disease has a Beta(19, 133) posterior distribution. Give the 95% equal-tail credible interval for the estimate of the probability of contracting a disease. qbeta(p = c(0.025,0.975), shape1 = 19, shape2 = 133) HPD Credible Intervals For the same example above hpd(posterior.icdf = qbeta,shape1 = 19, shape2 = 133, conf = 0.95) Note: hpd function is under the TeachingDemos package 7
8 Monte Carlo Integration Syntax mc.integration<-function(lower,upper,incr,n,term){ temp.upper<-lower+incr sum=0 while(temp.upper<=upper){ ran.values<-runif(n = n, min = lower, max = temp.upper) sum=sum + mean(term(ran.values))*incr lower<-lower + incr temp.upper<-temp.upper+incr } return(sum) } When to use? If the prior is not a conjugate distribution of the process or if the posterior is not a special distribution Monte Carlo HPD Credible Interval R function: HPDinterval(as.mcmc(),prob) HPDinterval function is under coda package Requires a random sample from the posterior distribution. A vector containing the random sample is placed inside the parentheses of the as.mcmc() argument prob specifies level of confidence Monte Carlo Equal Tail Credible Interval quantile(x,probs) x specifies random sample probs specifies tail probabilities Example: Approximate credible intervals for a Beta(19, 133) posterior distribution set.seed(22) sample<-rbeta(n = 1000,shape1 = 19,shape2 = 133) #Equal tail quantile(x = sample,probs = c(0.025,0.975) #HPD HPDinterval(as.mcmc(sample),prob = 0.95)
The GSM Package. August 9, 2007
The GSM Package August 9, 2007 Title Gamma Shape Mixture This package implements a Bayesian approach for estimation of a mixture of gamma distributions in which the mixing occurs over the shape parameter.
More information1. Start WinBUGS by double clicking on the WinBUGS icon (or double click on the file WinBUGS14.exe in the WinBUGS14 directory in C:\Program Files).
Hints on using WinBUGS 1 Running a model in WinBUGS 1. Start WinBUGS by double clicking on the WinBUGS icon (or double click on the file WinBUGS14.exe in the WinBUGS14 directory in C:\Program Files). 2.
More informationPackage HDInterval. June 9, 2018
Type Package Title Highest (Posterior) Density Intervals Version 0.2.0 Date 2018-06-09 Suggests coda Author Mike Meredith and John Kruschke Package HDInterval June 9, 2018 Maintainer Mike Meredith
More informationPackage simed. November 27, 2017
Version 1.0.3 Title Simulation Education Author Barry Lawson, Larry Leemis Package simed November 27, 2017 Maintainer Barry Lawson Imports graphics, grdevices, methods, stats, utils
More informationBayesian Methods. David Rosenberg. April 11, New York University. David Rosenberg (New York University) DS-GA 1003 April 11, / 19
Bayesian Methods David Rosenberg New York University April 11, 2017 David Rosenberg (New York University) DS-GA 1003 April 11, 2017 1 / 19 Classical Statistics Classical Statistics David Rosenberg (New
More information1 Pencil and Paper stuff
Spring 2008 - Stat C141/ Bioeng C141 - Statistics for Bioinformatics Course Website: http://www.stat.berkeley.edu/users/hhuang/141c-2008.html Section Website: http://www.stat.berkeley.edu/users/mgoldman
More informationLinear Modeling with Bayesian Statistics
Linear Modeling with Bayesian Statistics Bayesian Approach I I I I I Estimate probability of a parameter State degree of believe in specific parameter values Evaluate probability of hypothesis given the
More informationMarkov Chain Monte Carlo (part 1)
Markov Chain Monte Carlo (part 1) Edps 590BAY Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Spring 2018 Depending on the book that you select for
More informationR is a programming language of a higher-level Constantly increasing amount of packages (new research) Free of charge Website:
Introduction to R R R is a programming language of a higher-level Constantly increasing amount of packages (new research) Free of charge Website: http://www.r-project.org/ Code Editor: http://rstudio.org/
More informationSTATISTICAL LABORATORY, April 30th, 2010 BIVARIATE PROBABILITY DISTRIBUTIONS
STATISTICAL LABORATORY, April 3th, 21 BIVARIATE PROBABILITY DISTRIBUTIONS Mario Romanazzi 1 MULTINOMIAL DISTRIBUTION Ex1 Three players play 1 independent rounds of a game, and each player has probability
More informationLecture 3 - Object-oriented programming and statistical programming examples
Lecture 3 - Object-oriented programming and statistical programming examples Björn Andersson (w/ Ronnie Pingel) Department of Statistics, Uppsala University February 1, 2013 Table of Contents 1 Some notes
More informationPackage SHELF. August 16, 2018
Type Package Package SHELF August 16, 2018 Title Tools to Support the Sheffield Elicitation Framework Version 1.4.0 Date 2018-08-15 Author Jeremy Oakley Maintainer Implements various methods for eliciting
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 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 informationThe nor1mix Package. June 12, 2007
The nor1mix Package June 12, 2007 Title Normal (1-d) Mixture Models (S3 Classes and Methods) Version 1.0-7 Date 2007-03-15 Author Martin Mächler Maintainer Martin Maechler
More informationLab 4: Distributions of random variables
Lab 4: Distributions of random variables In this lab we ll investigate the probability distribution that is most central to statistics: the normal distribution If we are confident that our data are nearly
More informationPackage DPBBM. September 29, 2016
Type Package Title Dirichlet Process Beta-Binomial Mixture Version 0.2.5 Date 2016-09-21 Author Lin Zhang Package DPBBM September 29, 2016 Maintainer Lin Zhang Depends R (>= 3.1.0)
More informationPackage sspse. August 26, 2018
Type Package Version 0.6 Date 2018-08-24 Package sspse August 26, 2018 Title Estimating Hidden Population Size using Respondent Driven Sampling Data Maintainer Mark S. Handcock
More informationCommon Sta 101 Commands for R. 1 One quantitative variable. 2 One categorical variable. 3 Two categorical variables. Summary statistics
Common Sta 101 Commands for R 1 One quantitative variable summary(x) # most summary statitstics at once mean(x) median(x) sd(x) hist(x) boxplot(x) # horizontal = TRUE for horizontal plot qqnorm(x) qqline(x)
More informationBayesian Analysis of Extended Lomax Distribution
Bayesian Analysis of Extended Lomax Distribution Shankar Kumar Shrestha and Vijay Kumar 2 Public Youth Campus, Tribhuvan University, Nepal 2 Department of Mathematics and Statistics DDU Gorakhpur University,
More informationPackage ihs. February 25, 2015
Version 1.0 Type Package Title Inverse Hyperbolic Sine Distribution Date 2015-02-24 Author Carter Davis Package ihs February 25, 2015 Maintainer Carter Davis Depends R (>= 2.4.0),
More informationThe nor1mix Package. August 3, 2006
The nor1mix Package August 3, 2006 Title Normal (1-d) Mixture Models (S3 Classes and Methods) Version 1.0-6 Date 2006-08-02 Author: Martin Mächler Maintainer Martin Maechler
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 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 informationLaboratorio di Problemi Inversi Esercitazione 4: metodi Bayesiani e importance sampling
Laboratorio di Problemi Inversi Esercitazione 4: metodi Bayesiani e importance sampling Luca Calatroni Dipartimento di Matematica, Universitá degli studi di Genova May 19, 2016. Luca Calatroni (DIMA, Unige)
More informationLogistic Regression. (Dichotomous predicted variable) Tim Frasier
Logistic Regression (Dichotomous predicted variable) Tim Frasier Copyright Tim Frasier This work is licensed under the Creative Commons Attribution 4.0 International license. Click here for more information.
More informationCommon R commands used in Data Analysis and Statistical Inference
Common R commands used in Data Analysis and Statistical Inference 1 One numerical variable summary(x) # most summary statitstics at once mean(x) median(x) sd(x) hist(x) boxplot(x) # horizontal = TRUE for
More informationUnivariate Data - 2. Numeric Summaries
Univariate Data - 2. Numeric Summaries Young W. Lim 2018-08-01 Mon Young W. Lim Univariate Data - 2. Numeric Summaries 2018-08-01 Mon 1 / 36 Outline 1 Univariate Data Based on Numerical Summaries R Numeric
More informationPackage visualize. April 28, 2017
Type Package Package visualize April 28, 2017 Title Graph Probability Distributions with User Supplied Parameters and Statistics Version 4.3.0 Date 2017-04-27 Depends R (>= 3.0.0) Graphs the pdf or pmf
More informationTutorial 3: Probability & Distributions Johannes Karreth RPOS 517, Day 3
Tutorial 3: Probability & Distributions Johannes Karreth RPOS 517, Day 3 This tutorial shows you: how to simulate a random process how to plot the distribution of a variable how to assess the distribution
More informationJournal of Statistical Software
JSS Journal of Statistical Software December 2007, Volume 23, Issue 9. http://www.jstatsoft.org/ WinBUGSio: A SAS Macro for the Remote Execution of WinBUGS Michael K. Smith Pfizer Global Research and Development
More informationPackage OrthoPanels. November 11, 2016
Package OrthoPanels November 11, 2016 Title Dynamic Panel Models with Orthogonal Reparameterization of Fixed Effects Version 1.1-0 Implements the orthogonal reparameterization approach recommended by Lancaster
More informationBayesian Statistics Group 8th March Slice samplers. (A very brief introduction) The basic idea
Bayesian Statistics Group 8th March 2000 Slice samplers (A very brief introduction) The basic idea lacements To sample from a distribution, simply sample uniformly from the region under the density function
More informationPackage TBSSurvival. July 1, 2012
Package TBSSurvival July 1, 2012 Version 1.0 Date 2012-06-30 Title TBS Model R package Author Adriano Polpo , Cassio de Campos , D. Sinha , Stuart
More informationBayesian Estimation for Skew Normal Distributions Using Data Augmentation
The Korean Communications in Statistics Vol. 12 No. 2, 2005 pp. 323-333 Bayesian Estimation for Skew Normal Distributions Using Data Augmentation Hea-Jung Kim 1) Abstract In this paper, we develop a MCMC
More informationISyE8843A, Brani Vidakovic Handout 14
ISyE8843A, Brani Vidakovic Handout 4 BUGS BUGS is freely available software for constructing Bayesian statistical models and evaluating them using MCMC methodology. BUGS and WINBUGS are distributed freely
More informationA noninformative Bayesian approach to small area estimation
A noninformative Bayesian approach to small area estimation Glen Meeden School of Statistics University of Minnesota Minneapolis, MN 55455 glen@stat.umn.edu September 2001 Revised May 2002 Research supported
More informationPackage SCBmeanfd. December 27, 2016
Type Package Package SCBmeanfd December 27, 2016 Title Simultaneous Confidence Bands for the Mean of Functional Data Version 1.2.2 Date 2016-12-22 Author David Degras Maintainer David Degras
More informationPackage TBSSurvival. January 5, 2017
Version 1.3 Date 2017-01-05 Package TBSSurvival January 5, 2017 Title Survival Analysis using a Transform-Both-Sides Model Author Adriano Polpo , Cassio de Campos , D.
More informationPackage gaga. February 22, 2018
Version 2.24.0 Date 2015-06-14 Package gaga February 22, 2018 Title GaGa hierarchical model for high-throughput data analysis Author . Maintainer Depends
More informationBayesian Computation with JAGS
JAGS is Just Another Gibbs Sampler Cross-platform Accessible from within R Bayesian Computation with JAGS What I did Downloaded and installed JAGS. In the R package installer, downloaded rjags and dependencies.
More informationSimulation Extrapolation
Simulation Extrapolation mk:@msitstore:c:\progra~1\r\r-26~1.0\library\simex\chtml\simex.chm::/simex.html Page 1 of 3 simex(simex) Simulation Extrapolation Implementation of the SIMEX Algorithm for measurement
More informationPackage mmeta. R topics documented: March 28, 2017
Package mmeta March 28, 2017 Type Package Title Multivariate Meta-Analysis Version 2.3 Date 2017-3-26 Author Sheng Luo, Yong Chen, Xiao Su, Haitao Chu Maintainer Xiao Su Description
More informationPackage basictrendline
Version 2.0.3 Date 2018-07-26 Package basictrendline July 26, 2018 Title Add Trendline and Confidence Interval of Basic Regression Models to Plot Maintainer Weiping Mei Plot, draw
More informationUnivariate Data - 2. Numeric Summaries
Univariate Data - 2. Numeric Summaries Young W. Lim 2018-02-05 Mon Young W. Lim Univariate Data - 2. Numeric Summaries 2018-02-05 Mon 1 / 31 Outline 1 Univariate Data Based on Numerical Summaries Young
More informationA Handbook of Statistical Analyses Using R. Brian S. Everitt and Torsten Hothorn
A Handbook of Statistical Analyses Using R Brian S. Everitt and Torsten Hothorn CHAPTER 7 Density Estimation: Erupting Geysers and Star Clusters 7.1 Introduction 7.2 Density Estimation The three kernel
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 informationGenerating pseudo Normal Numbers : Ziggurat, Box Muller and R default rnorm method
Generating pseudo Normal Numbers : Ziggurat, Box Muller and R default rnorm method 1 Introduction Mriganka Aulia ISI There is many methods for generating pseudo Normal Random Numbers. Here we focus on
More informationMatrix algebra. Basics
Matrix.1 Matrix algebra Matrix algebra is very prevalently used in Statistics because it provides representations of models and computations in a much simpler manner than without its use. The purpose of
More informationThe R statistical computing environment
The R statistical computing environment Luke Tierney Department of Statistics & Actuarial Science University of Iowa June 17, 2011 Luke Tierney (U. of Iowa) R June 17, 2011 1 / 27 Introduction R is a language
More informationComputation for the Introduction to MCMC Chapter of Handbook of Markov Chain Monte Carlo By Charles J. Geyer Technical Report No.
Computation for the Introduction to MCMC Chapter of Handbook of Markov Chain Monte Carlo By Charles J. Geyer Technical Report No. 679 School of Statistics University of Minnesota July 29, 2010 Abstract
More informationContinuous-time stochastic simulation of epidemics in R
Continuous-time stochastic simulation of epidemics in R Ben Bolker May 16, 2005 1 Introduction/basic code Using the Gillespie algorithm, which assumes that all the possible events that can occur (death
More informationConfidence Intervals: Estimators
Confidence Intervals: Estimators Point Estimate: a specific value at estimates a parameter e.g., best estimator of e population mean ( ) is a sample mean problem is at ere is no way to determine how close
More informationSimulating from the Polya posterior by Glen Meeden, March 06
1 Introduction Simulating from the Polya posterior by Glen Meeden, glen@stat.umn.edu March 06 The Polya posterior is an objective Bayesian approach to finite population sampling. In its simplest form it
More informationPackage samplesizelogisticcasecontrol
Package samplesizelogisticcasecontrol February 4, 2017 Title Sample Size Calculations for Case-Control Studies Version 0.0.6 Date 2017-01-31 Author Mitchell H. Gail To determine sample size for case-control
More informationJAGS for Rats. Jonathan Rougier School of Mathematics University of Bristol UK. Version 2.0, compiled February 23, 2017
JAGS for Rats Jonathan Rougier School of Mathematics University of Bristol UK Version 2.0, compiled February 23, 2017 1 Introduction This is a worksheet about JAGS, using rjags. These must both be installed
More informationSemiparametric Mixed Effecs with Hierarchical DP Mixture
Semiparametric Mixed Effecs with Hierarchical DP Mixture R topics documented: April 21, 2007 hdpm-package........................................ 1 hdpm............................................ 2 hdpmfitsetup........................................
More informationIntroduction to the Practice of Statistics using R: Chapter 6
Introduction to the Practice of Statistics using R: Chapter 6 Ben Baumer Nicholas J. Horton March 10, 2013 Contents 1 Estimating with Confidence 2 1.1 Beyond the Basics.....................................
More informationGraphs and transformations 4G
Graphs and transformations 4G a f(x + ) is a translation by one unit to the left. d A (0, ), B ( ),0, C (, 4), D (, 0) A (, ), B (0, 0), C (, 4), D (5, 0) e f(x) is a stretch with scale factor b f(x) 4
More informationChapter 5. Data:Distribution
Chapter 5 Data:Distribution What you will learn in this chapter: How to create histograms and other graphics of sample distribution How to examine various distributions How to test for the normal distribution
More informationIST 3108 Data Analysis and Graphics Using R Week 9
IST 3108 Data Analysis and Graphics Using R Week 9 Engin YILDIZTEPE, Ph.D 2017-Spring Introduction to Graphics >y plot (y) In R, pictures are presented in the active graphical device or window.
More informationComputational statistics Jamie Griffin. Semester B 2018 Lecture 1
Computational statistics Jamie Griffin Semester B 2018 Lecture 1 Course overview This course is not: Statistical computing Programming This course is: Computational statistics Statistical methods that
More informationBUGS: Language, engines, and interfaces
BUGS: Language, engines, and interfaces Patrick Breheny January 17 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/18 The BUGS framework The BUGS project (Bayesian inference using Gibbs Sampling)
More informationSaratoga Springs City School District/Office of Continuing Education Introduction to Microsoft Excel 04 Charts. 1. Chart Types and Dimensions
1949 1954 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 Saratoga Springs City School District/Office of Continuing Education Introduction to Microsoft Excel 04 Charts 1. Chart Types and Dimensions
More informationPackage truncgof. R topics documented: February 20, 2015
Type Package Title GoF tests allowing for left truncated data Version 0.6-0 Date 2012-12-24 Author Thomas Wolter Package truncgof February 20, 2015 Maintainer Renato Vitolo
More informationA Basic Example of ANOVA in JAGS Joel S Steele
A Basic Example of ANOVA in JAGS Joel S Steele The purpose This demonstration is intended to show how a simple one-way ANOVA can be coded and run in the JAGS framework. This is by no means an exhaustive
More information6. More Loops, Control Structures, and Bootstrapping
6. More Loops, Control Structures, and Bootstrapping Ken Rice Timothy Thornotn University of Washington Seattle, July 2013 In this session We will introduce additional looping procedures as well as control
More informationPackage fitur. March 11, 2018
Title Fit Univariate Distributions Version 0.5.25 Package fitur March 11, 2018 Wrapper for computing parameters for univariate distributions using MLE. It creates an object that stores d, p, q, r functions
More informationDPP: Reference documentation. version Luis M. Avila, Mike R. May and Jeffrey Ross-Ibarra 17th November 2017
DPP: Reference documentation version 0.1.2 Luis M. Avila, Mike R. May and Jeffrey Ross-Ibarra 17th November 2017 1 Contents 1 DPP: Introduction 3 2 DPP: Classes and methods 3 2.1 Class: NormalModel.........................
More informationA GENERAL GIBBS SAMPLING ALGORITHM FOR ANALYZING LINEAR MODELS USING THE SAS SYSTEM
A GENERAL GIBBS SAMPLING ALGORITHM FOR ANALYZING LINEAR MODELS USING THE SAS SYSTEM Jayawant Mandrekar, Daniel J. Sargent, Paul J. Novotny, Jeff A. Sloan Mayo Clinic, Rochester, MN 55905 ABSTRACT A general
More informationLab 7: Bayesian analysis of a dice toss problem using C++ instead of python
Lab 7: Bayesian analysis of a dice toss problem using C++ instead of python Due date: Monday March 27, 11:59pm Short version of the assignment Take your python file from lab 6 and convert it into lab7
More informationPackage beast. March 16, 2018
Type Package Package beast March 16, 2018 Title Bayesian Estimation of Change-Points in the Slope of Multivariate Time-Series Version 1.1 Date 2018-03-16 Author Maintainer Assume that
More informationPackage bacon. October 31, 2018
Type Package Package October 31, 2018 Title Controlling bias and inflation in association studies using the empirical null distribution Version 1.10.0 Author Maarten van Iterson [aut, cre], Erik van Zwet
More informationUnit #13 : Integration to Find Areas and Volumes, Volumes of Revolution
Unit #13 : Integration to Find Areas and Volumes, Volumes of Revolution Goals: Beabletoapplyaslicingapproachtoconstructintegralsforareasandvolumes. Be able to visualize surfaces generated by rotating functions
More informationPackage CoImp. August 8, 2016
Title Copula Based Imputation Method Date 2016-07-31 Version 0.3-1 Package CoImp August 8, 2016 Author Francesca Marta Lilja Di Lascio, Simone Giannerini Depends R (>= 3.1.0), methods, copula Imports nnet,
More informationPackage robets. March 6, Type Package
Type Package Package robets March 6, 2018 Title Forecasting Time Series with Robust Exponential Smoothing Version 1.4 Date 2018-03-06 We provide an outlier robust alternative of the function ets() in the
More informationStatistics 202: Data Mining. c Jonathan Taylor. Outliers Based in part on slides from textbook, slides of Susan Holmes.
Outliers Based in part on slides from textbook, slides of Susan Holmes December 2, 2012 1 / 1 Concepts What is an outlier? The set of data points that are considerably different than the remainder of the
More informationUnivariate Extreme Value Analysis. 1 Block Maxima. Practice problems using the extremes ( 2.0 5) package. 1. Pearson Type III distribution
Univariate Extreme Value Analysis Practice problems using the extremes ( 2.0 5) package. 1 Block Maxima 1. Pearson Type III distribution (a) Simulate 100 maxima from samples of size 1000 from the gamma
More informationOne Factor Experiments
One Factor Experiments 20-1 Overview Computation of Effects Estimating Experimental Errors Allocation of Variation ANOVA Table and F-Test Visual Diagnostic Tests Confidence Intervals For Effects Unequal
More informationChapter 5: Joint Probability Distributions and Random
Chapter 5: Joint Probability Distributions and Random Samples Curtis Miller 2018-06-13 Introduction We may naturally inquire about collections of random variables that are related to each other in some
More informationPackage SmoothHazard
Package SmoothHazard September 19, 2014 Title Fitting illness-death model for interval-censored data Version 1.2.3 Author Celia Touraine, Pierre Joly, Thomas A. Gerds SmoothHazard is a package for fitting
More informationPackage evmix. February 9, 2018
Package evmix February 9, 2018 Title Extreme Value Mixture Modelling, Threshold Estimation and Boundary Corrected Kernel Density Estimation Version 2.9 Date 2018-02-08 Author Carl Scarrott, Yang Hu and
More informationPackage logitnorm. R topics documented: February 20, 2015
Title Functions for the al distribution. Version 0.8.29 Date 2012-08-24 Author Package February 20, 2015 Maintainer Density, distribution, quantile and random generation function
More informationS Basics. Statistics 135. Autumn Copyright c 2005 by Mark E. Irwin
S Basics Statistics 135 Autumn 2005 Copyright c 2005 by Mark E. Irwin S Basics When discussing the S environment, I will, or at least try to, make the following distinctions. S will be used when what is
More informationAz R adatelemzési nyelv
Az R adatelemzési nyelv alapjai II. Egészségügyi informatika és biostatisztika Gézsi András gezsi@mit.bme.hu Functions Functions Functions do things with data Input : function arguments (0,1,2, ) Output
More informationGCSE CCEA GCSE EXCEL 2010 USER GUIDE. Business and Communication Systems
GCSE CCEA GCSE EXCEL 2010 USER GUIDE Business and Communication Systems For first teaching from September 2017 Contents Page Define the purpose and uses of a spreadsheet... 3 Define a column, row, and
More informationCREATING THE DISTRIBUTION ANALYSIS
Chapter 12 Examining Distributions Chapter Table of Contents CREATING THE DISTRIBUTION ANALYSIS...176 BoxPlot...178 Histogram...180 Moments and Quantiles Tables...... 183 ADDING DENSITY ESTIMATES...184
More informationMachine Learning
Machine Learning 10-701 Tom M. Mitchell Machine Learning Department Carnegie Mellon University February 15, 2011 Today: Graphical models Inference Conditional independence and D-separation Learning from
More informationCluster Randomization Create Cluster Means Dataset
Chapter 270 Cluster Randomization Create Cluster Means Dataset Introduction A cluster randomization trial occurs when whole groups or clusters of individuals are treated together. Examples of such clusters
More informationR Command Summary. Steve Ambler Département des sciences économiques École des sciences de la gestion. April 2018
R Command Summary Steve Ambler Département des sciences économiques École des sciences de la gestion Université du Québec à Montréal c 2018 : Steve Ambler April 2018 This document describes some of the
More informationPackage Kernelheaping
Type Package Package Kernelheaping December 7, 2015 Title Kernel Density Estimation for Heaped and Rounded Data Version 1.2 Date 2015-12-01 Depends R (>= 2.15.0), evmix, MASS, ks, sparr Author Marcus Gross
More informationIntroduction to RStudio
First, take class through processes of: Signing in Changing password: Tools -> Shell, then use passwd command Installing packages Check that at least these are installed: MASS, ISLR, car, class, boot,
More informationSurvival Analysis: Exercises (Part 3)
Survival Analysis: Exercises (Part 3) This is the third part of a series surrounding survival analysis. For part 1, click here. For part 2, click here. This particular post follows the final part of fitting
More informationYou are free to use this program, for non-commercial purposes only, under two conditions:
Bayesian sample size determination for prevalence and diagnostic studies in the absence of a gold standard Sample size calculations and asymptotic results (Version 5.10, June 2016) 1. Introduction The
More informationComputational Methods for Finding Probabilities of Intervals
Computational Methods for Finding Probabilities of Intervals Often the standard normal density function is denoted ϕ(z) = (2π) 1/2 exp( z 2 /2) and its CDF by Φ(z). Because Φ cannot be expressed in closed
More informationThe boa Package. May 3, 2005
The boa Package May 3, 2005 Version 1.1.5-2 Date 2005-05-02 Title Bayesian Output Analysis Program (BOA) for MCMC Author Maintainer Depends R (>= 1.7) A menu-driven program and
More informationStatistical Programming with R
Statistical Programming with R Lecture 9: Basic graphics in R Part 2 Bisher M. Iqelan biqelan@iugaza.edu.ps Department of Mathematics, Faculty of Science, The Islamic University of Gaza 2017-2018, Semester
More informationPlotting: An Iterative Process
Plotting: An Iterative Process Plotting is an iterative process. First we find a way to represent the data that focusses on the important aspects of the data. What is considered an important aspect may
More informationCSSS 510: Lab 2. Introduction to Maximum Likelihood Estimation
CSSS 510: Lab 2 Introduction to Maximum Likelihood Estimation 2018-10-12 0. Agenda 1. Housekeeping: simcf, tile 2. Questions about Homework 1 or lecture 3. Simulating heteroskedastic normal data 4. Fitting
More informationPackage GLDreg. February 28, 2017
Type Package Package GLDreg February 28, 2017 Title Fit GLD Regression Model and GLD Quantile Regression Model to Empirical Data Version 1.0.7 Date 2017-03-15 Author Steve Su, with contributions from:
More information