[go: up one dir, main page]

Next Issue
Volume 15, April
Previous Issue
Volume 15, February
 
 
entropy-logo

Journal Browser

Journal Browser

Entropy, Volume 15, Issue 3 (March 2013) – 18 articles , Pages 721-1151

  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
250 KiB  
Article
Holographic Dark Information Energy: Predicted Dark Energy Measurement
by Michael Paul Gough
Entropy 2013, 15(3), 1135-1151; https://doi.org/10.3390/e15031135 - 22 Mar 2013
Cited by 5 | Viewed by 6099
Abstract
Several models have been proposed to explain the dark energy that is causing universe expansion to accelerate. Here the acceleration predicted by the Holographic Dark Information Energy (HDIE) model is compared to the acceleration that would be produced by a cosmological constant. While [...] Read more.
Several models have been proposed to explain the dark energy that is causing universe expansion to accelerate. Here the acceleration predicted by the Holographic Dark Information Energy (HDIE) model is compared to the acceleration that would be produced by a cosmological constant. While identical to a cosmological constant at low redshifts, z < 1, the HDIE model results in smaller Hubble parameter values at higher redshifts, z > 1, reaching a maximum difference of 2.6 ± 0.5% around z ~ 1.7. The next generation of dark energy measurements, both those scheduled to be made in space (ESA’s Euclid and NASA’s WFIRST missions) and those to be made on the ground (BigBOSS, LSST and Dark Energy Survey), should be capable of determining whether such a difference exists or not. In addition a computer simulation thought experiment is used to show that the algorithmic entropy of the universe always increases because the extra states produced by the accelerating expansion compensate for the loss of entropy from star formation. Full article
Show Figures

Figure 1

Figure 1
<p>Plotted against log of universe scale size, <span class="html-italic">a</span>, and redshift, <span class="html-italic">z</span>, are: (<b>a</b>) <span class="html-italic">Upper Panel</span>: Log plot of measured stellar densities (various symbols and blue line: see text for measurement sources) and resultant average baryon temperature, <span class="html-italic">T</span>, and the fraction of all baryons in stars, <span class="html-italic">f</span>. Red lines- best power law fits to data points are <span class="html-italic">a</span><sup>+0.98±0.1</sup> for <span class="html-italic">z</span> &lt; 1, and <span class="html-italic">a</span><sup>+2.8±0.3</sup> for <span class="html-italic">z</span> &gt; 1. (<b>b</b>) <span class="html-italic">Middle Panel</span>: Log plot of energy density contributions: red continuous line, HDIE energy density corresponding to the red line fit in the upper panel; dashed red line, cosmological constant; blue line, mass; solid black line, total for HDIE case; dashed black line, total for the case of a cosmological constant; data symbols, energy densities derived from recent Hubble parameter measurements normalised to Hubble constant measurement at <span class="html-italic">a</span>=1; grey continuous line, the gedanken experiment considered in <a href="#sec3dot2-entropy-15-01135" class="html-sec">Section 3.2</a>. (<b>c</b>) <span class="html-italic">Lower Panel</span>: Linear plot of relative differences in total energy, and in Hubble parameter, between the HDIE model and a cosmological constant. The resolving thresholds of three next generation space and ground-based measurements are shown for comparison in green together with the error bar of the existing measurement at z = 2.3 (see <a href="#sec2dot4-entropy-15-01135" class="html-sec">Section 2.4</a>).</p>
Full article ">
1228 KiB  
Article
Evidence of Large-Scale Quantization in Space Plasmas
by George Livadiotis and David J. McComas
Entropy 2013, 15(3), 1118-1134; https://doi.org/10.3390/e15031118 - 22 Mar 2013
Cited by 47 | Viewed by 7842
Abstract
In plasmas, Debye screening structures the possible correlations between particles. We identify a phase space minimum h* in non-equilibrium space plasmas that connects the energy of particles in a Debye sphere to an equivalent wave frequency. In particular, while there is no [...] Read more.
In plasmas, Debye screening structures the possible correlations between particles. We identify a phase space minimum h* in non-equilibrium space plasmas that connects the energy of particles in a Debye sphere to an equivalent wave frequency. In particular, while there is no a priori reason to expect a single value of h* across plasmas, we find a very similar value of h* ≈ (7.5 ± 2.4)×10−22 J·s using four independent methods: (1) Ulysses solar wind measurements, (2) space plasmas that typically reside in stationary states out of thermal equilibrium and spanning a broad range of physical properties, (3) an entropic limit emerging from statistical mechanics, (4) waiting-time distributions of explosive events in space plasmas. Finding a quasi-constant value for the phase space minimum in a variety of different plasmas, similar to the classical Planck constant but 12 orders of magnitude larger may be revealing a new type of quantization in many plasmas and correlated systems more generally. Full article
Show Figures

Figure 1

