Application software in Python 3 to apply the statistical proposal of Jelínek (1978) for a sample of several second-order tensors in order to obtain the mean tensor of the sample, its principal values with their confidence intervals, and the principal directions with their confidence regions.
This application program is able to plot the summary of the statistical model
described above in a stereographic projection for a better understanding of the
outcomes. Provided that, the next picture represents the aim of jelinekstat.
- Documentation
- PyPI
- GitHub
- Open source and free software: BSD-2-Clause.
The code was written in Python 3. The packages numpy,
scipy, matplotlib
and mplstereonet are
required for using jelinekstat. All of them are
downloadable from the PyPI repository by opening a terminal and typing the
following code lines:
pip install numpy pip install scipy pip install matplotlib pip install mplstereonet
To install jelinekstat open a terminal and type:
pip install jelinekstat
To produce the plot shown above execute the following script
from jelinekstat.jelinekstat import tensorStat
# Input data
sample = [[1.02327, 1.02946, 0.94727, -0.01495, -0.03599, -0.05574],
[1.02315, 1.01803, 0.95882, -0.00924, -0.02058, -0.03151],
[1.02801, 1.03572, 0.93627, -0.03029, -0.03491, -0.06088],
[1.02775, 1.00633, 0.96591, -0.01635, -0.04148, -0.02006],
[1.02143, 1.01775, 0.96082, -0.02798, -0.04727, -0.02384],
[1.01823, 1.01203, 0.96975, -0.01126, -0.02833, -0.03649],
[1.01486, 1.02067, 0.96446, -0.01046, -0.01913, -0.03864],
[1.04596, 1.01133, 0.94271, -0.01660, -0.04711, -0.03636]]
confLevel = 0.95
# Performing the calculation all in one function.
jelinekStatsSummary, stereonetPlot = tensorStat(
sample, confLevel=0.95, want2plot=True, plotName='tes', ext='pdf')
stereonetPlot.show()
Jelínek, V (1978). Statistical processing of anisotropy of magnetic susceptibility measured on group of specimens. Studia Geophysica et Geodaetica, 22 (1), pp. 50-62.