Category: Statistics

Phys Rev E. 2023 May;107(5-1):054128. doi: 10.1103/PhysRevE.107.054128.

ABSTRACT

Embedded random matrix ensembles with k-body interactions are well established to be appropriate for many quantum systems. For these ensembles the two point correlation function is not yet derived, though these ensembles are introduced 50 years back. Two-point correlation function in eigenvalues of a random matrix ensemble is the ensemble average of the product of the density of eigenvalues at two eigenvalues, say E and E^{‘}. Fluctuation measures such as the number variance and Dyson-Mehta Δ_{3} statistic are defined by the two-point function and so also the variance of the level motion in the ensemble. Recently, it is recognized that for the embedded ensembles with k-body interactions the one-point function (ensemble averaged density of eigenvalues) follows the so called q-normal distribution. With this, the eigenvalue density can be expanded by starting with the q-normal form and using the associated q-Hermite polynomials He_{ζ}(x|q). Covariances S_{ζ}S_{ζ^{‘}}[over ¯] (overline representing ensemble average) of the expansion coefficients S_{ζ} with ζ≥1 here determine the two-point function as they are a linear combination of the bivariate moments Σ_{PQ} of the two-point function. Besides describing all these, in this paper formulas are derived for the bivariate moments Σ_{PQ} with P+Q≤8, of the two-point correlation function, for the embedded Gaussian unitary ensembles with k-body interactions [EGUE(k)] as appropriate for systems with m fermions in N single particle states. Used for obtaining the formulas is the SU(N) Wigner-Racah algebra. These formulas with finite N corrections are used to derive formulas for the covariances S_{ζ}S_{ζ^{‘}}[over ¯] in the asymptotic limit. These show that the present work extends to all k values, the results known in the past in the two extreme limits with k/m→0 (same as q→1) and k=m (same as q=0).

PMID:37329068 | DOI:10.1103/PhysRevE.107.054128

Phys Rev E. 2023 May;107(5-2):055304. doi: 10.1103/PhysRevE.107.055304.

ABSTRACT

We present a general and numerically efficient method for calculation of collision integrals for interacting quantum gases on a discrete momentum lattice. Here we employ the original analytical approach based on Fourier transform covering a wide spectrum of solid-state problems with various particle statistics and arbitrary interaction models, including the case of momentum-dependent interaction. The comprehensive set of the transformation principles is given in detail and realized as a computer Fortran 90 library FLBE (Fast Library for Boltzmann Equation).

PMID:37329067 | DOI:10.1103/PhysRevE.107.055304

Phys Rev E. 2023 May;107(5-2):055212. doi: 10.1103/PhysRevE.107.055212.

ABSTRACT

A recent numerical treatment of data obtained by the Parker Solar Probe spacecraft describes the electron concentration in solar wind as a function of the heliocentric distance based on a Kappa distribution with spectral index κ=5. In this work, we derive and, subsequently, solve an entirely different class of nonlinear partial differential equations describing the one-dimensional diffusion of a suprathermal gas. The theory is applied to describe the aforementioned data and we find a spectral index κ≳1.5 providing the widely acknowledged identification of Kappa electrons in solar wind. We also find that suprathermal effects increase the length scale of classical diffusion by one order of magnitude. Such a result does not depend on the microscopic details of the diffusion coefficient since our theory is based on a macroscopic formulation. Forthcoming extensions of our theory by including magnetic fields and relating our formulation to nonextensive statistics are briefly addressed.

PMID:37329056 | DOI:10.1103/PhysRevE.107.055212

Phys Rev E. 2023 May;107(5-1):054139. doi: 10.1103/PhysRevE.107.054139.

ABSTRACT