Figure 1
<p>Phase space portion <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> calculated for the solar wind ion-electron plasma measurements and using Equation (1). (<b>a</b>) Diagram <math display="inline"> <semantics> <mrow> <mo stretchy="false">(</mo> <msub> <mi>ε</mi> <mi>C</mi> </msub> <mo>;</mo> <msub> <mi>t</mi> <mi>C</mi> </msub> <mo stretchy="false">)</mo> </mrow> </semantics> </math> (on a log-log scale), constructed from Ulysses daily measurements. (<b>b</b>) The product <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo>=</mo> <mn>2</mn> <msub> <mi>ε</mi> <mi>C</mi> </msub> <msub> <mi>t</mi> <mi>C</mi> </msub> </mrow> </semantics> </math> is depicted as a function of heliocentric distance <span class="html-italic">r</span>. (a) and (b) are two-dimensional normalized histograms. (<b>c</b>) Normalized histogram of the values of the values of <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math>. The fitted line in (a) has slope -1 and intercept <math display="inline"> <semantics> <mrow> <mo>≅</mo> <mo>−</mo> <mn>22.22</mn> </mrow> </semantics> </math> (= <math display="inline"> <semantics> <mrow> <mi>log</mi> <mo stretchy="false">(</mo> <mstyle scriptlevel="+1"> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </mstyle> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo stretchy="false">)</mo> </mrow> </semantics> </math>). The weighted mean of <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> values in (b) is found to be <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo stretchy="false">(</mo> <mtext>J</mtext> <mo>⋅</mo> <mtext>s</mtext> <mo stretchy="false">)</mo> <mo>≅</mo> <mo>−</mo> <mn>21.92</mn> <mtext> </mtext> <mo>±</mo> <mn>0.15</mn> </mrow> </semantics> </math>. (For details on the statistical method, see <a href="#app1-entropy-15-01118" class="html-app">Appendix A</a>).</p>
Full article ">Figure 2
<p>Phase space portion <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> in non-equilibrium space plasmas. (<b>a</b>) Types of space plasmas that are typically out of thermal equilibrium, across a broad range of electron density and temperature. (<b>b</b>) These plasmas produce a linear relation in <math display="inline"> <semantics> <mrow> <mo stretchy="false">(</mo> <msub> <mi>ε</mi> <mi>C</mi> </msub> <mo>;</mo> <msub> <mi>t</mi> <mi>C</mi> </msub> <mo stretchy="false">)</mo> </mrow> </semantics> </math> diagram (log-log scale) with slope -1, implying a relation <math display="inline"> <semantics> <mrow> <msub> <mi>ε</mi> <mi>C</mi> </msub> <msub> <mi>t</mi> <mi>C</mi> </msub> <mo>~</mo> </mrow> </semantics> </math> <span class="html-italic">constant</span>. (<b>c</b>) The constancy of <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> is also indicated by the respective log <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> values. (For details on the statistical method, see <a href="#app2-entropy-15-01118" class="html-app">Appendix B</a>. Details of plasma parameters used are provided in the <a href="#app3-entropy-15-01118" class="html-app">Appendix C</a>. Acronyms: CH: Corona Holes; CO: Corona; IH: Inner Heliosheath; IS: Interstellar; MA: Magnetosphere; MS: Magnetosheath; PS: Plasma Sheet; RC: Ring Current; TL: Tail Lobe; WH: solar Wind - Helios; WU: solar Wind-Ulysses).</p>
Full article ">Figure 3
<p>Waiting-time distributions of (<b>a–c</b>) Solar Flares [<a href="#B23-entropy-15-01118" class="html-bibr">23</a>,<a href="#B26-entropy-15-01118" class="html-bibr">26</a>,<a href="#B28-entropy-15-01118" class="html-bibr">28</a>], and (<b>d</b>) CMEs [<a href="#B26-entropy-15-01118" class="html-bibr">26</a>] (red data). The modeled distribution (blue lines) that is derived from (3) using one-particle kappa distribution of energy is well-fitted to the data (over six orders of magnitude), leading to an estimation of <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> consistent with other methods.</p>
Full article ">Figure 4
<p>Histograms of the relative errors of density (<b>a</b>,<b>c</b>) and magnetic field (<b>b</b>,<b>d</b>) for both their modeled (upper panels) and observational (lower panels) values (for the Ulysses data separated in 37 intervals of Δ<span class="html-italic">r</span> = 0.1 AU, where <span class="html-italic">ϑ</span> ≥10°). Note that the range of the top panels represent only ~3% that of the bottom ones.</p>
Full article ">Figure 5
<p>Parker relations modeled (black solid line) and observational values of <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> </mrow> </semantics> </math> (red points) plotted in terms of <span class="html-italic">r</span> (for the Ulysses data separated into 37 intervals). The Parker model was produced by fitting <span class="html-italic">B</span> and <span class="html-italic">n</span> separately to the actual data and using <math display="inline"> <semantics> <mrow> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo>≅</mo> <mn>7.053</mn> <mo>⋅</mo> <msup> <mrow> <mn>10</mn> </mrow> <mrow> <mo>−</mo> <mn>24</mn> </mrow> </msup> <mo>⋅</mo> <msup> <mi>B</mi> <mn>2</mn> </msup> <mo>/</mo> <msup> <mi>n</mi> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mrow> </semantics> </math>(see text).</p>
Full article ">Figure 6
<p>(<b>a</b>) Plot of the 11 values of <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo stretchy="false">(</mo> <mtext>J</mtext> <mo>⋅</mo> <mtext>s</mtext> <mo stretchy="false">)</mo> </mrow> </semantics> </math>, estimated for the 11 types of space plasmas (<a href="#app3-entropy-15-01118" class="html-app">Appendix C</a>), and their weighted mean <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo stretchy="false">(</mo> <mtext>J</mtext> <mo>⋅</mo> <mtext>s</mtext> <mo stretchy="false">)</mo> <mo>≅</mo> <mo>−</mo> <mn>21.87</mn> <mtext> </mtext> <mo>±</mo> <mn>0.18</mn> </mrow> </semantics> </math>. (<b>b</b>) Chi-square minimization (fitting). (<b>c</b>) Chi-Square distribution and p-value.</p>
Full article ">Figure 7
<p>(<b>a</b>) Plot of the four values of <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo stretchy="false">(</mo> <mtext>J</mtext> <mo>⋅</mo> <mtext>s</mtext> <mo stretchy="false">)</mo> </mrow> </semantics> </math>, estimated by the four independent methods that are examined in the main text and shown in <a href="#entropy-15-01118-t003" class="html-table">Table 3</a>; their weighted mean <math display="inline"> <semantics> <mrow> <mi>log</mi> <msub> <mi>ℏ</mi> <mo>*</mo> </msub> <mo stretchy="false">(</mo> <mtext>J</mtext> <mo>⋅</mo> <mtext>s</mtext> <mo stretchy="false">)</mo> <mo>≅</mo> <mo>−</mo> <mn>21.93</mn> <mtext> </mtext> <mo>±</mo> <mn>0.14</mn> </mrow> </semantics> </math> is also shown. (<b>b</b>) Chi-square minimization (fitting). (<b>c</b>) Chi-Square distribution and p-value.</p>
Full article ">
276 KiB  
Article
A Maximum Entropy Approach to Loss Distribution Analysis
by Marco Bee
Entropy 2013, 15(3), 1100-1117; https://doi.org/10.3390/e15031100 - 22 Mar 2013
Cited by 6 | Viewed by 5731
Abstract
In this paper we propose an approach to the estimation and simulation of loss distributions based on Maximum Entropy (ME), a non-parametric technique that maximizes the Shannon entropy of the data under moment constraints. Special cases of the ME density correspond to standard [...] Read more.
In this paper we propose an approach to the estimation and simulation of loss distributions based on Maximum Entropy (ME), a non-parametric technique that maximizes the Shannon entropy of the data under moment constraints. Special cases of the ME density correspond to standard distributions; therefore, this methodology is very general as it nests most classical parametric approaches. Sampling the ME distribution is essential in many contexts, such as loss models constructed via compound distributions. Given the difficulties in carrying out exact simulation,we propose an innovative algorithm, obtained by means of an extension of Adaptive Importance Sampling (AIS), for the approximate simulation of the ME distribution. Several numerical experiments confirm that the AIS-based simulation technique works well, and an application to insurance data gives further insights in the usefulness of the method for modelling, estimating and simulating loss distributions. Full article
Show Figures

Figure 1

