1 / 50

Synthetic plots: history and examples David R. Brillinger Statistics Department Universit

Synthetic plots: history and examples David R. Brillinger Statistics Department University of California, Berkeley www.stat.berkeley.edu/~brill brill@stat.berkeley.edu. SECTIONS. Model appraisal methods II. Synthetic plots

munin
Download Presentation

Synthetic plots: history and examples David R. Brillinger Statistics Department Universit

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. Synthetic plots: history and examples David R. Brillinger Statistics Department University of California, Berkeley www.stat.berkeley.edu/~brill brill@stat.berkeley.edu

  2. SECTIONS • Model appraisal methods • II. Synthetic plots • III. Spatial p.p. – galaxies • IV.Timeseries – river flow • V. Spatial-temporal p.p. - wildfires • VI. Trajectories – seals, elk • VII. Summary and discussion • ?. Explanatories

  3. MODEL APPRAISAL • Science needs appraisal methods • Cycle: • Model construction <--> model appraisal • Is model compatable with the data? • Classical chi-squared (df correction) • The method of synthetics • Neyman et al

  4. II. SYNTHETIC PLOTS Simulate realization of fitted model Put real and synthetic side by side Assessment Turing test? Compute same quantity for each?

  5. III. SPATIAL P. P. - galaxies Neyman, Scott and Shane (1953) On the spatial distribution of galaxies …Astr J, 117, 92-133

  6. “… Figure 1 was constructed assuming …, the Poisson law ….”“… it was decided to produce a synthetic plot …”“When the calculated scheme of distribution was compared with the actual distribution of galaxies recorded in Shane’s photographs of the sky [see page 192], it became apparent that the simple mechanism could not produce a distribution resembling the one we see. In the real universe there is a much more pronounced tendency for galaxies to be grouped in clusters.”

  7. Neyman, Scott & Shane (1954) On the index of clumpiness …Astr. J. Suppl. 1, 269-294. “In the third paper …, it was shone that the visual appearance of a ‘synthetic’ photographic plate, obtained by means of a large scale sampling experiment conforming exactly with the assumptions of the theory, is very similar to the actual plate. The only difference noticed between the two is concerned with the small-case clumpiness of images of galaxies.”

  8. In summary: Data {(xj,yj)}, spatial point process 1.Poisson rejected (visually) 2. Clustering (Neyman-Scott process) rejected (visually) 3. More clustering Detail: counting error, variation in limiting magnitude

  9. Results – Scientific American 1956 Turing test?

  10. Comparison – Scott et al (1953)

  11. IV. TIME SERIES – Saugeen River Average monthly flow 1915 – 1976 Walkerton

  12. Data: {y(t)}, time series Hippel and McLeod Periodic autoregression (PAR) Stack years - 62 by 12 matrix (Buys-Ballot) AR(1): Xij = ρjXi,j-1+εij Nonstationary Fit, generate synthetic series

  13. Turing test?

  14. Comparison Spectral ratio

  15. V. SPATIAL-TEMPORAL P. P. - wildfires

  16. Risk analysis Pixel model, {(xj,yj,tj)} logit P{Nxyt = 1} = g1(x,y)+g2(<t>)+h[t] (x,y): location, <t>: day, [t]: year, g1, g2 smooth Sampled 0’s

  17. Original and Bernoulli simulation Turing test?

  18. Comparisons Nearest neighbour distribution

  19. VI. TRAJECTORIES - Hawaiian monk seal, endangered

  20. Foraging, resting, …

  21. DEs. Newtonian motion Scalar potential function, H Planar case, location r = (x,y)’, time t dr(t) = v(t)dt dv(t) = - βv(t)dt – β grad H(r(t),t)dt v: velocity β: coefficient of friction dr = - grad H(r,t)dt = μ(r,t)dt, β >> 0

  22. Examples of H. Point of attraction H(r) = .5*σ2 log r – δr Point of repulsion H(r) = C/r Attraction/repulsion H(r) = α(1/r12 – 1/r6) General parametric H(r) = β10x + β01y + β20x2 + β11xy + β02y2 Nonparametric spline expansion

  23. SDEs. dr(t) = μ(r(t),t)dt + σ(r(t),t)dB(t) μ: drift, –grad H σ: diffusion {B(t)}: bivariate Brownian

  24. Data: {(x(tj),y(tj)),tj)} Solution/approximation. (r(ti+1)-r(ti))/(ti+1-ti) = μ(r(ti),ti) + σ(r(ti),ti) Zi+1/√(ti+1-ti) Euler scheme Approximate likelihood Boundary, startup effects

  25. Fitted potential: general parametric attraction & repulsion

  26. Synthetic Turing test?

  27. Comparison Bagplot cp. Boxplot centre is bivariate median “bag” contains 50% with greatest depth fence – inflate bag by 3 Rousseuw, Rutts, Tukey (1999)

  28. TRAJECTORIES – elk/wapiti Rocky Mountain elk (Cervus elaphus) Banff Starkey Reserve, Oregon Joint usage possible?

  29. Data: {(x(tj),y(tj)), tj)} 8 animals, Δt = 2hr Foraging, resting, hiding, … Model. dr = μ(r)dt + σdB(t) μ smooth – geography velocity field

  30. Boundary (NZ fence) dr= μ(r)dt + σ(r)dB(t) +dA(r) A, support on boundary, keeps particle constrained Synthetic paths. If generated point outside, keep pulling back by half til inside

  31. Turing test?

  32. Comparisons Distribution of distances to “centre”

  33. VII. SUMMARY & DISCUSSION Synthetic plots: method for appraising complex data-based models via Monte Carlo Criteria: EDA, formal Four examples: time series, spatial-temporal p.p., trajectories Found inadequacies in each case

  34. Corrections like Pearson’s for chi-squared Difficulties: land mask irregular times large time differences simulations based on same data

  35. Acknowledgements. Aager, Littman, Preisler, Stewart NSF, FS/USDA

  36. REFERENCES E. Nelson (1967). Dynamical Theories of Brownian Motion. Princeton H. S. Niwa (1996). Newtonian dynamical approach to fish schooling. J. theor. Biol.

  37. Part B. Experiment with explanatory Same 8 animals ATV days, Δt = 5min

  38. Next project

More Related