1 / 67

SiZer Analysis for Statistical Smoothing and Inference

SiZer analysis is a powerful tool for statistical smoothing, density estimation, and nonparametric regression. It helps identify significant patterns and modes in data, allowing for better visualization and understanding of trends. This analysis can be applied to various fields such as finance, internet traffic data analysis, and gene expression studies.

dmccombs
Download Presentation

SiZer Analysis for Statistical Smoothing and Inference

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Object Orie’d Data Analysis, Last Time • Statistical Smoothing • Histograms – Density Estimation • Scatterplot Smoothing – Nonpar. Regression • SiZer Analysis • Replaces bandwidth selection • Scale Space • Statistical Inference: Which bumps are “really there”? • Visualization

  2. Kernel Density Estimation Choice of bandwidth (window width)? Very important to performance Fundamental Issue: Which modes are “really there”?

  3. SiZer Background Fun Scale Spaces Views (Incomes Data) Surface View

  4. SiZer Background SiZer analysis of British Incomes data:

  5. SiZer Background Finance "tick data": (time, price) of single stock transactions Idea: "on line" version of SiZer for viewing and understanding trends

  6. SiZer Background Finance "tick data": (time, price) of single stock transactions Idea: "on line" version of SiZer for viewing and understanding trends Notes: trends depend heavily on scale double points and more background color transition (flop over at top)

  7. SiZer Background Internet traffic data analysis: SiZer analysis of time series of packet times at internet hub (UNC) Hannig, Marron, and Riedi (2001)

  8. SiZer Background Internet traffic data analysis: SiZer analysis of time series of packet times at internet hub (UNC) across very wide range of scales needs more pixels than screen allows thus do zooming view (zoom in over time) zoom in to yellow bd’ry in next frame readjust vertical axis

  9. SiZer Background Internet traffic data analysis (cont.) Insights from SiZer analysis: Coarse scales: amazing amount of significant structure Evidence of self-similar fractal type process? Fewer significant features at small scales But they exist, so not Poisson process Poisson approximation OK at small scale??? Smooths (top part) stable at large scales?

  10. DependentSiZer Rondonotti, Marron, and Park (2007) SiZer compares data with white noise Inappropriate in time series Dependent SiZer compares data with an assumed model Visual Goodness of Fit test

  11. Dep’entSiZer:2002 Apr 13 Sat 1 pm – 3 pm

  12. Zoomed view (tored region, i.e. “flat top”)

  13. Further Zoom: finds very periodic behavior!

  14. Possible Physical Explanation IP “Port Scan” Common device of hackers Searching for “break in points” Send query to every possible (within UNC domain): IP address Port Number Replies can indicate system weaknesses Internet Traffic is hard to model

  15. SiZer Overview Would you like to try a SiZer analysis? • Matlab software: http://www.unc.edu/depts/statistics/postscript/papers/marron/Matlab6Software/Smoothing/ • JAVA version (demo, beta): Follow the SiZer link from the Wagner Associates home page: http://www.wagner.com/www.wagner.com/SiZer/ • More details, examples and discussions: http://www.stat.unc.edu/faculty/marron/DataAnalyses/SiZer_Intro.html

  16. PCA to find clusters Return to PCA of Mass Flux Data:

  17. PCA to find clusters SiZer analysis of Mass Flux, PC1

  18. PCA to find clusters SiZer analysis of Mass Flux, PC1 Conclusion: • Found 3 significant clusters! • Correspond to 3 known “cloud types” • Worth deeper investigation

  19. Recall Yeast Cell Cycle Data • “Gene Expression”– Micro-array data • Data (after major preprocessing): Expression “level” of: • thousands of genes (d ~ 1,000s) • but only dozens of “cases” (n ~ 10s) • Interesting statistical issue: High Dimension Low Sample Size data (HDLSS)

  20. Yeast Cell Cycle Data, FDA View Central question: Which genes are “periodic” over 2 cell cycles?

  21. Yeast Cell Cycle Data, FDA View Periodic genes? Naïve approach: Simple PCA

  22. Yeast Cell Cycle Data, FDA View • Central question: which genes are “periodic” over 2 cell cycles? • Naïve approach: Simple PCA • No apparent (2 cycle) periodic structure? • Eigenvalues suggest large amount of “variation” • PCA finds “directions of maximal variation” • Often, but not always, same as “interesting directions” • Here need better approach to study periodicities

  23. Yeast Cell Cycles, Freq. 2 Proj. PCA on Freq. 2 Periodic Component Of Data

  24. Frequency 2 Analysis

  25. Frequency 2 Analysis • Project data onto 2-dim space of sin and cos (freq. 2) • Useful view: scatterplot • Angle (in polar coordinates) shows phase • Colors: Spellman’s cell cycle phase classification • Black was labeled “not periodic” • Within class phases approx’ly same, but notable differences • Now try to improve “phase classification”

  26. Yeast Cell Cycle Revisit “phase classification”, approach: • Use outer 200 genes (other numbers tried, less resolution) • Study distribution of angles • Use SiZer analysis (finds significant bumps, etc., in histogram) • Carefully redrew boundaries • Check by studying k.d.e. angles

  27. SiZer Study of Dist’n of Angles

  28. Reclassification of Major Genes

  29. Compare to Previous Classif’n

  30. New Subpopulation View

  31. OODA in Image Analysis First Generation Problems: • Denoising • Segmentation (find object boundaries) • Registration (align objects) (all about single images)

  32. OODA in Image Analysis Second Generation Problems: • Populations of Images • Understanding Population Variation • Discrimination (a.k.a. Classification) • Complex Data Structures (& Spaces) • HDLSS Statistics

  33. HDLSS Data in Image Analysis Why HDLSS (High Dim, Low Sample Size)? • Complex 3-d Objects Hard to Represent • Often need d = 100’s of parameters • Complex 3-d Objects Costly to Segment • Often have n = 10’s of cases

  34. Image Object Representation Major Approaches for Images: • Landmark Representations • Boundary Representations • Medial Representations

  35. Landmark Representations Main Idea: • On each object find important points • Treat point locations as features • I.e. represent objects by vectors of point locations (in 2-d or 3-d) (Fits in OODA framework)

  36. Landmark Representations Basis of Field of Statistical Shape Analysis: (important precursor of FDA & OODA) Main References: • Kendall (1981, 1984) • Bookstein (1984) • Dryden and Mardia (1998) (most readable and comprehnsive)

  37. Landmark Representations Nice Example: • Fly Wing Data (Drosophila fruit flies) • From George Gilchrist, W. & M. U. http://gwgilc.people.wm.edu/ • Graphic Illustrating Landmarks (next page) • Same veins appear in all flies • And always have same relationship • I.e. all landmarks always identifiable

  38. Landmark Representations Landmarks for fly wing data:

  39. Landmark Representations Important issue for landmark approaches: Location, i. e. Registration Illustration with Fly Wing Data (next slide) Problem: • coordinates are “locations in photo” • & unclear where wing is positioned…

  40. Landmark Representations Illustration of Registration, with Fly Wing Data

  41. Landmark Representations Standard Approach to Registration Problem: Procrustes Analysis Idea: mod out location • Can also mod out rotation • Can also mod out size Recommended reference: Dryden and Mardia (1988)

  42. Landmark Representations Procustes Results for Fly Wing Data

  43. Landmark Representations Effect of Procrustes Analysis: Study Difference Between Continents • Flies from Europe & South America • Look for important differences • Project onto mean difference direction • Visualize with movie • Equal time spacing • Through range of data

  44. Landmark Representations No Procrustes Adjustment: Movies on Difference Between Continents

  45. Landmark Representations Effect of Procrustes Analysis: Movies on Difference Between Continents • Raw Data • Driven by location effects • Strongly feels size • Hard to understand shape

  46. Landmark Representations Location, Rotation, Scale Procrustes: Movies on Difference Between Continents

  47. Landmark Representations Effect of Procrustes Analysis: Movies on Difference Between Continents • Raw Data • Driven by location effects • Strongly feels size • Hard to understand shape • Full Procrustes • Mods out location, size, rotation • Allows clear focus on shape

  48. Landmark Representations Major Drawback of Landmarks: • Need to always find each landmark • Need same relationship • I.e. Landmarks need to correspond • Often fails for medical images • E.g. How many corresponding landmarks on a set of kidneys, livers or brains???

  49. Landmark Representations Landmarks for brains??? (thanks to Liz Bullit) Very hard to identify

  50. Landmark Representations Look across people: Some structure in common But “folds” are different Consistent Landmarks???

More Related