Density estimation II

Similar documents
CS 2750 Machine Learning Lecture 5. Density estimation. Density estimation

Stats & Summary

Centre for Efficiency and Productivity Analysis

Introduction to mathematical Statistics

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

Sequences and summations

Mathematics HL and further mathematics HL formula booklet

Differential Entropy 吳家麟教授

Mathematics HL and further mathematics HL formula booklet

Chapter Gauss-Seidel Method

Solutions Manual for Polymer Science and Technology Third Edition

Some Unbiased Classes of Estimators of Finite Population Mean

CS 3710 Advanced Topics in AI Lecture 17. Density estimation. CS 3710 Probabilistic graphical models. Administration

Outline. Outline. Outline. Questions 2010/9/30. Introduction The Multivariate Normal Density and Its Properties

Chapter 7. Bounds for weighted sums of Random Variables

Math 1313 Final Exam Review

Problem Set 4 Solutions

Mathematics HL and further mathematics HL formula booklet

Section 7.2 Two-way ANOVA with random effect(s)

Outline. Outline. Outline. Questions. Multivariate Normal Distribution. Multivariate Normal Distribution

ICS141: Discrete Mathematics for Computer Science I

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

Advanced Algorithmic Problem Solving Le 3 Arithmetic. Fredrik Heintz Dept of Computer and Information Science Linköping University

PubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS

Mathematical models for computer systems behaviour

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]

An Extended Mixture Inverse Gaussian Distribution

MATRIX AND VECTOR NORMS

Modeling uncertainty using probabilities

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

STA261H1.doc. i 1 X n be a random sample. The sample mean is defined by i= 1 X 1 + ( ) X has a N ( σ ) 1 n. N distribution. Then n. distribution.

β (cf Khan, 2006). In this model, p independent

Chapter 2 Intro to Math Techniques for Quantum Mechanics

under the curve in the first quadrant.

Model Fitting and Robust Regression Methods

Random variables and sampling theory

Department of Statistics, Dibrugarh University, Dibrugarh, Assam, India. Department of Statistics, G. C. College, Silchar, Assam, India.

A Series Illustrating Innovative Forms of the Organization & Exposition of Mathematics by Walter Gottschalk

Xidian University Liu Congfeng Page 1 of 22

Mathacle. PSet Stats, Concepts In Statistics Level Number Name: Date:

Preliminary Examinations: Upper V Mathematics Paper 1

2. Independence and Bernoulli Trials

this is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ]

COMPLEX NUMBERS AND DE MOIVRE S THEOREM

Predicting Survival Outcomes Based on Compound Covariate Method under Cox Proportional Hazard Models with Microarrays

University of California at Berkeley College of Engineering Dept. of Electrical Engineering and Computer Sciences.

Expectation and Variance

Inference on Parameters of a Shifted and Rescaled Wiener Process Based on Nonidentically Distributed Observations

19 22 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y Write the given sum or difference as a single integral in

20 23 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y 2 2

Cooper and McGillem Chapter 4: Moments Linear Regression

Taylor Polynomials. The Tangent Line. (a, f (a)) and has the same slope as the curve y = f (x) at that point. It is the best

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

Fast Fourier Transform 1) Legendre s Interpolation 2) Vandermonde Matrix 3) Roots of Unity 4) Polynomial Evaluation

MTH 146 Class 7 Notes

Lecture 3 Naïve Bayes, Maximum Entropy and Text Classification COSI 134

GCE AS and A Level MATHEMATICS FORMULA BOOKLET. From September Issued WJEC CBAC Ltd.

GCE AS/A Level MATHEMATICS GCE AS/A Level FURTHER MATHEMATICS

Parametric Yield Maximization using Gate Sizing based on Efficient Statistical Power and Delay Gradient Computation

Inference on One Population Mean Hypothesis Testing

Lecture 3-4 Solutions of System of Linear Equations