Figure 1
<p>The simulated Lognormal–Pareto observations, the true Lognormal–Pareto mixture density and the fitted ME(7) density.</p>
Full article ">Figure 2
<p>The histogram of the General Liability Claims data with the estimated ME(3), ME(4) and ME(5) densities.</p>
Full article ">Figure 3
<p>The estimated ME(4) density superimposed on the normalized histogram of the simulated observations.</p>
Full article ">Figure 4
<p>The estimated ME(4) and lognormal densities superimposed on the histogram of the General Liability Claims observations.</p>
Full article ">
2378 KiB  
Article
Substrate Effect on Catalytic Loop and Global Dynamics of Triosephosphate Isomerase
by Zeynep Kurkcuoglu and Pemra Doruker
Entropy 2013, 15(3), 1085-1099; https://doi.org/10.3390/e15031085 - 18 Mar 2013
Cited by 4 | Viewed by 6007
Abstract
The opening/closure of the catalytic loop 6 over the active site in apo triosephosphate isomerase (TIM) has been previously shown to be driven by the global motions of the enzyme, specifically the counter-clockwise rotation of the subunits. In this work, the effect of [...] Read more.
The opening/closure of the catalytic loop 6 over the active site in apo triosephosphate isomerase (TIM) has been previously shown to be driven by the global motions of the enzyme, specifically the counter-clockwise rotation of the subunits. In this work, the effect of the substrate dihydroxyacetone phosphate (DHAP) on TIM dynamics is assessed using two apo and two DHAP-bound molecular dynamics (MD) trajectories (each 60 ns long). Multiple events of catalytic loop opening/closure take place during 60 ns runs for both apo TIM and its DHAP-complex. However, counter-clockwise rotation observed in apo TIM is suppressed and bending-type motions are linked to loop dynamics in the presence of DHAP. Bound DHAP molecules also reduce the overall mobility of the enzyme and change the pattern of orientational cross-correlations, mostly those within each subunit. The fluctuations of pseudodihedral angles of the loop 6 residues are enhanced towards the C-terminus, when DHAP is bound at the active site. Full article
(This article belongs to the Special Issue Loop Entropy)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Residue MSF for apo and complex simulations. Inset figure shows the complex residues with reduced (blue) and increased (red) MSF by 20% compared to apo. Green regions have similar MSF in both apo and complex simulations. DHAP is shown with stick representation.</p>
Full article ">Figure 2
<p>The distance between C<sup>α</sup> atoms of loop 6 tip residue (I170) and a reference residue (Y208) as a function of time. Left and right panels summarize the profiles for apo and complex simulations. Solid and dashed lines indicate the respective closed and open states of the loop based on crystal structures. Loop 6 samples open conformations even in the presence of DHAP. It opens/closes multiple times during apo and complex simulations.</p>
Full article ">Figure 3
<p>Distance histograms for I170-Y208, G171-Y208 and T172-Y208 for (<b>a</b>, <b>c</b>, <b>e</b>) apo and (<b>b</b>, <b>d</b>, <b>f</b>) complex simulations, respectively. Distance distribution is obtained by dividing the number of occurrence of the distance to the total number of snapshots for 10–60 ns range.</p>
Full article ">Figure 4
<p>Most dominant dynamics from PC1 of (<b>a</b>) Apo1, (<b>b</b>) Apo2, (<b>c</b>) Complex1 and (<b>d</b>) Complex2. Subunit A (orange), subunit B (blue) with loop 6 (red) are shown in ribbon and DHAP in magenta-orange stick representation. In (<b>a</b>) and (<b>b</b>), PC1 corresponds to counter-clockwise rotation (CCR) of two subunits, coupled with loop 6 opening/closing. In (<b>c</b>), PC1 of Complex1 simulation is a bending type motion rather than counter-clockwise rotation but loop 6 opening/closing is observed in this mode. In (<b>d</b>), subunit A (without DHAP) rotates like in the Apo1 PC1 coupled with loop opening/closing, whereas a mixed-type of motion is observed in subunit B (with DHAP), where the helix 5 undergoes bending without loop 6 opening/closing.</p>
Full article ">Figure 5
<p>Other PCs for Apo1, Apo2, Complex1 and Complex2 simulations.</p>
Full article ">Figure 6
<p>Normalized orientational cross-correlation maps for apo and complex1. The results for apo are averaged over all the subunits in Apo1 and Apo2 simulations and for complex, over Complex1 A and B subunits. The correlated regions with loop 6 (L6, represented by the average of three tip residues I170, G171, T172) are given in cartoon representations between intra and inter-subunit maps and colored using blue (negative correlation)-white (no correlation)-red (positive correlation) coloring scheme. In the complex intra-subunit map, two blocks- Regions 1 and 2- are negatively correlated with each other. Region 1 is also positively correlated with the same region of the other subunit (named as Region 3) in the inter-subunit map.</p>
Full article ">Figure 7
<p>Apo1 and Complex1 most dominant PCs from loop 6 A PCA. PC1 corresponds to loop opening/closing and PC2 lateral movement of the loop for both simulations.</p>
Full article ">Figure 8
<p>Average pseudodihedral RMS fluctuations for loop 6 residues in apo and complex simulations.</p>
Full article ">
579 KiB  
Article
Time Series Analysis Using Composite Multiscale Entropy
by Shuen-De Wu, Chiu-Wen Wu, Shiou-Gwo Lin, Chun-Chieh Wang and Kung-Yen Lee
Entropy 2013, 15(3), 1069-1084; https://doi.org/10.3390/e15031069 - 18 Mar 2013
Cited by 278 | Viewed by 16336
Abstract
Multiscale entropy (MSE) was recently developed to evaluate the complexity of time series over different time scales. Although the MSE algorithm has been successfully applied in a number of different fields, it encounters a problem in that the statistical reliability of the sample [...] Read more.
Multiscale entropy (MSE) was recently developed to evaluate the complexity of time series over different time scales. Although the MSE algorithm has been successfully applied in a number of different fields, it encounters a problem in that the statistical reliability of the sample entropy (SampEn) of a coarse-grained series is reduced as a time scale factor is increased. Therefore, in this paper, the concept of a composite multiscale entropy (CMSE) is introduced to overcome this difficulty. Simulation results on both white noise and 1/f noise show that the CMSE provides higher entropy reliablity than the MSE approach for large time scale factors. On real data analysis, both the MSE and CMSE are applied to extract features from fault bearing vibration signals. Experimental results demonstrate that the proposed CMSE-based feature extractor provides higher separability than the MSE-based feature extractor. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematic illustration of the coarse-grained procedure. Modified from reference [<a href="#B3-entropy-15-01069" class="html-bibr">3</a>].</p>
Full article ">Figure 2
<p>Brute force method. Modified from reference [<a href="#B18-entropy-15-01069" class="html-bibr">18</a>].</p>
Full article ">Figure 3
<p>Schematic illustration of the CMSE procedure.</p>
Full article ">Figure 4
<p>Flow charts of MSE and CMSE algorithms.</p>
Full article ">Figure 5
<p>MSE results of (<b>a</b>) white noise and (<b>b</b>) 1/<span class="html-italic">f</span> noise with different data lengths.</p>
Full article ">Figure 6
<p>CMSE results of (<b>a</b>) white noise and (<b>b</b>) 1/<span class="html-italic">f</span> noise with different data lengths.</p>
Full article ">Figure 7
<p>MSE and CMSE results of white noise with data lengths (<b>a</b>) <span class="html-italic">N</span> = 2,000 and (<b>b</b>) <span class="html-italic">N</span> = 10,000.</p>
Full article ">Figure 8
<p>MSE and CMSE results of 1/<span class="html-italic">f</span> noise with data lengths (<b>a</b>) <span class="html-italic">N</span> = 2,000 and (<b>b</b>) <span class="html-italic">N</span> = 10,000</p>
Full article ">Figure 9
<p>Measured acceleration signals of vibrations in the time domain of six different bearing conditions (<b>a</b>) normal state, ball fault and inner race fault; (<b>b</b>) outer race faults at 3, 6, and 12 o’clock positions.</p>
Full article ">Figure 10
<p>MSE and CMSE results on bearing vibration data (1,730 rpm, 7 mils). (<b>a</b>) Normal state. (<b>b</b>) Ball fault. (<b>c</b>) Inner race fault. (<b>d</b>) Outer race fault (3 o’clock position). (<b>e</b>) Outer race fault (6 o’clock position). (<b>f</b>) Outer race fault (12 o’clock position)</p>
Full article ">
229 KiB  
Article
Hawking and Unruh Effects of a 5-Dimensional Minimal Gauged Supergravity Black Hole by a Global Embedding Approach
by Hui-Hua Zhao, Li-Chun Zhang and Guang-Liang Li
Entropy 2013, 15(3), 1057-1068; https://doi.org/10.3390/e15031057 - 18 Mar 2013
Viewed by 5027
Abstract
Using the new global embedding approach we investigate Unruh/Hawking temperature of the 5-dimensional minimal gauged supergravity black hole with double rotating parameters in a general (1 + 1) space-time. Our results verify that views of Banerjee and Majhi, and extend this approach to [...] Read more.
Using the new global embedding approach we investigate Unruh/Hawking temperature of the 5-dimensional minimal gauged supergravity black hole with double rotating parameters in a general (1 + 1) space-time. Our results verify that views of Banerjee and Majhi, and extend this approach to a higher dimension situation. Full article
152 KiB  
Article
Entropy Principle and Galilean Relativity for Dense Gases, the General Solution without Approximations
by Maria Cristina Carrisi, Stefania Montisci and Sebastiano Pennisi
Entropy 2013, 15(3), 1035-1056; https://doi.org/10.3390/e15031035 - 11 Mar 2013
Cited by 10 | Viewed by 4252
Abstract
The many moments model for dense gases and macromolecular fluids is considered here, where the upper order moment is chosen in accordance to the suggestions of the non-relativistic limit of the corresponding relativistic model. The solutions of the restrictions imposed by the entropy [...] Read more.
The many moments model for dense gases and macromolecular fluids is considered here, where the upper order moment is chosen in accordance to the suggestions of the non-relativistic limit of the corresponding relativistic model. The solutions of the restrictions imposed by the entropy principle and that of Galilean relativity were, until now, obtained in the literature by using Taylor expansions around equilibrium and without proving convergence. Here, an exact solution without using expansions is found. The particular case with only 14 moments has already been treated in the literature in a completely different way. Here, it is proven that this particular closure is included in the presently more general one. Full article
731 KiB  
Article
Exergetic and Parametric Study of a Solar Aided Coal-Fired Power Plant
by Rongrong Zhai, Yong Zhu, Yongping Yang, Kaiyu Tan and Eric Hu
Entropy 2013, 15(3), 1014-1034; https://doi.org/10.3390/e15031014 - 11 Mar 2013
Cited by 33 | Viewed by 7883
Abstract
A solar-aided coal-fired power plant realizes the integration of a fossil fuel (coal or gas) and clean energy (solar). In this paper, a conventional 600 MW coal-fired power plant and a 600 MW solar-aided coal-fired power plant have been taken as the study [...] Read more.
A solar-aided coal-fired power plant realizes the integration of a fossil fuel (coal or gas) and clean energy (solar). In this paper, a conventional 600 MW coal-fired power plant and a 600 MW solar-aided coal-fired power plant have been taken as the study case to understand the merits of solar-aided power generation (SAPG) technology. The plants in the case study have been analyzed by using the First and Second Laws of Thermodynamics principles. The solar irradiation and load ratio have been considered in the analysis. We conclude that if the solar irradiation was 925 W/m2 and load ratio of the SAPG plant was 100%, the exergy efficiency would be 44.54% and the energy efficiency of the plant (46.35%). It was found that in the SAPG plant the largest exergy loss was from the boiler, which accounted for about 76.74% of the total loss. When the load ratio of the unit remains at 100%, and the solar irradiation varies from 500 W/m2 to 1,100 W/m2, the coal savings would be in the range of 8.6 g/kWh to 15.8 g/kWh. If the solar irradiation were kept at 925 W/m2 while the load ratio of the plant changed from 30% to 100%, the coal savings could be in the range of 11.99 g/kWh to 13.75 g/kWh. Full article
(This article belongs to the Special Issue Exergy: Analysis and Applications)
Show Figures

Figure 1

