390 research outputs found
JointAI: Joint Analysis and Imputation of Incomplete Data in R
Missing data occur in many types of studies and typically complicate the
analysis. Multiple imputation, either using joint modelling or the more
flexible fully conditional specification approach, are popular and work well in
standard settings. In settings involving non-linear associations or
interactions, however, incompatibility of the imputation model with the
analysis model is an issue often resulting in bias. Similarly, complex outcomes
such as longitudinal or survival outcomes cannot be adequately handled by
standard implementations. In this paper, we introduce the R package JointAI,
which utilizes the Bayesian framework to perform simultaneous analysis and
imputation in regression models with incomplete covariates. Using a fully
Bayesian joint modelling approach it overcomes the issue of uncongeniality
while retaining the attractive flexibility of fully conditional specification
multiple imputation by specifying the joint distribution of analysis and
imputation models as a sequence of univariate models that can be adapted to the
type of variable. JointAI provides functions for Bayesian inference with
generalized linear and generalized linear mixed models and extensions thereof
as well as survival models and joint models for longitudinal and survival data,
that take arguments analogous to corresponding well known functions for the
analysis of complete data from base R and other packages. Usage and features of
JointAI are described and illustrated using various examples and the
theoretical background is outlined.Comment: imputation, Bayesian, missing covariates, non-linear, interaction,
multi-level, survival, joint model R, JAG
Carbon-poor stellar cores as supernova progenitors
Exploring stellar models which ignite carbon off-center (in the mass range of
about 1.05 - 1.25 Msun, depending on the carbon mass fraction) we find that
they may present an interesting SN I progenitor scenario, since whereas in the
standard scenario runaway always takes place at the same density of about 2 X
10^9 gr/cm^3, in our case, due to the small amount of carbon ignited, we get a
whole range of densities from 1 X 10^9 up to 6 X 10^9 gr/cm^3. These results
could contribute in resolving the emerging recognition that at least some
diversity among SNe I exists, since runaway at various central densities is
expected to yield various outcomes in terms of the velocities and composition
of the ejecta, which should be modeled and compared to observations.Comment: 49 pages, 20 figure
Summarising salient information on historical controls: A structured assessment of validity and comparability across studies
BACKGROUND: While placebo-controlled randomised controlled trials remain the standard way to evaluate drugs for efficacy, historical data are used extensively across the development cycle. This ranges from supplementing contemporary data to increase the power of trials to cross-trial comparisons in estimating comparative efficacy. In many cases, these approaches are performed without in-depth review of the context of data, which may lead to bias and incorrect conclusions. METHODS: We discuss the original 'Pocock' criteria for the use of historical data and how the use of historical data has evolved over time. Based on these factors and personal experience, we created a series of questions that may be asked of historical data, prior to their use. Based on the answers to these questions, various statistical approaches are recommended. The strategy is illustrated with a case study in colorectal cancer. RESULTS: A number of areas need to be considered with historical data, which we split into three categories: outcome measurement, study/patient characteristics (including setting and inclusion/exclusion criteria), and disease process/intervention effects. Each of these areas may introduce issues if not appropriately handled, while some may preclude the use of historical data entirely. We present a tool (in the form of a table) for highlighting any such issues. Application of the tool to a colorectal cancer data set demonstrates under what conditions historical data could be used and what the limitations of such an analysis would be. CONCLUSION: Historical data can be a powerful tool to augment or compare with contemporary trial data, though caution is required. We present some of the issues that may be considered when involving historical data and what (if any) statistical approaches may account for differences between studies. We recommend that, where historical data are to be used in analyses, potential differences between studies are addressed explicitly
Sulphur-bearing molecules in diffuse molecular clouds: new results from SOFIA/GREAT and the IRAM 30 m telescope
We have observed five sulphur-bearing molecules in foreground diffuse
molecular clouds lying along the sight-lines to five bright continuum sources.
We have used the GREAT instrument on SOFIA to observe the 1383 GHz transitions of SH towards the star-forming regions W31C,
G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards
all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico
Veleta to detect the HS 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions.
In nine foreground absorption components detected towards these sources, the
inferred column densities of the four detected molecules showed relatively
constant ratios, with N(SH)/N(HS) in the range 1.1 - 3.0, N(CS)/N(HS)
in the range 0.32 - 0.61, and N(SO)/N(HS) in the range 0.08 - 0.30. The
observed SH/H ratios - in the range (0.5-2.6) - indicate
that SH (and other sulphur-bearing molecules) account for << 1% of the
gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules,
however, greatly exceed those predicted by standard models of cold diffuse
molecular clouds, providing further evidence for the enhancement of endothermic
reaction rates by elevated temperatures or ion-neutral drift. We have
considered the observed abundance ratios in the context of shock and turbulent
dissipation region (TDR) models. Using the TDR model, we find that the
turbulent energy available at large scale in the diffuse ISM is sufficient to
explain the observed column densities of SH and CS. Standard shock and TDR
models, however, fail to reproduce the column densities of HS and SO by a
factor of about 10; more elaborate shock models - in which account is taken of
the velocity drift, relative to H, of SH molecules produced by the
dissociative recombination of HS - reduce this discrepancy to a factor
~ 3.Comment: 30 pages, accepted for publication in A&
Low-velocity shocks: signatures of turbulent dissipation in diffuse irradiated gas
Context. Large-scale motions in galaxies (supernovae explosions, galaxy collisions, galactic shear etc.) generate turbulence, which allows a fraction of the available kinetic energy to cascade down to small scales before it is dissipated.
Aims. We establish and quantify the diagnostics of turbulent dissipation in mildly irradiated diffuse gas in the specific context of shock structures.
Methods. We incorporated the basic physics of photon-dominated regions into a state-of-the-art steady-state shock code. We examined the chemical and emission properties of mildly irradiated (G_0 = 1) magnetised shocks in diffuse media (n_H = 10^2 to 10^4 cm^(-3)) at low- to moderate velocities (from 3 to 40 km s^(-1)).
Results. The formation of some molecules relies on endoergic reactions. Their abundances in J-type shocks are enhanced by several orders of magnitude for shock velocities as low as 7 km s^(-1). Otherwise most chemical properties of J-type shocks vary over less than an order of magnitude between velocities from about 7 to about 30 km s^(-1), where H_2 dissociation sets in. C-type shocks display a more gradual molecular enhancement with increasing shock velocity.
We quantified the energy flux budget (fluxes of kinetic, radiated and magnetic energies) with emphasis on the main cooling lines of the cold interstellar medium. Their sensitivity to shock velocity is such that it allows observations to constrain statistical distributions of shock velocities.
We fitted various probability distribution functions (PDFs) of shock velocities to spectroscopic observations of the galaxy-wide shock in Stephan’s Quintet and of a Galactic line of sight which samples diffuse molecular gas in Chamaeleon. In both cases, low velocities bear the greatest statistical weight and the PDF is consistent with a bimodal distribution. In the very low velocity shocks (below 5 km s^(-1)), dissipation is due to ion-neutral friction and it powers H_2 low-energy transitions and atomic lines. In moderate velocity shocks (20 km s^(-1) and above), the dissipation is due to viscous heating and accounts for most of the molecular emission. In our interpretation a significant fraction of the gas in the line of sight is shocked (from 4% to 66%). For example, C^+ emission may trace shocks in UV irradiated gas where C^+ is the dominant carbon species.
Conclusions. Low- and moderate velocity shocks are important in shaping the chemical composition and excitation state of the interstellar gas. This allows one to probe the statistical distribution of shock velocities in interstellar turbulence
Towards quantitative tissue absorption imaging by combining photoacoustics and acousto-optics
We propose a strategy for quantitative photoacoustic mapping of chromophore
concentrations that can be performed purely experimentally. We exploit the
possibility of acousto-optic modulation using focused ultrasound, and the
principle that photons follow trajectories through a turbid medium in two
directions with equal probability. A theory is presented that expresses the
local absorption coefficient inside a medium in terms of noninvasively measured
quantities and experimental parameters. Proof of the validity of the theory is
given with Monte Carlo simulations.Comment: 14 pages, 5 figure
Detection of the tagged or untagged photons in acousto-optic imaging of thick highly scattering media by photorefractive adaptive holography
We propose an original adaptive wavefront holographic setup based on the
photorefractive effect (PR), to make real-time measurements of acousto-optic
signals in thick scattering media, with a high flux collection at high rates
for breast tumor detection. We describe here our present state of art and
understanding on the problem of breast imaging with PR detection of the
acousto-optic signal
- …