A generalized Fellner-Schall method for smoothing parameter estimation with application to Tweedie location, scale and shape models

We consider the estimation of smoothing parameters and variance components in models with a regular log likelihood subject to quadratic penalization of the model coefficients, via a generalization of the method of Fellner (1986) and Schall (1991). In particular: (i) we generalize the original method to the case of penalties that are linear in several smoothing parameters, thereby covering the important cases of tensor product and adaptive smoothers; (ii) we show why the method's steps increase the restricted marginal likelihood of the model, that it tends to converge faster than the EM algorithm, or obvious accelerations of this, and investigate its relation to Newton optimization; (iii) we generalize the method to any Fisher regular likelihood. The method represents a considerable simplification over existing methods of estimating smoothing parameters in the context of regular likelihoods, without sacrificing generality: for example, it is only necessary to compute with the same first and second derivatives of the log-likelihood required for coefficient estimation, and not with the third or fourth order derivatives required by alternative approaches. Examples are provided which would have been impossible or impractical with pre-existing Fellner-Schall methods, along with an example of a Tweedie location, scale and shape model which would be a challenge for alternative methods

Comparison of Material Properties and Microstructure of Specimens Built Using the 3D Systems Vanguard HS and Vanguard HiQ+HSSLS Systems

The HiQ upgrade to the 3D Systems Vanguard selective laser sintering (SLS) machine incorporates a revised thermal calibration system and new control software. The paper compares the tensile modulus, tensile strength, elongation at break, flexural modulus, Izod impact resistance and microstructure of two batteries of standard specimens built from recycled Duraform PA (Nylon 12). The first set is built on a Vanguard HS system and the second on the same system with the HiQ upgrade installed. The upgrade reduces user intervention, decreases total build time and improves surface finish. However, using the default processing parameters, tensile, flexure and impact properties are all found to decline after the upgrade is installed.Mechanical Engineerin

The Effects of Negative Legacies on the Adjustment of Parentally Bereaved Children and Adolescents

This is a report of a qualitative analysis of a sample of bereaved families in which one parent died and in which children scored in the clinical range on the Child Behavior Check List. The purpose of this analysis was to learn more about the lives of these children. They were considered to be at risk of developing emotional and behavioral problems associated with the death. We discovered that many of these “high risk” children had a continuing bond with the deceased that was primarily negative and troubling for them in contrast to a comparison group of children not at risk from the same study. Five types of legacies, not mutually exclusive, were identified: health related, role related, personal qualities, legacy of blame, and an emotional legacy. Coping behavior on the part of the surviving parent seemed to make a difference in whether or not a legacy was experienced as negative

Oral and Pharyngeal Cancer Control and Early Detection

Sixty-four standardized continuing education courses were given for dentists throughout the ten public health districts of the USA to determine if certain behaviors regarding oral and pharyngeal cancer (OPC) control could be modified. Questionnaires were obtained at baseline and at 6 months along with matched control groups. One thousand eight hundred two general dentists participated at baseline and 988 at a 6-month questionnaire follow-up. Analysis of the data indicated that continuing education courses had a positive influence on participants’ oral cancer attitudes, knowledge, and behavior that potentially could make a difference on prevention, early detection, and ultimately OPC control

Image Coaddition with Temporally Varying Kernels

Large, multi-frequency imaging surveys, such as the Large Synaptic Survey Telescope (LSST), need to do near-real time analysis of very large datasets. This raises a host of statistical and computational problems where standard methods do not work. In this paper, we study a proposed method for combining stacks of images into a single summary image, sometimes referred to as a template. This task is commonly referred to as image coaddition. In part, we focus on a method proposed in previous work, which outlines a procedure for combining stacks of images in an online fashion in the Fourier domain. We evaluate this method by comparing it to two straightforward methods through the use of various criteria and simulations. Note that the goal is not to propose these comparison methods for use in their own right, but to ensure that additional complexity also provides substantially improved performance

The Distribution of Nearby Stars in Velocity Space Inferred from Hipparcos Data

(abridged) The velocity distribution f(v) of nearby stars is estimated, via a maximum- likelihood algorithm, from the positions and tangential velocities of a kinematically unbiased sample of 14369 stars observed by the HIPPARCOS satellite. f(v) shows rich structure in the radial and azimuthal motions, v_R and v_phi, but not in the vertical velocity, v_z: there are four prominent and many smaller maxima, many of which correspond to well known moving groups. While samples of early-type stars are dominated by these maxima, also up to 25% of red main-sequence stars are associated with them. These moving groups are responsible for the vertex deviation measured even for samples of late-type stars; they appear more frequently for ever redder samples; and as a whole they follow an asymmetric-drift relation, in the sense that those only present in red samples predominantly have large |v_R| and lag in v_phi w.r.t. the local standard of rest (LSR). The question arise, how these old moving groups got on their eccentric orbits. A plausible mechanism, known from solar system dynamics, which is able to manage a shift in orbit space involves locking into an orbital resonance. Apart from these moving groups, there is a smooth background distribution, akin to Schwarzschild's ellipsoidal model, with axis ratio of about 1:0.6:0.35 in v_R, v_phi, and v_z. The contours are aligned with the $v_r$ direction, but not w.r.t. the v_phi and v_z axes: the mean v_z increases for stars rotating faster than the LSR. This effect can be explained by the stellar warp of the Galactic disk. If this explanation is correct, the warp's inner edge must not be within the solar circle, while its pattern rotates with frequency of about 13 km/s/kpc or more retrograde w.r.t. the stellar orbits.Comment: 16 pages LaTeX (aas2pp4.sty), 6 figures, accepted by A