A variety of transport processes in natural and man-made systems are intrinsically random. To model their stochasticity, lattice random walks have been employed for a long time, mainly by considering Cartesian lattices. However, in many applications in bounded space the geometry of the domain may have profound effects on the dynamics and ought to be accounted for. We consider here the cases of the six-neighbor (hexagonal) and three-neighbor (honeycomb) lattices, which are utilized in models ranging from adatoms diffusing in metals and excitations diffusing on single-walled carbon nanotubes to animal foraging strategy and the formation of territories in scent-marking organisms. In these and other examples, the main theoretical tool to study the dynamics of lattice random walks in hexagonal geometries has been via simulations. Analytic representations have in most cases been inaccessible, in particular in bounded hexagons, given the complicated “zigzag” boundary conditions that a walker is subject to. Here we generalize the method of images to hexagonal geometries and obtain closed-form expressions for the occupation probability, the so-called propagator, for lattice random walks both on hexagonal and honeycomb lattices with periodic, reflective, and absorbing boundary conditions. In the periodic case, we identify two possible choices of image placement and their corresponding propagators. Using them, we construct the exact propagators for the other boundary conditions, and we derive transport-related statistical quantities such as first-passage probabilities to one or multiple targets and their means, elucidating the effect of the boundary condition on transport properties.

PMID:37329046 | DOI:10.1103/PhysRevE.107.054139

Phys Rev E. 2023 May;107(5-2):055309. doi: 10.1103/PhysRevE.107.055309.

ABSTRACT

Digital cores can characterize the true internal structure of rocks at the pore scale. This method has become one of the most effective ways to quantitatively analyze the pore structure and other properties of digital cores in rock physics and petroleum science. Deep learning can precisely extract features from training images for a rapid reconstruction of digital cores. Usually, the reconstruction of three-dimensional (3D) digital cores is performed by optimization using generative adversarial networks. The training data required for the 3D reconstruction are 3D training images. In practice, two-dimensional (2D) imaging devices are widely used because they can achieve faster imaging, higher resolution, and easier identification of different rock phases, so replacing 3D images with 2D ones avoids the difficulty of acquiring 3D images. In this paper, we propose a method, named EWGAN-GP, for the reconstruction of 3D structures from a 2D image. Our proposed method includes an encoder, a generator, and three discriminators. The main purpose of the encoder is to extract statistical features of a 2D image. The generator extends the extracted features into 3D data structures. Meanwhile, the three discriminators have been designed to gauge the similarity of morphological characteristics between cross sections of the reconstructed 3D structure and the real image. The porosity loss function is used to control the distribution of each phase in general. In the entire optimization process, a strategy using Wasserstein distance with gradient penalty makes the convergence of the training process faster and the reconstruction result more stable; it also avoids the problems of gradient disappearance and mode collapse. Finally, the reconstructed 3D structure and the target 3D structure are visualized to ascertain their similar morphologies. The morphological parameter indicators of the reconstructed 3D structure were consistent with those of the target 3D structure. The microstructure parameters of the 3D structure were also compared and analyzed. The proposed method can achieve accurate and stable 3D reconstruction compared with classical stochastic methods of image reconstruction.

PMID:37329045 | DOI:10.1103/PhysRevE.107.055309

J Genet Couns. 2023 Jun 16. doi: 10.1002/jgc4.1738. Online ahead of print.

ABSTRACT

Pathogenic and likely pathogenic variants in BRCA1 and BRCA2 (BRCA1/2) are medically actionable and may inform hereditary breast and ovarian cancer (HBOC) treatment and prevention. However, rates of germline genetic testing (GT) in people with and without cancer are suboptimal. Individuals’ knowledge, attitudes, and beliefs may influence GT decisions. While genetic counseling (GC) provides decision support, the supply of genetic counselors is insufficient to meet demand. Accordingly, there is a need to explore the evidence on interventions that aim to support BRCA1/2 testing decisions. We conducted a scoping review of PubMed, CINAHL, Web of Science, and PsycINFO using search terms related to HBOC, GT, and decision making. First, we screened records to identify peer-reviewed reports that described interventions to support BRCA1/2 testing decisions. Next, we reviewed full-text reports and excluded studies that lacked statistical comparisons or enrolled previously tested individuals. Finally, we extracted study characteristics and findings into a table. All records and reports were reviewed independently by two authors; decisions were tracked in Rayyan, and discrepancies were resolved through discussion. Of 2116 unique citations, 25 met the eligibility criteria. Articles were published between 1997 and 2021 and described randomized trials and nonrandomized, quasi-experimental studies. Most studies tested technology-based (12/25, 48%) or written (9/25, 36%) interventions. Nearly half (12/25, 48%) of interventions were designed to complement traditional GC. Of the interventions compared to GC, 75% (6/8) increased or had a noninferior effect on knowledge, and 67% (4/6) decreased or had a noninferior effect on decisional conflict. Intervention effects on GT uptake were mixed, which may reflect evolving eligibility criteria for GT. Our findings suggest novel interventions may promote informed GT decision making, but many were developed to complement traditional GC. Trials that assess the effects of decision support interventions in diverse samples and evaluate implementation strategies for efficacious interventions are warranted.

