Extraction of diffuse correlation spectroscopy flow index by integration of Nth-order linear model with Monte Carlo simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Shang Y
PROVIDER: S-EPMC4032444 | biostudies-other | 2014 May
REPOSITORIES: biostudies-other
Similar Datasets
Project description:The Monte Carlo simulation algorithm of photon trajectory computation is implemented in object oriented R code. Diffuse reflectance, also called backscattering, is modeled in semi-infinite homogeneous media. Spatial photon flux leaving the surface of the media is collected. The profile of intensity along radii relative to the incident point is used to simulate measurement of computer vision systems. Four optical parameters of the media are used: absorption coefficient, scattering coefficient, anisotropy factor and refractive index. Five parameters are used to describe configuration of the vision system: number of photons, radius of circular light beam, limiting energy level of photons, radius of observed area, spatial resolution of the vision system.•The incident angle of the light beam is included in the photon launch procedure. Initial direction is typically assumed to be normal with x,y,z coordinates of 0,0,1. In the proposed modification, initial move vector is calculated based on the incident angle and refractive index of the media. Additionally, elliptic distortion of the circular light beam on the surface is calculated based on the incident angle.•Photon flux leaving media through the surface is corrected with Lambertian method to measure intensity captured by an imaging device in normal position.•The software implementing the method is written in R language, the R code is available as standard package.
| S-EPMC7327276 | biostudies-literature
Project description:BACKGROUND:Self-contained tests estimate and test the association between a phenotype and mean expression level in a gene set defined a priori. Many self-contained gene set analysis methods have been developed but the performance of these methods for phenotypes that are continuous rather than discrete and with multiple nuisance covariates has not been well studied. Here, I use Monte Carlo simulation to evaluate the performance of both novel and previously published (and readily available via R) methods for inferring effects of a continuous predictor on mean expression in the presence of nuisance covariates. The motivating data are a high-profile dataset which was used to show opposing effects of hedonic and eudaimonic well-being (or happiness) on the mean expression level of a set of genes that has been correlated with social adversity (the CTRA gene set). The original analysis of these data used a linear model (GLS) of fixed effects with correlated error to infer effects of Hedonia and Eudaimonia on mean CTRA expression. METHODS:The standardized effects of Hedonia and Eudaimonia on CTRA gene set expression estimated by GLS were compared to estimates using multivariate (OLS) linear models and generalized estimating equation (GEE) models. The OLS estimates were tested using O'Brien's OLS test, Anderson's permutation [Formula: see text]-test, two permutation F-tests (including GlobalAncova), and a rotation z-test (Roast). The GEE estimates were tested using a Wald test with robust standard errors. The performance (Type I, II, S, and M errors) of all tests was investigated using a Monte Carlo simulation of data explicitly modeled on the re-analyzed dataset. RESULTS:GLS estimates are inconsistent between data sets, and, in each dataset, at least one coefficient is large and highly statistically significant. By contrast, effects estimated by OLS or GEE are very small, especially relative to the standard errors. Bootstrap and permutation GLS distributions suggest that the GLS results in downward biased standard errors and inflated coefficients. The Monte Carlo simulation of error rates shows highly inflated Type I error from the GLS test and slightly inflated Type I error from the GEE test. By contrast, Type I error for all OLS tests are at the nominal level. The permutation F-tests have ?1.9X the power of the other OLS tests. This increased power comes at a cost of high sign error (?10%) if tested on small effects. DISCUSSION:The apparently replicated pattern of well-being effects on gene expression is most parsimoniously explained as "correlated noise" due to the geometry of multiple regression. The GLS for fixed effects with correlated error, or any linear mixed model for estimating fixed effects in designs with many repeated measures or outcomes, should be used cautiously because of the inflated Type I and M error. By contrast, all OLS tests perform well, and the permutation F-tests have superior performance, including moderate power for very small effects.
| S-EPMC5068350 | biostudies-literature
Project description:Conventional mammography can suffer from poor contrast between healthy and cancerous tissues due to the small difference in attenuation properties. Coherent scatter slot scan imaging is an imaging technique which provides additional information and is compatible with conventional mammography. A Monte Carlo simulation of coherent scatter slot scan imaging was performed to assess its performance and provide system optimization. Coherent scatter could be exploited using a system similar to conventional slot scan mammography system with antiscatter grids tilted at the characteristic angle of cancerous tissues. System optimization was performed across several parameters, including source voltage, tilt angle, grid distances, grid ratio, and shielding geometry. The simulated carcinomas were detectable for tumors as small as 5 mm in diameter, so coherent scatter analysis using a wide-slot setup could be promising as an enhancement for screening mammography. Employing coherent scatter information simultaneously with conventional mammography could yield a conventional high spatial resolution image with additional coherent scatter information.
| S-EPMC4999523 | biostudies-other
Project description:Traditional multiple-group confirmatory factor analysis (multiple-group CFA) is usually criticized for having too restrictive model assumption, namely the scalar measurement invariance. The new multiple-group analysis methodology, alignment, has become an effective alternative. The alignment evaluates measurement invariance and more importantly, permits factor mean comparisons without requiring scalar invariance which is usually required in traditional multiple-group CFA. Some simulation studies and empirical studies have investigated the applicability of alignment under different conditions, but some areas remain unexplored. Based on the simulation studies of Asparouhov and Muthén and of Flake and McCoach, this current simulation study is broken into two sections. The first study investigates the minimal group sizes required for alignment in three-factor models. The second study compares the performance of multiple-group CFA, multiple-group exploratory structural equation model (multiple-group ESEM), and alignment by including different proportions and magnitudes of cross-loadings in the models. Study 1 shows that when the model has no noninvariant parameters, the alignment requires relatively lower group sizes. Explicitly, the minimal group size required for alignment was 250 when the amount of groups was three, the minimal group size was 150 when the amount of groups was nine, and 200 when the amount of groups was 15. When there are noninvariant parameters in the model and the amount of groups is low, a group size of 350 is a safe rule of thumb. When there are noninvariant parameters in the model and the amount of groups is high, a group size of 250 is required for trustworthy results. The magnitude of noninvariance and the noninvariance rate do not affect the minimal group size required for alignment. Study 2 shows that multiple-group CFA provides accurate factor mean estimates when each factor had 20% factor loading (1 factor loading) with small-sized cross-loading. Multiple-group ESEM provides accurate factor mean estimates when the magnitude of cross-loading is small or when each factor had 20% factor loading (1 factor loading) with medium-sized cross-loading. Alignment provides accurate factor mean estimates when there are only small-sized cross-loadings in the model. The parameter estimates, coverage rates and ratios of average standard error to standard deviation for each methodology are not influenced by the amount of groups. Recommendations are concluded for using multiple-group CFA, multiple-group ESEM, traditional alignment and aligned ESEM (AESEM) based on the results. Multiple-group CFA is more suitable for use when scalar invariance is established. Multiple-group ESEM works best when there are small-sized or only a few medium-sized cross-loadings in the model. Traditional alignment allows for small-sized cross-loadings and a few noninvariant parameters in the model. AESEM integrates the advantages of alignment and ESEM, can provide accurate estimates when noninvariant parameters and cross-loadings both exist in the model. Compared to multiple-group CFA, multiple-group ESEM, the alignment methodology performs well in more situations.
| S-EPMC9263979 | biostudies-literature
Project description:A computational method is developed to carry out explicit solvent simulations of complex molecular systems under conditions of constant pH. In constant-pH simulations, preidentified ionizable sites are allowed to spontaneously protonate and deprotonate as a function of time in response to the environment and the imposed pH. The method, based on a hybrid scheme originally proposed by H. A. Stern (J. Chem. Phys. 2007, 126, 164112), consists of carrying out short nonequilibrium molecular dynamics (neMD) switching trajectories to generate physically plausible configurations with changed protonation states that are subsequently accepted or rejected according to a Metropolis Monte Carlo (MC) criterion. To ensure microscopic detailed balance arising from such nonequilibrium switches, the atomic momenta are altered according to the symmetric two-ends momentum reversal prescription. To achieve higher efficiency, the original neMD-MC scheme is separated into two steps, reducing the need for generating a large number of unproductive and costly nonequilibrium trajectories. In the first step, the protonation state of a site is randomly attributed via a Metropolis MC process on the basis of an intrinsic pKa; an attempted nonequilibrium switch is generated only if this change in protonation state is accepted. This hybrid two-step inherent pKa neMD-MC simulation method is tested with single amino acids in solution (Asp, Glu, and His) and then applied to turkey ovomucoid third domain and hen egg-white lysozyme. Because of the simple linear increase in the computational cost relative to the number of titratable sites, the present method is naturally able to treat extremely large systems.
| S-EPMC4535364 | biostudies-literature
Project description:BACKGROUND: In the biological sciences the TCID50 (median tissue culture infective dose) assay is often used to determine the strength of a virus. METHODS: When the so-called Spearman-Kaerber calculation is used, the ratio between the pfu (the number of plaque forming units, the effective number of virus particles) and the TCID50, theoretically approaches a simple function of Eulers constant. Further, the standard deviation of the logarithm of the TCID50 approaches a simple function of the dilution factor and the number of wells used for determining the ratios in the assay. However, these theoretical calculations assume that the dilutions of the assay are independent, and in practice this is not completely correct. The assay was simulated using Monte Carlo techniques. RESULTS: Our simulation studies show that the theoretical results actually hold true for practical implementations of the assay. Furthermore, the simulation studies show that the distribution of the (the log of) TCID50, although discrete in nature, has a close relationship to the normal distribution. CONCLUSION: The pfu is proportional to the TCID50 titre with a factor of about 0.56 when using the Spearman-Kaerber calculation method. The normal distribution can be used for statistical inferences and ANOVA on the (the log of) TCID50 values is meaningful with group sizes of 5 and above.
| S-EPMC3331834 | biostudies-other
Project description:Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10(7) - 10(9) detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS.The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels.On the order of only 10-100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10(7) reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10(6) CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU.The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.
| S-EPMC3978362 | biostudies-other
Project description:Two new statistical models based on Monte Carlo Simulation (MCS) have been developed to score peptide matches in shotgun proteomic data and incorporated in a database search program, MassMatrix (www.massmatrix.net). The first model evaluates peptide matches based on the total abundance of matched peaks in the experimental spectra. The second model evaluates amino acid residue tags within MS/MS spectra. The two models provide complementary scores for peptide matches that result in higher confidence in peptide identification when significant scores are returned from both models. The MCS-based models use a variance reduction technique that improves estimation precision. Due to the high computational expense of MCS-based models, peptide matches were prefiltered by other statistical models before further evaluation by the MCS-based models. Receiver operating characteristic analysis of the data sets confirmed that MCS-based models improved the overall performance of the MassMatrix search software, especially for low-mass accuracy data sets.
| S-EPMC2749500 | biostudies-literature
Project description:Molecular simulations based on classical force fields are computationally efficient but lack accuracy due to the empirical formulation of non-bonded interactions. Quantum mechanical (QM) methods, albeit accurate, have inhibitory computational costs for large molecules and clusters. Hence, to overcome the bottleneck, machine learning (ML)-based methods have been employed in the recent years. We had earlier reported a combined scheme of many-body expansion (MBE) and ML to predict the interaction energies of rigid water clusters. In this work, we proceed toward building a flexible water model using the ML-MBE scheme. This ML-MBE scheme has an error of <1% for interaction energy prediction in comparison to the parent QM method for flexible water decamers. Machine learning-based Monte Carlo simulations (MLMC) are performed with this water model, and the structural properties of these configurations are compared with those obtained from ab initio molecular dynamics (AIMD) and the TIP3P classical force field. The radial distribution functions, tetrahedral order parameters, and number of hydrogen bonds in AIMD and MLMC have a similar qualitative and quantitative trend, whereas the classical force fields show a significant deviation.
| S-EPMC7143414 | biostudies-literature
Project description:Protein-ligand docking is an important approach for virtual screening and protein function annotation. Although many docking methods have been developed, most require a high-resolution crystal structure of the receptor and a user-specified binding site to start. This information is, however, not available for the majority of unknown proteins, including many pharmaceutically important targets. Developing blind docking methods without predefined binding sites and working with low-resolution receptor models from protein structure prediction is thus essential. In this manuscript, we propose a novel Monte Carlo based method, EDock, for blind protein-ligand docking. For a given protein, binding sites are first predicted by sequence-profile and substructure-based comparison searches with initial ligand poses generated by graph matching. Next, replica-exchange Monte Carlo (REMC) simulations are performed for ligand conformation refinement under the guidance of a physical force field coupled with binding-site distance constraints. The method was tested on two large-scale datasets containing 535 protein-ligand pairs. Without specifying binding pockets on the experimental receptor structures, EDock achieves on average a ligand RMSD of 2.03 Å, which compares favorably with state-of-the-art docking methods including DOCK6 (2.68 Å) and AutoDock Vina (3.92 Å). When starting with predicted models from I-TASSER, EDock still generates reasonable docking models, with a success rate 159% and 67% higher than DOCK6 and AutoDock Vina, respectively. Detailed data analyses show that the major advantage of EDock lies in reliable ligand binding site predictions and extensive REMC sampling, which allows for the implementation of multiple van der Waals weightings to accommodate different levels of steric clashes and cavity distortions and therefore enhances the robustness of low-resolution docking with predicted protein structures.
| S-EPMC7251717 | biostudies-literature