Matrix. Definition 1... a1 ... (i) where a. are real numbers. for i 1, 2,, m and j = 1, 2,, n (iii) A is called a square matrix if m n.

Students must always use correct mathematical notation, not calculator notation. the set of positive integers and zero, {0,1, 2, 3,...

Prior distributions. July 29, 2002

Machine learning: Density estimation

On Several Inequalities Deduced Using a Power Series Approach

Measures of Explained Variation and the Base-Rate Problem for Logistic Regression

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini

: : 8.2. Test About a Population Mean. STT 351 Hypotheses Testing Case I: A Normal Population with Known. - null hypothesis states 0

Chapter Unary Matrix Operations

CS473-Algorithms I. Lecture 3. Solving Recurrences. Cevdet Aykanat - Bilkent University Computer Engineering Department

Z = = = = X np n. n n. npq. npq pq

Chapter 4: Distributions

ME 501A Seminar in Engineering Analysis Page 1

3. REVIEW OF PROPERTIES OF EIGENVALUES AND EIGENVECTORS

Multiplicative Versions of Infinitesimal Calculus

Classification : Logistic regression. Generative classification model.

5. Lighting & Shading

Chapter Linear Regression

St John s College. UPPER V Mathematics: Paper 1 Learning Outcome 1 and 2. Examiner: GE Marks: 150 Moderator: BT / SLS INSTRUCTIONS AND INFORMATION

Optimality of Strategies for Collapsing Expanded Random Variables In a Simple Random Sample Ed Stanek

BASIC PRINCIPLES OF STATISTICS

CS 310 (sec 20) - Winter Final Exam (solutions) SOLUTIONS

Part I: Background on the Binomial Distribution

CURVE FITTING LEAST SQUARES METHOD

Generative classification models

10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n

The z-transform. LTI System description. Prof. Siripong Potisuk

Roberto s Notes on Integral Calculus Chapter 4: Definite integrals and the FTC Section 2. Riemann sums

Numerical Analysis Topic 4: Least Squares Curve Fitting

Homework Assignment 6 Solution Set

A Technique for Constructing Odd-order Magic Squares Using Basic Latin Squares

Level-2 BLAS. Matrix-Vector operations with O(n 2 ) operations (sequentially) BLAS-Notation: S --- single precision G E general matrix M V --- vector

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files)

Unsupervised Learning and Other Neural Networks

Variable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia

COMPLEX ANALYSIS AND PROBABILITY DISTRIBUTION

Modern Physics. Unit 6: Hydrogen Atom - Radiation Lecture 6.1: The Radial Probability Density. Ron Reifenberger Professor of Physics Purdue University

Transcription:

