[go: up one dir, main page]

Next Issue
Volume 22, February
Previous Issue
Volume 21, December
 
 
entropy-logo

Journal Browser

Journal Browser

Entropy, Volume 22, Issue 1 (January 2020) – 127 articles

Cover Story (view full-size image): Optimal lossy data compression minimizes the storage cost of a data set X while retaining a given amount of information as possible about something (Y) that you care about. For example, what aspects of an image X contain the most information about whether it depicts a cat or not?
We present a method for efficiently solving this problem and apply it to the CIFAR-10, MNIST, and Fashion-MNIST datasets, illustrating how it can be interpreted as an information-theoretically optimal image clustering algorithm. For merely a handful of clusters, compressing the image into a single integer specifying its cluster number retains almost all the information about its class label Y. View this paper.
  • 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:
13 pages, 4382 KiB  
Article
Eigenvalues of Two-State Quantum Walks Induced by the Hadamard Walk
by Shimpei Endo, Takako Endo, Takashi Komatsu and Norio Konno
Entropy 2020, 22(1), 127; https://doi.org/10.3390/e22010127 - 20 Jan 2020
Cited by 7 | Viewed by 4316
Abstract
Existence of the eigenvalues of the discrete-time quantum walks is deeply related to localization of the walks. We revealed, for the first time, the distributions of the eigenvalues given by the splitted generating function method (the SGF method) of the space-inhomogeneous quantum walks [...] Read more.
Existence of the eigenvalues of the discrete-time quantum walks is deeply related to localization of the walks. We revealed, for the first time, the distributions of the eigenvalues given by the splitted generating function method (the SGF method) of the space-inhomogeneous quantum walks in one dimension we had treated in our previous studies. Especially, we clarified the characteristic parameter dependence for the distributions of the eigenvalues with the aid of numerical simulation. Full article
(This article belongs to the Special Issue Quantum Walks and Related Issues)
Show Figures

Figure 1