Figure 1
<p>The layout of the coal-fired power plant.</p>
Full article ">Figure 2
<p>The layout of the solar aided coal-fired power plant.</p>
Full article ">Figure 3
<p>The diagram of oil-water heat exchanger.</p>
Full article ">Figure 4
<p>The distribution of exergy destructions in the base plant.</p>
Full article ">Figure 5
<p>The distribution of energy destructions in the base plant.</p>
Full article ">Figure 6
<p>The distribution of exergy destructions in the new plant (SAPG).</p>
Full article ">Figure 7
<p>The distribution of energy destructions in the new plant (SAPG).</p>
Full article ">Figure 8
<p>(<b>a</b>) The energy and exergy efficiency of the main facilities in the base plant; (<b>b</b>) The energy and exergy efficiency of the main facilities in the new plant.</p>
Full article ">Figure 9
<p>(<b>a</b>) The exergy efficiency of the solar aided oil-water heat exchanger; (<b>b</b>) The energy efficiency of the solar aided oil-water heat exchanger.</p>
Full article ">Figure 10
<p>The energy efficiency and exergy efficiency with the change of solar irradiation (500 W/m<sup>2</sup>–1,100 W/m<sup>2</sup>) in SAPG.</p>
Full article ">Figure 11
<p>The coal consumption rate with the change of solar irradiation in SAPG.</p>
Full article ">Figure 12
<p>(<b>a</b>) The energy and exergy efficiency with the change of load ratio in the base plant (30%–100%); (<b>b</b>) The energy and exergy efficiency with the change of load ratio in the new plant (30%–100%); (<b>c</b>) The coal consumption rate with the change of load ratio.</p>
Full article ">Figure 12 Cont.
<p>(<b>a</b>) The energy and exergy efficiency with the change of load ratio in the base plant (30%–100%); (<b>b</b>) The energy and exergy efficiency with the change of load ratio in the new plant (30%–100%); (<b>c</b>) The coal consumption rate with the change of load ratio.</p>
Full article ">
5973 KiB  
Article
New Climatic Indicators for Improving Urban Sprawl: A Case Study of Tehran City
by Abdolazim Ghanghermeh, Gholamreza Roshan, José A. Orosa, José L. Calvo-Rolle and Ángel M. Costa
Entropy 2013, 15(3), 999-1013; https://doi.org/10.3390/e15030999 - 7 Mar 2013
Cited by 37 | Viewed by 7045
Abstract
In the modern world, the fine balance and delicate relationship between human society and the environment in which we exist has been affected by the phenomena of urbanisation and urban development. Today, various environmental factors give rise to horizontal dispersion, spread and growth [...] Read more.
In the modern world, the fine balance and delicate relationship between human society and the environment in which we exist has been affected by the phenomena of urbanisation and urban development. Today, various environmental factors give rise to horizontal dispersion, spread and growth of cities. One of the most important results of this is climatic change which is directly affected by the urban sprawl of every metropolis. The aim of this study is to identify the relationship between the various horizontally distributed components of Tehran city and changes in essential microclimate clusters, by means of the humidex index. Results showed that, when the humidex was calculated for each of the obtained clusters, it was evident that it had increased with time, in parallel with Shannon’s entropy, as a consequence of the average temperature and relative humidity of each cluster. At the same time, results have shown that both temperature and relative humidity of the study area are related with urban sprawl, urbanisation and development, as defined by Shannon’s entropy and, in consequence, with humidex. In consequence, this new concept must be considered in future research works to predict and control urban sprawl and microclimate conditions in cities. Full article
(This article belongs to the Special Issue Entropy and Urban Sprawl)
Show Figures

Figure 1

Figure 1
<p>Beyond the scope of Tehran, to reconstruct climate data.</p>
Full article ">Figure 2
<p>The physical development of Tehran [<a href="#B18-entropy-15-00999" class="html-bibr">18</a>,<a href="#B26-entropy-15-00999" class="html-bibr">26</a>,<a href="#B27-entropy-15-00999" class="html-bibr">27</a>].</p>
Full article ">Figure 3
<p>Comparison of daily minimum temperature reconstructed data using Kriging method, with experimental data at Mehrabad station.</p>
Full article ">Figure 4
<p>Changes of essential microclimate clusters in the study area for different decades. (a) 1966-1975. (b) 1976-1985 (c) 1986-1995 (d) 1996-2005.</p>
Full article ">Figure 5
<p>Shannon entropy coefficient variation process for studied decades.</p>
Full article ">Figure 6
<p>Humidex variation for different study periods.</p>
Full article ">
540 KiB  
Article
Experimental Assessment of a 2-D Entropy-Based Model for Velocity Distribution in Open Channel Flow
by Nicola Fontana, Gustavo Marini and Francesco De Paola
Entropy 2013, 15(3), 988-998; https://doi.org/10.3390/e15030988 - 6 Mar 2013
Cited by 22 | Viewed by 6536
Abstract
Velocity distribution in an open channel flow can be very useful to model many hydraulic phenomena. Among the others, several 1D models based on the concept of entropy are available in the literature, which allow estimating the velocity distribution by measuring velocities only [...] Read more.
Velocity distribution in an open channel flow can be very useful to model many hydraulic phenomena. Among the others, several 1D models based on the concept of entropy are available in the literature, which allow estimating the velocity distribution by measuring velocities only in a few points. Nevertheless, since 1D models have often a limited practical use, a 2D entropy based model was recently developed. The model provides a reliable estimation of the velocity distribution for open channel flow with a rectangular cross section, if the maximum velocity and the average velocity are known. In this paper results from the proposed model were compared with measured velocities carried out from laboratory experiments. Calculated values were also compared with results inferred from a 2D model available in the literature, resulting in a greater ease of use and a more reliable estimate of the velocity profile. Full article
(This article belongs to the Special Issue Maximum Entropy and Bayes Theorem)
Show Figures

Figure 1

Figure 1
<p>Rectangular symmetrical domain in (<b>a</b>) dimensional and (<b>b</b>) non-dimensional coordinates.</p>
Full article ">Figure 2
<p>Contour sketch of CDF <span class="html-italic">F</span>(<span class="html-italic">u</span>) (<span class="html-italic">H</span>/<span class="html-italic">B</span> = 0.5 and <span class="html-italic">ψ<sub>0</sub></span> = 0.8).</p>
Full article ">Figure 3
<p>Flume used for experimental measurements.</p>
Full article ">Figure 4
<p>Comparison between experimental and theoretical data calculated by 2D proposed model and Chiu’s model (experiment n. 3).</p>
Full article ">Figure 5
<p>Comparison between experimental and theoretical data calculated by 2D proposed model and Chiu’s model (experiment n. 17).</p>
Full article ">Figure 6
<p>Comparison between experimental and theoretical data; triangle corresponds to <span class="html-italic">c</span> = 0, circle to <span class="html-italic">c</span> = 1.</p>
Full article ">
245 KiB  
Article
Information Properties of Boundary Line Models for N2O Emissions from Agricultural Soils
by Cairistiona F.E. Topp, Weijin Wang, Joanna M. Cloy, Robert M. Rees and Gareth Hughes
Entropy 2013, 15(3), 972-987; https://doi.org/10.3390/e15030972 - 5 Mar 2013
Cited by 9 | Viewed by 6045
Abstract
Boundary line models for N2O emissions from agricultural soils provide a means of estimating emissions within defined ranges. Boundary line models partition a two-dimensional region of parameter space into sub-regions by means of thresholds based on relationships between N2O [...] Read more.
Boundary line models for N2O emissions from agricultural soils provide a means of estimating emissions within defined ranges. Boundary line models partition a two-dimensional region of parameter space into sub-regions by means of thresholds based on relationships between N2O emissions and explanatory variables, typically using soil data available from laboratory or field studies. Such models are intermediate in complexity between the use of IPCC emission factors and complex process-based models. Model calibration involves characterizing the extent to which observed data are correctly forecast. Writing the numerical results from graphical two-threshold boundary line models as 3×3 prediction-realization tables facilitates calculation of expected mutual information, a measure of the amount of information about the observations contained in the forecasts. Whereas mutual information characterizes the performance of a forecaster averaged over all forecast categories, specific information and relative entropy both characterize aspects of the amount of information contained in particular forecasts. We calculate and interpret these information quantities for experimental N2O emissions data. Full article
(This article belongs to the Special Issue Applications of Information Theory in the Geosciences)
Show Figures

Figure 1