Dissipation enhanced vibrational sensing in an olfactory molecular switch

Motivated by a proposed olfactory mechanism based on a vibrationally-activated molecular switch, we study electron transport within a donor-acceptor pair that is coupled to a vibrational mode and embedded in a surrounding environment. We derive a polaron master equation with which we study the dynamics of both the electronic and vibrational degrees of freedom beyond previously employed semiclassical (Marcus-Jortner) rate analyses. We show: (i) that in the absence of explicit dissipation of the vibrational mode, the semiclassical approach is generally unable to capture the dynamics predicted by our master equation due to both its assumption of one-way (exponential) electron transfer from donor to acceptor and its neglect of the spectral details of the environment; (ii) that by additionally allowing strong dissipation to act on the odorant vibrational mode we can recover exponential electron transfer, though typically at a rate that differs from that given by the Marcus-Jortner expression; (iii) that the ability of the molecular switch to discriminate between the presence and absence of the odorant, and its sensitivity to the odorant vibrational frequency, are enhanced significantly in this strong dissipation regime, when compared to the case without mode dissipation; and (iv) that details of the environment absent from previous Marcus-Jortner analyses can also dramatically alter the sensitivity of the molecular switch, in particular allowing its frequency resolution to be improved. Our results thus demonstrate the constructive role dissipation can play in facilitating sensitive and selective operation in molecular switch devices, as well as the inadequacy of semiclassical rate equations in analysing such behaviour over a wide range of parameters.Comment: 12 pages, 6 figures, close to published version, comments welcom

Nonlocal Phases of Local Quantum Mechanical Wavefunctions in Static and Time-Dependent Aharonov-Bohm Experiments

We show that the standard Dirac phase factor is not the only solution of the gauge transformation equations. The full form of a general gauge function (that connects systems that move in different sets of scalar and vector potentials), apart from Dirac phases also contains terms of classical fields that act nonlocally (in spacetime) on the local solutions of the time-dependent Schr\"odinger equation: the phases of wavefunctions in the Schr\"odinger picture are affected nonlocally by spatially and temporally remote magnetic and electric fields, in ways that are fully explored. These contributions go beyond the usual Aharonov-Bohm effects (magnetic or electric). (i) Application to cases of particles passing through static magnetic or electric fields leads to cancellations of Aharonov-Bohm phases at the observation point; these are linked to behaviors at the semiclassical level (to the old Werner & Brill experimental observations, or their "electric analogs" - or to recent reports of Batelaan & Tonomura) but are shown to be far more general (true not only for narrow wavepackets but also for completely delocalized quantum states). By using these cancellations, certain previously unnoticed sign-errors in the literature are corrected. (ii) Application to time-dependent situations provides a remedy for erroneous results in the literature (on improper uses of Dirac phase factors) and leads to phases that contain an Aharonov-Bohm part and a field-nonlocal part: their competition is shown to recover Relativistic Causality in earlier "paradoxes" (such as the van Kampen thought-experiment), while a more general consideration indicates that the temporal nonlocalities found here demonstrate in part a causal propagation of phases of quantum mechanical wavefunctions in the Schr\"odinger picture. This may open a direct way to address time-dependent double-slit experiments and the associated causal issuesComment: 49 pages, 1 figure, presented in Conferences "50 years of the Aharonov-Bohm effect and 25 years of the Berry's phase" (Tel Aviv and Bristol), published in Journ. Phys. A. Compared to the published paper, this version has 17 additional lines after eqn.(14) for maximum clarity, and the Abstract has been slightly modified and reduced from the published 2035 characters to the required 1920 character

Granger causality and transfer entropy are equivalent for Gaussian variables

Granger causality is a statistical notion of causal influence based on prediction via vector autoregression. Developed originally in the field of econometrics, it has since found application in a broader arena, particularly in neuroscience. More recently transfer entropy, an information-theoretic measure of time-directed information transfer between jointly dependent processes, has gained traction in a similarly wide field. While it has been recognized that the two concepts must be related, the exact relationship has until now not been formally described. Here we show that for Gaussian variables, Granger causality and transfer entropy are entirely equivalent, thus bridging autoregressive and information-theoretic approaches to data-driven causal inference.Comment: In review, Phys. Rev. Lett., Nov. 200