Figure 1
<p>Pattern diagram of the QW.</p>
Full article ">Figure 2
<p>Arrangement of the unitary operators of the Wojcik model.</p>
Full article ">Figure 3
<p>Arrangement of the unitary operators of the Hadamard walk with one defect.</p>
Full article ">Figure 4
<p>The illustrations of the eigenvalues movements of the Hadamard walk with one defect <math display="inline"><semantics> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> <mo> </mo> <mi>β</mi> <mo>=</mo> <mo>−</mo> <mi>i</mi> <mi>α</mi> </mrow> </semantics></math> case. <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ξ</mi> <mo>)</mo> </mrow> <mo>,</mo> <mo> </mo> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ξ</mi> <mo>)</mo> </mrow> <mo> </mo> <mrow> <mo>(</mo> <mi>ξ</mi> <mo>∈</mo> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mi>π</mi> <mn>4</mn> </mfrac> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math> <math display="inline"><semantics> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> <mo> </mo> <mi>β</mi> <mo>=</mo> <mi>i</mi> <mi>α</mi> </mrow> </semantics></math> case. <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ξ</mi> <mo>)</mo> </mrow> <mo>,</mo> <mo> </mo> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ξ</mi> <mo>)</mo> </mrow> <mo> </mo> <mrow> <mo>(</mo> <mi>ξ</mi> <mo>∈</mo> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mi>π</mi> <mn>4</mn> </mfrac> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math> (The red part is the region of the continuous spectrum of the Hadamard walk, i.e., <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>(</mo> <mi>H</mi> <mo>)</mo> </mrow> </semantics></math>).</p>
Full article ">Figure 5
<p>Arrangement of the unitary operators of the two-phase QW with one defect.</p>
Full article ">Figure 6
<p>The illustrations of the eigenvalues movements of the Wojcik model <math display="inline"><semantics> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> <mo> </mo> <mi>β</mi> <mo>=</mo> <mi>i</mi> <mi>α</mi> </mrow> </semantics></math> case. <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ϕ</mi> <mo>)</mo> </mrow> <mo>,</mo> <mo> </mo> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ϕ</mi> <mo>)</mo> </mrow> <mo> </mo> <mrow> <mo>(</mo> <mi>ϕ</mi> <mo>∈</mo> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mn>4</mn> </mfrac> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math> <math display="inline"><semantics> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> <mo> </mo> <mi>β</mi> <mo>=</mo> <mo>−</mo> <mi>i</mi> <mi>α</mi> </mrow> </semantics></math> case. <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ϕ</mi> <mo>)</mo> </mrow> <mo>,</mo> <mo> </mo> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>ϕ</mi> <mo>)</mo> </mrow> <mo> </mo> <mrow> <mo>(</mo> <mi>ϕ</mi> <mo>∈</mo> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mn>3</mn> <mn>4</mn> </mfrac> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math> (The red part is the region of the continuous spectrum of the Hadamard walk, i.e., <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>(</mo> <mi>H</mi> <mo>)</mo> </mrow> </semantics></math>).</p>
Full article ">Figure 7
<p>The illustration of the eigenvalues movements of the two-phase quantum walk with one defect For <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </semantics></math> we have <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>∈</mo> <mrow> <mo>[</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mn>5</mn> <mn>4</mn> </mfrac> <mi>π</mi> <mo>)</mo> </mrow> <mo>∪</mo> <mrow> <mo>(</mo> <mfrac> <mn>7</mn> <mn>4</mn> </mfrac> <mi>π</mi> <mo>,</mo> <mn>2</mn> <mi>π</mi> <mo>]</mo> </mrow> <mo>.</mo> </mrow> </semantics></math> For <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </semantics></math> we have <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>∈</mo> <mrow> <mo>[</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mn>1</mn> <mn>4</mn> </mfrac> <mi>π</mi> <mo>)</mo> </mrow> <mo>∪</mo> <mrow> <mo>(</mo> <mfrac> <mn>3</mn> <mn>4</mn> </mfrac> <mi>π</mi> <mo>,</mo> <mn>2</mn> <mi>π</mi> <mo>]</mo> </mrow> <mo>.</mo> </mrow> </semantics></math> (The red part is the region of the continuous spectrum of the Hadamard walk, i.e., <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>(</mo> <mi>H</mi> <mo>)</mo> </mrow> </semantics></math>).</p>
Full article ">Figure 8
<p>Arrangement of the unitary operators of the two-phase QW with one defect.</p>
Full article ">Figure 9
<p>The illustration of the eigenvalues movements of the complete two-phase quantum walk For <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </semantics></math> we have <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>∈</mo> <mrow> <mo>[</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>π</mi> <mo>,</mo> <mi>π</mi> <mo>)</mo> </mrow> <mo>∪</mo> <mrow> <mo>(</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <mi>π</mi> <mo>,</mo> <mn>2</mn> <mi>π</mi> <mo>]</mo> </mrow> <mo>.</mo> </mrow> </semantics></math> For <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msup> <mi>λ</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <mi>σ</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </semantics></math> we have <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>∈</mo> <mrow> <mo>[</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>π</mi> <mo>)</mo> </mrow> <mo>∪</mo> <mrow> <mo>(</mo> <mi>π</mi> <mo>,</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <mi>π</mi> <mo>]</mo> </mrow> <mo>.</mo> </mrow> </semantics></math> (The red part is the region of the continuous spectrum of the Hadamard walk, i.e., <math display="inline"><semantics> <mrow> <mi>σ</mi> <mo>(</mo> <mi>H</mi> <mo>)</mo> </mrow> </semantics></math>).</p>
Full article ">
14 pages, 636 KiB  
Article
A Standardized Project Gutenberg Corpus for Statistical Analysis of Natural Language and Quantitative Linguistics
by Martin Gerlach and Francesc Font-Clos
Entropy 2020, 22(1), 126; https://doi.org/10.3390/e22010126 - 20 Jan 2020
Cited by 43 | Viewed by 9453
Abstract
The use of Project Gutenberg (PG) as a text corpus has been extremely popular in statistical analysis of language for more than 25 years. However, in contrast to other major linguistic datasets of similar importance, no consensual full version of PG exists to [...] Read more.
The use of Project Gutenberg (PG) as a text corpus has been extremely popular in statistical analysis of language for more than 25 years. However, in contrast to other major linguistic datasets of similar importance, no consensual full version of PG exists to date. In fact, most PG studies so far either consider only a small number of manually selected books, leading to potential biased subsets, or employ vastly different pre-processing strategies (often specified in insufficient details), raising concerns regarding the reproducibility of published results. In order to address these shortcomings, here we present the Standardized Project Gutenberg Corpus (SPGC), an open science approach to a curated version of the complete PG data containing more than 50,000 books and more than 3 × 10 9 word-tokens. Using different sources of annotated metadata, we not only provide a broad characterization of the content of PG, but also show different examples highlighting the potential of SPGC for investigating language variability across time, subjects, and authors. We publish our methodology in detail, the code to download and process the data, as well as the obtained corpus itself on three different levels of granularity (raw text, timeseries of word tokens, and counts of words). In this way, we provide a reproducible, pre-processed, full-size version of Project Gutenberg as a new scientific resource for corpus linguistics, natural language processing, and information retrieval. Full article
(This article belongs to the Special Issue Information Theory and Language)
Show Figures

Figure 1

Figure 1
<p>Sketch of the pre-processing pipeline of the Project Gutenberg (PG) data. The folder structure (<b>left</b>) organizes each PG book on four different levels of granularity, see example books (<b>middle</b>): raw, text, tokens, and counts. On the right we show the basic python commands used in the pre-processing.</p>
Full article ">Figure 2
<p>Basic summary statistics from the processed PG data. (<b>a</b>) Number of books with a text length larger than <span class="html-italic">m</span>; (<b>b</b>) Number of books which are compatible with being published in year <span class="html-italic">t</span>, i.e., year of author’s birth is 20 years prior and year of author’s death is after <span class="html-italic">t</span>; (<b>c</b>) Number of books (left axis) and number of tokens (right axis) which are assigned to a given language based on the metadata. en: English, fr: French, fi: Finnish, de: German, nl: Dutch, it: Italian, es: Spanish, pt: Portuguese, zh: Chinese, el: Greek, sv: Swedish, hu: Hungarian, eo: Esperanto, la: Latin, da: Danish, tl: Tagalog, ca: Catalan, pl: Polish, ja: Japanese, no: Norwegian, cy: Welsh, cs: Czech.</p>
Full article ">Figure 3
<p>Comparison between bookshelf labels (top, green) and subject labels (bottom, red). (<b>a</b>,<b>c</b>) Number of labels with a given number of books; (<b>b</b>,<b>d</b>) Fraction of books with a given number of labels.</p>
Full article ">Figure 4
<p>2-dimensional embedding shows clustering of books from the same bookshelf. Approximate visualization of the pair-wise distances between 370 PG books using Uniform Manifold Approximation and Projection (UMAP) (for details, see <a href="#sec4-entropy-22-00126" class="html-sec">Section 4</a>). Each dot corresponds to one book colored according to the bookshelf membership.</p>
Full article ">Figure 5
<p>Distance between books from the same author is significantly smaller than distance between books from different authors. For each author, the boxplots shows the <math display="inline"><semantics> <mrow> <mn>5</mn> <mo>,</mo> <mn>25</mn> <mo>,</mo> <mn>50</mn> <mo>,</mo> <mn>75</mn> <mo>,</mo> <mn>95</mn> </mrow> </semantics></math>-percentile of the distribution of distances from <math display="inline"><semantics> <mrow> <mn>1000</mn> </mrow> </semantics></math> pairs of books from the same author (<b>green</b>) and to a different author (<b>gray</b>).</p>
Full article ">Figure 6
<p>Distance between books increases with their time separation. Average and standard error of the distance between 1000 pairs of books, where the two books in each pair is drawn from two different 20-year time intervals. We fix the first interval and continuously increase the second time interval.</p>
Full article ">
19 pages, 2160 KiB  
Article
Low-Complexity Rate-Distortion Optimization of Sampling Rate and Bit-Depth for Compressed Sensing of Images
by Qunlin Chen, Derong Chen, Jiulu Gong and Jie Ruan
Entropy 2020, 22(1), 125; https://doi.org/10.3390/e22010125 - 20 Jan 2020
Cited by 5 | Viewed by 3432
Abstract
Compressed sensing (CS) offers a framework for image acquisition, which has excellent potential in image sampling and compression applications due to the sub-Nyquist sampling rate and low complexity. In engineering practices, the resulting CS samples are quantized by finite bits for transmission. In [...] Read more.
Compressed sensing (CS) offers a framework for image acquisition, which has excellent potential in image sampling and compression applications due to the sub-Nyquist sampling rate and low complexity. In engineering practices, the resulting CS samples are quantized by finite bits for transmission. In circumstances where the bit budget for image transmission is constrained, knowing how to choose the sampling rate and the number of bits per measurement (bit-depth) is essential for the quality of CS reconstruction. In this paper, we first present a bit-rate model that considers the compression performance of CS, quantification, and entropy coder. The bit-rate model reveals the relationship between bit rate, sampling rate, and bit-depth. Then, we propose a relative peak signal-to-noise ratio (PSNR) model for evaluating distortion, which reveals the relationship between relative PSNR, sampling rate, and bit-depth. Finally, the optimal sampling rate and bit-depth are determined based on the rate-distortion (RD) criteria with the bit-rate model and the relative PSNR model. The experimental results show that the actual bit rate obtained by the optimized sampling rate and bit-depth is very close to the target bit rate. Compared with the traditional CS coding method with a fixed sampling rate, the proposed method provides better rate-distortion performance, and the additional calculation amount amounts to less than 1%. Full article
(This article belongs to the Section Signal and Data Analysis)
Show Figures

Figure 1

Figure 1
<p>Proposed adaptive compressive sampling framework with rate-distortion optimization.</p>
Full article ">Figure 2
<p>Four-layer feedforward neural network model for the relative peak signal-to-noise ratio (PSNR).</p>
Full article ">Figure 3
<p>Four testing images. (<b>a</b>) Monarch; (<b>b</b>) Cameraman; (<b>c</b>) Peppers; (<b>d</b>) Lena.</p>
Full article ">Figure 3 Cont.
<p>Four testing images. (<b>a</b>) Monarch; (<b>b</b>) Cameraman; (<b>c</b>) Peppers; (<b>d</b>) Lena.</p>
Full article ">Figure 4
<p>Comparison of rate-distortion (RD) performances. (<b>a</b>) Monarch; (<b>b</b>) Cameraman; (<b>c</b>) Peppers; (<b>d</b>) Lena.</p>
Full article ">
4 pages, 1717 KiB  
Editorial
Entropy 2019 Best Paper Award
by Entropy Editorial Office
Entropy 2020, 22(1), 124; https://doi.org/10.3390/e22010124 - 20 Jan 2020
Viewed by 5280
Abstract
On behalf of the Editor-in-Chief, Prof [...] Full article
12 pages, 2083 KiB  
Article
Analyzing Uncertainty in Complex Socio-Ecological Networks
by Ana D. Maldonado, María Morales, Pedro A. Aguilera and Antonio Salmerón
Entropy 2020, 22(1), 123; https://doi.org/10.3390/e22010123 - 19 Jan 2020
Cited by 12 | Viewed by 3622
Abstract
Socio-ecological systems are recognized as complex adaptive systems whose multiple interactions might change as a response to external or internal changes. Due to its complexity, the behavior of the system is often uncertain. Bayesian networks provide a sound approach for handling complex domains [...] Read more.
Socio-ecological systems are recognized as complex adaptive systems whose multiple interactions might change as a response to external or internal changes. Due to its complexity, the behavior of the system is often uncertain. Bayesian networks provide a sound approach for handling complex domains endowed with uncertainty. The aim of this paper is to analyze the impact of the Bayesian network structure on the uncertainty of the model, expressed as the Shannon entropy. In particular, three strategies for model structure have been followed: naive Bayes (NB), tree augmented network (TAN) and network with unrestricted structure (GSS). Using these network structures, two experiments are carried out: (1) the impact of the Bayesian network structure on the entropy of the model is assessed and (2) the entropy of the posterior distribution of the class variable obtained from the different structures is compared. The results show that GSS constantly outperforms both NB and TAN when it comes to evaluating the uncertainty of the entire model. On the other hand, NB and TAN yielded lower entropy values of the posterior distribution of the class variable, which makes them preferable when the goal is to carry out predictions. Full article
(This article belongs to the Special Issue Computation in Complex Networks)
Show Figures

Figure 1

Figure 1
<p>An example of a Bayesian network structure with 5 variables.</p>
Full article ">Figure 2
<p>Structure of a naive Bayes model with <span class="html-italic">n</span> features (<b>a</b>) and a tree augmented network (TAN) model with 4 features (<b>b</b>).</p>
Full article ">Figure 3
<p>Structure of the Bayesian network used as reference in the experiments.</p>
Full article ">Figure 4
<p>Shannon entropy vs. sample size for the Bayesian networks used in Experiment 1.</p>
Full article ">Figure 5
<p>Shannon entropy of the class posterior distribution vs. sample size for scenario 1 in Experiment 2.</p>
Full article ">Figure 6
<p>Shannon entropy of the class posterior distribution vs. sample size for scenario 2 in Experiment 2.</p>
Full article ">Figure 7
<p>Shannon entropy of the class posterior distribution vs. sample size for scenario 3 in Experiment 2.</p>
Full article ">
15 pages, 826 KiB  
Article
Stabilization of Port Hamiltonian Chaotic Systems with Hidden Attractors by Adaptive Terminal Sliding Mode Control
by Ahmad Taher Azar and Fernando E. Serrano
Entropy 2020, 22(1), 122; https://doi.org/10.3390/e22010122 - 19 Jan 2020
Cited by 25 | Viewed by 3690
Abstract
In this study, the design of an adaptive terminal sliding mode controller for the stabilization of port Hamiltonian chaotic systems with hidden attractors is proposed. This study begins with the design methodology of a chaotic oscillator with a hidden attractor implementing the topological [...] Read more.
In this study, the design of an adaptive terminal sliding mode controller for the stabilization of port Hamiltonian chaotic systems with hidden attractors is proposed. This study begins with the design methodology of a chaotic oscillator with a hidden attractor implementing the topological framework for its respective design. With this technique it is possible to design a 2-D chaotic oscillator, which is then converted into port-Hamiltonia to track and analyze these models for the stabilization of the hidden chaotic attractors created by this analysis. Adaptive terminal sliding mode controllers (ATSMC) are built when a Hamiltonian system has a chaotic behavior and a hidden attractor is detected. A Lyapunov approach is used to formulate the adaptive device controller by creating a control law and the adaptive law, which are used online to make the system states stable while at the same time suppressing its chaotic behavior. The empirical tests obtaining the discussion and conclusions of this thesis should verify the theoretical findings. Full article
Show Figures

Figure 1

Figure 1
<p>Phase portrait of <math display="inline"><semantics> <msub> <mi>p</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>p</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 2
<p>Phase portrait of <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Phase portrait of <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Phase portrait of <math display="inline"><semantics> <msub> <mi>p</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>p</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 5
<p>Sliding surface <math display="inline"><semantics> <msub> <mi>σ</mi> <mn>3</mn> </msub> </semantics></math> fo the experiment 1.</p>
Full article ">Figure 6
<p>Evolution in time of the variables <math display="inline"><semantics> <msub> <mi>q</mi> <mn>2</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>p</mi> <mn>1</mn> </msub> </semantics></math>.</p>
Full article ">Figure 7
<p>Input variables <math display="inline"><semantics> <msub> <mi>u</mi> <mn>2</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>u</mi> <mn>4</mn> </msub> </semantics></math>.</p>
Full article ">Figure 8
<p>Evolution in time of the gain variable <math display="inline"><semantics> <msub> <mi>k</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 9
<p>Evolution in time of the variables <math display="inline"><semantics> <msub> <mover accent="true"> <mi>q</mi> <mo>˙</mo> </mover> <mn>3</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mover accent="true"> <mi>p</mi> <mo>˙</mo> </mover> <mn>3</mn> </msub> </semantics></math>.</p>
Full article ">Figure 10
<p>Evolution in time of the input variables <math display="inline"><semantics> <msub> <mi>u</mi> <mn>3</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>u</mi> <mn>6</mn> </msub> </semantics></math>.</p>
Full article ">Figure 11
<p>Evolution in time of the gain variable <math display="inline"><semantics> <msub> <mi>k</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 12
<p>Evolution in time of the variable <math display="inline"><semantics> <msub> <mi>σ</mi> <mn>3</mn> </msub> </semantics></math>.</p>
Full article ">
27 pages, 1087 KiB  
Article
Linear Programming and Fuzzy Optimization to Substantiate Investment Decisions in Tangible Assets
by Marcel-Ioan Boloș, Ioana-Alexandra Bradea and Camelia Delcea
Entropy 2020, 22(1), 121; https://doi.org/10.3390/e22010121 - 19 Jan 2020
Cited by 4 | Viewed by 4338
Abstract
This paper studies the problem of tangible assets acquisition within the company by proposing a new hybrid model that uses linear programming and fuzzy numbers. Regarding linear programming, two methods were implemented in the model, namely: the graphical method and the primal simplex [...] Read more.
This paper studies the problem of tangible assets acquisition within the company by proposing a new hybrid model that uses linear programming and fuzzy numbers. Regarding linear programming, two methods were implemented in the model, namely: the graphical method and the primal simplex algorithm. This hybrid model is proposed for solving investment decision problems, based on decision variables, objective function coefficients, and a matrix of constraints, all of them presented in the form of triangular fuzzy numbers. Solving the primal simplex algorithm using fuzzy numbers and coefficients, allowed the results of the linear programming problem to also be in the form of fuzzy variables. The fuzzy variables compared to the crisp variables allow the determination of optimal intervals for which the objective function has values depending on the fuzzy variables. The major advantage of this model is that the results are presented as value ranges that intervene in the decision-making process. Thus, the company’s decision makers can select any of the result values as they satisfy two basic requirements namely: minimizing/maximizing the objective function and satisfying the basic requirements regarding the constraints resulting from the company’s activity. The paper is accompanied by a practical example. Full article
Show Figures

Figure 1

Figure 1
<p>The triangular fuzzy number <span class="html-italic">C</span> used in fuzzy modeling.</p>
Full article ">Figure 2
<p>The graphical solution of the linear programming method.</p>
Full article ">Figure 3
<p>The flow chart of problem solving using simplex algorithms with fuzzy coefficients.</p>
Full article ">
12 pages, 385 KiB  
Article
Generalized Independence in the q-Voter Model: How Do Parameters Influence the Phase Transition?
by Angelika Abramiuk and Katarzyna Sznajd-Weron
Entropy 2020, 22(1), 120; https://doi.org/10.3390/e22010120 - 19 Jan 2020
Cited by 10 | Viewed by 3259
Abstract
We study the q-voter model with flexibility, which allows for describing a broad spectrum of independence from zealots, inflexibility, or stubbornness through noisy voters to self-anticonformity. Analyzing the model within the pair approximation allows us to derive the analytical formula for the [...] Read more.
We study the q-voter model with flexibility, which allows for describing a broad spectrum of independence from zealots, inflexibility, or stubbornness through noisy voters to self-anticonformity. Analyzing the model within the pair approximation allows us to derive the analytical formula for the critical point, below which an ordered (agreement) phase is stable. We determine the role of flexibility, which can be understood as an amount of variability associated with an independent behavior, as well as the role of the average network degree in shaping the character of the phase transition. We check the existence of the scaling relation, which previously was derived for the Sznajd model. We show that the scaling is universal, in a sense that it does not depend neither on the size of the group of influence nor on the average network degree. Analyzing the model in terms of the rescaled parameter, we determine the critical point, the jump of the order parameter, as well as the width of the hysteresis as a function of the average network degree k and the size of the group of influence q. Full article
(This article belongs to the Special Issue Statistical Mechanics of Complex Systems)
Show Figures

Figure 1

Figure 1
<p>Dependence between the stationary value of the ratio of positive votes <span class="html-italic">c</span> and the probability of independence <span class="html-italic">p</span> for different values of flexibility <span class="html-italic">f</span> for <math display="inline"><semantics> <mrow> <mi>q</mi> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math> (upper panels) and <math display="inline"><semantics> <mrow> <mi>q</mi> <mo>=</mo> <mn>6</mn> </mrow> </semantics></math> (bottom panels). Original, unscaled results are shown in the left panels, whereas rescaled results, in terms of <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mi>x</mi> <mo>(</mo> <mi>p</mi> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> </semantics></math>, defined by Equation (<a href="#FD9-entropy-22-00120" class="html-disp-formula">9</a>), are presented in the right ones. In this example, the average degree of a graph <math display="inline"><semantics> <mrow> <mo>〈</mo> <mi>k</mi> <mo>〉</mo> <mo>=</mo> <mn>50</mn> </mrow> </semantics></math>. Solid and dashed lines correspond to stable and unstable steady states, respectively.</p>
Full article ">Figure 2
<p>The stationary value of the ratio of positive votes <span class="html-italic">c</span> as a function of <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>=</mo> <mi>x</mi> <mo>(</mo> <mi>p</mi> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> </semantics></math>, defined by Equation (<a href="#FD9-entropy-22-00120" class="html-disp-formula">9</a>), for: different sizes of the group of influence <span class="html-italic">q</span> and the fixed average degree <math display="inline"><semantics> <mrow> <mo>〈</mo> <mi>k</mi> <mo>〉</mo> <mo>=</mo> <mn>50</mn> </mrow> </semantics></math> (left panel), and different average degree <math display="inline"><semantics> <mrow> <mo>〈</mo> <mi>k</mi> <mo>〉</mo> </mrow> </semantics></math> and the fixed size of the group of influence <math display="inline"><semantics> <mrow> <mi>q</mi> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math> (right panel). Solid and dashed lines correspond to stable and unstable steady states, respectively.</p>
Full article ">Figure 3
<p>The critical value of the rescaled independence <math display="inline"><semantics> <msubsup> <mi>x</mi> <mn>1</mn> <mo>*</mo> </msubsup> </semantics></math> (lower spinodal), below which disagreement is unstable, as a function of the size of the group of influence <span class="html-italic">q</span> and the average degree of a graph <math display="inline"><semantics> <mrow> <mo>〈</mo> <mi>k</mi> <mo>〉</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The size of the hysteresis as a function of the size of the group of influence <span class="html-italic">q</span> and the average degree of a graph <math display="inline"><semantics> <mrow> <mo>〈</mo> <mi>k</mi> <mo>〉</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>The jump of the public opinion at upper spinodal <math display="inline"><semantics> <msubsup> <mi>x</mi> <mn>2</mn> <mo>*</mo> </msubsup> </semantics></math> as a function of the size of the group of influence <span class="html-italic">q</span> and the average degree of a graph <math display="inline"><semantics> <mrow> <mo>〈</mo> <mi>k</mi> <mo>〉</mo> </mrow> </semantics></math>.</p>
Full article ">
20 pages, 1055 KiB  
Article
Activeness and Loyalty Analysis in Event-Based Social Networks
by Thanh Trinh, Dingming Wu, Joshua Zhexue Huang and Muhammad Azhar
Entropy 2020, 22(1), 119; https://doi.org/10.3390/e22010119 - 18 Jan 2020
Cited by 11 | Viewed by 3866
Abstract
Event-based social networks (EBSNs) are widely used to create online social groups and organize offline events for users. Activeness and loyalty are crucial characteristics of these online social groups in terms of determining the growth or inactiveness of the social groups in a [...] Read more.
Event-based social networks (EBSNs) are widely used to create online social groups and organize offline events for users. Activeness and loyalty are crucial characteristics of these online social groups in terms of determining the growth or inactiveness of the social groups in a specific time frame. However, there is less research on these concepts to clarify the existence of groups in event-based social networks. In this paper, we study the problem of group activeness and user loyalty to provide a novel insight into online social networks. First, we analyze the structure of EBSNs and generate features from the crawled datasets. Second, we define the concepts of group activeness and user loyalty based on a series of time windows, and propose a method to measure the group activeness. In this proposed method, we first compute a ratio of a number of events between two consecutive time windows. We then develop an association matrix to assign the activeness label for each group after several consecutive time windows. Similarly, we measure the user loyalty in terms of attended events gathered in time windows and treat loyalty as a contributive feature of the group activeness. Finally, three well-known machine learning techniques are used to verify the activeness label and to generate features for each group. As a consequence, we also find a small group of features that are highly correlated and result in higher accuracy as compared to the whole features. Full article
(This article belongs to the Special Issue Social Networks and Information Diffusion II)
Show Figures

Figure 1

Figure 1
<p>Example of an event-based social network (EBSN).</p>
Full article ">Figure 2
<p>Example of the time frame of group <span class="html-italic">G</span>.</p>
Full article ">Figure 3
<p>Category type distribution in groups created in four cities in the year 2014.</p>
Full article ">Figure 4
<p>Distribution of average days between two consecutive events in each group.</p>
Full article ">Figure 4 Cont.
<p>Distribution of average days between two consecutive events in each group.</p>
Full article ">Figure 5
<p>Example of the time window <span class="html-italic">T</span> in <math display="inline"><semantics> <msup> <mi>G</mi> <mn>1</mn> </msup> </semantics></math> and <math display="inline"><semantics> <msup> <mi>G</mi> <mn>2</mn> </msup> </semantics></math>.</p>
Full article ">Figure 6
<p>New York—The numbers of groups and events in 8 time windows.</p>
Full article ">Figure 7
<p>San Francisco—The numbers of groups and events in 8 time windows.</p>
Full article ">Figure 8
<p>London—The numbers of groups and events in 8 time windows.</p>
Full article ">Figure 9
<p>Sydney—The numbers of groups and events in 8 time windows.</p>
Full article ">Figure 10
<p>Distributions in terms of percentage and the number of attended events that users participate in among the total events of their groups in the two-year period. (<b>a</b>) New York—Percentage of attended events in total events for users. (<b>b</b>) New York—Number of attended events for users. (<b>c</b>) San Francisco—Percentage of attended events in total events for users. (<b>d</b>) San Francisco—Number of attended events for users. (<b>e</b>) London—Percentage of attended events in total events for users. (<b>f</b>) London—Number of attended events for users. (<b>g</b>) Sydney—Percentage of attended events in total events for users. (<b>h</b>) Sydney—Number of attended events for users.</p>
Full article ">Figure 11
<p>Numbers of loyal users and disloyal users varying in different <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>s</mi> </mrow> </semantics></math> for the four cities. (<b>a</b>) New York—Number of loyal users and disloyal users. (<b>b</b>) San Francisco—Number of loyal users and disloyal users. (<b>c</b>) London—Number of loyal users and disloyal users. (<b>d</b>) Sydney—Number of loyal users and disloyal users.</p>
Full article ">Figure 12
<p>Number of loyal users in groups varying in different <math display="inline"><semantics> <mrow> <mi>T</mi> <mi>s</mi> </mrow> </semantics></math>. Groups are obtained from <a href="#entropy-22-00119-t005" class="html-table">Table 5</a>.</p>
Full article ">
22 pages, 11340 KiB  
Article
Entropy-Based Image Fusion with Joint Sparse Representation and Rolling Guidance Filter
by Yudan Liu, Xiaomin Yang, Rongzhu Zhang, Marcelo Keese Albertini, Turgay Celik and Gwanggil Jeon
Entropy 2020, 22(1), 118; https://doi.org/10.3390/e22010118 - 18 Jan 2020
Cited by 22 | Viewed by 5104
Abstract
Image fusion is a very practical technology that can be applied in many fields, such as medicine, remote sensing and surveillance. An image fusion method using multi-scale decomposition and joint sparse representation is introduced in this paper. First, joint sparse representation is applied [...] Read more.
Image fusion is a very practical technology that can be applied in many fields, such as medicine, remote sensing and surveillance. An image fusion method using multi-scale decomposition and joint sparse representation is introduced in this paper. First, joint sparse representation is applied to decompose two source images into a common image and two innovation images. Second, two initial weight maps are generated by filtering the two source images separately. Final weight maps are obtained by joint bilateral filtering according to the initial weight maps. Then, the multi-scale decomposition of the innovation images is performed through the rolling guide filter. Finally, the final weight maps are used to generate the fused innovation image. The fused innovation image and the common image are combined to generate the ultimate fused image. The experimental results show that our method’s average metrics are: mutual information ( M I )—5.3377, feature mutual information ( F M I )—0.5600, normalized weighted edge preservation value ( Q A B / F )—0.6978 and nonlinear correlation information entropy ( N C I E )—0.8226. Our method can achieve better performance compared to the state-of-the-art methods in visual perception and objective quantification. Full article
(This article belongs to the Special Issue Entropy-Based Algorithms for Signal Processing)
Show Figures

Figure 1

Figure 1
<p>Rolling guidance filtering.</p>
Full article ">Figure 2
<p>The schematic diagram of our proposed fusion method.</p>
Full article ">Figure 3
<p>An example of JSR decomposition. (<b>a</b>,<b>b</b>) Source images; (<b>c</b>) The common image; (<b>d</b>) The innovation image of (<b>a</b>); (<b>e</b>) The innovation image of (<b>b</b>).</p>
Full article ">Figure 4
<p>Infrared-visible image sets. (<b>a</b>–<b>d</b>) Four sets of infrared–visible source images.</p>
Full article ">Figure 5
<p>Medical image sets. (<b>a</b>–<b>d</b>) Four sets of medical source images.</p>
Full article ">Figure 6
<p>Multi-focus image sets. (<b>a</b>–<b>d</b>) Four sets of multi-focus source images.</p>
Full article ">Figure 7
<p>Remote sensing image sets. (<b>a</b>–<b>d</b>) Four sets of remote sensing source images.</p>
Full article ">Figure 8
<p>Objective evaluation of different decomposition level <span class="html-italic">K</span>. (<b>a</b>–<b>d</b>) The values of <math display="inline"><semantics> <mrow> <mi>M</mi> <mi>I</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>F</mi> <mi>M</mi> <mi>I</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <msup> <mi>Q</mi> <mrow> <mi>A</mi> <mi>B</mi> <mo>/</mo> <mi>F</mi> </mrow> </msup> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>N</mi> <mi>C</mi> <mi>I</mi> <mi>E</mi> </mrow> </semantics></math>, respectively.</p>
Full article ">Figure 9
<p>Some fused images of JSR, RGF and our proposed method. (<b>a</b>,<b>b</b>) Source images; (<b>c</b>) The fused results of JSR; (<b>d</b>) The fused results of RGF; (<b>e</b>) The fused results of our proposed method.</p>
Full article ">Figure 10
<p>Examples of the fusion results of infrared-visible images. (<b>a</b>,<b>b</b>) Source images; (<b>c</b>–<b>p</b>) The fused results of ASR, CSR, CVT, DTCWT, GTF, H-MSD, CNN, LP, MSSR, MSVD, NSCT, WLS, FFIF, and our proposed method, respectively.</p>
Full article ">Figure 11
<p>Examples of the fusion results of medical images. (<b>a</b>,<b>b</b>) Source images; (<b>c</b>–<b>p</b>) The fused results of ASR, CSR, CVT, DTCWT, GTF, H-MSD, CNN, LP, MSSR, MSVD, NSCT, WLS, FFIF, and our proposed method, respectively.</p>
Full article ">Figure 12
<p>Examples of the fusion results of multi-focus images. (<b>a</b>,<b>b</b>) Source images; (<b>c</b>–<b>p</b>) The fused results of ASR, CSR, CVT, DTCWT, GTF, H-MSD, CNN, LP, MSSR, MSVD, NSCT, WLS, FFIF, and our proposed method, respectively.</p>
Full article ">Figure 13
<p>Examples of the fusion results of remote sensing images. (<b>a</b>,<b>b</b>) Source images; (<b>c</b>–<b>p</b>) The fused results of ASR, CSR, CVT, DTCWT, GTF, H-MSD, CNN, LP, MSSR, MSVD, NSCT, WLS, FFIF, and our proposed method, respectively.</p>
Full article ">
15 pages, 2398 KiB  
Article
Chemical Reaction Networks Possess Intrinsic, Temperature-Dependent Functionality
by Stephan O. Adler and Edda Klipp
Entropy 2020, 22(1), 117; https://doi.org/10.3390/e22010117 - 18 Jan 2020
Cited by 6 | Viewed by 4322
Abstract
Temperature influences the life of many organisms in various ways. A great number of organisms live under conditions where their ability to adapt to changes in temperature can be vital and largely determines their fitness. Understanding the mechanisms and principles underlying this ability [...] Read more.
Temperature influences the life of many organisms in various ways. A great number of organisms live under conditions where their ability to adapt to changes in temperature can be vital and largely determines their fitness. Understanding the mechanisms and principles underlying this ability to adapt can be of great advantage, for example, to improve growth conditions for crops and increase their yield. In times of imminent, increasing climate change, this becomes even more important in order to find strategies and help crops cope with these fundamental changes. There is intense research in the field of acclimation that comprises fluctuations of various environmental conditions, but most acclimation research focuses on regulatory effects and the observation of gene expression changes within the examined organism. As thermodynamic effects are a direct consequence of temperature changes, these should necessarily be considered in this field of research but are often neglected. Additionally, compensated effects might be missed even though they are equally important for the organism, since they do not cause observable changes, but rather counteract them. In this work, using a systems biology approach, we demonstrate that even simple network motifs can exhibit temperature-dependent functional features resulting from the interplay of network structure and the distribution of activation energies over the involved reactions. The demonstrated functional features are (i) the reversal of fluxes within a linear pathway, (ii) a thermo-selective branched pathway with different flux modes and (iii) the increased flux towards carbohydrates in a minimal Calvin cycle that was designed to demonstrate temperature compensation within reaction networks. Comparing a system’s response to either temperature changes or changes in enzyme activity we also dissect the influence of thermodynamic changes versus genetic regulation. By this, we expand the scope of thermodynamic modelling of biochemical processes by addressing further possibilities and effects, following established mathematical descriptions of biophysical properties. Full article
(This article belongs to the Special Issue Information Flow and Entropy Production in Biomolecular Networks)
Show Figures

Figure 1

Figure 1
<p>Energy diagram for reversible Michaelis-Menten kinetics. The figure sketches the different levels of Gibbs free energy <span class="html-italic">G</span><sub>1</sub> to <span class="html-italic">G</span><sub>5</sub> for a reversible enzymatic reaction with a substrate <span class="html-italic">S</span>, an enzyme <span class="html-italic">E</span> and a product <span class="html-italic">P</span>. The relation of these energy levels is a determining factor for the overall reaction rate and direction.</p>
Full article ">Figure 2
<p>Illustration of exponential temperature response curves. Panel (<b>A</b>) shows a typical exponential temperature response curve of a reaction flux, following the Arrhenius equation. Panel (<b>B</b>) illustrates a case where the forward flux exceeds the backward flux of a reversible reaction at temperature <span class="html-italic">T<sub>J</sub></span><sub>=0</sub>. The resulting net flux at this temperature is zero and changes its direction when passing it.</p>
Full article ">Figure 3
<p>Flux reversal in a linear pathway. The upper part of (<b>A</b>) shows the structure of the reversible linear pathway, comprising the two species A and B. The rate constants of the forward reactions are indicated by <span class="html-italic">k</span><sub>1</sub> to <span class="html-italic">k</span><sub>3</sub>, while <span class="html-italic">k</span><sub>−1</sub> to <span class="html-italic">k</span><sub>−3</sub> are the rate constants of the backward reactions. The lower part of (<b>A</b>) shows the corresponding net fluxes for a temperature range from 273 to 308 K. The three panels of (<b>B</b>) illustrate the system’s response to continuous cyclical temperature changes. The upper panel shows the time courses for three different temperatures oscillating by ± 5 K around 290 K (blue), 300 K (green) and 310 K (orange). The middle panel shows the resulting temporal changes in concentrations of species A and B. The flux changes are depicted in the lower panel. The colors of the concentration and flux curves match the colors of their corresponding temperature curves.</p>
Full article ">Figure 4
<p>Energy diagram for a single reaction within a linear, reversible pathway. Panel (<b>A</b>) shows the Gibbs free energies of activation for the forward (<span class="html-italic">ΔG<sub>+i</sub></span>) and backward (<span class="html-italic">ΔG<sub>−i</sub></span>) direction of reaction i of the linear, reversible pathway shown in panel (<b>B</b>). Here, <span class="html-italic">ΔG<sub>i</sub></span> is the difference in Gibbs free energy of the substrate S<sub>i−1</sub> and the product S<sub>i</sub>.</p>
Full article ">Figure 5
<p>Thermo-selective, branched pathway. The figure shows the structure of a reversible branched pathway with three species (upper part) and the net fluxes (considering forward and backward direction) for all involved reactions (lower part). The numbers in top left corner in each of the panels denote the corresponding reaction. Here we see three specifiable temperature regions with distinctive flux directions (colored areas). While the flux comes in from the upper branch and flows out through the lower branch within the blue temperature region, it can change to flowing outwards at both branches (red region) and eventually invert the initial flux direction with rising temperature (green region).</p>
Full article ">Figure 6
<p>Comparison of the effect of gene expression changes versus temperature changes. Considering a minimal branched network, we changed either the concentrations of the three enzymes individually or the temperature that affects all three reactions at the same time. The panels in the top row show the flux changes in response to changes in temperature (left) and enzyme concentration (right) separately. Notice the reversal of flux direction of reaction 1 (blue). The bottom row shows three heatmaps illustrating the combined effects of enzyme 2 and temperature <span class="html-italic">T</span> on the concentration of species <span class="html-italic">S</span> (left). The relative change of <span class="html-italic">S</span> compared to the case that only enzyme 2 is varied and <span class="html-italic">T</span> is kept at 290 K (middle), and the relative change of S compared to the case that temperature is varied and enzyme 2 is kept at concentration 1.</p>
Full article ">Figure 7
<p>Increased flux towards carbohydrates in a minimal Calvin cycle. The figure shows the structure of a Calvin cycle model designed by P. Ruoff et al. (2007) [<a href="#B10-entropy-22-00117" class="html-bibr">10</a>] which comprises only its most essential processes as lumped single reactions. The structure of this network has the potential to exhibit an increased flux towards carbohydrates at lower temperatures that is based on the distribution of activation energies and pre-exponential factors of the reactions within the network. An accumulation of carbohydrates at lower temperatures has been observed and reported frequently.</p>
Full article ">Figure 8
<p>Temperature response of <span class="html-italic">σ</span> for a linear reversible reaction system. The figure shows the temperature response of the entropy production density <span class="html-italic">σ</span> of all internal reactions and reaction 2 from the linear pathway illustrated in <a href="#entropy-22-00117-f003" class="html-fig">Figure 3</a>, separately. The entropy production vanishes at <span class="html-italic">T<sub>J=</sub></span><sub>0</sub>, the temperature at which the net flux is zero. It increases for both higher and lower temperatures.</p>
Full article ">
20 pages, 6819 KiB  
Article
Permutation Entropy and Statistical Complexity in Mild Cognitive Impairment and Alzheimer’s Disease: An Analysis Based on Frequency Bands
by Ignacio Echegoyen, David López-Sanz, Johann H. Martínez, Fernando Maestú and Javier M. Buldú
Entropy 2020, 22(1), 116; https://doi.org/10.3390/e22010116 - 18 Jan 2020
Cited by 21 | Viewed by 4580
Abstract
We present one of the first applications of Permutation Entropy (PE) and Statistical Complexity (SC) (measured as the product of PE and Jensen-Shanon Divergence) on Magnetoencephalography (MEG) recordings of 46 subjects suffering from Mild Cognitive Impairment (MCI), 17 individuals diagnosed with Alzheimer’s Disease [...] Read more.
We present one of the first applications of Permutation Entropy (PE) and Statistical Complexity (SC) (measured as the product of PE and Jensen-Shanon Divergence) on Magnetoencephalography (MEG) recordings of 46 subjects suffering from Mild Cognitive Impairment (MCI), 17 individuals diagnosed with Alzheimer’s Disease (AD) and 48 healthy controls. We studied the differences in PE and SC in broadband signals and their decomposition into frequency bands ( δ , θ , α and β ), considering two modalities: (i) raw time series obtained from the magnetometers and (ii) a reconstruction into cortical sources or regions of interest (ROIs). We conducted our analyses at three levels: (i) at the group level we compared SC in each frequency band and modality between groups; (ii) at the individual level we compared how the [PE, SC] plane differs in each modality; and (iii) at the local level we explored differences in scalp and cortical space. We recovered classical results that considered only broadband signals and found a nontrivial pattern of alterations in each frequency band, showing that SC does not necessarily decrease in AD or MCI. Full article
(This article belongs to the Special Issue Permutation Entropy: Theory and Applications)
Show Figures

Figure 1

Figure 1
<p>Example of a spectrogram from a signal in sensors (<b>a</b>) and regions of interest (ROIs) (<b>b</b>). Selected from a randomly chosen channel/ROI of a healthy subject.</p>
Full article ">Figure 2
<p>Entropy-complexity [PE,SC] planes for sensors (<b>a</b>) and ROIs (<b>b</b>), for the three groups (control group, CG; mild cognitive impairment, MCI; and Alzheimer’s disease, AD). Both modalities show the same qualitative results. There are clear differences between bands, especially in the case of <math display="inline"><semantics> <mi>α</mi> </semantics></math>, <math display="inline"><semantics> <mi>δ</mi> </semantics></math> and <math display="inline"><semantics> <mi>θ</mi> </semantics></math>. Each band seems homogeneous across groups and modalities.</p>
Full article ">Figure 3
<p>Average complexity (SC) for broadband, <math display="inline"><semantics> <mi>δ</mi> </semantics></math>, <math display="inline"><semantics> <mi>θ</mi> </semantics></math>, <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>β</mi> </semantics></math> bands, for sensors (<b>a</b>) and ROIs (<b>b</b>). One asterisk: <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>&lt;</mo> <mn>0.05</mn> </mrow> </semantics></math>, two asterisks: <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>&lt;</mo> <mn>0.01</mn> </mrow> </semantics></math>, three asterisks: <math display="inline"><semantics> <mrow> <mi>p</mi> <mo>&lt;</mo> <mn>0.001</mn> </mrow> </semantics></math>. We find differences between the three groups in broadband, in both modalities. Also, <math display="inline"><semantics> <mi>θ</mi> </semantics></math> and <math display="inline"><semantics> <mi>α</mi> </semantics></math> bands show clear differences between CG and AD and MCI and AD, a pattern only preserved in <math display="inline"><semantics> <mi>α</mi> </semantics></math> band in ROIs.</p>
Full article ">Figure 4
<p>Topographic map of complexity over the scalp for each band and group. Each dot represents a sensor. First row (<b>a</b>): values in complexity as the difference between MCI and Controls divided by controls (<math display="inline"><semantics> <mrow> <mo>(</mo> <mi>M</mi> <mi>C</mi> <mi>I</mi> <mo>−</mo> <mi>C</mi> <mi>G</mi> <mo>)</mo> <mo>/</mo> <mi>C</mi> <mi>G</mi> </mrow> </semantics></math>). Second row (<b>b</b>): difference between AD and Controls divided by controls (<math display="inline"><semantics> <mrow> <mo>(</mo> <mi>A</mi> <mi>D</mi> <mo>−</mo> <mi>C</mi> <mi>G</mi> <mo>)</mo> <mo>/</mo> <mi>C</mi> <mi>G</mi> </mrow> </semantics></math>). Third row (<b>c</b>): difference between AD and MCI divided by controls (<math display="inline"><semantics> <mrow> <mo>(</mo> <mi>A</mi> <mi>D</mi> <mo>−</mo> <mi>M</mi> <mi>C</mi> <mi>I</mi> <mo>)</mo> <mo>/</mo> <mi>C</mi> <mi>G</mi> </mrow> </semantics></math>).</p>
Full article ">Figure 5
<p>Distribution of complexity over the cortical surface for each band and group (estimated sources). First row (<b>a</b>): values in complexity as the difference between MCI and Controls divided by controls (<math display="inline"><semantics> <mrow> <mo>(</mo> <mi>M</mi> <mi>C</mi> <mi>I</mi> <mo>−</mo> <mi>C</mi> <mi>G</mi> <mo>)</mo> <mo>/</mo> <mi>C</mi> <mi>G</mi> </mrow> </semantics></math>). Second row (<b>b</b>): difference between AD and Controls divided by controls (<math display="inline"><semantics> <mrow> <mo>(</mo> <mi>A</mi> <mi>D</mi> <mo>−</mo> <mi>C</mi> <mi>G</mi> <mo>)</mo> <mo>/</mo> <mi>C</mi> <mi>G</mi> </mrow> </semantics></math>). Third row (<b>c</b>): difference between AD and MCI divided by controls (<math display="inline"><semantics> <mrow> <mo>(</mo> <mi>A</mi> <mi>D</mi> <mo>−</mo> <mi>M</mi> <mi>C</mi> <mi>I</mi> <mo>)</mo> <mo>/</mo> <mi>C</mi> <mi>G</mi> </mrow> </semantics></math>).</p>
Full article ">Figure 6
<p>Schematic representation of the pipeline followed to conduct this research. (<b>1</b>) Data acquisition from MEG recordings (sensors). (<b>2</b>) Data curation. (<b>3</b>) Estimated signals in Regions of Interest. (<b>4</b>) Frequency decomposition into four different frequency bands: <math display="inline"><semantics> <mi>δ</mi> </semantics></math>, <math display="inline"><semantics> <mi>θ</mi> </semantics></math>, <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>β</mi> </semantics></math>. (<b>5</b>) Ordinal patterns extraction and probability density estimation. (<b>6</b>) Entropy, disequilibrium and complexity calculation.</p>
Full article ">
17 pages, 1176 KiB  
Article
Association between Mean Heart Rate and Recurrence Quantification Analysis of Heart Rate Variability in End-Stage Renal Disease
by Martín Calderón-Juárez, Gertrudis Hortensia González-Gómez, Juan C. Echeverría, Héctor Pérez-Grovas and Claudia Lerma
Entropy 2020, 22(1), 114; https://doi.org/10.3390/e22010114 - 18 Jan 2020
Cited by 11 | Viewed by 4424
Abstract
Linear heart rate variability (HRV) indices are dependent on the mean heart rate, which has been demonstrated in different models (from sinoatrial cells to humans). The association between nonlinear HRV indices, including those provided by recurrence plot quantitative analysis (RQA), and the mean [...] Read more.
Linear heart rate variability (HRV) indices are dependent on the mean heart rate, which has been demonstrated in different models (from sinoatrial cells to humans). The association between nonlinear HRV indices, including those provided by recurrence plot quantitative analysis (RQA), and the mean heart rate (or the mean cardiac period, also called meanNN) has been scarcely studied. For this purpose, we analyzed RQA indices of five minute-long HRV time series obtained in the supine position and during active standing from 30 healthy subjects and 29 end-stage renal disease (ESRD) patients (before and after hemodialysis). In the supine position, ESRD patients showed shorter meanNN (i.e., faster heart rate) and decreased variability compared to healthy subjects. The healthy subjects responded to active standing by shortening the meanNN and decreasing HRV indices to reach similar values of ESRD patients. Bivariate correlations between all RQA indices and meanNN were significant in healthy subjects and ESRD after hemodialysis and for most RQA indices in ESRD patients before hemodialysis. Multiple linear regression analyses showed that RQA indices were also dependent on the position and the ESRD condition. Then, future studies should consider the association among RQA indices, meanNN, and these other factors for a correct interpretation of HRV. Full article
Show Figures

Figure 1

Figure 1
<p>Assessment of the autocorrelation function (<b>upper panels</b>) and false nearest neighbors method (<b>lower panels</b>) for each recording (thin lines) and averaged values per group (thick lines). Black lines correspond to healthy subjects, green lines correspond to ESRD patients before hemodialysis, and magenta lines correspond to ESRD patients after hemodialysis.</p>
Full article ">Figure 2
<p>Examples of RR time series and recurrence plots from a healthy subject (<b>upper panels</b>) and an end-stage renal disease (ESRD) patient (<b>lower panels</b>). The recording from the healthy subject shows a reduction in RR intervals and their variability, as a consequence of changing to an active standing position. The recurrence plot shows extended global diagonal lines formed by short diagonals covering a regularly spaced pattern of the plot surface. In contrast, during the standing position, there is an increase in white zones (i.e., there are less accessible dynamic states in the central area). In the patient before hemodialysis, there is a marked reduction in the RR intervals and their associated variability, even in the supine position. Its recurrence plot quantitative analysis (RQA) plot has less regularly extended diagonals with more concentration around the identity line and its parallel lines. During the active standing condition, the plot is almost empty; there are very fewer accessible dynamical states. In the same patient after hemodialysis, there is a slight improvement in the RR intervals during the supine condition compared to those before hemodialysis, which is also reflected in the RQA plot. Again, during the active standing challenge, there is a reduction in RR intervals and reachable dynamical states.</p>
Full article ">
19 pages, 4838 KiB  
Article
Energy and Exergy Evaluation of a Two-Stage Axial Vapour Compressor on the LNG Carrier
by Igor Poljak, Josip Orović, Vedran Mrzljak and Dean Bernečić
Entropy 2020, 22(1), 115; https://doi.org/10.3390/e22010115 - 17 Jan 2020
Cited by 7 | Viewed by 4971
Abstract
Data from a two-stage axial vapor cryogenic compressor on the dual-fuel diesel–electric (DFDE) liquefied natural gas (LNG) carrier were measured and analyzed to investigate compressor energy and exergy efficiency in real exploitation conditions. The running parameters of the two-stage compressor were collected while [...] Read more.
Data from a two-stage axial vapor cryogenic compressor on the dual-fuel diesel–electric (DFDE) liquefied natural gas (LNG) carrier were measured and analyzed to investigate compressor energy and exergy efficiency in real exploitation conditions. The running parameters of the two-stage compressor were collected while changing the main propeller shafts rpm. As the compressor supply of vaporized gas to the main engines increases, so does the load and rpm in propulsion electric motors, and vice versa. The results show that when the main engine load varied from 46 to 56 rpm at main propulsion shafts increased mass flow rate of vaporized LNG at a two-stage compressor has an influence on compressor performance. Compressor average energy efficiency is around 50%, while the exergy efficiency of the compressor is significantly lower in all measured ranges and on average is around 34%. The change in the ambient temperature from 0 to 50 °C also influences the compressor’s exergy efficiency. Higher exergy efficiency is achieved at lower ambient temperatures. As temperature increases, overall compressor exergy efficiency decreases by about 7% on average over the whole analyzed range. The proposed new concept of energy-saving and increasing the compressor efficiency based on pre-cooling of the compressor second stage is also analyzed. The temperature at the second stage was varied in the range from 0 to −50 °C, which results in power savings up to 26 kW for optimal running regimes. Full article
(This article belongs to the Special Issue Carnot Cycle and Heat Engine Fundamentals and Applications)
Show Figures

Figure 1

Figure 1
<p>Two-stage compressor outline with given flows.</p>
Full article ">Figure 2
<p>The <span class="html-italic">h</span>–<span class="html-italic">s</span> diagram of the compression process in the two-stage compressor.</p>
Full article ">Figure 3
<p>Two-stage compressor energy loss with main propulsion shafts speed variation.</p>
Full article ">Figure 4
<p>Two-stage compressor energy efficiency with main propulsion shafts speed variation.</p>
Full article ">Figure 5
<p>Two-stage compressor exergy destruction with main propulsion shafts speed variation.</p>
Full article ">Figure 6
<p>Two-stage compressor exergy efficiency with main propulsion shafts speed variation.</p>
Full article ">Figure 7
<p>Two-stage compressor energy and exergy efficiency with main propulsion shafts speed variation.</p>
Full article ">Figure 8
<p>First stage compressor exergy efficiency with surrounding temperature change.</p>
Full article ">Figure 9
<p>Compressor’s second stage exergy efficiency with surrounding temperature change.</p>
Full article ">Figure 10
<p>Two-stage compressor exergy efficiency with surrounding temperature change.</p>
Full article ">Figure 11
<p>Two-stage compressor cumulative power optimum saving with second stage inlet temperature variation.</p>
Full article ">
17 pages, 7722 KiB  
Article
Entropy-Based Effect Evaluation of Delineators in Tunnels on Drivers’ Gaze Behavior
by Xueyan Han, Yang Shao, Shaowei Yang and Peng Yu
Entropy 2020, 22(1), 113; https://doi.org/10.3390/e22010113 - 17 Jan 2020
Cited by 18 | Viewed by 3647
Abstract
Driving safety in tunnels has always been an issue of great concern. Establishing delineators to improve drivers’ instantaneous cognition of the surrounding environment in tunnels can effectively enhance driver safety. Through a simulation study, this paper explored how delineators affect drivers’ gaze behavior [...] Read more.
Driving safety in tunnels has always been an issue of great concern. Establishing delineators to improve drivers’ instantaneous cognition of the surrounding environment in tunnels can effectively enhance driver safety. Through a simulation study, this paper explored how delineators affect drivers’ gaze behavior (including fixation and scanpath) in tunnels. In addition to analyzing typical parameters, such as fixation position and fixation duration in areas of interest (AOIs), by modeling drivers’ switching process as Markov chains and calculating Shannon’s entropy of the fit Markov model, this paper quantified the complexity of individual switching patterns between AOIs under different delineator configurations and with different road alignments. A total of 25 subjects participated in this research. The results show that setting delineators in tunnels can attract drivers’ attention and make them focus on the pavement. When driving in tunnels equipped with delineators, especially tunnels with both wall delineators and pavement delineators, the participants exhibited a smaller transition entropy H t and stationary entropy H s , which can greatly reduce drivers’ visual fatigue. Compared with left curve and right curve, participants obtained higher H t and H s values in the straight section. Full article
(This article belongs to the Section Multidisciplinary Applications)
Show Figures

Figure 1

Figure 1
<p>Test route alignment design (Qinling Mountain No. 1, No. 2 and No. 3 tunnels of the G5 Expressway).</p>
Full article ">Figure 2
<p>The location of delineators in tunnel.</p>
Full article ">Figure 3
<p>Simulation view in the simulator. This picture is the simulation shown on the middle liquid crystal display when participants drove in scenario A.</p>
Full article ">Figure 4
<p>(<b>a</b>) Straight; (<b>b</b>) Left curve; (<b>c</b>) Right curve. Participants’ fixation positions in three different scenarios (In scenario A, both the wall delineators and pavement delineators were set. In scenario B, only the pavement delineators were set. In scenario C, there were no delineators.) combined with three tunnel alignments (straight, left curve and right curve).</p>
Full article ">Figure 5
<p>Participant’s visual area division.</p>
Full article ">Figure 6
<p>(<b>a</b>) straight line in scenario A; (<b>b</b>) straight line in scenario B; (<b>c</b>) straight line in scenario C; (<b>d</b>) left curve in scenario A; (<b>e</b>) left curve in scenario B; (<b>f</b>) left curve in scenario C; (<b>g</b>) right curve in scenario A; (<b>h</b>) right curve in scenario B; (<b>i</b>) right curve in scenario C. The percentage of fixation duration in three different scenarios (In scenario A, both the wall delineators and pavement delineators were set. In scenario B, only the pavement delineators were set. In scenario C, there was no delineator.) combined with three tunnel alignments (straight, left curve and right curve) for the 21 participants. There are 9 subgraphs, in each subgraph, there are 21 concentric circles. Each concentric circle represents one participant. The six different colors in each concentric circle represent the fixation duration percentage of participant for six different areas, the black, red, blue, green, purple, and yellow parts represent the proportion of gaze time spent in the pavement area (PA), the central area (CA), the top wall (TW), the left wall (LW), the right wall (RW) and the white space (WS), respectively.</p>
Full article ">Figure 7
<p>(<b>a</b>) <math display="inline"><semantics> <msub> <mi>H</mi> <mi>t</mi> </msub> </semantics></math> in straight line; (<b>b</b>) <math display="inline"><semantics> <msub> <mi>H</mi> <mi>s</mi> </msub> </semantics></math> in straight line; (<b>c</b>) <math display="inline"><semantics> <msub> <mi>H</mi> <mi>t</mi> </msub> </semantics></math> on left curve; (<b>d</b>) <math display="inline"><semantics> <msub> <mi>H</mi> <mi>s</mi> </msub> </semantics></math> on left curve; (<b>e</b>) <math display="inline"><semantics> <msub> <mi>H</mi> <mi>t</mi> </msub> </semantics></math> on right curve; (<b>f</b>) <math display="inline"><semantics> <msub> <mi>H</mi> <mi>s</mi> </msub> </semantics></math> on right curve. The 21 participants’ entropy values in three different scenarios (In scenario A, both the wall delineators and pavement delineators were set. In scenario B, only the pavement delineators were set. In scenario C, there were no delineators.) combined with three tunnel alignments (straight, left curve and right curve).</p>
Full article ">
28 pages, 6517 KiB  
Article
Fast and Efficient Image Encryption Algorithm Based on Modular Addition and SPD
by Khushbu Khalid Butt, Guohui Li, Sajid Khan and Sohaib Manzoor
Entropy 2020, 22(1), 112; https://doi.org/10.3390/e22010112 - 16 Jan 2020
Cited by 10 | Viewed by 5093
Abstract
Bit-level and pixel-level methods are two classifications for image encryption, which describe the smallest processing elements manipulated in diffusion and permutation respectively. Most pixel-level permutation methods merely alter the positions of pixels, resulting in similar histograms for the original and permuted images. Bit-level [...] Read more.
Bit-level and pixel-level methods are two classifications for image encryption, which describe the smallest processing elements manipulated in diffusion and permutation respectively. Most pixel-level permutation methods merely alter the positions of pixels, resulting in similar histograms for the original and permuted images. Bit-level permutation methods, however, have the ability to change the histogram of the image, but are usually not preferred due to their time-consuming nature, which is owed to bit-level computation, unlike that of other permutation techniques. In this paper, we introduce a new image encryption algorithm which uses binary bit-plane scrambling and an SPD diffusion technique for the bit-planes of a plain image, based on a card game trick. Integer values of the hexadecimal key SHA-512 are also used, along with the adaptive block-based modular addition of pixels to encrypt the images. To prove the first-rate encryption performance of our proposed algorithm, security analyses are provided in this paper. Simulations and other results confirmed the robustness of the proposed image encryption algorithm against many well-known attacks; in particular, brute-force attacks, known/chosen plain text attacks, occlusion attacks, differential attacks, and gray value difference attacks, among others. Full article
(This article belongs to the Section Multidisciplinary Applications)
Show Figures

Figure 1

Figure 1
<p>General block diagram.</p>
Full article ">Figure 2
<p>Encryption and decryption results. (<b>a</b>) Plain image of fruits; (<b>b</b>) cipher image of fruits; (<b>c</b>) decrypted image of fruits.</p>
Full article ">Figure 3
<p>All test images.</p>
Full article ">Figure 4
<p>Key sensitivity analysis. (<b>a</b>) Number of bit change rate (NBCR) between C1 and C2; (<b>b</b>) NBCR between D1 and D2.</p>
Full article ">Figure 5
<p>Key sensitivity test: (<b>a</b>) encrypted with K1; (<b>b</b>) encrypted with K2; (<b>c</b>) subtraction of (<b>a</b>,<b>b</b>); (<b>d</b>) histogram of (<b>a</b>); (<b>e</b>) histogram of (<b>b</b>); (<b>f</b>) histogram of (<b>c</b>); (<b>g</b>) decrypted with K1; (<b>h</b>) decrypted with K2; (<b>i</b>) subtraction of (<b>g</b>,<b>h</b>); (<b>j</b>) histogram of (<b>g</b>); (<b>k</b>) histogram of (<b>h</b>); (<b>l</b>) histogram of (<b>i</b>).</p>
Full article ">Figure 6
<p>Histograms of plain test images.</p>
Full article ">Figure 7
<p>Histograms of cipher images.</p>
Full article ">Figure 8
<p>Histograms of RGB channels of fruits. (<b>a</b>) Plain red channel; (<b>b</b>) plain green channel; (<b>c</b>) plain blue channel; (<b>d</b>) ciphered red channel; (<b>e</b>) ciphered green channel; (<b>f</b>) ciphered blue channel.</p>
Full article ">Figure 8 Cont.
<p>Histograms of RGB channels of fruits. (<b>a</b>) Plain red channel; (<b>b</b>) plain green channel; (<b>c</b>) plain blue channel; (<b>d</b>) ciphered red channel; (<b>e</b>) ciphered green channel; (<b>f</b>) ciphered blue channel.</p>
Full article ">Figure 9
<p>Pixel correlation of test image tree. (<b>a</b>) Plain red channel pixels, horizontal; (<b>b</b>) plain green channel pixels, vertical; (<b>c</b>) plain blue channel, diagonal; (<b>d</b>) ciphered red channel pixels, horizontal; (<b>e</b>) ciphered green channel pixels, vertical; (<b>f</b>) ciphered blue channel pixels, diagonal.</p>
Full article ">Figure 10
<p>Encryption result of all white and full black images. (<b>a</b>) All white image; (<b>b</b>) full black image; (<b>c</b>) cipher image of (<b>a</b>); (<b>d</b>) cipher image of (<b>b</b>); (<b>e</b>) histogram of (<b>c</b>); (<b>f</b>) histogram of (<b>d</b>); (<b>g</b>) cipher of one pixel change of (<b>a</b>); (<b>h</b>) subtraction of (<b>b</b>,<b>g</b>); (<b>i</b>) histogram of (<b>g</b>).</p>
Full article ">Figure 11
<p>Occlusion attack test. (<b>a</b>) Cropped image of Lena; (<b>b</b>) cropped image of a panda; (<b>c</b>) cropped image of a house; (<b>d</b>) cropped image of baboon; (<b>e</b>) cropped image of a panda; (<b>f</b>) cropped image of Lena; (<b>g</b>) retrieved image of (<b>a</b>); (<b>h</b>) retrieved image of (<b>b</b>); (<b>i</b>) retrieved image of (<b>c</b>); (<b>j</b>) retrieved image of (<b>d</b>); (<b>k</b>) retrieved image of (<b>e</b>); (<b>l</b>) retrieved image of (<b>f</b>).</p>
Full article ">Figure 12
<p>Encryption time in seconds.</p>
Full article ">
29 pages, 970 KiB  
Article
Quantifying Athermality and Quantum Induced Deviations from Classical Fluctuation Relations
by Zoë Holmes, Erick Hinds Mingo, Calvin Y.-R. Chen and Florian Mintert
Entropy 2020, 22(1), 111; https://doi.org/10.3390/e22010111 - 16 Jan 2020
Cited by 2 | Viewed by 4304
Abstract
In recent years, a quantum information theoretic framework has emerged for incorporating non-classical phenomena into fluctuation relations. Here, we elucidate this framework by exploring deviations from classical fluctuation relations resulting from the athermality of the initial thermal system and quantum coherence of the [...] Read more.
In recent years, a quantum information theoretic framework has emerged for incorporating non-classical phenomena into fluctuation relations. Here, we elucidate this framework by exploring deviations from classical fluctuation relations resulting from the athermality of the initial thermal system and quantum coherence of the system’s energy supply. In particular, we develop Crooks-like equalities for an oscillator system which is prepared either in photon added or photon subtracted thermal states and derive a Jarzynski-like equality for average work extraction. We use these equalities to discuss the extent to which adding or subtracting a photon increases the informational content of a state, thereby amplifying the suppression of free energy increasing process. We go on to derive a Crooks-like equality for an energy supply that is prepared in a pure binomial state, leading to a non-trivial contribution from energy and coherence on the resultant irreversibility. We show how the binomial state equality fits in relation to a previously derived coherent state equality and offers a richer feature-set. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Relation between prepared states and measurements. In the forwards (reverse) process, the state <math display="inline"><semantics> <mrow> <msubsup> <mi>ρ</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>ρ</mi> <mi>S</mi> <mi>i</mi> </msubsup> <mo>⊗</mo> <msubsup> <mi>ρ</mi> <mi>B</mi> <mi>i</mi> </msubsup> </mrow> </semantics></math><math display="inline"><semantics> <mfenced separators="" open="(" close=")"> <msubsup> <mi>ρ</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>f</mi> </msubsup> <mo>=</mo> <msubsup> <mi>ρ</mi> <mi>S</mi> <mi>f</mi> </msubsup> <mo>⊗</mo> <msubsup> <mi>ρ</mi> <mi>B</mi> <mi>f</mi> </msubsup> </mfenced> </semantics></math> is prepared, it evolves under <span class="html-italic">U</span> as indicated by the wiggly arrow, and then the measurement <math display="inline"><semantics> <mrow> <msubsup> <mi>X</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>f</mi> </msubsup> <mo>=</mo> <msubsup> <mi>X</mi> <mi>S</mi> <mi>f</mi> </msubsup> <mo>⊗</mo> <msubsup> <mi>X</mi> <mi>B</mi> <mi>f</mi> </msubsup> </mrow> </semantics></math><math display="inline"><semantics> <mfenced separators="" open="(" close=")"> <msubsup> <mi>X</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <msubsup> <mi>X</mi> <mi>S</mi> <mi>i</mi> </msubsup> <mo>⊗</mo> <msubsup> <mi>X</mi> <mi>B</mi> <mi>i</mi> </msubsup> </mfenced> </semantics></math> is performed. As indicated by the solid lines, the measurements <math display="inline"><semantics> <msubsup> <mi>X</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>i</mi> </msubsup> </semantics></math> and <math display="inline"><semantics> <msubsup> <mi>X</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>f</mi> </msubsup> </semantics></math> are related to the states <math display="inline"><semantics> <msubsup> <mi>ρ</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>i</mi> </msubsup> </semantics></math> and <math display="inline"><semantics> <msubsup> <mi>ρ</mi> <mrow> <mi>S</mi> <mi>B</mi> </mrow> <mi>f</mi> </msubsup> </semantics></math>, respectively, by the mapping <math display="inline"><semantics> <mi mathvariant="script">M</mi> </semantics></math>, defined in Equation (<a href="#FD7-entropy-22-00111" class="html-disp-formula">7</a>).</p>
Full article ">Figure 2
<p>Generalised Free Energies. The solid red and dark blue lines show the generalised free energy, <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msup> <mi>F</mi> <mo>+</mo> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msup> <mi>F</mi> <mo>−</mo> </msup> </mrow> </semantics></math>, of the oscillator system for the photon added and photon subtracted equalities, respectively. These are plotted as a function of <math display="inline"><semantics> <mrow> <mi>χ</mi> <mo>=</mo> <mi>β</mi> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math>, the ratio between the initial vacuum fluctuations, <math display="inline"><semantics> <mrow> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math>, and the thermal fluctuations, <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>B</mi> </msub> <mi>T</mi> </mrow> </semantics></math>, a measure which quantifies the temperature and thus effectively delineates the classical and quantum regimes. The grey dashed line is the usual change in energy <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>F</mi> </mrow> </semantics></math>. The dotted lines indicate the contribution of <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msub> <mi>E</mi> <mrow> <mi>vac</mi> </mrow> </msub> </mrow> </semantics></math> (purple) and <math display="inline"><semantics> <mrow> <mn>2</mn> <mo>Δ</mo> <mi>F</mi> </mrow> </semantics></math> (light blue) to <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msup> <mi>F</mi> <mo>+</mo> </msup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>Δ</mo> <msup> <mi>F</mi> <mo>+</mo> </msup> </mrow> </semantics></math>. In this plot, we suppose <math display="inline"><semantics> <mrow> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>f</mi> </msub> <mo>=</mo> <mn>1.5</mn> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>i</mi> </msub> </mrow> </semantics></math> and energies are given in units of <math display="inline"><semantics> <mrow> <msub> <mi>k</mi> <mi>B</mi> </msub> <mi>T</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Predicted ratio and <math display="inline"><semantics> <mi mathvariant="script">R</mi> </semantics></math> prefactor. The left figure plots the predicted ratio of the forwards and reverse transition probabilities, i.e., the right hand side of Equation (<a href="#FD25-entropy-22-00111" class="html-disp-formula">25</a>), for the photon added (subtracted) Crooks equality as a function of <math display="inline"><semantics> <mrow> <mi>χ</mi> <mo>=</mo> <mi>β</mi> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>i</mi> </msub> <mo>/</mo> <mn>2</mn> </mrow> </semantics></math>. The right figure plots <math display="inline"><semantics> <mi mathvariant="script">R</mi> </semantics></math> as a function of <math display="inline"><semantics> <mi>χ</mi> </semantics></math>. The red (blue) lines indicates the photon added (subtracted) case and the grey lines indicate the equivalent classical limit. That is, in the left plot the grey line is the right hand side of the classical Crooks equality, Equation (<a href="#FD1-entropy-22-00111" class="html-disp-formula">1</a>), and in the right plot the grey line is <math display="inline"><semantics> <mrow> <mi mathvariant="script">R</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>. The solid lines plot the case <math display="inline"><semantics> <mrow> <mi>W</mi> <mo>=</mo> <mn>2</mn> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>i</mi> </msub> </mrow> </semantics></math> and the dashed lines, <math display="inline"><semantics> <mrow> <mi>W</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>. Here, we suppose <math display="inline"><semantics> <mrow> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>f</mi> </msub> <mo>=</mo> <mn>5</mn> <mi>ℏ</mi> <msub> <mi>ω</mi> <mi>i</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Quantum distortions of fluctuation relations due to binomial battery states: The left and right plots correspond to <math display="inline"><semantics> <msub> <mi>q</mi> <mrow> <mi>align</mi> </mrow> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>q</mi> <mrow> <mi>size</mi> </mrow> </msub> </semantics></math>, respectively. The left plot is evaluated for a fixed value <math display="inline"><semantics> <mrow> <msub> <mi>p</mi> <mi>f</mi> </msub> <mo>=</mo> <mn>0.8</mn> </mrow> </semantics></math>. Both functions are plotted against the quantum-thermodynamic ratio <math display="inline"><semantics> <mrow> <mi>χ</mi> <mo>=</mo> <mfrac> <mrow> <mi>β</mi> <mi>ℏ</mi> <mi>ω</mi> </mrow> <mn>2</mn> </mfrac> </mrow> </semantics></math>. The plots show that the distortion due to quantum features can both enhance and suppress irreversibility in a process as compared to a “classical equivalent” solely involving energy exchanges. In both cases, we typically find suppressed irreversibility as quantum features dominate for large values of <math display="inline"><semantics> <mi>χ</mi> </semantics></math>. However, when thermodynamic and quantum energy scales are of similar magnitude, we observe unexpected behaviour.</p>
Full article ">Figure 5
<p>Linear optic implementation schematic. A photon added (or subtracted) thermal state is sent into one input arm of a linear optical set up and a coherent state the other. The linear optical set up, consisting of a series of linear optical elements, such as beamsplitters, phase-shifters and mirrors (the particular sequence sketched here is chosen arbitrarily), drives the photonic system and battery with an energy conserving and time reversal invariant operation. Finally, a coherent state measurement is performed on one output arm of the optical setup using a homodyne detection and the number of photons out put is measured in the other arm.</p>
Full article ">
22 pages, 4286 KiB  
Article
Gravity Wave Activity in the Stratosphere before the 2011 Tohoku Earthquake as the Mechanism of Lithosphere-atmosphere-ionosphere Coupling
by Shih-Sian Yang and Masashi Hayakawa
Entropy 2020, 22(1), 110; https://doi.org/10.3390/e22010110 - 16 Jan 2020
Cited by 28 | Viewed by 4192
Abstract
The precursory atmospheric gravity wave (AGW) activity in the stratosphere has been investigated in our previous paper by studying an inland Kumamoto earthquake (EQ). We are interested in whether the same phenomenon occurs or not before another major EQ, especially an oceanic EQ. [...] Read more.
The precursory atmospheric gravity wave (AGW) activity in the stratosphere has been investigated in our previous paper by studying an inland Kumamoto earthquake (EQ). We are interested in whether the same phenomenon occurs or not before another major EQ, especially an oceanic EQ. In this study, we have examined the stratospheric AGW activity before the oceanic 2011 Tohoku EQ (Mw 9.0), while using the temperature profiles that were retrieved from ERA5. The potential energy (EP) of AGW has enhanced from 3 to 7 March, 4–8 days before the EQ. The active region of the precursory AGW first appeared around the EQ epicenter, and then expanded omnidirectionally, but mainly toward the east, covering a wide area of 2500 km (in longitude) by 1500 km (in latitude). We also found the influence of the present AGW activity on some stratospheric parameters. The stratopause was heated and descended; the ozone concentration was also reduced and the zonal wind was reversed at the stratopause altitude before the EQ. These abnormalities of the stratospheric AGW and physical/chemical parameters are most significant on 5–6 March, which are found to be consistent in time and spatial distribution with the lower ionospheric perturbation, as detected by our VLF network observations. We have excluded the other probabilities by the processes of elimination and finally concluded that the abnormal phenomena observed in the present study are EQ precursors, although several potential sources can generate AGW activities and chemical variations in the stratosphere. The present paper shows that the abnormal stratospheric AGW activity has also been detected even before an oceanic EQ, and the AGW activity has obliquely propagated upward and further disturbed the lower ionosphere. This case study has provided further support to the AGW hypothesis of the lithosphere-atmosphere-ionosphere coupling process. Full article
Show Figures

Figure 1

Figure 1
<p>The time-altitude intensity of potential energy (E<sub>P</sub>). The magenta line indicates the occurrence time of the earthquake (EQ) (main shock). Thin/bold black lines circle the regions that the E<sub>P</sub> values exceed the high-E<sub>P</sub>/extremely-high-E<sub>P</sub> threshold (quadruple/octuple of the reference value), respectively.</p>
Full article ">Figure 2
<p>The E<sub>P</sub> maps from 1 to 8 March. Thin/bold black lines circle the areas that the E<sub>P</sub> values exceed the high-E<sub>P</sub>/extremely-high-E<sub>P</sub> threshold, respectively. The location of the EQ epicenter is marked by a magenta dot.</p>
Full article ">Figure 3
<p>The compared plots of (<b>a</b>) E<sub>P</sub> as shown in <a href="#entropy-22-00110-f001" class="html-fig">Figure 1</a>, (<b>b</b>) temperature, and (<b>c</b>) ozone mass mixing ratio. The stratopause height is also displayed by a bold curve in both panels (<b>b</b>,<b>c</b>). The EQ occurrence is indicated by the vertical magenta line.</p>
Full article ">Figure 4
<p>Same as <a href="#entropy-22-00110-f003" class="html-fig">Figure 3</a>, but observed by TIMED/SABER.</p>
Full article ">Figure 5
<p>The relative locations of two Japanese very low frequency/low frequency (VLF/LF) transmitters (JJY (Fukushima) and JJI (Miyazaki) indicated by blue diamonds) and VLF receiving stations (Moshiri (MSR), Chofu (CHF), Kasugai (KSG), and Kochi (KCH) shown with red stars). Three red lines are the propagation paths associated with the American transmitter, NLK, together with their corresponding wave sensitive areas in thin black lines defined by Fresnel zones (elliptic zones). The most important one is the NLK-CHF path. For comparison, the great-circle path for JJY-MSR is also plotted. The epicenters of the main shock and its foreshock were nearly the same, as indicated by a red star in the sea. Taken from [<a href="#B52-entropy-22-00110" class="html-bibr">52</a>].</p>
Full article ">Figure 6
<p>Temporal evolutions of the propagation characteristics for the NLK-CHF path. The top panel refers to the average nighttime amplitude (called trend), and the bottom, to the dispersion (in amplitude). All of these values are normalized by their corresponding standard deviation (σ). A clear anomaly as extraordinary depletion in trend is observed on 5 and 6 March. Taken from [<a href="#B52-entropy-22-00110" class="html-bibr">52</a>].</p>
Full article ">Figure 7
<p>Illustration of the possible suggested region of the precursory seismo-ionospheric perturbation for the Tohoku EQ as inferred from a combination of propagation characteristics of all propagation paths we have studied. However, the north and east extension is quite uncertain due to the limitation of our VLF network. Taken from [<a href="#B52-entropy-22-00110" class="html-bibr">52</a>].</p>
Full article ">Figure 8
<p>The convective available potential energy (CAPE) at the EQ epicenter.</p>
Full article ">Figure 9
<p>(<b>a</b>) The topographical elevation along 39.85° N, which is the latitude of Mount Iwate. The vertical magenta line further indicates the longitude (141° E) of Mount Iwate. (<b>b</b>) The longitude-altitude intensity of E<sub>P</sub> along the 39.85°N latitude. The vertical magenta line again indicates the longitude of Mount Iwate. (<b>c</b>) The E<sub>P</sub> map at 17.6 km altitude. Magenta triangles sketch the position of Ou Mountains, and the solid one is Mount Iwate. The magenta dot marks the EQ epicenter. The interval of contour lines is 1 J/kg in both (<b>b</b>,<b>c</b>).</p>
Full article ">
15 pages, 3180 KiB  
Article
Visual Analysis on Information Theory and Science of Complexity Approaches in Healthcare Research
by Xiaoyu Wang, Wang Zhao, Yongzhong Wang, Qin Zhao, Xuejie Yang, Kaixiang Su and Dongxiao Gu
Entropy 2020, 22(1), 109; https://doi.org/10.3390/e22010109 - 16 Jan 2020
Cited by 6 | Viewed by 4901
Abstract
In order to explore the knowledge base, research hotspot, development status, and future research direction of healthcare research based on information theory and complex science, a total of 3031 literature data samples from the core collection of Web of Science from 2003 to [...] Read more.
In order to explore the knowledge base, research hotspot, development status, and future research direction of healthcare research based on information theory and complex science, a total of 3031 literature data samples from the core collection of Web of Science from 2003 to 2019 were selected for bibliometric analysis. HistCite, CiteSpace, Excel, and other analytical tools were used to deeply analyze and visualize the temporal distribution, spatial distribution, knowledge evolution, literature co-citation, and research hotspots of this field. This paper reveals the current development of healthcare research field based on information theory and science of complexity, analyzes and discusses the research hotspots and future development that trends in this field, and provides important knowledge support for researchers in this field for further relevant research. Full article
Show Figures

Figure 1

Figure 1
<p>Annual number of published articles.</p>
Full article ">Figure 2
<p>Annual number of authors input.</p>
Full article ">Figure 3
<p>Annual author input–output ratio.</p>
Full article ">Figure 4
<p>Countries’ collaboration network.</p>
Full article ">Figure 5
<p>Institution collaboration network.</p>
Full article ">Figure 6
<p>Author collaboration network.</p>
Full article ">Figure 7
<p>Articles in the co-citation network.</p>
Full article ">Figure 8
<p>Keyword co-occurrence network.</p>
Full article ">
32 pages, 1107 KiB  
Article
Generalizing Information to the Evolution of Rational Belief
by Jed A. Duersch and Thomas A. Catanach
Entropy 2020, 22(1), 108; https://doi.org/10.3390/e22010108 - 16 Jan 2020
Cited by 7 | Viewed by 3441
Abstract
Information theory provides a mathematical foundation to measure uncertainty in belief. Belief is represented by a probability distribution that captures our understanding of an outcome’s plausibility. Information measures based on Shannon’s concept of entropy include realization information, Kullback–Leibler divergence, Lindley’s information in experiment, [...] Read more.
Information theory provides a mathematical foundation to measure uncertainty in belief. Belief is represented by a probability distribution that captures our understanding of an outcome’s plausibility. Information measures based on Shannon’s concept of entropy include realization information, Kullback–Leibler divergence, Lindley’s information in experiment, cross entropy, and mutual information. We derive a general theory of information from first principles that accounts for evolving belief and recovers all of these measures. Rather than simply gauging uncertainty, information is understood in this theory to measure change in belief. We may then regard entropy as the information we expect to gain upon realization of a discrete latent random variable. This theory of information is compatible with the Bayesian paradigm in which rational belief is updated as evidence becomes available. Furthermore, this theory admits novel measures of information with well-defined properties, which we explored in both analysis and experiment. This view of information illuminates the study of machine learning by allowing us to quantify information captured by a predictive model and distinguish it from residual information contained in training data. We gain related insights regarding feature selection, anomaly detection, and novel Bayesian approaches. Full article
Show Figures

Figure 1

Figure 1
<p>Example of the evolution of plausible labels for an image. Without evidence, the probability distribution <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="bold">q</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </semantics></math> assigns equal plausibility to all outcomes. A machine learning model processes the image and produces predictions <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="bold">q</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>y</mi> <mo>|</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </semantics></math>. The incorrect label <b>3</b> is represented by <math display="inline"><semantics> <mrow> <mi mathvariant="bold">r</mi> <mo>(</mo> <mi>y</mi> <mo>|</mo> <mover accent="true"> <mi>y</mi> <mo>ˇ</mo> </mover> <mo>)</mo> </mrow> </semantics></math>. After observing the image, shown on the right, the label is corrected to <b>0</b> in <math display="inline"><semantics> <mrow> <mi mathvariant="bold">r</mi> <mo>(</mo> <mi>y</mi> <mo>|</mo> <mover accent="true"> <mi>y</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>Information measurements before and after label correction. Neither construction of prediction information allows the computation to account for claimed labels. Residual information, however, decreases in the KL construction when the label is corrected. The Lindley forms are totally unaffected by relabeling.</p>
Full article ">Figure 3
<p>Information measurements using our proposed framework. Total information is a conserved quantity, and when our belief changes, so do the information measurements. Negative prediction information forewarns either potential mislabeling or a poor prediction.</p>
Full article ">Figure 4
<p>Prior distribution of <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> and prior predictive distribution of individual <span class="html-italic">y</span> samples. The domain of plausible <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> values is large before any observations are made.</p>
Full article ">Figure 5
<p>Typical inference of <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> from observation of 10 samples of <span class="html-italic">y</span> (left) followed by 10, 20, and 40 additional samples, respectively. Both prior belief and first inference deciles of <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> are shown in gray. As observations accumulate, the domain of plausible <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> values tightens.</p>
Full article ">Figure 6
<p>Histogram of the earliest inference information in the observation sequence. The vertical line at 5.36 bits is mutual information. Information is positive after first inference, but may drop with additional observations. The limiting view of infinite samples (realization) is shown on the right.</p>
Full article ">Figure 7
<p>Minimum first inference information out of one million independent experiments. This particularly rare case shows how first samples can mislead inference, which is later corrected by additional observations. The fourth inference (right) bears remarkably little overlap with the first.</p>
Full article ">Figure 8
<p>Maximum first inference information out of one million independent experiments. The true value of <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> has taken an extremely rare value. As evidence accumulates, plausible ranges of <math display="inline"><semantics> <mi mathvariant="bold-italic">θ</mi> </semantics></math> confirm the first inference.</p>
Full article ">Figure 9
<p>Inconsistent inference. The first 10 samples are drawn from a different ground truth than subsequent samples, but inference proceeds as usual. As additional data become available, first inference information becomes markedly negative.</p>
Full article ">Figure 10
<p>Histogram of information outcomes for mismatched and original labels. Correct label information is highly concentrated at 3.2 bits, which is <math display="inline"><semantics> <mrow> <mn>95.9</mn> <mo>%</mo> </mrow> </semantics></math> of the total information contained in labels. Mislabeled cases have mean information at –18 bits, and information is negative for <math display="inline"><semantics> <mrow> <mn>99.15</mn> <mo>%</mo> </mrow> </semantics></math> of mislabeled cases.</p>
Full article ">Figure 11
<p>Original MNIST labels. The top row shows lowest predictive label information among original labels. Notably, the two leading images appear to be genuinely mislabeled in the original dataset. Subsequent predictions are poor. The bottom row shows the highest information among original labels. Labels and predictions are consistent in these cases.</p>
Full article ">Figure 12
<p>Mislabeled digits. The top row shows the lowest predictive label information among mislabeled cases. In each case, the claimed label is implausible and the prediction is correct. The bottom row shows the highest prediction information among mislabeled cases. Although claimed labels are incorrect, most images share identifiable features with the claim.</p>
Full article ">
12 pages, 325 KiB  
Article
On the Composability of Statistically Secure Random Oblivious Transfer
by Rafael Dowsley, Jörn Müller-Quade and Anderson C. A. Nascimento
Entropy 2020, 22(1), 107; https://doi.org/10.3390/e22010107 - 16 Jan 2020
Cited by 2 | Viewed by 2891
Abstract
We show that random oblivious transfer protocols that are statistically secure according to a definition based on a list of information-theoretical properties are also statistically universally composable. That is, they are simulatable secure with an unlimited adversary, an unlimited simulator, and an unlimited [...] Read more.
We show that random oblivious transfer protocols that are statistically secure according to a definition based on a list of information-theoretical properties are also statistically universally composable. That is, they are simulatable secure with an unlimited adversary, an unlimited simulator, and an unlimited environment machine. Our result implies that several previous oblivious transfer protocols in the literature that were proven secure under weaker, non-composable definitions of security can actually be used in arbitrary statistically secure applications without lowering the security. Full article
Show Figures

Figure 1

Figure 1
<p>The one-out-of-two bit random oblivious transfer functionality.</p>
Full article ">Figure 2
<p>The functionality, given valid inputs, samples outputs according to the conditional probability distribution and delivers the outputs to <math display="inline"><semantics> <mi mathvariant="sans-serif">Alice</mi> </semantics></math> and <math display="inline"><semantics> <mi mathvariant="sans-serif">Bob</mi> </semantics></math>.</p>
Full article ">
15 pages, 4537 KiB  
Article
Spatio-Temporal Evolution Analysis of Drought Based on Cloud Transformation Algorithm over Northern Anhui Province
by Xia Bai, Yimin Wang, Juliang Jin, Shaowei Ning, Yanfang Wang and Chengguo Wu
Entropy 2020, 22(1), 106; https://doi.org/10.3390/e22010106 - 16 Jan 2020
Cited by 5 | Viewed by 2799
Abstract
Drought is one of the most typical and serious natural disasters, which occurs frequently in most of mainland China, and it is crucial to explore the evolution characteristics of drought for developing effective schemes and strategies of drought disaster risk management. Based on [...] Read more.
Drought is one of the most typical and serious natural disasters, which occurs frequently in most of mainland China, and it is crucial to explore the evolution characteristics of drought for developing effective schemes and strategies of drought disaster risk management. Based on the application of Cloud theory in the drought evolution research field, the cloud transformation algorithm, and the conception zooming coupling model was proposed to re-fit the distribution pattern of SPI instead of the Pearson-III distribution. Then the spatio-temporal evolution features of drought were further summarized utilizing the cloud characteristics, average, entropy, and hyper-entropy. Lastly, the application results in Northern Anhui province revealed that the drought condition was the most serious during the period from 1957 to 1970 with the SPI12 index in 49 months being less than −0.5 and 12 months with an extreme drought level. The overall drought intensity varied with the highest certainty level but lowest stability level in winter, but this was opposite in the summer. Moreover, drought hazard would be more significantly intensified along the elevation of latitude in Northern Anhui province. The overall drought hazard in Suzhou and Huaibei were the most serious, which is followed by Bozhou, Bengbu, and Fuyang. Drought intensity in Huainan was the lightest. The exploration results of drought evolution analysis were reasonable and reliable, which would supply an effective decision-making basis for establishing drought risk management strategies. Full article
(This article belongs to the Special Issue Spatial Information Theory)
Show Figures

Figure 1

Figure 1
<p>Geographical location of the Northern Anhui province in China.</p>
Full article ">Figure 2
<p>Distribution of a qualitative concept cloud <span class="html-italic">C</span>(1.5, 0.5, 0.1).</p>
Full article ">Figure 3
<p>Framework of drought evolution characteristic analysis of this study.</p>
Full article ">Figure 4
<p>Variation of SPI series for different time scales in Northern Anhui province, 1957-2010.</p>
Full article ">Figure 5
<p>Cloud distribution of SPI3 in Huaibei, (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 6
<p>Cloud distribution of SPI3 in Bozhou, (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 7
<p>Cloud distribution of SPI3 in Suzhou, (<b>a</b>) spring, (<b>b</b>) summer; (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 8
<p>Cloud distribution of SPI3 in Bengbu, (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 9
<p>Cloud distribution of SPI3 in Fuyang, (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 10
<p>Cloud distribution of SPI3 in Huainan, (<b>a</b>) spring, (<b>b</b>) summer, (<b>c</b>) autumn, and (<b>d</b>) winter.</p>
Full article ">Figure 11
<p>Determination process of cloud distribution of drought by CTA in Huaibei city, (<b>a</b>) probability density histogram, (<b>b</b>) initial conceptual cloud distribution, and (<b>c</b>) final cloud distribution of drought.</p>
Full article ">Figure 12
<p>Conceptual cloud distribution of drought in Northern Anhui province, (<b>a</b>) Bozhou, (<b>b</b>) Suzhou, (<b>c</b>) Bengbu, (<b>d</b>) Fuyang, and (<b>e</b>) Huainan.</p>
Full article ">
26 pages, 2607 KiB  
Article
Statistical Complexity Analysis of Turing Machine tapes with Fixed Algorithmic Complexity Using the Best-Order Markov Model
by Jorge M. Silva, Eduardo Pinho, Sérgio Matos and Diogo Pratas
Entropy 2020, 22(1), 105; https://doi.org/10.3390/e22010105 - 16 Jan 2020
Cited by 3 | Viewed by 4537
Abstract
Sources that generate symbolic sequences with algorithmic nature may differ in statistical complexity because they create structures that follow algorithmic schemes, rather than generating symbols from a probabilistic function assuming independence. In the case of Turing machines, this means that machines with the [...] Read more.
Sources that generate symbolic sequences with algorithmic nature may differ in statistical complexity because they create structures that follow algorithmic schemes, rather than generating symbols from a probabilistic function assuming independence. In the case of Turing machines, this means that machines with the same algorithmic complexity can create tapes with different statistical complexity. In this paper, we use a compression-based approach to measure global and local statistical complexity of specific Turing machine tapes with the same number of states and alphabet. Both measures are estimated using the best-order Markov model. For the global measure, we use the Normalized Compression (NC), while, for the local measures, we define and use normal and dynamic complexity profiles to quantify and localize lower and higher regions of statistical complexity. We assessed the validity of our methodology on synthetic and real genomic data showing that it is tolerant to increasing rates of editions and block permutations. Regarding the analysis of the tapes, we localize patterns of higher statistical complexity in two regions, for a different number of machine states. We show that these patterns are generated by a decrease of the tape’s amplitude, given the setting of small rule cycles. Additionally, we performed a comparison with a measure that uses both algorithmic and statistical approaches (BDM) for analysis of the tapes. Naturally, BDM is efficient given the algorithmic nature of the tapes. However, for a higher number of states, BDM is progressively approximated by our methodology. Finally, we provide a simple algorithm to increase the statistical complexity of a Turing machine tape while retaining the same algorithmic complexity. We supply a publicly available implementation of the algorithm in C++ language under the GPLv3 license. All results can be reproduced in full with scripts provided at the repository. Full article
(This article belongs to the Special Issue Shannon Information and Kolmogorov Complexity)
Show Figures

Figure 1

Figure 1
<p>Heat map of Normalized Compression with an increase in permutation and edition rate. Generated string starting with 500 zeros followed by 500 ones (<b>top</b>); NC_007044.1 <span class="html-italic">Microplitis demolitor bracovirus segment O</span>, complete genome (<b>bottom-left</b>); and MH201455.1 <span class="html-italic">Human parvovirus B19 isolate BX1</span>, complete genome (<b>bottom-right</b>).</p>
Full article ">Figure 2
<p>Plot of all TMs in <a href="#entropy-22-00105-t002" class="html-table">Table 2</a>. NC value is in blue and the tape’s normalized amplitude size is in yellow. The x-axes of the plots represent the index of the Turing machine computed according to Algorithm A1. The <b>blue background</b> is the plot that corresponds to the group of TMs with <math display="inline"><semantics> <mrow> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>; all other plots have <math display="inline"><semantics> <mrow> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>The average value for the amplitude of TM’s tape (<b>top-left</b>); average required bits to perform compression of the tape (<b>top-middle</b>); and average NC value (<b>top-right</b>), inside and outside the regions marked with circles in <a href="#entropy-22-00105-f002" class="html-fig">Figure 2</a>. The average bits required (<b>bottom-left</b>); and the average NC value obtained for the rules used by the TM (<b>bottom-right</b>), inside and outside the regions marked with circles in <a href="#entropy-22-00105-f002" class="html-fig">Figure 2</a></p>
Full article ">Figure 4
<p>Regional capture of average rule complexity profiles obtained from pseudo-randomly selected TMs with <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>∈</mo> <mo>{</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>5</mn> <mo>,</mo> <mn>6</mn> <mo>}</mo> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> up to 1000 iterations.</p>
Full article ">Figure 5
<p>Normal complexity profiles (<b>left</b>); and dynamic complexity profiles (<b>right</b>) obtained for some of the filtered TMs. Each TM has a different cardinality of states or alphabet.</p>
Full article ">Figure 6
<p>Comparison between the NC and BDM for 10,000 TM with <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> that ran over 50,000 iterations: (<b>Left</b>) BDM scaled by a factor of <math display="inline"><semantics> <msup> <mn>10</mn> <mn>2</mn> </msup> </semantics></math>; and (<b>Right</b>) same example but with non-scaled BDM.</p>
Full article ">Figure 7
<p>Comparison of: Method I (<b>left</b>); and Method II (<b>right</b>). (<b>Top</b>) Plots show the amplitude of the tapes, bits required to represent the sequence, and the NC obtained for 200 TMs after a low-pass filter was applied. (<b>Bottom</b>) Plots show the average tape amplitude (<b>bottom-left</b>); average bits required (<b>bottom-middle</b>); and NC (<b>bottom-right</b>). Green and red colors represent TMs before and after the method was applied, respectively. For Method I, the average corresponds to 200 instances and for Method II to 2000.</p>
Full article ">Figure 8
<p>First 59 characters of TMs’ tapes before and after the Method II was applied.</p>
Full article ">Figure 9
<p>Average final amplitude of the tape (<b>top-left</b>); variation of the bits required to represent the string (<b>top-right</b>); and variation of the NC (<b>bottom</b>), with the increase in number of rule iterations and tape iterations.</p>
Full article ">Figure A1
<p>Average rule complexity profiles obtained from pseudo-randomly selected TMs with <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>∈</mo> <mo>{</mo> <mn>2</mn> <mo>,</mo> <mo>…</mo> <mo>,</mo> <mn>6</mn> <mo>}</mo> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mo>{</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> <mo>}</mo> </mrow> </semantics></math> up to 1000 iterations.</p>
Full article ">Figure A2
<p>Comparison between the NC and BDM for 10,000 TM that have run over 50,000 iterations. (<b>top-left</b>) TMs with <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>=</mo> <mn>6</mn> <mo>,</mo> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>; (<b>top-right</b>) TMs with <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>=</mo> <mn>8</mn> <mo>,</mo> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>; (<b>bottom-left</b>) TMs with <math display="inline"><semantics> <mrow> <mo>#</mo> <mi>Q</mi> <mo>=</mo> <mn>10</mn> <mo>,</mo> <mo>#</mo> <mo>Σ</mo> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>; and (<b>bottom-right</b>) an example with non-scaled BDM results.</p>
Full article ">
12 pages, 474 KiB  
Article
Complexity of Cardiotocographic Signals as A Predictor of Labor
by João Monteiro-Santos, Teresa Henriques, Inês Nunes, Célia Amorim-Costa, João Bernardes and Cristina Costa-Santos
Entropy 2020, 22(1), 104; https://doi.org/10.3390/e22010104 - 16 Jan 2020
Cited by 5 | Viewed by 3370
Abstract
Prediction of labor is of extreme importance in obstetric care to allow for preventive measures, assuring that both baby and mother have the best possible care. In this work, the authors studied how important nonlinear parameters (entropy and compression) can be as labor [...] Read more.
Prediction of labor is of extreme importance in obstetric care to allow for preventive measures, assuring that both baby and mother have the best possible care. In this work, the authors studied how important nonlinear parameters (entropy and compression) can be as labor predictors. Linear features retrieved from the SisPorto system for cardiotocogram analysis and nonlinear measures were used to predict labor in a dataset of 1072 antepartum tracings, at between 30 and 35 weeks of gestation. Two groups were defined: Group A—fetuses whose traces date was less than one or two weeks before labor, and Group B—fetuses whose traces date was at least one or two weeks before labor. Results suggest that, compared with linear features such as decelerations and variability indices, compression improves labor prediction both within one (C-Statistics of 0.728) and two weeks (C-Statistics of 0.704). Moreover, the correlation between compression and long-term variability was significantly different in groups A and B, denoting that compression and heart rate variability look at different information associated with whether the fetus is closer to or further from labor onset. Nonlinear measures, compression in particular, may be useful in improving labor prediction as a complement to other fetal heart rate features. Full article
Show Figures

Figure 1

Figure 1
<p>Example of a fetal heart rate (FHR) time series.</p>
Full article ">
14 pages, 3812 KiB  
Article
Determining the Bulk Parameters of Plasma Electrons from Pitch-Angle Distribution Measurements
by Georgios Nicolaou, Robert Wicks, George Livadiotis, Daniel Verscharen, Christopher Owen and Dhiren Kataria
Entropy 2020, 22(1), 103; https://doi.org/10.3390/e22010103 - 16 Jan 2020
Cited by 15 | Viewed by 4878
Abstract
Electrostatic analysers measure the flux of plasma particles in velocity space and determine their velocity distribution function. There are occasions when science objectives require high time-resolution measurements, and the instrument operates in short measurement cycles, sampling only a portion of the velocity distribution [...] Read more.
Electrostatic analysers measure the flux of plasma particles in velocity space and determine their velocity distribution function. There are occasions when science objectives require high time-resolution measurements, and the instrument operates in short measurement cycles, sampling only a portion of the velocity distribution function. One such high-resolution measurement strategy consists of sampling the two-dimensional pitch-angle distributions of the plasma particles, which describes the velocities of the particles with respect to the local magnetic field direction. Here, we investigate the accuracy of plasma bulk parameters from such high-resolution measurements. We simulate electron observations from the Solar Wind Analyser’s (SWA) Electron Analyser System (EAS) on board Solar Orbiter. We show that fitting analysis of the synthetic datasets determines the plasma temperature and kappa index of the distribution within 10% of their actual values, even at large heliocentric distances where the expected solar wind flux is very low. Interestingly, we show that although measurement points with zero counts are not statistically significant, they provide information about the particle distribution function which becomes important when the particle flux is low. We also examine the convergence of the fitting algorithm for expected plasma conditions and discuss the sources of statistical and systematic uncertainties. Full article
(This article belongs to the Special Issue Theoretical Aspects of Kappa Distributions)
Show Figures

Figure 1

Figure 1
<p>Schematic of a Solar Wind Analyser’s Electron Analyser System (SWA-EAS) top-hat analyser head and its angular field of view. (<b>Left</b>) The elevation angle is defined as the complement of the angle between the particle velocity vector and the <span class="html-italic">z</span>-axis, perpendicular to the top-hat plane. The elevation angle of the electrons is resolved in 16 electrostatic uniform steps. (<b>Right</b>) The azimuth angle is the angle within the projection of the velocity vector on the top-hat plane and the <span class="html-italic">x</span>-axis. Both SWA-EAS analyser heads resolve the azimuth direction on MCP detectors using 32 sectors.</p>
Full article ">Figure 2
<p>Modelled counts as a function energy and azimuth direction on the analyser’s head frame for (<b>left</b>) plasma density <span class="html-italic">n</span> = 5 cm<sup>−3</sup> and (<b>right</b>) <span class="html-italic">n</span> = 50 cm<sup>−3</sup>. For both examples, the magnetic field vector (magenta) is in the top-hat plane (<span class="html-italic">Θ</span> = <span class="html-italic">θ</span><sub>B</sub> = 0°) in azimuth direction <span class="html-italic">Φ</span> = 45°. The bulk flow of the electrons <span class="html-italic">u</span><sub>0</sub> = 500 kms<sup>−1</sup> along the <span class="html-italic">x</span>-axis (<span class="html-italic">Θ</span> = <span class="html-italic">Φ</span> = 0°). The parallel temperature <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math> = 10 eV, the perpendicular temperature <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math> = 20 eV, and the kappa index <span class="html-italic">κ</span> = 3.</p>
Full article ">Figure 3
<p>(<b>Left</b>) Modelled counts as a function of energy and azimuth direction (instrument frame), using <span class="html-italic">n</span> = 20 cm<sup>−3</sup>, <span class="html-italic">u</span><sub>0</sub> = 500 kms<sup>−1</sup> towards the <span class="html-italic">x</span>-axis (<span class="html-italic">Θ</span> = <span class="html-italic">Φ</span> = 0°), <span class="html-italic">κ</span> = 3, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math> = 10 eV, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math> = 20 eV, and a magnetic-field direction (magenta) in the top-hat plane (<span class="html-italic">Θ</span> = 0° and <span class="html-italic">Φ</span> = 45°). (<b>Right</b>) Result of our fit to the modelled observations. The model finds the optimal combination of <span class="html-italic">n</span>, <span class="html-italic">κ</span>, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math> that minimizes the <span class="html-italic">χ</span><sup>2</sup> value (see text for more).</p>
Full article ">Figure 4
<p>Histograms of (<b>top left</b>) density <span class="html-italic">n</span><sub>out</sub>, (<b>top right</b>) kappa index <span class="html-italic">κ</span><sub>out</sub>, (<b>bottom left</b>) parallel temperature <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mrow> <mo>∥</mo> <mo>,</mo> <mi>out</mi> </mrow> </msub> </mrow> </semantics></math> and (<b>bottom right</b>) perpendicular temperature <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mrow> <mo>⊥</mo> <mo>,</mo> <mi>out</mi> </mrow> </msub> </mrow> </semantics></math>, as determined from the analysis of 200 measurement samples of plasma with <span class="html-italic">n</span> = 7 cm<sup>−3</sup>, <span class="html-italic">u</span><sub>0</sub> = 500 kms<sup>−1</sup> pointing along the <span class="html-italic">x</span>-axis, <span class="html-italic">κ</span> = 3, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math> = 10 eV and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math> = 20 eV. The blue histograms correspond to values derived by a fit that includes points with <span class="html-italic">C<sub>i</sub></span> = 0, while the red histograms represent values derived by a fit that excludes points with <span class="html-italic">C<sub>i</sub></span> = 0.</p>
Full article ">Figure 5
<p>(<b>From top to bottom</b>) The derived electron density over input density, kappa index, parallel and perpendicular temperature as functions of the input plasma density. The red points represent the mean values (over 200 samples) of the parameters derived by fitting only the measurements with <span class="html-italic">C<sub>i</sub></span> ≥ 1. The blue points represent the mean values of the parameters derived by fitting to all measurements including those with <span class="html-italic">C<sub>i</sub></span> = 0. The shadowed regions represent the standard deviations of the derived parameters.</p>
Full article ">Figure 6
<p>2D histograms of the <span class="html-italic">χ</span><sup>2</sup> value as a function of (<b>top</b>) the modelled <span class="html-italic">κ</span> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math> and (<b>bottom</b>) the modelled <span class="html-italic">κ</span> and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math>, as derived for plasma with two different input densities; (<b>left</b>) <span class="html-italic">n</span> = 10 cm<sup>−3</sup>, and (<b>right</b>) <span class="html-italic">n</span> = 50 cm<sup>−3</sup>. In both examples, we use <span class="html-italic">u</span><sub>0</sub> = 500 kms<sup>−1</sup> pointing along the x-axis, <span class="html-italic">κ</span> = 3, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math> = 10 eV, and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math> = 20 eV as input parameters.</p>
Full article ">Figure 7
<p>Number of counts as a function of energy for the pitch-angle with the maximum flux assuming a plasma with <span class="html-italic">n</span> = 5 cm<sup>−3</sup>, <span class="html-italic">κ</span> = 3, <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>∥</mo> </msub> </mrow> </semantics></math> = 10 eV and <math display="inline"><semantics> <mrow> <msub> <mi>T</mi> <mo>⊥</mo> </msub> </mrow> </semantics></math>= 20 eV. The blue line is the fitted model to the observations by (<b>left</b>) excluding points with <span class="html-italic">C<sub>i</sub></span> = 0 which are shown with red colour, and (<b>right</b>) including points with <span class="html-italic">C<sub>i</sub></span> = 0. The magenta line is the expected counts <span class="html-italic">C</span><sub>exp</sub>, given by Equation (2). The labels in each panel show the parameters as derived by the corresponding fit.</p>
Full article ">Figure 8
<p>Poisson distribution with average value (<b>top</b>) <span class="html-italic">C</span><sub>exp</sub> = 1, (<b>middle</b>) <span class="html-italic">C</span><sub>exp</sub> = 3, and (<b>bottom</b>) <span class="html-italic">C</span><sub>exp</sub> = 5. The vertical lines indicate the two modes of the distribution, <span class="html-italic">C</span><sub>exp</sub> (blue) and <span class="html-italic">C</span><sub>exp</sub> − 1 (orange) respectively. For small average values, the Poisson distribution is asymmetric, and the probability to measure number of counts lower than the average value is significant. This can bias the results to lower densities.</p>
Full article ">Figure 9
<p>Number of the expected average counts <span class="html-italic">C</span><sub>exp</sub> as a function of energy in the pitch-angle bin with the maximum particle flux, considering the same plasma conditions as in the example shown in <a href="#entropy-22-00103-f007" class="html-fig">Figure 7</a>. The blue curve is the model fitted to the observations by (<b>left</b>) excluding points with <span class="html-italic">C<sub>i</sub></span> = 0 which are shown with red colour, and (<b>right</b>) including points with <span class="html-italic">C<sub>i</sub></span> =0. The orange curve is the mode <span class="html-italic">C</span><sub>exp</sub> − 1. In each panel, we show the parameters as derived by the corresponding fit. In the absence of statistical fluctuations, both fitting strategies derive identical bulk parameters.</p>
Full article ">
19 pages, 2322 KiB  
Article
Learning in Feedforward Neural Networks Accelerated by Transfer Entropy
by Adrian Moldovan, Angel Caţaron and Răzvan Andonie
Entropy 2020, 22(1), 102; https://doi.org/10.3390/e22010102 - 16 Jan 2020
Cited by 18 | Viewed by 5464
Abstract
Current neural networks architectures are many times harder to train because of the increasing size and complexity of the used datasets. Our objective is to design more efficient training algorithms utilizing causal relationships inferred from neural networks. The transfer entropy (TE) was initially [...] Read more.
Current neural networks architectures are many times harder to train because of the increasing size and complexity of the used datasets. Our objective is to design more efficient training algorithms utilizing causal relationships inferred from neural networks. The transfer entropy (TE) was initially introduced as an information transfer measure used to quantify the statistical coherence between events (time series). Later, it was related to causality, even if they are not the same. There are only few papers reporting applications of causality or TE in neural networks. Our contribution is an information-theoretical method for analyzing information transfer between the nodes of feedforward neural networks. The information transfer is measured by the TE of feedback neural connections. Intuitively, TE measures the relevance of a connection in the network and the feedback amplifies this connection. We introduce a backpropagation type training algorithm that uses TE feedback connections to improve its performance. Full article
(This article belongs to the Section Information Theory, Probability and Statistics)
Show Figures

Figure 1

Figure 1
<p>This illustrates how the two neurons with indices <span class="html-italic">i</span> and <span class="html-italic">j</span> from a network produce a series of activations. The <span class="html-italic">g</span> threshold is the red line that splits these activations into two groups: the ones above the threshold (blue) and the ones below the threshold (red). They correspond to the <math display="inline"><semantics> <msubsup> <mi>o</mi> <mrow> <mi>i</mi> </mrow> <mrow> <mi>r</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> </semantics></math> and <math display="inline"><semantics> <msubsup> <mi>o</mi> <mrow> <mi>j</mi> </mrow> <mrow> <mi>r</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> </semantics></math> time series, which produce the time series of binary values <math display="inline"><semantics> <msubsup> <mi>s</mi> <mrow> <mi>i</mi> </mrow> <mrow> <mi>r</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> </semantics></math> and <math display="inline"><semantics> <msubsup> <mi>s</mi> <mrow> <mi>j</mi> </mrow> <mrow> <mi>r</mi> <mo>,</mo> <mi>n</mi> </mrow> </msubsup> </semantics></math> used to calculate the TE. The process is applied to all pairs of connected neurons.</p>
Full article ">Figure 2
<p>The <span class="html-italic">FF+FB</span> architecture for the XOR problem. The <math display="inline"><semantics> <mrow> <mi>t</mi> <msubsup> <mi>e</mi> <mrow> <mi>j</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>l</mi> </msubsup> </mrow> </semantics></math> values are calculated between neurons <span class="html-italic">j</span> and <span class="html-italic">i</span>. Neuron <span class="html-italic">i</span> is in layer <math display="inline"><semantics> <mrow> <mi>l</mi> <mo>−</mo> <mn>1</mn> </mrow> </semantics></math>, neuron <span class="html-italic">j</span> is in layer <span class="html-italic">l</span>. The colored arrows from the neurons show how the outputs from the neurons are used to calculate the <math display="inline"><semantics> <mrow> <mi>t</mi> <mi>e</mi> </mrow> </semantics></math> values. The same color in a dotted line arrow shows to which weight the <math display="inline"><semantics> <mrow> <mi>t</mi> <mi>e</mi> </mrow> </semantics></math> is applied (see Equation <a href="#FD3-entropy-22-00102" class="html-disp-formula">3</a>). The bias units are implemented but not shown here since they do not use the <math display="inline"><semantics> <mrow> <mi>t</mi> <mi>e</mi> </mrow> </semantics></math> values in the algorithm.</p>
Full article ">Figure 3
<p>Ten runs on XOR dataset. Each <span class="html-italic">x</span> axis finishes when the last of the <span class="html-italic">FF+FB</span> or <span class="html-italic">FF</span> reaches either the maximum number of epochs or 100% training accuracy (log scale).</p>
Full article ">Figure A1
<p>Ten runs of Abalone dataset. Each X axis finishes when the last of the <span class="html-italic">FF+FB</span> or <span class="html-italic">FF</span> reaches either the maximum number of epochs or the maximum set validation accuracy (log scale).</p>
Full article ">Figure A2
<p>We tried several <math display="inline"><semantics> <mi>η</mi> </semantics></math> values for this benchmark. To improve stability of both models, we selected a smaller value. It can be observed that <span class="html-italic">FF+FB</span> failed to converge on the 9th run and even to properly learn on other runs.</p>
Full article ">Figure A2 Cont.
<p>We tried several <math display="inline"><semantics> <mi>η</mi> </semantics></math> values for this benchmark. To improve stability of both models, we selected a smaller value. It can be observed that <span class="html-italic">FF+FB</span> failed to converge on the 9th run and even to properly learn on other runs.</p>
Full article ">Figure A2 Cont.
<p>We tried several <math display="inline"><semantics> <mi>η</mi> </semantics></math> values for this benchmark. To improve stability of both models, we selected a smaller value. It can be observed that <span class="html-italic">FF+FB</span> failed to converge on the 9th run and even to properly learn on other runs.</p>
Full article ">Figure A2 Cont.
<p>We tried several <math display="inline"><semantics> <mi>η</mi> </semantics></math> values for this benchmark. To improve stability of both models, we selected a smaller value. It can be observed that <span class="html-italic">FF+FB</span> failed to converge on the 9th run and even to properly learn on other runs.</p>
Full article ">Figure A3
<p>The <span class="html-italic">divorce</span> dataset constantly requires only a few epochs to reach the target accuracy.</p>
Full article ">
7 pages, 276 KiB  
Article
A Geometric Interpretation of Stochastic Gradient Descent Using Diffusion Metrics
by Rita Fioresi, Pratik Chaudhari and Stefano Soatto
Entropy 2020, 22(1), 101; https://doi.org/10.3390/e22010101 - 15 Jan 2020
Cited by 3 | Viewed by 5483
Abstract
This paper is a step towards developing a geometric understanding of a popular algorithm for training deep neural networks named stochastic gradient descent (SGD). We built upon a recent result which observed that the noise in SGD while training typical networks is highly [...] Read more.
This paper is a step towards developing a geometric understanding of a popular algorithm for training deep neural networks named stochastic gradient descent (SGD). We built upon a recent result which observed that the noise in SGD while training typical networks is highly non-isotropic. That motivated a deterministic model in which the trajectories of our dynamical systems are described via geodesics of a family of metrics arising from a certain diffusion matrix; namely, the covariance of the stochastic gradients in SGD. Our model is analogous to models in general relativity: the role of the electromagnetic field in the latter is played by the gradient of the loss function of a deep network in the former. Full article
(This article belongs to the Special Issue The Information Bottleneck in Deep Learning)
13 pages, 5513 KiB  
Article
On Heat Transfer Performance of Cooling Systems Using Nanofluid for Electric Motor Applications
by Ali Deriszadeh and Filippo de Monte
Entropy 2020, 22(1), 99; https://doi.org/10.3390/e22010099 - 14 Jan 2020
Cited by 29 | Viewed by 7804
Abstract
This paper studies the fluid flow and heat transfer characteristics of nanofluids as advance coolants for the cooling system of electric motors. Investigations are carried out using numerical analysis for a cooling system with spiral channels. To solve the governing equations, computational fluid [...] Read more.
This paper studies the fluid flow and heat transfer characteristics of nanofluids as advance coolants for the cooling system of electric motors. Investigations are carried out using numerical analysis for a cooling system with spiral channels. To solve the governing equations, computational fluid dynamics and 3D fluid motion analysis are used. The base fluid is water with a laminar flow. The fluid Reynolds number and turn-number of spiral channels are evaluation parameters. The effect of nanoparticles volume fraction in the base fluid on the heat transfer performance of the cooling system is studied. Increasing the volume fraction of nanoparticles leads to improving the heat transfer performance of the cooling system. On the other hand, a high-volume fraction of the nanofluid increases the pressure drop of the coolant fluid and increases the required pumping power. This paper aims at finding a trade-off between effective parameters by studying both fluid flow and heat transfer characteristics of the nanofluid. Full article
Show Figures

Figure 1

Figure 1
<p>Thermal limit for transient operation of an electric motor.</p>
Full article ">Figure 2
<p>Schematic of the cooling system with (<b>a</b>) 4 turns, (<b>b</b>) 6 turns and (<b>c</b>) 8 turns.</p>
Full article ">Figure 3
<p>Applied boundary conditions.</p>
Full article ">Figure 4
<p>Calculated heat transfer coefficient versus number of nodes.</p>
Full article ">Figure 5
<p>Average temperature of the motor versus Reynolds number for cooling systems using pure water coolant with different turns numbers.</p>
Full article ">Figure 6
<p>Effect of Reynolds number and turns number on heat transfer performance of the cooling system using pure water as a coolant.</p>
Full article ">Figure 7
<p>Pressure drop versus the Reynolds number for cooling systems using pure water coolant with three different turn numbers.</p>
Full article ">Figure 8
<p>The effect of nanoparticle volume fraction on the heat transfer.</p>
Full article ">Figure 9
<p>The effect of nanoparticles volume fraction on the pressure drop.</p>
Full article ">Figure 10
<p>Performance coefficient versus Reynolds number.</p>
Full article ">Figure 11
<p>Temperature distribution contours of the cooling system using the nanofluid at Reynolds number of 2000 with (<b>a</b>) 4 turns, (<b>b</b>) 6 turns and (<b>c</b>) 8 turns.</p>
Full article ">Figure 12
<p>Pressure distribution contours of the cooling system using the nanofluid with different turns numbers (<b>a</b>) 4 turns, (<b>b</b>) 6 turns and (<b>c</b>) 8 turns.</p>
Full article ">Figure 13
<p>Temperature distribution contours at Reynolds number of 2000 with volume fractions of (<b>a</b>) 2% and (<b>b</b>) 4%.</p>
Full article ">
27 pages, 1642 KiB  
Article
The Convex Information Bottleneck Lagrangian
by Borja Rodríguez Gálvez, Ragnar Thobaben and Mikael Skoglund
Entropy 2020, 22(1), 98; https://doi.org/10.3390/e22010098 - 14 Jan 2020
Cited by 20 | Viewed by 5201
Abstract
The information bottleneck (IB) problem tackles the issue of obtaining relevant compressed representations T of some random variable X for the task of predicting Y. It is defined as a constrained optimization problem that maximizes the information the representation has about the [...] Read more.
The information bottleneck (IB) problem tackles the issue of obtaining relevant compressed representations T of some random variable X for the task of predicting Y. It is defined as a constrained optimization problem that maximizes the information the representation has about the task, I ( T ; Y ) , while ensuring that a certain level of compression r is achieved (i.e., I ( X ; T ) r ). For practical reasons, the problem is usually solved by maximizing the IB Lagrangian (i.e., L IB ( T ; β ) = I ( T ; Y ) β I ( X ; T ) ) for many values of β [ 0 , 1 ] . Then, the curve of maximal I ( T ; Y ) for a given I ( X ; T ) is drawn and a representation with the desired predictability and compression is selected. It is known when Y is a deterministic function of X, the IB curve cannot be explored and another Lagrangian has been proposed to tackle this problem: the squared IB Lagrangian: L sq IB ( T ; β sq ) = I ( T ; Y ) β sq I ( X ; T ) 2 . In this paper, we (i) present a general family of Lagrangians which allow for the exploration of the IB curve in all scenarios; (ii) provide the exact one-to-one mapping between the Lagrange multiplier and the desired compression rate r for known IB curve shapes; and (iii) show we can approximately obtain a specific compression level with the convex IB Lagrangian for both known and unknown IB curve shapes. This eliminates the burden of solving the optimization problem for many values of the Lagrange multiplier. That is, we prove that we can solve the original constrained problem with a single optimization. Full article
(This article belongs to the Special Issue Information Bottleneck: Theory and Applications in Deep Learning)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The top row shows the results for the power information bottleneck (IB) Lagrangian with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>, and the bottom row for the exponential IB Lagrangian with <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>, both in the MNIST dataset. In each row, from left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set, the blue dots the values computed with the validation set and the green stars the theoretical values computed as dictated by Proposition 3. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mo form="prefix">log</mo> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </semantics></math> in a dashed, orange line. All values are shown in bits.</p>
Full article ">Figure 2
<p>Depiction of the clusterization behavior of the bottleneck variable for the power IB Lagrangian in the MNIST dataset with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>. The clusters were obtained using the DBSCAN algorithm [<a href="#B44-entropy-22-00098" class="html-bibr">44</a>,<a href="#B45-entropy-22-00098" class="html-bibr">45</a>].</p>
Full article ">Figure 3
<p>Example of value convergence with the exponential IB Lagrangian with <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mn>3</mn> </mrow> </semantics></math>. We show the intersection of the isolines of <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="script">L</mi> <mrow> <mi>IB</mi> <mo>,</mo> <mi>exp</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>T</mi> <mo>;</mo> <msub> <mi>β</mi> <mrow> <mi>exp</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </semantics></math> for different <math display="inline"><semantics> <mrow> <msub> <mi>β</mi> <mrow> <mi>exp</mi> </mrow> </msub> <mo>∈</mo> <msub> <mi>B</mi> <mrow> <mi>exp</mi> </mrow> </msub> <mo>≈</mo> <mrow> <mo>[</mo> <mn>1.56</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>5</mn> </mrow> </msup> <mo>,</mo> <msup> <mn>3</mn> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> <mo>]</mo> </mrow> </mrow> </semantics></math> using Corollary 2.</p>
Full article ">Figure 4
<p>Example of value convergence exploitation with the shifted exponential Lagrangian with <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mn>200</mn> </mrow> </semantics></math>. In the top row, for the MNIST dataset aiming for a compression level <math display="inline"><semantics> <mrow> <msup> <mi>r</mi> <mo>∗</mo> </msup> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> and in the bottom row, for the TREC-6 dataset aiming for a compression level of <math display="inline"><semantics> <mrow> <msup> <mi>r</mi> <mo>∗</mo> </msup> <mo>=</mo> <mn>16</mn> </mrow> </semantics></math>. In each row, from left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set, the blue dots the values computed with the validation set and the green stars the theoretical values computed as dictated by Proposition 3. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>H</mi> <mo>(</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> in a dashed, orange line. All values are shown in bits.</p>
Full article ">Figure A1
<p>Graphical representation of the IB curve in the information plane. Dashed lines in orange represent tight bounds confining the region (in light orange) of possible IB curves (delimited by the red line, also known as the Pareto frontier). Black dotted lines are informative values. In blue we show an example of a possible IB curve confining a region (in darker orange) of an IB curve that does not achieve the Pareto frontier. Finally, the yellow star represents the point where the representation keeps the same information about the input and the output.</p>
Full article ">Figure A2
<p>Results for the power IB Lagrangian in the MNIST dataset with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mo>{</mo> <mn>0.5</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>}</mo> </mrow> </semantics></math>, from top to bottom. In each row, from left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set, the blue dots the values computed with the validation set and the green stars the theoretical values computed as dictated by Proposition 3. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mo form="prefix">log</mo> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </semantics></math> in a dashed, orange line. All values are shown in bits.</p>
Full article ">Figure A3
<p>Results for the exponential IB Lagrangian in the MNIST dataset with <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mo>{</mo> <mo form="prefix">log</mo> <mo>(</mo> <mn>2</mn> <mo>)</mo> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>1.5</mn> <mo>}</mo> </mrow> </semantics></math>, from top to bottom. In each row, from left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set, the blue dots the values computed with the validation set and the gren stars the theoretical values computed as dictated by Proposition 3. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mo form="prefix">log</mo> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </semantics></math> in a dashed, orange line. All values are shown in bits.</p>
Full article ">Figure A4
<p>Depiction of the clusterization behavior of the bottleneck variable in the MNIST dataset. In the first row, from left to right, the power IB Lagrangian with different values of <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mo>{</mo> <mn>0.5</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>}</mo> </mrow> </semantics></math>. In the second row, from left to right, the exponential IB Lagrangian with different values of <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mo>{</mo> <mo form="prefix">log</mo> <mo>(</mo> <mn>2</mn> <mo>)</mo> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mn>1.5</mn> <mo>}</mo> </mrow> </semantics></math>.</p>
Full article ">Figure A5
<p>Results for the exponential IB Lagrangian in the Fashion MNIST dataset with <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>. From left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set and the blue dots the values computed with the validation set. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mo form="prefix">log</mo> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </semantics></math>. All values are shown in bits.</p>
Full article ">Figure A6
<p>The top row shows the results for the normal IB Lagrangian, and the bottom row for the exponential IB Lagrangian with <math display="inline"><semantics> <mrow> <mi>η</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>, both in the California housing dataset. In each row, from left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set and the blue dots the values computed with the validation set. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as the empirical value obtained maximizing <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> without compression limitations as in [<a href="#B26-entropy-22-00098" class="html-bibr">26</a>]. All values are shown in bits.</p>
Full article ">Figure A7
<p>The top row shows the results for the normal IB Lagrangian, and the bottom row for the power IB Lagrangian with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>0.1</mn> </mrow> </semantics></math>, both in the TREC-6 dataset. In each row, from left to right it is shown (i) the information plane, where the region of possible solutions of the IB problem is shadowed in light orange and the information-theoretic limits are the dashed orange line; (ii) <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>T</mi> <mo>;</mo> <mi>Y</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>; and (iii) the compression <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> as a function of <math display="inline"><semantics> <msub> <mi>β</mi> <mi>u</mi> </msub> </semantics></math>. In all plots, the red crosses joined by a dotted line represent the values computed with the training set and the blue dots the values computed with the validation set. Moreover, in all plots, it is indicated <math display="inline"><semantics> <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mo form="prefix">log</mo> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </semantics></math>. All values are shown in bits.</p>
Full article ">Figure A8
<p>Theoretical bijection between <math display="inline"><semantics> <mrow> <mi>I</mi> <mo>(</mo> <mi>X</mi> <mo>;</mo> <mi>T</mi> <mo>)</mo> </mrow> </semantics></math> and different <math display="inline"><semantics> <mi>α</mi> </semantics></math> from <math display="inline"><semantics> <msub> <mi>β</mi> <mrow> <mi>u</mi> <mo>,</mo> <mi>min</mi> </mrow> </msub> </semantics></math> to 1.5 in the power IB Lagrangian (<b>top</b>), and different <math display="inline"><semantics> <mi>η</mi> </semantics></math> in the domain <math display="inline"><semantics> <msub> <mi>B</mi> <mi>u</mi> </msub> </semantics></math> in the exponential IB Lagrangian (<b>bottom</b>).</p>
Full article ">
20 pages, 1424 KiB  
Review
A Review of the Application of Information Theory to Clinical Diagnostic Testing
by William A. Benish
Entropy 2020, 22(1), 97; https://doi.org/10.3390/e22010097 - 14 Jan 2020
Cited by 14 | Viewed by 5418
Abstract
The fundamental information theory functions of entropy, relative entropy, and mutual information are directly applicable to clinical diagnostic testing. This is a consequence of the fact that an individual’s disease state and diagnostic test result are random variables. In this paper, we review [...] Read more.
The fundamental information theory functions of entropy, relative entropy, and mutual information are directly applicable to clinical diagnostic testing. This is a consequence of the fact that an individual’s disease state and diagnostic test result are random variables. In this paper, we review the application of information theory to the quantification of diagnostic uncertainty, diagnostic information, and diagnostic test performance. An advantage of information theory functions over more established test performance measures is that they can be used when multiple disease states are under consideration as well as when the diagnostic test can yield multiple or continuous results. Since more than one diagnostic test is often required to help determine a patient’s disease state, we also discuss the application of the theory to situations in which more than one diagnostic test is used. The total diagnostic information provided by two or more tests can be partitioned into meaningful components. Full article
(This article belongs to the Special Issue Applications of Information Theory to Epidemiology)
Show Figures

Figure 1

Figure 1
<p>Surprisal (in bits) as a function of probability.</p>
Full article ">Figure 2
<p>The area under the receiver operating characteristic curve (AUC), a transformation of the AUC (AUC*), and the mutual information between the disease state and the test result (<math display="inline"><semantics> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mrow> <mi>D</mi> <mo>;</mo> <mi>R</mi> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math> in units of bits), as a function of the distance between the means of the distributions of test results for both healthy (d<math display="inline"><semantics> <mo>−</mo> </semantics></math>) and diseased (d<math display="inline"><semantics> <mo>+</mo> </semantics></math>) individuals (see text). <math display="inline"><semantics> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mrow> <mi>D</mi> <mo>;</mo> <mi>R</mi> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math> is plotted for three pretest probabilities of disease: 0.1, 0.2, and 0.5.</p>
Full article ">Figure 3
<p>Mutual information as a function of pretest probability of significant coronary heart disease (CHD) for the exercise electrocardiogram. The plot compares the performance of a theoretical ideal test with the actual performance when the results are either (1) dichotomized using the criterion of ST segment depression of <math display="inline"><semantics> <mo>≥</mo> </semantics></math> 1 mm or (2) not dichotomized. This plot has been reconstructed with permission from the paper by Diamond et al. [<a href="#B6-entropy-22-00097" class="html-bibr">6</a>].</p>
Full article ">Figure 4
<p>Diagnostic information (relative entropy) provided by the findings of an ST segment depression (STdep) &lt; 0.5 mm and an ST depression (STdep) ≥ 2.5 mm as a function of pretest probability of significant coronary heart disease (CHD). Also shown are relative entropy plots for a theoretical ideal test when significant CHD is present (d+) and when significant CHD is absent (d<math display="inline"><semantics> <mo>−</mo> </semantics></math>). For the theoretical ideal test when significant CHD is present, relative entropy increases indefinitely as pretest probability of significant CHD approaches zero; and for the theoretical ideal test when significant CHD is absent, relative entropy increases indefinitely as pretest probability of significant CHD approaches one.</p>
Full article ">Figure 5
<p>Contour plots showing the probability densities of the results of tests A and B for the healthy and diseased populations. The test results are expressed as standard deviations from the means of the distribution for healthy patients <math display="inline"><semantics> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>μ</mi> <mi>A</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Contour plot showing the information (relative entropy in units of bits) provided by specific combinations of results of tests A and B. The test results are expressed as standard deviations from the means of the distribution for healthy patients <math display="inline"><semantics> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>μ</mi> <mi>A</mi> </msub> <mo>=</mo> <msub> <mi>μ</mi> <mi>B</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> <mo>)</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop