# What Could We Do better? Alternative Statistical Methods Jim Crooks and Xingye Qiao.

## Presentation on theme: "What Could We Do better? Alternative Statistical Methods Jim Crooks and Xingye Qiao."— Presentation transcript:

What Could We Do better? Alternative Statistical Methods Jim Crooks and Xingye Qiao

Review of the Statistical Model We take measurements of the beam displacement, y, at times t 1,…,t n What we actually observe is which is a noisy version of y, –or is the error resulting from imperfect measurement at time t j

A Statistical Model For Displacement Under the spring model it is assumed that the displacement over time is governed by the spring model: So we could write:

A Statistical Model for Displacement Remember the assumptions: –Data at different time points are independent –Residuals are normally distributed –Residuals variance is constant over time With these assumptions the model can be written:

What if we have Replicates? We may have repeated measurements of the same beam Notation: Let t ij be the j th time point; then i indexes the repeats at t j Denote the repeated measurement of the beam displacement at tj by i (t j )= ij. If we believe C and K are the same across replicates then we may write the model as: independent over time and replicate

The likelihood Because of our independence assumption, the likelihood of the model (which we think of as a function of the parameters C and K, not the data) is the product of the individual density functions evaluated at the data : where N(x;, 2 ) denotes the normal density with mean and variance 2 evaluated at x.

Maximum Likelihood Estimates The Maximum Likelihood Estimates (MLEs) for C and K, denoted and, are the values of C and K that maximize the likelihood function Given 2 known, the MLEs are the same as what youd get with a Least-Squares procedure (the former tends to justify the use of the latter)

How Good is the Model? We can asses the goodness of fit by using the spring model to predict the observed measurements The predicted (AKA fitted) values, ŷ(t) are obtained by evaluating the spring model with and : We can compare the fitted values at the observed times ŷ(t ij ) = ŷ ij to the observed values ij. Run the MATLAB file inv_beam.m by typing: > inv_beam

Spring Model

Model Residuals We need to know the difference between beam displacement data and our models predictions for the beam displacement: These are called the model residuals The residuals are our best guess for the values of ij. Hence from our current model we would expect the e ij to look independent and normally distributed with constant variance

Spring Model Residuals Are the residuals normally distributed? Are they independent (i.e., is there correlation in time)? Is their variance constant? Run the MATLAB file plotresidual.m by typing: > plotresidual

Spring Model Residuals

Spring Model QQ-plot

Coefficient of Determination R 2 One criteria to use when judging a model is the fraction of the variability in the data it can explain: SSTot is the total variability in the data SSE is the variability left over after fitting the model (SSE SSTot) So R 2 represents the fraction of variability in the data that is explained by the model

Coefficient of Determination R 2 In the example shown above we can find that: This means the spring model accounts for about 52% of the variability in our displacement measurements Is 52% a lot?

Coefficient of Determination R 2 Brief aside: note that we can also get an estimate of 2 from SSE: where n is the number of data points and df is the number of degrees of freedom (AKA the number of unknown parameters)

How Good is this Model? Is R 2 = 52% any good? It depends. It can be useful (or even necessary) to set up a naïve straw man alternative against which to compare a physical model There are many possible alternatives and choosing between them is subjective To illustrate we will use a smoothing spline alternative

Smoothing Splines A cubic spline is a function that is a piecewise cubic polynomial: –Between each sequential pair of time points the function is a cubic polynomial –At each time point the function is continuous and has continuous first and second derivatives –The time points are called knots

Smoothing Splines A smoothing spline is a type of cubic spline where: –The time points are specifically those at which measurements are made, t j –Given (y j,t j ) for all j, is determined to be that cubic spline that minimizes –Here is called the smoothing parameter What happens to the smoothness as α?

Smoothing Splines The value of parameterizes the relative importance of smoothness to fit –Larger values of result in a bigger penalty for curvature and hence results in a smoother fit that may not fit the data –Smaller values of result in a wigglier spline that more closely follows the data Straight Line Exact Interpolator

Smoothing Splines But how do we choose the value of ? Another choice without an objectively correct answer! One useful answer is the value that minimizes the leave-one-out predictive error: –Fit the spline to all the displacement data except one point –Use the spline to predict the displacement at this time point –Repeat over all displacement points and sum the residual errors This is called leave-one-out cross-validation

Fitted Spline

Fitted Spline Residuals Are the residuals normally distributed? Are they independent (i.e., is there correlation in time)? Is their variance constant? You can make your own cross-validated spline using the MATLAB file splineplot.m Dont do it now!!!

Spline Residuals

Spline QQ-plot

Compare the Spring Model to Splines If you compare residuals, those for the spline are generally smaller (i.e., it fits the data better) Spline Coefficient of variation is Our spring model explains less of the variation than does a naive spline (52% < 88%)

Compare the Spring Model to Splines Is the difference big enough to reject the use of the spring model? Again, this is subjective, but can use statistical tests to answer the question as objectively as possible Such tests are beyond the scope of this workshop, but if you are interested in supercharging your group project using them please ask me

Good Luck!!!

Download ppt "What Could We Do better? Alternative Statistical Methods Jim Crooks and Xingye Qiao."

Similar presentations
{"found_math": false, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0}