CS 750 Mche Lerg Lecture 6 esty estmto II Mlos Husrecht mlos@tt.edu 539 Seott Squre t: esty estmto {.. } vector of ttrute vlues Ojectve: estmte the model of the uderlyg rolty dstruto over vrles X X usg emles true dstruto smles X.. } { estmte ˆ X

esty estmto yes of desty estmto: rmetrc the dstruto s modeled usg set of rmeters ˆ X X Emle: me d covrces of multvrte orml Estmto: fd rmeters descrg dt o-rmetrc he model of the dstruto utlzes ll emles As f ll emles were rmeters of the dstruto Emles: erest-eghor Model ML rmeter estmto ˆ X X Θ Mmum lelhood ML rg m ML ML.. t log log {.. } Log lelhood hs the sme mmum s lelhood Ideedet emles

Byes rmeter estmto Byes rmeter estmto Uses the osteror dstruto for rmeters osteror covers ll ossle rmeter vlues & ther weghts rmeter osteror t Lelhood rmeter ror ror o rmeter + t + = Cojugte choces: ror dstruto mtches the dt dstruto osteror s the sme tye rmeter estmto Other crter: Mmum osteror rolty MA mode of the osteror Model: MA ˆ X X Θ rg m Θ Eected vlue of the rmeter ˆ X X ΘEV Eectto te wth regrd to osteror me of the osteror EV E d MA 3

Beroull dstruto Model for rdom vrle wth two outcomes Rdom vrle: wo outcomes: 0 or struto: where s the rolty of = Emle: Co toss Outcomes: Hed = l =0 rolty of Hed ror dstruto Choce of ror: Bet dstruto Bet osteror dstruto s g Bet dstruto Bet Bet CS 750 Mche Lerg - Gmm fucto For teger vlues of Why to use Bet dstruto? Bet dstruto fts Beroull smle - cojugte choces 4

5 CS 750 Mche Lerg Bet dstruto Bet osteror dstruto * = Bet Bet Bet Bet

Boml dstruto = * + 3* Emle rolem: co fls where ech co fl c hve two results: hed or tl Outcome: - umer of heds see - umer of tls see trls Model: rolty of hed rolty of tl rolty of outcome: Boml dstruto Boml dstruto: models order deedet sequece of Beroull trls Boml dstruto: Boml dstruto Bm0 0.5 Mtchg ror: Bet dstruto 6

7 Multoml dstruto Emle: multle rolls of de wth 6 results Outcome: couts of occurreces of ossle outcomes of trls: Model rmeters: rolty dstruto: ML estmte: ML s.t. - umer of tmes outcome hs ee see - rolty of outcome Multoml dstruto osteror d MA estmte Choce of the ror: rchlet dstruto.... r r MA.. MA estmte: osteror desty.. r rchlet s the cojugte choce for the multoml smlg

rchlet dstruto: Assume: =3 rchlet dstruto r.. 3 Other dstrutos he sme des c e led to other dstrutos yclly we choose dstrutos tht ehve well so tht comuttos led to ce solutos Eoetl fmly of dstrutos Cojugte choces for some of the dstrutos from the eoetl fmly: Boml Bet Multoml - rchlet Eoetl Gmm osso Iverse Gmm Guss - Guss me d Wshrt covrce 8

Guss orml dstruto Guss: ~ rmeters: - me - stdrd devto esty fucto: e[ Emle: ] 0.4 0.35 0 0.3 0.5 0. 0.5 0. 0.05 0-4 -3 - - 0 3 4 rmeter estmtes Loglelhood l log ML estmtes of the me d vrce: ˆ ˆ ˆ ML vrce estmte s sed E E ˆ Used estmte: ˆ ˆ 9

Multvrte orml dstruto Multvrte orml: ~ rmeters: - me - covrce mtr esty fucto: e d / / Emle: rttoed Guss strutos Multvrte Guss: Emle: recso mtr Wht re the dstrutos for mrgls d codtols? 0

Codtol desty: Codtols d Mrgls Mrgl esty: Codtols d Mrgls

rmeter estmtes Loglelhood ML estmtes of the me d covrces: Covrce estmte s sed Used estmte: ˆ ˆ ˆ ˆ log l ˆ ˆ ˆ E E ˆ ˆ ˆ osteror of the me of multvrte orml Assume ror o the me tht s ormlly dstruted: he the osteror of s ormlly dstruted d / / e e * / / d e / / d

3 CS 750 Mche Lerg osteror of the me of multvrte orml he the osteror of s ormlly dstruted e / / d CS 750 Mche Lerg Other dstrutos Gmm dstruto: Eoetl dstruto: A secl cse of Gmm for = osso dstruto: e e e for ] 0 [ for } 0 {

4 Other dstrutos Gmm dstruto: e for ] 0 [ E vr Sequetl Byes rmeter estmto Sequetl Byes roch Uder the d the estmtes of the osteror c e comuted cremetlly for sequece of dt ots If we use cojugte ror we get c the sme osteror Assume we slt the dt the lst elemet d the rest he: d Θ Θ Θ d Θ Θ Θ Θ A ew ror