Smoothing, penalized least squares and splines
|
|
- Laureen Rose
- 6 years ago
- Views:
Transcription
1 Smthing, penalized least squares and splines Duglas Nychka, Lcally weighted averages Penalized least squares smthers Prperties f smthers Splines and Reprducing Kernels The interplatin prf CV and the smthing parameter Supprted by the Natinal Science Fundatin L Ecle D Ete, Ovrnnaz Sep 2005
2 Estimating a curve r surface. The additive statistical mdel: Given n pairs f bservatins (x i, y i ), i = 1,..., n y i = g(x i ) + ɛ i ɛ i s are randm errrs and g is an unknwn smth functin. The gal is t estimate a functin g based n the bservatins
3 A 2-d example Predict surface zne where it is nt mnitred Ambient daily zne in PPB June 16, 1987, US Midwestern Regin
4 Lcal linear surface estimates Use lcal infrmatin t predict unbserved values Use a linear regressin based n clse by bservatins. find ˆβ by least squares. y i = β 1 + ln i β 2 + lat i β 3 + ɛ i The predictin at lcatin ( 88, 41) is just a weighted average f the bservatins. ĝ( 88, 41) = ˆβ ˆβ ˆβ 3 = i=1,n w i y i
5 A kernel estimatr Determine the weights based n the distance t the predictin pints w i (1/h)K((x x i )/h) (and nrmalize s that the weights sum t ne.) Kernel: K K is bump shaped e.g. a nrmal Bandwidth: h h cntrls the spread f K as h gets large the estimate is just the average.
6 Sme kernel estimates fr zne
7 Linear smthers Let ĝ = g(x 1 ),..., g(x n ) be the predictin vectr at the bserved pints. A smther matrix satisfies ĝ = Ay where A is an n n matrix eigenvalues f A are in the range [0,1]. Nte: Ay y Usually values in between the data are filled in by interplating the predictins at the bservatins.
8 Prblems with lcal regressin and kernels Hw large shuld the neighbrhd/bandwidth be? What is the uncertainty f the predictin? Predicting in between bservatins is ad hc and can get weird when the errr is small.
9 Prblems with lcal regressin and kernels Hw large shuld the neighbrhd/bandwidth be? What is the uncertainty f the predictin? Predicting in between bservatins is ad hc and can get weird when the errr is small. But the theretical prperties f kernel estimatrs are well understd... E [g(x) ĝ(x)] 2 = h 4 K 2 /4 + σ2 nh K 0 MSE= Bias 2 + Variance
10 Penalized least squares Ridge regressin Start with yur favrite n basis functins {ψ k } n k=1 estimate has the frm f(x) = n l=1 θ k ψ k (x) where θ = (θ 1,..., θ n ) are the cefficients. The Let W i,k = ψ k (x i ) s f = W θ
11 Estimate the cefficients by a penalized least squares. Sum f squares(θ) + penalty n θ min θ n i=1 (y [W θ] i ) 2 + λθ T Bθ with λ > 0 a hyperparameter and B a nnnegative definite matrix.
12 Estimate the cefficients by a penalized least squares. Sum f squares(θ) + penalty n θ min θ n i=1 (y [W θ] i ) 2 + λθ T Bθ with λ > 0 a hyperparameter and B a nnnegative definite matrix. r in general, - lg likelihd + λ penalty n θ In any case nce we have the parameter estimates these can be used t evaluate ĝ at any pint.
13 The frm f the smther matrix Just calculus... Take derivatives f the penalized likelihd w/r t θ, set equal t zer, slve fr θ
14 The frm f the smther matrix Just calculus... Take derivatives f the penalized likelihd w/r t θ, set equal t zer, slve fr θ The mnster... ˆθ = (W T W + λb) 1 W T y ĝ = W ˆθ = W (W T W + λb) 1 W T y = A(λ)y
15 Why is this a smthing matrix? If W is symmetric A(λ) = W (W T W + λb) 1 W T = (I + λw 1 BW 1 ) 1
16 Effective degrees f freedm in the smther Fr linear regressin trace A(λ) gives us the number f parameters. (Because it is a prjectin matrix) By analgy, tra(λ) is measure f the effective number f degrees f freedm attributed t the smth surface
17 Effective degrees f freedm in the smther Fr linear regressin trace A(λ) gives us the number f parameters. (Because it is a prjectin matrix) By analgy, tra(λ) is measure f the effective number f degrees f freedm attributed t the smth surface A useful decmpsitin Recall: A(λ) = (I + λw 1 BW 1 ) 1 One can always find an rthgnal matrix, U s that U T U = I and (W 1 BW 1 ) = UΓU T where Γ is diagnal. A(λ) = (I + λuγu T ) 1 = U(I + λγ) 1 U T
18 A simple frmula fr the trace S tra(λ) = tru(i + λγ) 1 U T = tr(i + λγ) 1 U T U = tr(i + λγ) 1 = n i= λγ ii The Gamma s are all nnnegative s this must be increasing as λ decreases. The relatinship is ne-t-ne with λ and independent f the data s we can always talk abut λ in terms f the effective degrees f freedm.
19 Splines One btains a spline estimate using a specific basis and a specific penalty matrix. Splines are cnfusing because the basis is a bit mysterius. The classic cubic smthing spline: Fr curve smthing in ne dimensin, min f n i=1 (y i f(x i )) 2 + λ (f (x)) 2 dx The secnd derivative measures the rughness f the fitted curve.
20 Splines One btains a spline estimate using a specific basis and a specific penalty matrix. Splines are cnfusing because the basis is a bit mysterius. The classic cubic smthing spline: Fr curve smthing in ne dimensin, min f n i=1 (y i f(x i )) 2 + λ (f (x)) 2 dx The secnd derivative measures the rughness f the fitted curve. The slutin, is cntinuus up t its secnd derivative and is a piecewise cubic plynmial in between the bservatin pints. Where des this cme frm?
21 Climate fr Clrad CO met statin lcatins elevatin July average min/max temp
22 Cubic splines with different λ s July average min/max temp elevatin
23 Sme abstractin Reprducing kernels Think f these as cvariance functins... k(x, x ) T get basis functins we will hld ne argument fixed at sme value f x and let the ther vary. A Space f functins Let H be the smallest Hilbert space with inner prduct, such that k(x, ) H fr all x and k reprduces itself! k(x, ), k(x, ) = k(x, x )
24 Sme abstractin Reprducing kernels Think f these as cvariance functins... k(x, x ) T get basis functins we will hld ne argument fixed at sme value f x and let the ther vary. A Space f functins Let H be the smallest Hilbert space with inner prduct, such that k(x, ) H fr all x and k reprduces itself! k(x, ), k(x, ) = k(x, x ) Fr f in H, k(x, ), f = f(x)!
25 A variatinal prblem With H and its inner prduct, min f H n i=1 (y i f(x i )) 2 + λ f, f
26 A variatinal prblem With H and its inner prduct, min f H n i=1 (y i f(x i )) 2 + λ f, f The slutin Chse basis functins k(., x i ) fr 1 i n ˆf(x) = n i=1 ˆθ i k(x, x i ) Chse rughness penalty matrix B = W
27 A variatinal prblem With H and its inner prduct, min f H n i=1 (y i f(x i )) 2 + λ f, f The slutin Chse basis functins k(., x i ) fr 1 i n ˆf(x) = n i=1 ˆθ i k(x, x i ) Chse rughness penalty matrix B = W this gives ˆθ = (W W T + λw ) 1 W T y = (W + λi) 1 y ĝ = W ˆθ = W (W + λi) 1 y = (I + λw 1 ) 1 y
28 Cvariances Here W ij = k(x i, x j ) this lks like a cvariance matrix desn t it. The Prf The main part is t use the minimum interplatin prperties f splines. It is easiest t prve this in mre generality fr any reprducing kernel Hilbert space. A simple prf is t guess at the minimizing slutin and use the ptimal interplatin results t shw that there can nt be anther slutin. mre n this later.
29 A 1-d cubic smthing spline The key step is a decmpsitin: f(x) = β 1 + β 2 x + h(x) becmes min β,h min f n i=1 n i=1 (y i f(x i )) 2 + λ (y i β 1 + β 2 x i + h(x i )) 2 + λ H: All functins with 2 derivatives Inner prduct: related t the integral Usual (Wahba style) Reprducing Kernel: k(x, x ) = x x 3 + linear terms (f (x)) 2 dx (h (x)) 2 dx
30 A less scary reprducing kernel fr the same cubic spline prblem On the interval [0, 1] Let k(u, v) = u 2 v/2 u 3 /6 fr u < v and k(u, v) = v 2 u/2 v 3 /6 fr u v k(., v) fr fixed v, is a functin that is a cubic plynmial frm 0 t v and is then a linear functin fr u > v. It is twice differentiable. Let H be a Hilbert space f functins n [0, 1] where f(0) = 0 and f (0) = 0 and with inner prduct f, g = 1 0 f (u)g (u)du k is a reprducing kernel fr this space! Shw this just using integratin by parts.
31 Mre n cubic splines T be explicit the cubic smthing spline slutin has the frm: ˆf(x) = ˆβ 1 + ˆβ 2 x + n i=1 ˆθ i k(x, x i ) Besides verifying that the reprducing n the previus page wrks with the integrated secnd derivative inner prduct ne can als check that it is the cvariance functin fr integrated Brwnian mtin,b(t): with B(0)=0. i.e. X(u) = u 0 B(t)dt, E(X(u)X(v)) = k(u, v) S as we will see later, the Bayesian prir assciated with the cubic smthing spline is an integrated Brwnian mtin.
32 The last details... Given the reprducing kernel, the prblem is min β,θ y Xβ W θ 2 + λθ T W θ X is the regressin matrix with clumns 1 and {x i } First minimize ver θ with β fixed. Plug in slutin and then minimize ver β. Ω = (W + λi) ˆβ = (X T Ω 1 X) 1 X T Ω 1 y ˆθ = (W + λi) 1 (y X ˆβ)
33 A 2-d thin plate smthing spline min f n i=1 (y i f i ) 2 + λ R 2 ( 2 ) 2 ( f 2 f u u v ) 2 + ( 2 ) 2 f dudv 2 v Cllectin f secnd partials is invariant t a rtatin. Again, separate ff the linear part f f. f(x) = β 1 + β 2 x 1 + β 3 x 2 + h(x) Reprducing Kernel: k(x, x ) = x x 2 lg( x x ) + linear terms leading t basis functins that are bumps at the bservatin lcatins.
34 Sme thin plate splines fr the zne data
35 The interplatin prblem: a editrial Splines arse in numerical analysis with the interplatin prblem: Find a curve g that minimizes (g (x)) 2 dx, subject t g(x i ) = y i
36 The interplatin prblem: a editrial Splines arse in numerical analysis with the interplatin prblem: Find a curve g that minimizes (g (x)) 2 dx, subject t g(x i ) = y i We already knw hw t d this by letting λ 0 Statisticians are, f curse, suspicius f fitting data exactly 1-d splines have sme very fast cmputatinal prperties that d nt extend t higher dimensins! Assuming a reprducing kernel is equivalent t a mdel fr the unknwn functin. What is that mdel?
37 The general interplatin prblem Given {x i, y i } and a reprducing kernel Hilbert space. Find a functin g in the space s that minimizes g, g. subject t the cnstraints g(x i ) = y i, 1 i n The slutin g = n i=1 θ i k(., x i ) Where θ is determined by slving a system f n linear equatins t guarentee the interplatin cnstraints.
38 Interplatin Prf Suppse g = n i=1 θ i k(., x i ) interplates {x i, y i } and s des sme ther functin h. We will shw that h, h > g, g. Let δ = h g and s δ, δ > 0 h, h = g + δ, g + δ = g, g + 2 g, δ + δ, δ The cl part, using the reprducing prperty g, δ = n i=1 θ i k(., x i ), δ = n i=1 θ i k(., x i ), δ = because δ(x i ) = h(x i ) g(x i ) = y i y i = 0. S h, h = g, g + δ, δ > g, g Dne! n i=1 θ i δ(x i ) = 0
39 A prf f the smthing spline slutin The prf is by cntradictin and uses the interplatin result. Let ĝ be the smthing spline btained as a linear cmbinatin f the kernel basis functins and pssibly a linear r lw rder plynmial. This is fund as a penalized smther by plugging this frm int the penalized least squares criterin and minimizing by rdinary calculus. Suppse there is an h that has a smaller penalized least squares than ĝ. Cnstruct a g that interplates h at the x i and uses the same basis functins as ĝ. By the interplatin result we knw that h, h > g, g and since bth h and g have the same residuals sums f squares the penalized least squares criterin will be smaller fr g. Thus, h can nt be a minimizer. Because ĝ hwever, has the same frm as the interplant and minimizes the criterin. Thus it must in fact be the minimizer ver all functins in the Hilbert space.
40 Chsing λ by Crss-validatin Sequentially leave each bservatin ut and predict it using the rest f the data. Find the λ that gives the best ut f sample predictins. Refitting the spline when each data pint is mitted, and fr a grid f λ values is cmputatinally demanding. Frtunately there is a shrtcut. The magic frmula residual fr g(x i ) having mitted y i (y i ĝ i ) = (y i ĝ i )/(1 A(λ)) i,i This has a simple frm because adding a data pair (x i, ĝ 1 ) t the data des nt change the estimate.
41 CV and Generalized CV criterin CV (λ) (1/n) n (y i ĝ i ) 2 = (1/n) n i=1 i=1 (y i ĝ i ) 2 (1 A(λ)) i,i ) 2 GCV (λ) ni=1 (y i ĝ i ) 2 (1/n) (1 tra(λ)/n) 2 Minimize CV r GCV ver λ t determine a gd value
42 GCV fr the zne data GCV( eff. degrees f freedm), the estimated surface GCV functin GCV pints, slid GCV mdel, dashed GCV ne Eff. number f parameters
43 GCV fr the climate data GCV( eff. degrees f freedm), the estimated surface GCV July average min/max temp eff.df elevatin
44 Summary We have frmulated the curve/surface fitting prblem as penalized least squares. Splines treat estimating the entire curve but als have a finite basis related t a cvariance functin (reprducing kernel). One can use CV r GCV t find the smthing parameter.
45 Thank yu!
Distributions, spatial statistics and a Bayesian perspective
Distributins, spatial statistics and a Bayesian perspective Dug Nychka Natinal Center fr Atmspheric Research Distributins and densities Cnditinal distributins and Bayes Thm Bivariate nrmal Spatial statistics
More informationStatistics, Numerical Models and Ensembles
Statistics, Numerical Mdels and Ensembles Duglas Nychka, Reinhard Furrer,, Dan Cley Claudia Tebaldi, Linda Mearns, Jerry Meehl and Richard Smith (UNC). Spatial predictin and data assimilatin Precipitatin
More informationIAML: Support Vector Machines
1 / 22 IAML: Supprt Vectr Machines Charles Suttn and Victr Lavrenk Schl f Infrmatics Semester 1 2 / 22 Outline Separating hyperplane with maimum margin Nn-separable training data Epanding the input int
More informationWhat is Statistical Learning?
What is Statistical Learning? Sales 5 10 15 20 25 Sales 5 10 15 20 25 Sales 5 10 15 20 25 0 50 100 200 300 TV 0 10 20 30 40 50 Radi 0 20 40 60 80 100 Newspaper Shwn are Sales vs TV, Radi and Newspaper,
More informationMidwest Big Data Summer School: Machine Learning I: Introduction. Kris De Brabanter
Midwest Big Data Summer Schl: Machine Learning I: Intrductin Kris De Brabanter kbrabant@iastate.edu Iwa State University Department f Statistics Department f Cmputer Science June 24, 2016 1/24 Outline
More information3.4 Shrinkage Methods Prostate Cancer Data Example (Continued) Ridge Regression
3.3.4 Prstate Cancer Data Example (Cntinued) 3.4 Shrinkage Methds 61 Table 3.3 shws the cefficients frm a number f different selectin and shrinkage methds. They are best-subset selectin using an all-subsets
More informationLecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff
Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeff Reading: Chapter 2 STATS 202: Data mining and analysis September 27, 2017 1 / 20 Supervised vs. unsupervised learning In unsupervised
More informationIn SMV I. IAML: Support Vector Machines II. This Time. The SVM optimization problem. We saw:
In SMV I IAML: Supprt Vectr Machines II Nigel Gddard Schl f Infrmatics Semester 1 We sa: Ma margin trick Gemetry f the margin and h t cmpute it Finding the ma margin hyperplane using a cnstrained ptimizatin
More informationSimple Linear Regression (single variable)
Simple Linear Regressin (single variable) Intrductin t Machine Learning Marek Petrik January 31, 2017 Sme f the figures in this presentatin are taken frm An Intrductin t Statistical Learning, with applicatins
More informationLecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff
Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeff Reading: Chapter 2 STATS 202: Data mining and analysis September 27, 2017 1 / 20 Supervised vs. unsupervised learning In unsupervised
More informationPattern Recognition 2014 Support Vector Machines
Pattern Recgnitin 2014 Supprt Vectr Machines Ad Feelders Universiteit Utrecht Ad Feelders ( Universiteit Utrecht ) Pattern Recgnitin 1 / 55 Overview 1 Separable Case 2 Kernel Functins 3 Allwing Errrs (Sft
More informationSupport-Vector Machines
Supprt-Vectr Machines Intrductin Supprt vectr machine is a linear machine with sme very nice prperties. Haykin chapter 6. See Alpaydin chapter 13 fr similar cntent. Nte: Part f this lecture drew material
More informationInternal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.
Sectin 7 Mdel Assessment This sectin is based n Stck and Watsn s Chapter 9. Internal vs. external validity Internal validity refers t whether the analysis is valid fr the ppulatin and sample being studied.
More informationk-nearest Neighbor How to choose k Average of k points more reliable when: Large k: noise in attributes +o o noise in class labels
Mtivating Example Memry-Based Learning Instance-Based Learning K-earest eighbr Inductive Assumptin Similar inputs map t similar utputs If nt true => learning is impssible If true => learning reduces t
More informationResampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017
Resampling Methds Crss-validatin, Btstrapping Marek Petrik 2/21/2017 Sme f the figures in this presentatin are taken frm An Intrductin t Statistical Learning, with applicatins in R (Springer, 2013) with
More informationCOMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d)
COMP 551 Applied Machine Learning Lecture 9: Supprt Vectr Machines (cnt d) Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Class web page: www.cs.mcgill.ca/~hvanh2/cmp551 Unless therwise
More informationModule 3: Gaussian Process Parameter Estimation, Prediction Uncertainty, and Diagnostics
Mdule 3: Gaussian Prcess Parameter Estimatin, Predictin Uncertainty, and Diagnstics Jerme Sacks and William J Welch Natinal Institute f Statistical Sciences and University f British Clumbia Adapted frm
More informationA Matrix Representation of Panel Data
web Extensin 6 Appendix 6.A A Matrix Representatin f Panel Data Panel data mdels cme in tw brad varieties, distinct intercept DGPs and errr cmpnent DGPs. his appendix presents matrix algebra representatins
More informationResampling Methods. Chapter 5. Chapter 5 1 / 52
Resampling Methds Chapter 5 Chapter 5 1 / 52 1 51 Validatin set apprach 2 52 Crss validatin 3 53 Btstrap Chapter 5 2 / 52 Abut Resampling An imprtant statistical tl Pretending the data as ppulatin and
More informationT Algorithmic methods for data mining. Slide set 6: dimensionality reduction
T-61.5060 Algrithmic methds fr data mining Slide set 6: dimensinality reductin reading assignment LRU bk: 11.1 11.3 PCA tutrial in mycurses (ptinal) ptinal: An Elementary Prf f a Therem f Jhnsn and Lindenstrauss,
More informationDifferentiation Applications 1: Related Rates
Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm
More informationSupport Vector Machines and Flexible Discriminants
12 Supprt Vectr Machines and Flexible Discriminants This is page 417 Printer: Opaque this 12.1 Intrductin In this chapter we describe generalizatins f linear decisin bundaries fr classificatin. Optimal
More informationPrincipal Components
Principal Cmpnents Suppse we have N measurements n each f p variables X j, j = 1,..., p. There are several equivalent appraches t principal cmpnents: Given X = (X 1,... X p ), prduce a derived (and small)
More informationFall 2013 Physics 172 Recitation 3 Momentum and Springs
Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.
More informationCHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.
MATH 1342 Ch. 24 April 25 and 27, 2013 Page 1 f 5 CHAPTER 24: INFERENCE IN REGRESSION Chapters 4 and 5: Relatinships between tw quantitative variables. Be able t Make a graph (scatterplt) Summarize the
More informationCOMP 551 Applied Machine Learning Lecture 11: Support Vector Machines
COMP 551 Applied Machine Learning Lecture 11: Supprt Vectr Machines Instructr: (jpineau@cs.mcgill.ca) Class web page: www.cs.mcgill.ca/~jpineau/cmp551 Unless therwise nted, all material psted fr this curse
More information4th Indian Institute of Astrophysics - PennState Astrostatistics School July, 2013 Vainu Bappu Observatory, Kavalur. Correlation and Regression
4th Indian Institute f Astrphysics - PennState Astrstatistics Schl July, 2013 Vainu Bappu Observatry, Kavalur Crrelatin and Regressin Rahul Ry Indian Statistical Institute, Delhi. Crrelatin Cnsider a tw
More informationBuilding to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.
Building t Transfrmatins n Crdinate Axis Grade 5: Gemetry Graph pints n the crdinate plane t slve real-wrld and mathematical prblems. 5.G.1. Use a pair f perpendicular number lines, called axes, t define
More informationMisc. ArcMap Stuff Andrew Phay
Misc. ArcMap Stuff Andrew Phay aphay@whatcmcd.rg Prjectins Used t shw a spherical surface n a flat surface Distrtin Shape Distance True Directin Area Different Classes Thse that minimize distrtin in shape
More informationCN700 Additive Models and Trees Chapter 9: Hastie et al. (2001)
CN700 Additive Mdels and Trees Chapter 9: Hastie et al. (2001) Madhusudana Shashanka Department f Cgnitive and Neural Systems Bstn University CN700 - Additive Mdels and Trees March 02, 2004 p.1/34 Overview
More informationCOMP 551 Applied Machine Learning Lecture 5: Generative models for linear classification
COMP 551 Applied Machine Learning Lecture 5: Generative mdels fr linear classificatin Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Jelle Pineau Class web page: www.cs.mcgill.ca/~hvanh2/cmp551
More informationLinear programming III
Linear prgramming III Review 1/33 What have cvered in previus tw classes LP prblem setup: linear bjective functin, linear cnstraints. exist extreme pint ptimal slutin. Simplex methd: g thrugh extreme pint
More information, which yields. where z1. and z2
The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin
More informationCHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS
CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS 1 Influential bservatins are bservatins whse presence in the data can have a distrting effect n the parameter estimates and pssibly the entire analysis,
More informationAP Physics Kinematic Wrap Up
AP Physics Kinematic Wrap Up S what d yu need t knw abut this mtin in tw-dimensin stuff t get a gd scre n the ld AP Physics Test? First ff, here are the equatins that yu ll have t wrk with: v v at x x
More informationRigid Body Dynamics (continued)
Last time: Rigid dy Dynamics (cntinued) Discussin f pint mass, rigid bdy as useful abstractins f reality Many-particle apprach t rigid bdy mdeling: Newtn s Secnd Law, Euler s Law Cntinuus bdy apprach t
More information[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )
(Abut the final) [COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t m a k e s u r e y u a r e r e a d y ) The department writes the final exam s I dn't really knw what's n it and I can't very well
More informationCOMP 551 Applied Machine Learning Lecture 4: Linear classification
COMP 551 Applied Machine Learning Lecture 4: Linear classificatin Instructr: Jelle Pineau (jpineau@cs.mcgill.ca) Class web page: www.cs.mcgill.ca/~jpineau/cmp551 Unless therwise nted, all material psted
More informationNUMBERS, MATHEMATICS AND EQUATIONS
AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS NUMBERS, MATHEMATICS AND EQUATIONS An integral part t the understanding f ur physical wrld is the use f mathematical mdels which can be used t
More informationYou need to be able to define the following terms and answer basic questions about them:
CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f
More informationHomology groups of disks with holes
Hmlgy grups f disks with hles THEOREM. Let p 1,, p k } be a sequence f distinct pints in the interir unit disk D n where n 2, and suppse that fr all j the sets E j Int D n are clsed, pairwise disjint subdisks.
More informationThe general linear model and Statistical Parametric Mapping I: Introduction to the GLM
The general linear mdel and Statistical Parametric Mapping I: Intrductin t the GLM Alexa Mrcm and Stefan Kiebel, Rik Hensn, Andrew Hlmes & J-B J Pline Overview Intrductin Essential cncepts Mdelling Design
More informationPSU GISPOPSCI June 2011 Ordinary Least Squares & Spatial Linear Regression in GeoDa
There are tw parts t this lab. The first is intended t demnstrate hw t request and interpret the spatial diagnstics f a standard OLS regressin mdel using GeDa. The diagnstics prvide infrmatin abut the
More informationPreparation work for A2 Mathematics [2018]
Preparatin wrk fr A Mathematics [018] The wrk studied in Y1 will frm the fundatins n which will build upn in Year 13. It will nly be reviewed during Year 13, it will nt be retaught. This is t allw time
More informationMath Foundations 20 Work Plan
Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant
More informationBootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >
Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);
More informationSAMPLING DYNAMICAL SYSTEMS
SAMPLING DYNAMICAL SYSTEMS Melvin J. Hinich Applied Research Labratries The University f Texas at Austin Austin, TX 78713-8029, USA (512) 835-3278 (Vice) 835-3259 (Fax) hinich@mail.la.utexas.edu ABSTRACT
More informationThe blessing of dimensionality for kernel methods
fr kernel methds Building classifiers in high dimensinal space Pierre Dupnt Pierre.Dupnt@ucluvain.be Classifiers define decisin surfaces in sme feature space where the data is either initially represented
More informationKinetic Model Completeness
5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins
More informationTree Structured Classifier
Tree Structured Classifier Reference: Classificatin and Regressin Trees by L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stne, Chapman & Hall, 98. A Medical Eample (CART): Predict high risk patients
More informationSPH3U1 Lesson 06 Kinematics
PROJECTILE MOTION LEARNING GOALS Students will: Describe the mtin f an bject thrwn at arbitrary angles thrugh the air. Describe the hrizntal and vertical mtins f a prjectile. Slve prjectile mtin prblems.
More informationMath 302 Learning Objectives
Multivariable Calculus (Part I) 13.1 Vectrs in Three-Dimensinal Space Math 302 Learning Objectives Plt pints in three-dimensinal space. Find the distance between tw pints in three-dimensinal space. Write
More informationMATHEMATICS SYLLABUS SECONDARY 5th YEAR
Eurpean Schls Office f the Secretary-General Pedaggical Develpment Unit Ref. : 011-01-D-8-en- Orig. : EN MATHEMATICS SYLLABUS SECONDARY 5th YEAR 6 perid/week curse APPROVED BY THE JOINT TEACHING COMMITTEE
More informationExercise 3 Identification of parameters of the vibrating system with one degree of freedom
Exercise 3 Identificatin f parameters f the vibrating system with ne degree f freedm Gal T determine the value f the damping cefficient, the stiffness cefficient and the amplitude f the vibratin excitatin
More informationSlide04 (supplemental) Haykin Chapter 4 (both 2nd and 3rd ed): Multi-Layer Perceptrons
Slide04 supplemental) Haykin Chapter 4 bth 2nd and 3rd ed): Multi-Layer Perceptrns CPSC 636-600 Instructr: Ynsuck Che Heuristic fr Making Backprp Perfrm Better 1. Sequential vs. batch update: fr large
More informationLHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers
LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the
More informationAP Statistics Notes Unit Two: The Normal Distributions
AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).
More informationIntroduction to Regression
Intrductin t Regressin Administrivia Hmewrk 6 psted later tnight. Due Friday after Break. 2 Statistical Mdeling Thus far we ve talked abut Descriptive Statistics: This is the way my sample is Inferential
More informationData Analysis, Statistics, Machine Learning
Data Analysis, Statistics, Machine Learning Leland Wilkinsn Adjunct Prfessr UIC Cmputer Science Chief Scien
More informationStats Classification Ji Zhu, Michigan Statistics 1. Classification. Ji Zhu 445C West Hall
Stats 415 - Classificatin Ji Zhu, Michigan Statistics 1 Classificatin Ji Zhu 445C West Hall 734-936-2577 jizhu@umich.edu Stats 415 - Classificatin Ji Zhu, Michigan Statistics 2 Examples f Classificatin
More informationChapter 3: Cluster Analysis
Chapter 3: Cluster Analysis } 3.1 Basic Cncepts f Clustering 3.1.1 Cluster Analysis 3.1. Clustering Categries } 3. Partitining Methds 3..1 The principle 3.. K-Means Methd 3..3 K-Medids Methd 3..4 CLARA
More informationChapter Summary. Mathematical Induction Strong Induction Recursive Definitions Structural Induction Recursive Algorithms
Chapter 5 1 Chapter Summary Mathematical Inductin Strng Inductin Recursive Definitins Structural Inductin Recursive Algrithms Sectin 5.1 3 Sectin Summary Mathematical Inductin Examples f Prf by Mathematical
More informationCHM112 Lab Graphing with Excel Grading Rubric
Name CHM112 Lab Graphing with Excel Grading Rubric Criteria Pints pssible Pints earned Graphs crrectly pltted and adhere t all guidelines (including descriptive title, prperly frmatted axes, trendline
More informationExperiment #3. Graphing with Excel
Experiment #3. Graphing with Excel Study the "Graphing with Excel" instructins that have been prvided. Additinal help with learning t use Excel can be fund n several web sites, including http://www.ncsu.edu/labwrite/res/gt/gt-
More informationMath 105: Review for Exam I - Solutions
1. Let f(x) = 3 + x + 5. Math 105: Review fr Exam I - Slutins (a) What is the natural dmain f f? [ 5, ), which means all reals greater than r equal t 5 (b) What is the range f f? [3, ), which means all
More informationPart a: Writing the nodal equations and solving for v o gives the magnitude and phase response: tan ( 0.25 )
+ - Hmewrk 0 Slutin ) In the circuit belw: a. Find the magnitude and phase respnse. b. What kind f filter is it? c. At what frequency is the respnse 0.707 if the generatr has a ltage f? d. What is the
More informationHow do we solve it and what does the solution look like?
Hw d we slve it and what des the slutin l lie? KF/PFs ffer slutins t dynamical systems, nnlinear in general, using predictin and update as data becmes available. Tracing in time r space ffers an ideal
More informationLyapunov Stability Stability of Equilibrium Points
Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),
More informationPhysics 2010 Motion with Constant Acceleration Experiment 1
. Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin
More informationFive Whys How To Do It Better
Five Whys Definitin. As explained in the previus article, we define rt cause as simply the uncvering f hw the current prblem came int being. Fr a simple causal chain, it is the entire chain. Fr a cmplex
More informationx 1 Outline IAML: Logistic Regression Decision Boundaries Example Data
Outline IAML: Lgistic Regressin Charles Suttn and Victr Lavrenk Schl f Infrmatics Semester Lgistic functin Lgistic regressin Learning lgistic regressin Optimizatin The pwer f nn-linear basis functins Least-squares
More informationTrigonometric Ratios Unit 5 Tentative TEST date
1 U n i t 5 11U Date: Name: Trignmetric Ratis Unit 5 Tentative TEST date Big idea/learning Gals In this unit yu will extend yur knwledge f SOH CAH TOA t wrk with btuse and reflex angles. This extensin
More informationLeast Squares Optimal Filtering with Multirate Observations
Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical
More informationChapter 9 Vector Differential Calculus, Grad, Div, Curl
Chapter 9 Vectr Differential Calculus, Grad, Div, Curl 9.1 Vectrs in 2-Space and 3-Space 9.2 Inner Prduct (Dt Prduct) 9.3 Vectr Prduct (Crss Prduct, Outer Prduct) 9.4 Vectr and Scalar Functins and Fields
More informationMODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b
. REVIEW OF SOME BASIC ALGEBRA MODULE () Slving Equatins Yu shuld be able t slve fr x: a + b = c a d + e x + c and get x = e(ba +) b(c a) d(ba +) c Cmmn mistakes and strategies:. a b + c a b + a c, but
More informationAdmin. MDP Search Trees. Optimal Quantities. Reinforcement Learning
Admin Reinfrcement Learning Cntent adapted frm Berkeley CS188 MDP Search Trees Each MDP state prjects an expectimax-like search tree Optimal Quantities The value (utility) f a state s: V*(s) = expected
More informationLim f (x) e. Find the largest possible domain and its discontinuity points. Why is it discontinuous at those points (if any)?
THESE ARE SAMPLE QUESTIONS FOR EACH OF THE STUDENT LEARNING OUTCOMES (SLO) SET FOR THIS COURSE. SLO 1: Understand and use the cncept f the limit f a functin i. Use prperties f limits and ther techniques,
More information20 Faraday s Law and Maxwell s Extension to Ampere s Law
Chapter 20 Faraday s Law and Maxwell s Extensin t Ampere s Law 20 Faraday s Law and Maxwell s Extensin t Ampere s Law Cnsider the case f a charged particle that is ming in the icinity f a ming bar magnet
More informationMATCHING TECHNIQUES. Technical Track Session VI. Emanuela Galasso. The World Bank
MATCHING TECHNIQUES Technical Track Sessin VI Emanuela Galass The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Emanuela Galass fr the purpse f this wrkshp When can we use
More informationStatistical Learning. 2.1 What Is Statistical Learning?
2 Statistical Learning 2.1 What Is Statistical Learning? In rder t mtivate ur study f statistical learning, we begin with a simple example. Suppse that we are statistical cnsultants hired by a client t
More information1 The limitations of Hartree Fock approximation
Chapter: Pst-Hartree Fck Methds - I The limitatins f Hartree Fck apprximatin The n electrn single determinant Hartree Fck wave functin is the variatinal best amng all pssible n electrn single determinants
More informationL a) Calculate the maximum allowable midspan deflection (w o ) critical under which the beam will slide off its support.
ecture 6 Mderately arge Deflectin Thery f Beams Prblem 6-1: Part A: The department f Highways and Public Wrks f the state f Califrnia is in the prcess f imprving the design f bridge verpasses t meet earthquake
More informationAdmissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs
Admissibility Cnditins and Asympttic Behavir f Strngly Regular Graphs VASCO MOÇO MANO Department f Mathematics University f Prt Oprt PORTUGAL vascmcman@gmailcm LUÍS ANTÓNIO DE ALMEIDA VIEIRA Department
More informationPreparation work for A2 Mathematics [2017]
Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will
More information5 th grade Common Core Standards
5 th grade Cmmn Cre Standards In Grade 5, instructinal time shuld fcus n three critical areas: (1) develping fluency with additin and subtractin f fractins, and develping understanding f the multiplicatin
More informationInference in the Multiple-Regression
Sectin 5 Mdel Inference in the Multiple-Regressin Kinds f hypthesis tests in a multiple regressin There are several distinct kinds f hypthesis tests we can run in a multiple regressin. Suppse that amng
More informationPhysics 212. Lecture 12. Today's Concept: Magnetic Force on moving charges. Physics 212 Lecture 12, Slide 1
Physics 1 Lecture 1 Tday's Cncept: Magnetic Frce n mving charges F qv Physics 1 Lecture 1, Slide 1 Music Wh is the Artist? A) The Meters ) The Neville rthers C) Trmbne Shrty D) Michael Franti E) Radiatrs
More informationElements of Machine Intelligence - I
ECE-175A Elements f Machine Intelligence - I Ken Kreutz-Delgad Nun Vascncels ECE Department, UCSD Winter 2011 The curse The curse will cver basic, but imprtant, aspects f machine learning and pattern recgnitin
More informationExample 1. A robot has a mass of 60 kg. How much does that robot weigh sitting on the earth at sea level? Given: m. Find: Relationships: W
Eample 1 rbt has a mass f 60 kg. Hw much des that rbt weigh sitting n the earth at sea level? Given: m Rbt = 60 kg ind: Rbt Relatinships: Slutin: Rbt =589 N = mg, g = 9.81 m/s Rbt = mrbt g = 60 9. 81 =
More informationWork, Energy, and Power
rk, Energy, and Pwer Physics 1 There are many different TYPES f Energy. Energy is expressed in JOULES (J 419J 4.19 1 calrie Energy can be expressed mre specifically by using the term ORK( rk The Scalar
More informationPHYS 314 HOMEWORK #3
PHYS 34 HOMEWORK #3 Due : 8 Feb. 07. A unifrm chain f mass M, lenth L and density λ (measured in k/m) hans s that its bttm link is just tuchin a scale. The chain is drpped frm rest nt the scale. What des
More informationAP Statistics Practice Test Unit Three Exploring Relationships Between Variables. Name Period Date
AP Statistics Practice Test Unit Three Explring Relatinships Between Variables Name Perid Date True r False: 1. Crrelatin and regressin require explanatry and respnse variables. 1. 2. Every least squares
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationand the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:
Algrithm fr Estimating R and R - (David Sandwell, SIO, August 4, 2006) Azimith cmpressin invlves the alignment f successive eches t be fcused n a pint target Let s be the slw time alng the satellite track
More informationx x
Mdeling the Dynamics f Life: Calculus and Prbability fr Life Scientists Frederick R. Adler cfrederick R. Adler, Department f Mathematics and Department f Bilgy, University f Utah, Salt Lake City, Utah
More informationDataflow Analysis and Abstract Interpretation
Dataflw Analysis and Abstract Interpretatin Cmputer Science and Artificial Intelligence Labratry MIT Nvember 9, 2015 Recap Last time we develped frm first principles an algrithm t derive invariants. Key
More informationProfessional Development. Implementing the NGSS: High School Physics
Prfessinal Develpment Implementing the NGSS: High Schl Physics This is a dem. The 30-min vide webinar is available in the full PD. Get it here. Tday s Learning Objectives NGSS key cncepts why this is different
More informationCS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007
CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is
More informationLead/Lag Compensator Frequency Domain Properties and Design Methods
Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin
More informationContents. This is page i Printer: Opaque this
Cntents This is page i Printer: Opaque this Supprt Vectr Machines and Flexible Discriminants. Intrductin............. The Supprt Vectr Classifier.... Cmputing the Supprt Vectr Classifier........ Mixture
More information