Figure 1
<p>The parameter space delimited by observed ranges of water filled pore space (WFPS, %) and soil temperature (T, °C) in 2011-2012 at a grassland site in Dumfries, SW Scotland, receiving inorganic fertilizer, is the basis for a boundary line model. Observed N<sub>2</sub>O emissions were categorized as ‘low’ (&lt;10 g N<sub>2</sub>O-N ha<sup>−1</sup> day<sup>−1</sup>), ‘medium’ (10-100 g N<sub>2</sub>O-N ha<sup>−1</sup> day<sup>−1</sup>) or ‘high’ (&gt;100 g N<sub>2</sub>O-N ha<sup>−1</sup> day<sup>−1</sup>), as in Conen <span class="html-italic">et al</span>. [<a href="#B15-entropy-15-00972" class="html-bibr">15</a>]. There were 715 ‘low’ observations, 322 ‘medium’ observations and 19 ‘high’ observations (<span class="html-italic">N</span> = 1056), resulting in many overlapping data points on the graph. The boundary lines between forecast emission categories are WFPS(%) + 2∙T(°C) = 90 (low-medium) and WFPS(%) + 2∙T(°C) = 105 (medium-high), as described in Conen <span class="html-italic">et al</span>. [<a href="#B15-entropy-15-00972" class="html-bibr">15</a>].</p>
Full article ">Figure 2
<p>For each forecast category <span class="html-italic">i</span>, the bar comprises a red component <span class="html-italic">I<sub>S</sub></span>(<span class="html-italic">f<sub>i</sub></span>), and a blue component <span class="html-italic">H</span>(<span class="html-italic">O</span>|<span class="html-italic">f<sub>i</sub></span>) which together sum to <span class="html-italic">H</span>(<span class="html-italic">O</span>) in each case. The weighted average of red components is equal to <span class="html-italic">I<sub>M</sub></span>(<span class="html-italic">O</span>,<span class="html-italic">F</span>) (the Pr(<span class="html-italic">f<sub>i</sub></span>) provide the appropriate weights).</p>
Full article ">Figure 3
<p>For each forecast category <span class="html-italic">i</span>, the bar comprises a red component <span class="html-italic">I</span>(<span class="html-italic">f<sub>i</sub></span>), and a blue component <span class="html-italic">H</span>(<span class="html-italic">O</span>|<span class="html-italic">f<sub>i</sub></span>). The weighted average of the sums of the two components is equal to <span class="html-italic">H</span>(<span class="html-italic">O</span>). The weighted average of red components is equal to <span class="html-italic">I<sub>M</sub></span>(<span class="html-italic">O</span>,<span class="html-italic">F</span>). In each case the Pr(<span class="html-italic">f<sub>i</sub></span>) provide the appropriate weights.</p>
Full article ">
243 KiB  
Article
On Classical Ideal Gases
by Jacques Arnaud, Laurent Chusseau and Fabrice Philippe
Entropy 2013, 15(3), 960-971; https://doi.org/10.3390/e15030960 - 27 Feb 2013
Cited by 3 | Viewed by 6970
Abstract
We show that the thermodynamics of ideal gases may be derived solely from the Democritean concept of corpuscles moving in vacuum plus a principle of simplicity, namely that these laws are independent of the laws of motion, aside from the law of energy [...] Read more.
We show that the thermodynamics of ideal gases may be derived solely from the Democritean concept of corpuscles moving in vacuum plus a principle of simplicity, namely that these laws are independent of the laws of motion, aside from the law of energy conservation. Only a single corpuscle in contact with a heat bath submitted to a z and t-invariant force is considered. Most of the end results are known but the method appears to be novel. The mathematics being elementary, the present paper should facilitate the understanding of the ideal gas law and of classical thermodynamics even though not-usually-taught concepts are being introduced. Full article
Show Figures

Figure 1

Figure 1
<p>Space-time (<math display="inline"> <mrow> <mi>z</mi> <mo>,</mo> <mi>t</mi> </mrow> </math>) trajectory for a corpuscle of weight <span class="html-italic">w</span> bouncing off the ground (<span class="html-italic">z</span> = 0). The maximum altitude reached by the corpuscle is <math display="inline"> <mrow> <msub> <mi>z</mi> <mi>m</mi> </msub> <mo>=</mo> <mi>E</mi> <mo>/</mo> <mi>w</mi> </mrow> </math>, where <span class="html-italic">E</span> denotes the energy. The motion is periodic with period <math display="inline"> <mrow> <mi>τ</mi> <mo>(</mo> <msub> <mi>z</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> </math>, where <math display="inline"> <mrow> <mi>τ</mi> <mo>(</mo> <mi>Z</mi> <mo>)</mo> </mrow> </math> denotes the corpuscle round-trip time at a distance <span class="html-italic">Z</span> from the top of the trajectory. When the altitude is restricted to <span class="html-italic">h</span> by a plate (dashed horizontal line) the motion remains periodic with a period evidently equal to: <math display="inline"> <mrow> <mi>τ</mi> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> <mo>−</mo> <mi>τ</mi> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>m</mi> </msub> <mo>−</mo> <mi>h</mi> <mo>)</mo> </mrow> </mrow> </math>. Note that this expression holds even if the motion is not symmetric in time.</p>
Full article ">
284 KiB  
Article
Minimum-Information-Entropy-Based Control Performance Assessment
by Qing-Wei Meng, Fang Fang and Ji-Zhen Liu
Entropy 2013, 15(3), 943-959; https://doi.org/10.3390/e15030943 - 27 Feb 2013
Cited by 13 | Viewed by 5289
Abstract
Generally, the controller design should be performed to narrow the shape of the probability density function of the tracking error. A small information entropy value corresponds to a narrow distribution function, which means that the uncertainty of the related random variable is small. [...] Read more.
Generally, the controller design should be performed to narrow the shape of the probability density function of the tracking error. A small information entropy value corresponds to a narrow distribution function, which means that the uncertainty of the related random variable is small. In this paper, information entropy is introduced in the field of control performance assessment (CPA). For the unknown time delay case, the minimum information entropy (MIE) benchmark is presented, and a MIE-based performance index is defined. For the known time delay case, a tight upper bound of MIE is derived and adopted as a performance benchmark to assess the stochastic control performance. Based on these, the control performance assessment procedures are developed for both the steady and the transient processes. Simulation tests and an industrial case study of a main steam pressure system of a 1,000MW power unit are utilized to verify the effectiveness of the proposed procedures. Full article
Show Figures

Figure 1

Figure 1
<p>The procedure of the transient CPA.</p>
Full article ">Figure 2
<p>(<b>a</b>)The tracking error; (<b>b</b>) the PDF.</p>
Full article ">Figure 3
<p>Comparison of the entropy estimations of the tracking error.</p>
Full article ">Figure 4
<p>Comparison of upper bound of the MIE.</p>
Full article ">Figure 5
<p>Comparison of the MIE index.</p>
Full article ">Figure 6
<p>Estimation of the stochastic disturbance.</p>
Full article ">
1336 KiB  
Article
Mutual Information-Based Inputs Selection for Electric Load Time Series Forecasting
by Miloš Božić, Miloš Stojanović, Zoran Stajić and Nenad Floranović
Entropy 2013, 15(3), 926-942; https://doi.org/10.3390/e15030926 - 27 Feb 2013
Cited by 23 | Viewed by 5666
Abstract
Providing accurate load forecast to electric utility corporations is essential in order to reduce their operational costs and increase profits. Hence, training set selection is an important preprocessing step which has to be considered in practice in order to increase the accuracy of [...] Read more.
Providing accurate load forecast to electric utility corporations is essential in order to reduce their operational costs and increase profits. Hence, training set selection is an important preprocessing step which has to be considered in practice in order to increase the accuracy of load forecasts. The usage of mutual information (MI) has been recently proposed in regression tasks, mostly for feature selection and for identifying the real instances from training sets that contains noise and outliers. This paper proposes a methodology for the training set selection in a least squares support vector machines (LS-SVMs) load forecasting model. A new application of the concept of MI is presented for the selection of a training set based on MI computation between initial training set instances and testing set instances. Accordingly, several LS-SVMs models have been trained, based on the proposed methodology, for hourly prediction of electric load for one day ahead. The results obtained from a real-world data set indicate that the proposed method increases the accuracy of load forecasting as well as reduces the size of the initial training set needed for model training. Full article
Show Figures

Figure 1