PMID:37328917 | DOI:10.1002/jgc4.1738

Vet Res. 2023 Jun 16;54(1):49. doi: 10.1186/s13567-023-01180-y.

ABSTRACT

Swine dysentery (SD) caused by pathogenic Brachyspira spp. is an economic challenge for the swine industry. In research settings, experimental reproduction of swine dysentery typically relies on intragastric inoculation which has shown variable success. This project aimed to improve the consistency of the experimental inoculation protocol used for swine dysentery in our laboratory. Over six experiments, we evaluated the influence of group housing in inoculated pigs using a frozen-thawed broth culture of strongly hemolytic B. hyodysenteriae strain D19 (Trial A), compared the relative virulence of B. hyodysenteriae strains D19 and G44 (Trial B), compared inoculum volumes (50 mL vs 100 mL) for G44 and B. hampsonii 30446 (Trial C), and performed three independent trials evaluating intragastric inoculation using different oral inoculation methods: oral feed balls (Trial D), and oral syringe bolus of 100 mL (Trial E) or 300 mL (Trial F). Intragastric inoculation with a fresh broth culture of B. hyodysenteriae strain G44 resulted in a shorter incubation period and a higher proportionate duration of mucohemorrhagic diarrhea (MMHD) compared to D19. Intragastric inoculation with either 50 or 100 mL of B. hampsonii 30446 or B. hyodysenteriae (G44) were statistically equivalent. Oral inoculation with 100 mL or 300 mL also yielded similar results to intragastric inoculation but was more expensive due to the additional work and supplies associated with syringe training. Our future research will use intragastric inoculation with 100 mL of a fresh broth culture containing B. hyodysenteriae strain G44 as it yields a high incidence of mucohaemorrhagic diarrhea with a reasonable cost.

PMID:37328906 | DOI:10.1186/s13567-023-01180-y

BMC Oral Health. 2023 Jun 16;23(1):397. doi: 10.1186/s12903-023-03104-w.

ABSTRACT

BACKGROUND: Three-dimensional (3D) printing is a novel innovation in the field of craniomaxillofacial surgery, however, a lack of evidence exists related to the comparison of the trueness of skull models fabricated using different technology-based printers belonging to different cost segments.

METHODS: A study was performed to investigate the trueness of cone-beam computed tomography-derived skull models fabricated using different technology based on low-, medium-, and high-cost 3D printers. Following the segmentation of a patient’s skull, the model was printed by: (i) a low-cost fused filament fabrication printer; (ii) a medium-cost stereolithography printer; and (iii) a high-cost material jetting printer. The fabricated models were later scanned by industrial computed tomography and superimposed onto the original reference virtual model by applying surface-based registration. A part comparison color-coded analysis was conducted for assessing the difference between the reference and scanned models. A one-way analysis of variance (ANOVA) with Bonferroni correction was applied for statistical analysis.

RESULTS: The model printed with the low-cost fused filament fabrication printer showed the highest mean absolute error ([Formula: see text]), whereas both medium-cost stereolithography-based and the high-cost material jetting models had an overall similar dimensional error of [Formula: see text] and [Formula: see text], respectively. Overall, the models printed with medium- and high-cost printers showed a significantly ([Formula: see text]) lower error compared to the low-cost printer.

