Statsmodels: Econometric and Statistical Modeling with Python

Seabold, Skipper; Perktold, Josef

doi:10.25080/majora-92bf1922-011

Cited by 4,376 publications

(2,836 citation statements)

References 7 publications

Supporting

Mentioning

2,823

Contrasting

Unclassified

Order By: Relevance

“…This study was funded by the Research Council of Norway (project SEISMOGLAC, 213359/F20 and CalvingSEIS, 244196/E10; support from NORRUSS program, project 233973/H30). We use ObsPy [Beyreuther et al, 2010] and EP/DP [Fyen, 1989;Schweitzer et al, 2012] software for seismic data analysis, and Python packages "statsmodel" [Seabold and Perktold, 2010] for GLMs and "scikit-learn" [Pedregosa et al, 2011] for Bayes classifiers. Figures are produced using GMT [Wessel and Smith, 1998].…”

Section: Acknowledgmentsmentioning

confidence: 99%

A 15 year record of frontal glacier ablation rates estimated from seismic data

Köhler

Nuth

Kohler

et al. 2016

Geophysical Research Letters

View full text Add to dashboard Cite

We present a unique time series of continuous glacier frontal ablation rates with weekly resolution over 15 years estimated from seismic calving observations at Kronebreen, Svalbard. Using linear statistical models, we calibrate the seismic record with 7 years of satellite‐derived frontal ablation measurements. The two basic input parameters required for our models are the cumulative duration of individual seismic calving events and the incompleteness of the seismic record to correct for the effect of seismic background noise. Frontal ablation follows the seasonal glacier speedup, peaking 1–2 months after the melt season maximum. Short‐lived peaks are associated with melt and rain events. Cumulative frontal ablation of Kronebreen between 2001 and 2015 is about 4.0 km3 (3.7 Gt), with the greatest annual loss (0.45 km3) between 2013 and 2014 at the onset of the recent accelerated retreat of the glacier. Our approach provides a potential method for monitoring tidewater glaciers worldwide that have sufficiently close seismic instrumentation.

show abstract

Section: Acknowledgmentsmentioning

confidence: 99%

A 15 year record of frontal glacier ablation rates estimated from seismic data

Köhler

Nuth

Kohler

et al. 2016

Geophysical Research Letters

View full text Add to dashboard Cite

show abstract

“…13 Statistical testing used the "statsmodels" package. 14 Varian Linacs have been shown to have comparable dosimetric characteristics for machines of the same model and energy; beyond that, many different models have similar dosimetric properties. 1,5,8,14,15 This is understandable because different model names do not necessarily relate to differences in dosimetry-for example, the EX and iX models differ only in the inclusion of an OBI system.…”

Section: B Data Analysismentioning

confidence: 99%

Technical Report: Reference photon dosimetry data for Varian accelerators based on IROC-Houston site visit data

et al. 2016

View full text Add to dashboard Cite

Purpose: Accurate data regarding linear accelerator (Linac) radiation characteristics are important for treatment planning system modeling as well as regular quality assurance of the machine. The Imaging and Radiation Oncology Core-Houston (IROC-H) has measured the dosimetric characteristics of numerous machines through their on-site dosimetry review protocols. Photon data are presented and can be used as a secondary check of acquired values, as a means to verify commissioning a new machine, or in preparation for an IROC-H site visit. Methods: Photon data from IROC-H on-site reviews from 2000 to 2014 were compiled and analyzed. Specifically, data from approximately 500 Varian machines were analyzed. Each dataset consisted of point measurements of several dosimetric parameters at various locations in a water phantom to assess the percentage depth dose, jaw output factors, multileaf collimator small field output factors, off-axis factors, and wedge factors. The data were analyzed by energy and parameter, with similarly performing machine models being assimilated into classes. Common statistical metrics are presented for each machine class. Measurement data were compared against other reference data where applicable. Results: Distributions of the parameter data were shown to be robust and derive from a student's t distribution. Based on statistical and clinical criteria, all machine models were able to be classified into two or three classes for each energy, except for 6 MV for which there were eight classes. Quantitative analysis of the measurements for 6, 10, 15, and 18 MV photon beams is presented for each parameter; supplementary material has also been made available which contains further statistical information.

show abstract

“…Convolutions were computed using a 2-D fast Fourier transform. All calculations were implemented in Python (version 2.7.3) using IPython [Pérez and Granger, 2007], NumPy (version 1.7) [Oliphant, 2007], and Pandas (version 12) [McKinney, 2010]; additional statistical analysis was performed with StatsModels [Seabold and Perktold, 2010]. We created an open-source Python package to calculate topographic stresses in a reasonably automated way, which is available at https://github.com/cossatot/halfspace.…”

Section: Numerical Implementationmentioning

confidence: 99%

The weight of the mountains: Constraints on tectonic stress, friction, and fluid pressure in the 2008 Wenchuan earthquake from estimates of topographic loading

Styron

Hetland

2015

JGR Solid Earth

View full text Add to dashboard Cite

Citation:Styron, R. H., and E. A. Hetland (2015), The weight of the mountains: Constraints on tectonic stress, friction, and fluid pressure in the 2008 Wenchuan earthquake from estimates of topographic loading, J. Geophys. Res. Solid Earth, 120, 2697-2716, doi:10.1002 Abstract Though it is widely recognized that large mountain ranges produce significant stresses in the Earth's crust, these stresses are not commonly quantified. Nonetheless, near large mountains topography may affect fault activity by changing the stress balance on the faults. In this work, we calculate the stress field from topography in the Longmen Shan (Sichuan, China) and resolve those stresses on several models of the faults that ruptured in the 2008 M w 7.9 Wenchuan earthquake. We find that the topography results in shear stresses up to 20 MPa and normal stresses up to 80 MPa on the faults, with significant variability across the faults. Topographic stresses generally load the fault in a normal and left-lateral shear sense, opposite to the inferred coseismic slip sense, and thus inhibit the coseismic slip. We estimate the tectonic stress needed to overcome topographic and lithostatic stresses by assuming that the direction of maximum shear accumulated on the faults is roughly collinear with the inferred coseismic slip. We further estimate the static friction and pore fluid pressure assuming that the fault was, on average, at Mohr-Coulomb failure at the time of the Wenchuan earthquake. We use a Bayesian inversion strategy, yielding posterior probability distributions for the estimated parameters. We find most likely estimates of maximum tectonic compressive stress near 0.6 gz and oriented ∼E-W, and minimum tectonic stress near 0.2 gz. Static friction on the fault is near 0.2, and pore fluid pressure is between 0 and 0.4 of the total pressure.

show abstract

Statsmodels: Econometric and Statistical Modeling with Python

Cited by 4,376 publications

References 7 publications

A 15 year record of frontal glacier ablation rates estimated from seismic data

A 15 year record of frontal glacier ablation rates estimated from seismic data

Technical Report: Reference photon dosimetry data for Varian accelerators based on IROC-Houston site visit data

The weight of the mountains: Constraints on tectonic stress, friction, and fluid pressure in the 2008 Wenchuan earthquake from estimates of topographic loading

Contact Info

Product

Resources

About