Figure 1
<p>Hourly load for September 2010.</p>
Full article ">Figure 2
<p>Hourly load during the week.</p>
Full article ">Figure 3
<p>Hourly load during the day.</p>
Full article ">Figure 4
<p>Architecture of the proposed inputs selection algorithm and forecasting strategy.</p>
Full article ">Figure 5
<p>Daily MAPEs for all of the generated models. (<b>a</b>) “MI threshold” selection option with 0.5 for M1, 0.8 for M2 and 0.98 for M3 (TH<sub>1</sub>), (<b>b</b>) “MI threshold” selection option with 0.6 for M1, 0.9 for M2 and 0.99 for M3 (TH<sub>2</sub>), (<b>c</b>) “number of inputs” selection option with 50 vectors for M1–M3 (NI<sub>1</sub>), (<b>d</b>) “number of inputs” selection option with 100 vectors for M1–M3 (NI<sub>2</sub>), (<b>e</b>) models without input selection with M1 for comparison.</p>
Full article ">Figure 5 Cont.
<p>Daily MAPEs for all of the generated models. (<b>a</b>) “MI threshold” selection option with 0.5 for M1, 0.8 for M2 and 0.98 for M3 (TH<sub>1</sub>), (<b>b</b>) “MI threshold” selection option with 0.6 for M1, 0.9 for M2 and 0.99 for M3 (TH<sub>2</sub>), (<b>c</b>) “number of inputs” selection option with 50 vectors for M1–M3 (NI<sub>1</sub>), (<b>d</b>) “number of inputs” selection option with 100 vectors for M1–M3 (NI<sub>2</sub>), (<b>e</b>) models without input selection with M1 for comparison.</p>
Full article ">Figure 6
<p>Vector number distribution in the initial training for the first hour of each day. (<b>a</b>) by <span class="html-italic">k</span>NN MI. (<b>b</b>) by kernel MI. (<b>c</b>) by correlation coefficient.</p>
Full article ">Figure 7
<p>Real and predicted loads of models M0 and M1. <b>(a)</b> period from September 17 to 23, (<b>b)</b> period from September 24 to 30.</p>
Full article ">
592 KiB  
Review
Quantum Models of Classical World
by Petr Hájíček
Entropy 2013, 15(3), 789-925; https://doi.org/10.3390/e15030789 - 27 Feb 2013
Cited by 2 | Viewed by 5184
Abstract
This paper is a review of our recent work on three notorious problems of non-relativistic quantum mechanics: realist interpretation, quantum theory of classical properties, and the problem of quantum measurement. A considerable progress has been achieved, based on four distinct new ideas. First, [...] Read more.
This paper is a review of our recent work on three notorious problems of non-relativistic quantum mechanics: realist interpretation, quantum theory of classical properties, and the problem of quantum measurement. A considerable progress has been achieved, based on four distinct new ideas. First, objective properties are associated with states rather than with values of observables. Second, all classical properties are selected properties of certain high entropy quantum states of macroscopic systems. Third, registration of a quantum system is strongly disturbed by systems of the same type in the environment. Fourth, detectors must be distinguished from ancillas and the states of registered systems are partially dissipated and lost in the detectors. The paper has two aims: a clear explanation of all new results and a coherent and contradiction-free account of the whole quantum mechanics including all necessary changes of its current textbook version. Full article
(This article belongs to the Special Issue Maximum Entropy and Bayes Theorem)
1163 KiB  
Article
Transfer Entropy for Coupled Autoregressive Processes
by Daniel W. Hahs and Shawn D. Pethel
Entropy 2013, 15(3), 767-788; https://doi.org/10.3390/e15030767 - 25 Feb 2013
Cited by 25 | Viewed by 8764
Abstract
A method is shown for computing transfer entropy over multiple time lags for coupled autoregressive processes using formulas for the differential entropy of multivariate Gaussian processes. Two examples are provided: (1) a first-order filtered noise process whose state is measured with additive noise, [...] Read more.
A method is shown for computing transfer entropy over multiple time lags for coupled autoregressive processes using formulas for the differential entropy of multivariate Gaussian processes. Two examples are provided: (1) a first-order filtered noise process whose state is measured with additive noise, and (2) two first-order coupled processes each of which is driven by white process noise. We found that, for the first example, increasing the first-order AR coefficient while keeping the correlation coefficient between filtered and measured process fixed, transfer entropy increased since the entropy of the measured process was itself increased. For the second example, the minimum correlation coefficient occurs when the process noise variances match. It was seen that matching of these variances results in minimum information flow, expressed as the sum of transfer entropies in both directions. Without a match, the transfer entropy is larger in the direction away from the process having the larger process noise. Fixing the process noise variances, transfer entropies in both directions increase with the coupling strength. Finally, we note that the method can be generally employed to compute other information theoretic quantities as well. Full article
(This article belongs to the Special Issue Transfer Entropy)
Show Figures

Figure 1

Figure 1
<p>Example 1: Transfer entropy TE<sup>(k)</sup><sub>x→y</sub> <span class="html-italic">versus</span> correlation coefficient ρ for three values of parameter <span class="html-italic">a</span> (see legend). Solid trace: k = 10, dotted trace: k = 2.</p>
Full article ">Figure 2
<p>Example 1: Logarithm of R <span class="html-italic">versus</span> ρ for three values of parameter <span class="html-italic">a</span> (see legend).</p>
Full article ">Figure 3
<p>Example 1: Process entropies H<sub>X</sub> and H<sub>Y</sub> <span class="html-italic">versus</span> correlation coefficient ρ for three values of parameter <span class="html-italic">a</span> (see legend).</p>
Full article ">Figure 4
<p>Example 1: Information transfer IT<sup>(k)</sup><sub>x→y</sub> <span class="html-italic">versus</span> correlation coefficient ρ for three different values of parameter <span class="html-italic">a</span> (see legend) for k = 10 (solid trace) and k = 2 (dotted trace).</p>
Full article ">Figure 5
<p>Example 1: Information transfer IT<sup>(10)</sup><sub>x→y</sub> <span class="html-italic">versus</span> measurement error variance R for three different values of parameter a (see legend).</p>
Full article ">Figure 6
<p>Example 2: Process noise variance R <span class="html-italic">versus</span> correlation coefficient ρ for a set of ε parameter values (see figure legend).</p>
Full article ">Figure 7
<p>Example 2: Transfer entropy values <span class="html-italic">versus</span> correlation ρ for a set of ε parameter values (see figure legend). Arrows indicate direction of increasing R values.</p>
Full article ">Figure 8
<p>Example 2: Transfer entropies difference (TE<sub>x− &gt;y</sub> – TE<sub>y − &gt; x</sub>) and sum (TE<sub>x− &gt; y</sub> + TE<sub>y− &gt; x</sub>) <span class="html-italic">versus</span> correlation ρ for a set of ε parameter values (see figure legend). Arrow indicates direction of increasing R values.</p>
Full article ">Figure 9
<p><span class="html-italic">Example 2:</span> Transfer entropies TE<sub>x→y</sub> and TE<sub>y→x <span class="html-italic">versus</span></sub> process noise variance R for a set of ε parameter values (see figure legend).</p>
Full article ">Figure 10
<p>Example 2: Transfer entropy TE<sub>x− &gt; y</sub> plotted <span class="html-italic">versus</span> TE<sub>y− &gt; x</sub> for a set of ε parameter values (see figure legend). The black diagonal line indicates locations where equality obtains. Arrow indicates direction of increasing R values.</p>
Full article ">Figure 11
<p>Example 2: Correlation coefficient ρ vs coupling coefficient ε for a set of R values (see figure legend).</p>
Full article ">Figure 12
<p>Example 2: Transfer entropies TE<sub>x→y</sub> (solid lines) vs TE<sub>y→x</sub> (dashed lines) vs coupling coefficient ε for a set of R values (see figure legend).</p>
Full article ">
177 KiB  
Article
Statistical Analysis of Gait Maturation in Children Using Nonparametric Probability Density Function Modeling
by Ning Xiang, Suxian Cai, Shanshan Yang, Zhangting Zhong, Fang Zheng, Jia He and Yunfeng Wu
Entropy 2013, 15(3), 753-766; https://doi.org/10.3390/e15030753 - 25 Feb 2013
Cited by 11 | Viewed by 8038
Abstract
Analysis of gait dynamics in children may help understand the development of neuromuscular control and maturation of locomotor function. This paper applied the nonparametric Parzen-window estimation method to establish the probability density function (PDF) models for the stride interval time series of 50 [...] Read more.
Analysis of gait dynamics in children may help understand the development of neuromuscular control and maturation of locomotor function. This paper applied the nonparametric Parzen-window estimation method to establish the probability density function (PDF) models for the stride interval time series of 50 children (25 boys and 25 girls). Four statistical parameters, in terms of averaged stride interval (ASI), variation of stride interval (VSI), PDF skewness (SK), and PDF kurtosis (KU), were computed with the Parzen-window PDFs to study the maturation of stride interval in children. By analyzing the results of the children in three age groups (aged 3–5 years, 6–8 years, and 10–14 years), we summarize the key findings of the present study as follows. (1) The gait cycle duration, in terms of ASI, increases until 14 years of age. On the other hand, the gait variability, in terms of VSI, decreases rapidly until 8 years of age, and then continues to decrease at a slower rate. (2) The SK values of both the histograms and Parzen-window PDFs for all of the three age groups are positive, which indicates an imbalance in the stride interval distribution within an age group. However, such an imbalance would be meliorated when the children grow up. (3) The KU values of both the histograms and Parzen-window PDFs decrease with the body growth in children, which suggests that the musculoskeletal growth enables the children to modulate a gait cadence with ease. (4) The SK and KU results also demonstrate the superiority of the Parzen-window PDF estimation method to the Gaussian distribution modeling, for the study of gait maturation in children. Full article
Show Figures

