Estimating Nuisance Parameters in Inverse Problems

Many inverse problems include nuisance parameters which, while not of direct interest, are required to recover primary parameters. Structure present in these problems allows efficient optimization strategies - a well known example is variable projection, where nonlinear least squares problems which are linear in some parameters can be very efficiently optimized. In this paper, we extend the idea of projecting out a subset over the variables to a broad class of maximum likelihood (ML) and maximum a posteriori likelihood (MAP) problems with nuisance parameters, such as variance or degrees of freedom. As a result, we are able to incorporate nuisance parameter estimation into large-scale constrained and unconstrained inverse problem formulations. We apply the approach to a variety of problems, including estimation of unknown variance parameters in the Gaussian model, degree of freedom (d.o.f.) parameter estimation in the context of robust inverse problems, automatic calibration, and optimal experimental design. Using numerical examples, we demonstrate improvement in recovery of primary parameters for several large- scale inverse problems. The proposed approach is compatible with a wide variety of algorithms and formulations, and its implementation requires only minor modifications to existing algorithms.Comment: 16 pages, 5 figure

Chinese-style pragmatic objectivity in war reporting

China plays an increasing role in the wars and conflicts around the world with its expanding political and economic interests overseas, and its diplomatic role in international affairs. More and more Chinese journalists go to the frontlines overseas to cover distant conflicts for domestic audiences. Based on semi-structured in-depth interviews with sixteen Chinese correspondents who have covered conflicts outside China, this study examines Chinese journalists’ perceptions and reflections on objectivity in the war zones. The author adopts a term of Chinese-style pragmatic objectivity to mean that objectivity is a convenient approach for Chinese journalists to do war journalism in the field. At the level of objectivity-as-a- value, objectivity is defined as a pragmatic value and a practical ritual for Chinese journalists to do news within the scope they can reach, to protect themselves from criticisms, and to justify their version of the truth. It promotes allegiance and patriotism. At the level of objectivity-as-a-practice, objectivity in war coverage is compromised by China’s foreign policies, military constraints, the press’s political orientations and editorial polices, and journalists’ personal experiences and values. Chinese journalists use Chinese-style objectivity to negotiate their roles in the power struggle with the state, foreign militaries, the newsroom and journalists

The milieu of accountability of early companies in the Qing Dynasty: evidence from the Shànghai-based print media.

This study looks at evidence from the Shànghǎi (Shanghai)-based print media of the accountability of early Chinese companies from the middle period of the Qīng (Qing) Dynasty when the Opium War broke out in 1840 until the imperial monarchy's overthrow in the revolution of 1911. The Qīng Dynasty is known for its technical accounting stagnation. Yet, an examination of the Shànghǎi-based print media shows the existence of a strong sense of public reporting by early companies of the Qīng Dynasty. The findings of this study indicate that the print media displayed a rich milieu of accountability of these early companies by incorporating components of Western and Chinese benchmarks of accountability expressed in commercial metrics and key financial ratios

RNAalifold: improved consensus structure prediction for RNA alignments

<p>Abstract</p> <p>Background</p> <p>The prediction of a consensus structure for a set of related RNAs is an important first step for subsequent analyses. RNAalifold, which computes the minimum energy structure that is simultaneously formed by a set of aligned sequences, is one of the oldest and most widely used tools for this task. In recent years, several alternative approaches have been advocated, pointing to several shortcomings of the original RNAalifold approach.</p> <p>Results</p> <p>We show that the accuracy of RNAalifold predictions can be improved substantially by introducing a different, more rational handling of alignment gaps, and by replacing the rather simplistic model of covariance scoring with more sophisticated RIBOSUM-like scoring matrices. These improvements are achieved without compromising the computational efficiency of the algorithm. We show here that the new version of RNAalifold not only outperforms the old one, but also several other tools recently developed, on different datasets.</p> <p>Conclusion</p> <p>The new version of RNAalifold not only can replace the old one for almost any application but it is also competitive with other approaches including those based on SCFGs, maximum expected accuracy, or hierarchical nearest neighbor classifiers.</p

The VLA-COSMOS 3 GHz Large Project: continuum data and source catalog release

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys

Energetic eruptions leading to a peculiar hydrogen-rich explosion of a massive star

Every supernova so far observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower-moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining1. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the progenitor star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95 to 130 solar masses, which experience the pulsational pair instability2,3,4,5. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required

The Deepest Radio Observations of Nearby Type IA Supernovae: Constraining Progenitor Types and Optimizing Future Surveys

We report deep radio observations of nearby Type Ia Supernovae (SNe Ia) with the electronic Multi-Element Radio Linked Interferometer Net-work (e-MERLIN), and the Australia Telescope Compact Array (ATCA). No detections were made. With standard assumptions for the energy densities of relativistic electrons going into a power-law energy distribution, and the magnetic field strength (epsilon_e = epsilon_B = 0.1), we arrive at the upper limits on mass-loss rate for the progenitor system of SN 2013dy (2016coj, 2018gv, 2018pv, 2019np), to be less than 12 (2.8,1.3, 2.1, 1.7)EE(-8) solar masses per year (for a wind speed of 100 km/s). To SNe 2016coj, 2018gv, 2018pv and 2019np we add radio data for 17 other nearby SNe Ia, and model their non-detections. With the same model as described, all 21 SNe Ia have mass-loss rates less than 4EE(-8) solar masses per year (for a wind speed of 100 km/s). We compare those limits with the expected mass loss rates in different single-degenerate progenitor scenarios. We also discuss how information on epsilon_e and epsilon_B can be obtained from late observations of SNe Ia and the youngest SN Ia remnant detected in radio, G1.9+0.3, as well as stripped-envelope core-collapse SNe. We highlight SN 2011dh, and argue for epsilon_e approximately equal to 0.1 and epsilon_B approximately equal to 0.0033. Finally, we discuss strategies to observe at radio frequencies to maximize the chance of detection, given the time since explosion, the distance to the supernova and the telescope sensitivity