CONCLUSIONS: Both stereolithography and material jetting based printers, belonging to the medium- and high-cost market segment, were able to replicate the skeletal anatomy with optimal trueness, which might be suitable for patient-specific treatment planning tasks in craniomaxillofacial surgery. In contrast, the low-cost fused filament fabrication printer could serve as a cost-effective alternative for anatomical education, and/or patient communication.

PMID:37328901 | DOI:10.1186/s12903-023-03104-w

Algorithms Mol Biol. 2023 Jun 17;18(1):3. doi: 10.1186/s13015-023-00226-2.

ABSTRACT

We consider the problem of representing a set of [Formula: see text]-mers and their abundance counts, or weights, in compressed space so that assessing membership and retrieving the weight of a [Formula: see text]-mer is efficient. The representation is called a weighted dictionary of [Formula: see text]-mers and finds application in numerous tasks in Bioinformatics that usually count [Formula: see text]-mers as a pre-processing step. In fact, [Formula: see text]-mer counting tools produce very large outputs that may result in a severe bottleneck for subsequent processing. In this work we extend the recently introduced SSHash dictionary (Pibiri in Bioinformatics 38:185-194, 2022) to also store compactly the weights of the [Formula: see text]-mers. From a technical perspective, we exploit the order of the [Formula: see text]-mers represented in SSHash to encode runs of weights, hence allowing much better compression than the empirical entropy of the weights. We study the problem of reducing the number of runs in the weights to improve compression even further and give an optimal algorithm for this problem. Lastly, we corroborate our findings with experiments on real-world datasets and comparison with competitive alternatives. Up to date, SSHash is the only [Formula: see text]-mer dictionary that is exact, weighted, associative, fast, and small.

PMID:37328897 | DOI:10.1186/s13015-023-00226-2

BMC Musculoskelet Disord. 2023 Jun 16;24(1):495. doi: 10.1186/s12891-023-06632-2.

ABSTRACT

BACKGROUND: Data on the epidemiological characteristics and prognostic factors of patients with pelvic fractures are lacking, particularly in China. This study aimed to summarise the clinical and epidemiological characteristics of patients with pelvic fractures in eastern Zhejiang Province, China, and to identify risk factors for poor prognosis.

METHODS: The clinical data of 369 patients with pelvic fractures admitted to the Ningbo No. 6 Hospital between September 2020 and September 2021 were retrospectively analysed. Data on the demographic characteristics; fracture classification; injury time, cause, and site; treatment plan; and prognosis were collected using the Picture Archiving and Communication System and the Hospital Information System. Differences in constituent proportions were analysed using the chi-square test. Logistic regression analysis was used to identify factors affecting patient prognosis. Statistical significance was set at p ≤ 0.05.

RESULTS: Among the 369 patients, there were 206 men and 163 women, at a ratio of 1.26:1, and the average age was 53.64 ± 0.78 years. More than 50% of patients were aged 41-65 years. The average length of hospital stay was 18.88 ± 1.78 days. The three most common causes of pelvic fractures were traffic accidents (51.2%), falls from height (31.44%), and flat-ground falls (14.09%). There were significant differences in the distribution of the three causes of injury depending on age (p < 0.001), sex (p < 0.001), and occupation (p < 0.0001). Most patients were manual workers (48.8%). Furthermore, most patients (n = 262, 71.0%) underwent surgical treatment for pelvic fractures. Postoperative complications occurred in 26 patients (7.05%), and infection was the main complication (73.08%). Age (p = 0.013), occupation (p = 0.034), cause of injury (p = 0.022), treatment options (p = 0.001), and complications (p < 0.0001) were independent factors affecting the prognosis of patients with pelvic fractures. One death (0.027%) occurred, which was due to severe blood loss.

CONCLUSIONS: Age, occupation, cause of injury, treatment options and complications were factors affecting patient prognosis. In addition, changes in blood flow and prevention of infection warrant attention.

PMID:37328894 | DOI:10.1186/s12891-023-06632-2