Figure 1

Figure 1
<p>Original time series of stride interval of the children (<b>a</b>) aged 45 months (in the group aged 3–5 years); (<b>b</b>) aged 80 months (in the group aged 6–8 years); and (<b>c</b>) aged 129 months (in the group aged 10–14 years), respectively. Outliers detected, along with one stride before or after the outliers, are marked with asterisks; Subfigures (<b>d</b>)–(<b>f</b>) plot the corresponding outlier-free time series of stride interval. The first strides in (<b>a</b>)–(<b>f</b>) start after the start-up 60 s (1 min).</p>
Full article ">Figure 2
<p>Histograms and Parzen-window probability density functions (PDFs) estimated for the stride interval time series of the children (<b>a</b>) aged 45 months (in the group aged 3–5 years); (<b>b</b>) aged 80 months (in the group aged 6–8 years); and (<b>c</b>) aged 129 months (in the group aged 10–14 years), respectively.</p>
Full article ">Figure 3
<p>Bar graphics of (<b>a</b>) averaged stride interval (ASI) and (<b>b</b>) variation stride interval (VSI) computed from the Parzen-window probability density functions (PDFs) of the children in the groups aged 3–5 years, 6–8 years, and 10–14 years, respectively. Vertical lines on the tops of the bars denote the standard deviation (SD) values. The values (mean ± SD) of the bars for the children of 3–5 years old (ASI: <math display="inline"> <mrow> <mn>0.904</mn> <mo>±</mo> <mn>0.041</mn> </mrow> </math> s, VSI: <math display="inline"> <mrow> <mn>0.058</mn> <mo>±</mo> <mn>0.015</mn> </mrow> </math> s), of 6–8 years old (ASI: <math display="inline"> <mrow> <mn>0.96</mn> <mo>±</mo> <mn>0.056</mn> </mrow> </math> s, VSI: <math display="inline"> <mrow> <mn>0.035</mn> <mo>±</mo> <mn>0.009</mn> </mrow> </math> s), and of 10–14 years old (ASI: <math display="inline"> <mrow> <mn>1.059</mn> <mo>±</mo> <mn>0.063</mn> </mrow> </math> s, VSI: <math display="inline"> <mrow> <mn>0.027</mn> <mo>±</mo> <mn>0.006</mn> </mrow> </math> s).</p>
Full article ">Figure 4
<p>Illustration of the probability density functions (PDFs) with three types of skewness (SK): right-skewed PDF, <span class="html-italic">dash-dot curve</span>, SK <math display="inline"> <mrow> <mo>&gt;</mo> <mn>0</mn> </mrow> </math>; symmetric PDF, <span class="html-italic">solid curve</span>, SK <math display="inline"> <mrow> <mo>=</mo> <mn>0</mn> </mrow> </math> (in particular Gaussian distributions); left-skewed PDF, <span class="html-italic">dashed curve</span>, SK <math display="inline"> <mrow> <mo>&lt;</mo> <mn>0</mn> </mrow> </math>. <span class="html-italic">Dot line</span> represents a central axis up to the mean of the symmetric PDF. au: arbitrary units.</p>
Full article ">Figure 5
<p>Illustration of four typical probability density functions (PDFs) that possess different kurtosis (KU) values: uniform distribution, <span class="html-italic">dot curve</span>, KU <math display="inline"> <mrow> <mo>=</mo> <mn>1</mn> <mo>.</mo> <mn>8</mn> </mrow> </math>; raised cosine distribution, <span class="html-italic">dashed curve</span>, KU <math display="inline"> <mrow> <mo>=</mo> <mn>2</mn> <mo>.</mo> <mn>41</mn> </mrow> </math>; Gaussian distribution, <span class="html-italic">solid curve</span>, KU <math display="inline"> <mrow> <mo>=</mo> <mn>3</mn> </mrow> </math>; hyperbolic secant distribution, <span class="html-italic">dash-dot curve</span>, KU <math display="inline"> <mrow> <mo>=</mo> <mn>5</mn> </mrow> </math>. au: arbitrary units.</p>
Full article ">Figure 6
<p>Bar graphics of the mean values of (<b>a</b>) skewness (SK) and (<b>b</b>) kurtosis (KU) computed from the histograms and the Parzen-window probability density functions (PDFs) of the 3- to 5-year-old, 6- to 8-year-old, and 10- to 14-year-old age groups, respectively. Statistics of the skewness and the kurtosis for the three age groups are listed in <a href="#entropy-15-00753-t001" class="html-table">Table 1</a>.</p>
Full article ">
637 KiB  
Article
Minimum Mutual Information and Non-Gaussianity through the Maximum Entropy Method: Estimation from Finite Samples
by Carlos A. L. Pires and Rui A. P. Perdigão
Entropy 2013, 15(3), 721-752; https://doi.org/10.3390/e15030721 - 25 Feb 2013
Cited by 8 | Viewed by 6978
Abstract
The Minimum Mutual Information (MinMI) Principle provides the least committed, maximum-joint-entropy (ME) inferential law that is compatible with prescribed marginal distributions and empirical cross constraints. Here, we estimate MI bounds (the MinMI values) generated by constraining sets Tcr comprehended by mcr [...] Read more.
The Minimum Mutual Information (MinMI) Principle provides the least committed, maximum-joint-entropy (ME) inferential law that is compatible with prescribed marginal distributions and empirical cross constraints. Here, we estimate MI bounds (the MinMI values) generated by constraining sets Tcr comprehended by mcr linear and/or nonlinear joint expectations, computed from samples of N iid outcomes. Marginals (and their entropy) are imposed by single morphisms of the original random variables. N-asymptotic formulas are given both for the distribution of cross expectation’s estimation errors, the MinMI estimation bias, its variance and distribution. A growing Tcr leads to an increasing MinMI, converging eventually to the total MI. Under N-sized samples, the MinMI increment relative to two encapsulated sets Tcr1 Tcr2 (with numbers of constraints mcr1<mcr2 ) is the test-difference δH = Hmax 1, N - Hmax 2, N ≥ 0 between the two respective estimated MEs. Asymptotically, δH follows a Chi-Squared distribution 1/2NΧ2 (mcr2-mcr1) whose upper quantiles determine if constraints in Tcr2/Tcr1 explain significant extra MI. As an example, we have set marginals to being normally distributed (Gaussian) and have built a sequence of MI bounds, associated to successive non-linear correlations due to joint non-Gaussianity. Noting that in real-world situations available sample sizes can be rather low, the relationship between MinMI bias, probability density over-fitting and outliers is put in evidence for under-sampled data. Full article
(This article belongs to the Special Issue Estimating Information-Theoretic Quantities from Data)
Show Figures

Figure 1

