Results from the paper

These scripts are located in ./rft1d/examples/paper and reproduce the figures from:

Pataky TC (2015) RFT1D: Smooth One-Dimensional Random Field Upcrossing Probabilities in Python. Journal of Statistical Software (accepted for publication 18 April 2015).

Note

The dataset associated with Figures 13, 15, 16 and 17 is available in rft1d.data.weather

Script list:

  • all_results.py — Single script which replicates all paper results.

  1. fig01_fields.py — 1D Gaussian fields

  2. fig02_EC_HC.py — Description of the Euler and Hadwiger characteristics

  3. fig03_EC.py — Euler characteristic validation

  4. fig04_fields1d_broken.py — pieceswise continuous 1D Gaussian fields

  5. fig05_EC_broken.py — topological validation for broken fields

  6. fig06_sf.py — survival function for 1D Gaussian fields

  7. fig07_bonf_sf.py — RFT vs. Bonferroni survival functions for rough fields

  8. fig08_bonf_isf.py — RFT and Bonferroni inverse survival functions crossing points

  9. fig09_val_maxima.py — validation for 1D Gaussian field maxima

  10. fig10_upcrossing.py — anatomy of field upcrossings

  11. fig11_val_cluster.py — validation for 1D Gaussian cluster-level inference

  12. fig12_val_set.py — validation for 1D Gaussian set-level inference

  13. fig13_weather_0_data.py — dataset overview

  14. fig14_fwhm_estimation.py — validation of field smoothness estimation

  15. fig15_weather_1_rft.py — parametric inference using RFT

  16. fig16_weather_2_wrapped.py — parametric inference, assuming a circular field

  17. fig17_weather_3_nonparam.py — nonparametric inference using an RFT-like procedure