Figure 1
<p>Squared empirical bias: <inline-formula> <mml:math display="block" id="mm416"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mrow> <mml:mrow> <mml:mo>‖</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>b</mml:mi> </mml:mstyle> <mml:mo>‖</mml:mo> </mml:mrow> </mml:mrow> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (black lines) of <italic>N</italic>-based <inline-formula> <mml:math display="block" id="mm417"> <mml:semantics> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> </mml:semantics> </mml:math> </inline-formula>- expectations as function of <italic>N</italic>, empirical variances: <inline-formula> <mml:math display="block" id="mm418"> <mml:semantics> <mml:mrow> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:msub> <mml:mi>E</mml:mi> <mml:mi>N</mml:mi> </mml:msub> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">)</mml:mo> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (red lines), approximated variances: <inline-formula> <mml:math display="block" id="mm419"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mi>N</mml:mi> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">|</mml:mo> <mml:mi>l</mml:mi> <mml:mi>m</mml:mi> <mml:mi>s</mml:mi> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (blue lines) and variance for the case of <italic>N iid</italic> trials: <inline-formula> <mml:math display="block" id="mm420"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mi>N</mml:mi> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (green lines). <inline-formula> <mml:math display="block" id="mm421"> <mml:semantics> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> </mml:semantics> </mml:math> </inline-formula> stands for different bivariate monomials: <inline-formula> <mml:math display="block" id="mm422"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mi>X</mml:mi> <mml:mn>4</mml:mn> </mml:msup> <mml:msup> <mml:mi>Y</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (a), <inline-formula> <mml:math display="block" id="mm423"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mi>X</mml:mi> <mml:mn>6</mml:mn> </mml:msup> <mml:msup> <mml:mi>Y</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (b) and <inline-formula> <mml:math display="block" id="mm424"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mi>X</mml:mi> <mml:mn>8</mml:mn> </mml:msup> <mml:msup> <mml:mi>Y</mml:mi> <mml:mn>2</mml:mn> </mml:msup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (c).</p>
Full article ">Figure 2
<p><italic>N</italic> times Monte-Carlo variances: <inline-formula> <mml:math display="block" id="mm447"> <mml:semantics> <mml:mrow> <mml:mi>N</mml:mi> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:msub> <mml:mi>E</mml:mi> <mml:mi>N</mml:mi> </mml:msub> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">)</mml:mo> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> thick solid lines) and its theoretical analytical value <inline-formula> <mml:math display="block" id="mm448"> <mml:semantics> <mml:mrow> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">|</mml:mo> <mml:mi>l</mml:mi> <mml:mi>m</mml:mi> <mml:mi>s</mml:mi> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (thick dashed lines), both under imposed marginals (morphisms) and analytical value of <inline-formula> <mml:math display="block" id="mm449"> <mml:semantics> <mml:mrow> <mml:mi>N</mml:mi> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:msub> <mml:mi>E</mml:mi> <mml:mi>N</mml:mi> </mml:msub> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">)</mml:mo> <mml:mo stretchy="false">)</mml:mo> <mml:mo>=</mml:mo> <mml:mi>var</mml:mi> <mml:mo stretchy="false">(</mml:mo> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> <mml:mo stretchy="false">)</mml:mo> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> for <italic>iid</italic> data (<bold>thin solid lines</bold>). <inline-formula> <mml:math display="block" id="mm450"> <mml:semantics> <mml:mstyle mathvariant="bold" mathsize="normal"> <mml:mi>T</mml:mi> </mml:mstyle> </mml:semantics> </mml:math> </inline-formula> means different bivariate monomials: <inline-formula> <mml:math display="block" id="mm451"> <mml:semantics> <mml:mrow> <mml:mi>X</mml:mi> <mml:mi>Y</mml:mi> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<bold>black curves</bold>), <inline-formula> <mml:math display="block" id="mm452"> <mml:semantics> <mml:mrow> <mml:msup> <mml:mi>X</mml:mi> <mml:mn>2</mml:mn> </mml:msup> <mml:mi>Y</mml:mi> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<bold>red curves</bold>). <italic>N</italic> = 200.</p>
Full article ">Figure 3
<p>Test statistics: bias (black lines), standard deviation (red lines) and 95%-quantiles (green lines), provided by the Monte-Carlo approach (tick full lines), the semi-analytical approach (thin dashed lines) and the analytical approach (tick full lines). The tests are <inline-formula> <mml:math display="block" id="mm622"> <mml:semantics> <mml:mrow> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>g</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (a); <inline-formula> <mml:math display="block" id="mm623"> <mml:semantics> <mml:mrow> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>n</mml:mi> <mml:mi>g</mml:mi> <mml:mo>,</mml:mo> <mml:mi>j</mml:mi> <mml:mo>=</mml:mo> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (b); <inline-formula> <mml:math display="block" id="mm624"> <mml:semantics> <mml:mrow> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>n</mml:mi> <mml:mi>g</mml:mi> <mml:mo>,</mml:mo> <mml:mi>j</mml:mi> <mml:mo>=</mml:mo> <mml:mn>6</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (c) and <inline-formula> <mml:math display="block" id="mm625"> <mml:semantics> <mml:mrow> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>n</mml:mi> <mml:mi>g</mml:mi> <mml:mo>,</mml:mo> <mml:mi>j</mml:mi> <mml:mo>=</mml:mo> <mml:mn>8</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (d).</p>
Full article ">Figure 4
<p>Monte-Carlo empirical cumulative histogram (solid lines) and theoretical cumulative Chi-Squared fit (dashed lines) normalized by <italic>N</italic>: <inline-formula> <mml:math display="block" id="mm644"> <mml:semantics> <mml:mrow> <mml:mn>2</mml:mn> <mml:mi>N</mml:mi> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>g</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<inline-formula> <mml:math display="block" id="mm645"> <mml:semantics> <mml:mrow> <mml:msubsup> <mml:mi>χ</mml:mi> <mml:mn>1</mml:mn> <mml:mn>2</mml:mn> </mml:msubsup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula>) for <inline-formula> <mml:math display="block" id="mm646"> <mml:semantics> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>=</mml:mo> <mml:mn>50</mml:mn> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<bold>black curves</bold>); <inline-formula> <mml:math display="block" id="mm647"> <mml:semantics> <mml:mrow> <mml:mn>2</mml:mn> <mml:mi>N</mml:mi> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>n</mml:mi> <mml:mi>g</mml:mi> <mml:mo>,</mml:mo> <mml:mi>j</mml:mi> <mml:mo>=</mml:mo> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<inline-formula> <mml:math display="block" id="mm648"> <mml:semantics> <mml:mrow> <mml:msubsup> <mml:mi>χ</mml:mi> <mml:mn>5</mml:mn> <mml:mn>2</mml:mn> </mml:msubsup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula>) for <inline-formula> <mml:math display="block" id="mm649"> <mml:semantics> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>=</mml:mo> <mml:mn>400</mml:mn> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<bold>red curves</bold>); <inline-formula> <mml:math display="block" id="mm650"> <mml:semantics> <mml:mrow> <mml:mn>2</mml:mn> <mml:mi>N</mml:mi> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>n</mml:mi> <mml:mi>g</mml:mi> <mml:mo>,</mml:mo> <mml:mn>6</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<inline-formula> <mml:math display="block" id="mm651"> <mml:semantics> <mml:mrow> <mml:msubsup> <mml:mi>χ</mml:mi> <mml:mrow> <mml:mn>14</mml:mn> </mml:mrow> <mml:mn>2</mml:mn> </mml:msubsup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula>) for <inline-formula> <mml:math display="block" id="mm652"> <mml:semantics> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>=</mml:mo> <mml:mn>1600</mml:mn> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<bold>green curves</bold>) and <inline-formula> <mml:math display="block" id="mm653"> <mml:semantics> <mml:mrow> <mml:mn>2</mml:mn> <mml:mi>N</mml:mi> <mml:mi>δ</mml:mi> <mml:msub> <mml:mi>I</mml:mi> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>,</mml:mo> <mml:mi>n</mml:mi> <mml:mi>g</mml:mi> <mml:mo>,</mml:mo> <mml:mn>8</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<inline-formula> <mml:math display="block" id="mm654"> <mml:semantics> <mml:mrow> <mml:msubsup> <mml:mi>χ</mml:mi> <mml:mrow> <mml:mn>27</mml:mn> </mml:mrow> <mml:mn>2</mml:mn> </mml:msubsup> </mml:mrow> </mml:semantics> </mml:math> </inline-formula>) for <inline-formula> <mml:math display="block" id="mm655"> <mml:semantics> <mml:mrow> <mml:mi>N</mml:mi> <mml:mo>=</mml:mo> <mml:mn>3200</mml:mn> </mml:mrow> </mml:semantics> </mml:math> </inline-formula> (<bold>blue curves</bold>).</p>
Full article ">Figure 5
<p>Scatter-plot of the Gaussianized variables <italic>X</italic> (in abscissas) <italic>Y</italic> (in ordinates) (see text for details) along with ME-PDF fitting constrained by monomial bivariate moments up to order <italic>j</italic> = 2 (<bold>a</bold>), <italic>j</italic> = 4 (<bold>b</bold>), <italic>j</italic> = 6 (<bold>c</bold>) and <italic>j</italic> = 8 (<bold>d</bold>). Contour levels are set to 0.0005, 0.005, 0.05, 0.5, and 5.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop