-
Particle acceleration and magnetic field amplification by relativistic shocks in inhomogeneous media
Authors:
Kanji Morikawa,
Yutaka Ohira,
Takumi Ohmura
Abstract:
Particle acceleration and magnetic field amplification in relativistic shocks propagating in inhomogeneous media are investigated by three-dimensional magnetohydrodynamical (MHD) simulations and test-particle simulations. The MHD simulations show that the interaction between the relativistic shock and dense clumps amplifies the downstream magnetic field to the value expected from observations of t…
▽ More
Particle acceleration and magnetic field amplification in relativistic shocks propagating in inhomogeneous media are investigated by three-dimensional magnetohydrodynamical (MHD) simulations and test-particle simulations. The MHD simulations show that the interaction between the relativistic shock and dense clumps amplifies the downstream magnetic field to the value expected from observations of the gamma-ray burst. The test-particle simulations in the electromagnetic field given by the MHD simulation show that particles are accelerated by the downstream turbulence and the relativistic shock. We provide the injection energy to the shock acceleration in this system. If the amplitude of upstream density fluctuations is sufficiently large, low-energy particles are initially accelerated to the injection energy by the downstream turbulence and then rapidly accelerated to higher energies by the relativistic shock. Therefore, the density fluctuation significantly affects particle acceleration in the relativistic shock.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Efficient Multiple-Robust Estimation for Nonresponse Data Under Informative Sampling
Authors:
Kosuke Morikawa,
Kenji Beppu,
Wataru Aida
Abstract:
Nonresponse after probability sampling is a universal challenge in survey sampling, often necessitating adjustments to mitigate sampling and selection bias simultaneously. This study explored the removal of bias and effective utilization of available information, not just in nonresponse but also in the scenario of data integration, where summary statistics from other data sources are accessible. W…
▽ More
Nonresponse after probability sampling is a universal challenge in survey sampling, often necessitating adjustments to mitigate sampling and selection bias simultaneously. This study explored the removal of bias and effective utilization of available information, not just in nonresponse but also in the scenario of data integration, where summary statistics from other data sources are accessible. We reformulate these settings within a two-step monotone missing data framework, where the first step of missingness arises from sampling and the second originates from nonresponse. Subsequently, we derive the semiparametric efficiency bound for the target parameter. We also propose adaptive estimators utilizing methods of moments and empirical likelihood approaches to attain the lower bound. The proposed estimator exhibits both efficiency and double robustness. However, attaining efficiency with an adaptive estimator requires the correct specification of certain working models. To reinforce robustness against the misspecification of working models, we extend the property of double robustness to multiple robustness by proposing a two-step empirical likelihood method that effectively leverages empirical weights. A numerical study is undertaken to investigate the finite-sample performance of the proposed methods. We further applied our methods to a dataset from the National Health and Nutrition Examination Survey data by efficiently incorporating summary statistics from the National Health Interview Survey data.
△ Less
Submitted 11 November, 2023;
originally announced November 2023.
-
Verifiable identification condition for nonignorable nonresponse data with categorical instrumental variables
Authors:
Kenji Beppu,
Kosuke Morikawa
Abstract:
We consider a model identification problem in which an outcome variable contains nonignorable missing values. Statistical inference requires a guarantee of the model identifiability to obtain estimators enjoying theoretically reasonable properties such as consistency and asymptotic normality. Recently, instrumental or shadow variables, combined with the completeness condition in the outcome model,…
▽ More
We consider a model identification problem in which an outcome variable contains nonignorable missing values. Statistical inference requires a guarantee of the model identifiability to obtain estimators enjoying theoretically reasonable properties such as consistency and asymptotic normality. Recently, instrumental or shadow variables, combined with the completeness condition in the outcome model, have been highlighted to make a model identifiable. However, the completeness condition may not hold even for simple models when the instrument is categorical. We propose a sufficient condition for model identifiability, which is applicable to cases where establishing the completeness condition is difficult. Using observed data, we demonstrate that the proposed conditions are easy to check for many practical models and outline their usefulness in numerical experiments and real data analysis.
△ Less
Submitted 3 July, 2023; v1 submitted 5 April, 2023;
originally announced April 2023.
-
An empirical likelihood approach to reduce selection bias in voluntary samples
Authors:
Jae Kwang Kim,
Kosuke Morikawa
Abstract:
We address the weighting problem in voluntary samples under a nonignorable sample selection model. Under the assumption that the sample selection model is correctly specified, we can compute a consistent estimator of the model parameter and construct the propensity score estimator of the population mean. We use the empirical likelihood method to construct the final weights for voluntary samples by…
▽ More
We address the weighting problem in voluntary samples under a nonignorable sample selection model. Under the assumption that the sample selection model is correctly specified, we can compute a consistent estimator of the model parameter and construct the propensity score estimator of the population mean. We use the empirical likelihood method to construct the final weights for voluntary samples by incorporating the bias calibration constraints and the benchmarking constraints. Linearization variance estimation of the proposed method is developed. A limited simulation study is also performed to check the performance of the proposed methods.
△ Less
Submitted 11 May, 2023; v1 submitted 5 November, 2022;
originally announced November 2022.
-
Semiparametric imputation using latent sparse conditional Gaussian mixtures for multivariate mixed outcomes
Authors:
Shonosuke Sugasawa,
Jae Kwang Kim,
Kosuke Morikawa
Abstract:
This paper proposes a flexible Bayesian approach to multiple imputation using conditional Gaussian mixtures. We introduce novel shrinkage priors for covariate-dependent mixing proportions in the mixture models to automatically select the suitable number of components used in the imputation step. We develop an efficient sampling algorithm for posterior computation and multiple imputation via Markov…
▽ More
This paper proposes a flexible Bayesian approach to multiple imputation using conditional Gaussian mixtures. We introduce novel shrinkage priors for covariate-dependent mixing proportions in the mixture models to automatically select the suitable number of components used in the imputation step. We develop an efficient sampling algorithm for posterior computation and multiple imputation via Markov Chain Monte Carlo methods. The proposed method can be easily extended to the situation where the data contains not only continuous variables but also discrete variables such as binary and count values. We also propose approximate Bayesian inference for parameters defined by loss functions based on posterior predictive distributing of missing observations, by extending bootstrap-based Bayesian inference for complete data. The proposed method is demonstrated through numerical studies using simulated and real data.
△ Less
Submitted 16 August, 2022;
originally announced August 2022.
-
Semiparametric adaptive estimation under informative sampling
Authors:
Kosuke Morikawa,
Yoshikazu Terada,
Jae Kwang Kim
Abstract:
In survey sampling, survey data do not necessarily represent the target population, and the samples are often biased. However, information on the survey weights aids in the elimination of selection bias. The Horvitz-Thompson estimator is a well-known unbiased, consistent, and asymptotically normal estimator; however, it is not efficient. Thus, this study derives the semiparametric efficiency bound…
▽ More
In survey sampling, survey data do not necessarily represent the target population, and the samples are often biased. However, information on the survey weights aids in the elimination of selection bias. The Horvitz-Thompson estimator is a well-known unbiased, consistent, and asymptotically normal estimator; however, it is not efficient. Thus, this study derives the semiparametric efficiency bound for various target parameters by considering the survey weight as a random variable and consequently proposes a semiparametric optimal estimator with certain working models on the survey weights. The proposed estimator is consistent, asymptotically normal, and efficient in a class of the regular and asymptotically linear estimators. Further, a limited simulation study is conducted to investigate the finite sample performance of the proposed method. The proposed method is applied to the 1999 Canadian Workplace and Employee Survey data.
△ Less
Submitted 3 April, 2024; v1 submitted 11 August, 2022;
originally announced August 2022.
-
Identification enhanced generalised linear model estimation with nonignorable missing outcomes
Authors:
Kenji Beppu,
Jinung Choi,
Kosuke Morikawa,
Jongho Im
Abstract:
Missing data often result in undesirable bias and loss of efficiency. These become substantial problems when the response mechanism is nonignorable, such that the response model depends on unobserved variables. It is necessary to estimate the joint distribution of unobserved variables and response indicators to manage nonignorable nonresponse. However, model misspecification and identification iss…
▽ More
Missing data often result in undesirable bias and loss of efficiency. These become substantial problems when the response mechanism is nonignorable, such that the response model depends on unobserved variables. It is necessary to estimate the joint distribution of unobserved variables and response indicators to manage nonignorable nonresponse. However, model misspecification and identification issues prevent robust estimates despite careful estimation of the target joint distribution. In this study, we modelled the distribution of the observed parts and derived sufficient conditions for model identifiability, assuming a logistic regression model as the response mechanism and generalised linear models as the main outcome model of interest. More importantly, the derived sufficient conditions are testable with the observed data and do not require any instrumental variables, which are often assumed to guarantee model identifiability but cannot be practically determined beforehand. To analyse missing data, we propose a new imputation method which incorporates verifiable identifiability using only observed data. Furthermore, we present the performance of the proposed estimators in numerical studies and apply the proposed method to two sets of real data: exit polls for the 19th South Korean election data and public data collected from the Korean Survey of Household Finances and Living Conditions.
△ Less
Submitted 20 February, 2024; v1 submitted 22 April, 2022;
originally announced April 2022.
-
Adjusting for publication bias in meta-analysis via inverse probability weighting using clinical trial registries
Authors:
Ao Huang,
Kosuke Morikawa,
Tim Friede,
Satoshi Hattori
Abstract:
Publication bias is a major concern in conducting systematic reviews and meta-analyses. Various sensitivity analysis or bias-correction methods have been developed based on selection models and they have some advantages over the widely used bias-correction method of the trim-and-fill method. However, likelihood methods based on selection models may have difficulty in obtaining precise estimates an…
▽ More
Publication bias is a major concern in conducting systematic reviews and meta-analyses. Various sensitivity analysis or bias-correction methods have been developed based on selection models and they have some advantages over the widely used bias-correction method of the trim-and-fill method. However, likelihood methods based on selection models may have difficulty in obtaining precise estimates and reasonable confidence intervals or require a complicated sensitivity analysis process. In this paper, we develop a simple publication bias adjustment method utilizing information on conducted but still unpublished trials from clinical trial registries. We introduce an estimating equation for parameter estimation in the selection function by regarding the publication bias issue as a missing data problem under missing not at random. With the estimated selection function, we introduce the inverse probability weighting (IPW) method to estimate the overall mean across studies. Furthermore, the IPW versions of heterogeneity measures such as the between-study variance and the I2 measure are proposed. We propose methods to construct asymptotic confidence intervals and suggest intervals based on parametric bootstrapping as an alternative. Through numerical experiments, we observed that the estimators successfully eliminate biases and the confidence intervals had empirical coverage probabilities close to the nominal level. On the other hand, the asymptotic confidence interval is much wider in some scenarios than the bootstrap confidence interval. Therefore, the latter is recommended for practical use.
△ Less
Submitted 26 September, 2021;
originally announced September 2021.
-
Precise Neutron Lifetime Measurement Using Pulsed Neutron Beams at J-PARC
Authors:
N. Sumi,
K. Hirota,
G. Ichikawa,
T. Ino,
Y. Iwashita,
S. Kajiwara,
Y. Kato,
M. Kitaguchi,
K. Mishima,
K. Morikawa,
T. Mogi,
H. Oide,
H. Okabe,
H. Otono,
T. Shima,
H. M. Shimizu,
Y. Sugisawa,
T. Tanabe,
S. Yamashita,
K. Yano,
T. Yoshioka
Abstract:
A neutron decays into a proton, an electron, and an anti-neutrino through the beta-decay process. The decay lifetime ($\sim$880 s) is an important parameter in the weak interaction. For example, the neutron lifetime is a parameter used to determine the |$V_{\rm ud}$| parameter of the CKM quark mixing matrix. The lifetime is also one of the input parameters for the Big Bang Nucleosynthesis, which p…
▽ More
A neutron decays into a proton, an electron, and an anti-neutrino through the beta-decay process. The decay lifetime ($\sim$880 s) is an important parameter in the weak interaction. For example, the neutron lifetime is a parameter used to determine the |$V_{\rm ud}$| parameter of the CKM quark mixing matrix. The lifetime is also one of the input parameters for the Big Bang Nucleosynthesis, which predicts light element synthesis in the early universe. However, experimental measurements of the neutron lifetime today are significantly different (8.4 s or 4.0$σ$) depending on the methods. One is a bottle method measuring surviving neutron in the neutron storage bottle. The other is a beam method measuring neutron beam flux and neutron decay rate in the detector. There is a discussion that the discrepancy comes from unconsidered systematic error or undetectable decay mode, such as dark decay. A new type of beam experiment is performed at the BL05 MLF J-PARC. This experiment measured neutron flux and decay rate simultaneously with a time projection chamber using a pulsed neutron beam. We will present the world situation of neutron lifetime and the latest results at J-PARC.
△ Less
Submitted 19 February, 2021;
originally announced February 2021.
-
A Linked Aggregate Code for Processing Faces (Revised Version)
Authors:
Michael Lyons,
Kazunori Morikawa
Abstract:
A model of face representation, inspired by the biology of the visual system, is compared to experimental data on the perception of facial similarity. The face representation model uses aggregate primary visual cortex (V1) cell responses topographically linked to a grid covering the face, allowing comparison of shape and texture at corresponding points in two facial images. When a set of relativel…
▽ More
A model of face representation, inspired by the biology of the visual system, is compared to experimental data on the perception of facial similarity. The face representation model uses aggregate primary visual cortex (V1) cell responses topographically linked to a grid covering the face, allowing comparison of shape and texture at corresponding points in two facial images. When a set of relatively similar faces was used as stimuli, this Linked Aggregate Code (LAC) predicted human performance in similarity judgment experiments. When faces of perceivable categories were used, dimensions such as apparent sex and race emerged from the LAC model without training. The dimensional structure of the LAC similarity measure for the mixed category task displayed some psychologically plausible features but also highlighted differences between the model and the human similarity judgements. The human judgements exhibited a racial perceptual bias that was not shared by the LAC model. The results suggest that the LAC based similarity measure may offer a fertile starting point for further modelling studies of face representation in higher visual areas, including studies of the development of biases in face perception.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
Foundations of Temperature Theory
Authors:
Katsuyuki Bando,
Eitetsu Ken,
Kota Morikawa
Abstract:
Temperature of combinatorial games have been long studied since when Conway established the modern combinatorial game theory, and there are several variations of the concepts. In this article, we focus on one of the classical versions of temperature, and give precise proofs to the fundamental claims on temperature, namely, the existence, order-preservation, and homomorphism. Besides, a general rel…
▽ More
Temperature of combinatorial games have been long studied since when Conway established the modern combinatorial game theory, and there are several variations of the concepts. In this article, we focus on one of the classical versions of temperature, and give precise proofs to the fundamental claims on temperature, namely, the existence, order-preservation, and homomorphism. Besides, a general relationship between the value and the thermograph of a game is investigated.
△ Less
Submitted 22 November, 2020; v1 submitted 6 September, 2020;
originally announced September 2020.
-
Neutron lifetime measurement with pulsed cold neutrons
Authors:
K. Hirota,
G. Ichikawa,
S. Ieki,
T. Ino,
Y. Iwashita,
M. Kitaguchi,
R. Kitahara,
J. Koga,
K. Mishima,
T. Mogi,
K. Morikawa,
A. Morishita,
N. Nagakura,
H. Oide,
H. Okabe,
H. Otono,
Y. Seki,
D. Sekiba,
T. Shima,
H. M. Shimizu,
N. Sumi,
H. Sumino,
T. Tomita,
H. Uehara,
T. Yamada
, et al. (4 additional authors not shown)
Abstract:
The neutron lifetime has been measured by comparing the decay rate with the reaction rate of $^3$He nuclei of a pulsed neutron beam from the spallation neutron source at the Japan Proton Accelerator Research Complex (J-PARC). The decay rate and the reaction rate were determined by simultaneously detecting electrons from the neutron decay and protons from the $^3$He(n,p)$^3$H reaction using a gas c…
▽ More
The neutron lifetime has been measured by comparing the decay rate with the reaction rate of $^3$He nuclei of a pulsed neutron beam from the spallation neutron source at the Japan Proton Accelerator Research Complex (J-PARC). The decay rate and the reaction rate were determined by simultaneously detecting electrons from the neutron decay and protons from the $^3$He(n,p)$^3$H reaction using a gas chamber of which working gas contains diluted $^3$He. The measured neutron lifetime was $898\,\pm\,10\,_{\rm stat}\,^{+15}_{-18}\,_{\rm sys}\,$s.
△ Less
Submitted 25 November, 2020; v1 submitted 22 July, 2020;
originally announced July 2020.
-
Forecasting temporal variation of aftershocks immediately after a main shock using Gaussian process regression
Authors:
Kosuke Morikawa,
Hiromichi Nagao,
Shin-ichi Ito,
Yoshikazu Terada,
Shin'ichi Sakai,
Naoshi Hirata
Abstract:
Uncovering the distribution of magnitudes and arrival times of aftershocks is a key to comprehend the characteristics of the sequence of earthquakes, which enables us to predict seismic activities and hazard assessments. However, identifying the number of aftershocks immediately after the main shock is practically difficult due to contaminations of arriving seismic waves. To overcome the difficult…
▽ More
Uncovering the distribution of magnitudes and arrival times of aftershocks is a key to comprehend the characteristics of the sequence of earthquakes, which enables us to predict seismic activities and hazard assessments. However, identifying the number of aftershocks immediately after the main shock is practically difficult due to contaminations of arriving seismic waves. To overcome the difficulty, we construct a likelihood based on the detected data incorporating a detection function to which the Gaussian process regression (GPR) is applied. The GPR is capable of estimating not only the parameters of the distribution of aftershocks together with the detection function but also credible intervals for both of the parameters and the detection function. A property that distributions of both the Gaussian process and aftershocks are exponential functions leads to an efficient Bayesian computational algorithm to estimate the hyperparameters. After the validations through numerical tests, the proposed method is retrospectively applied to the catalog data related to the 2004 Chuetsu earthquake towards early forecasting of the aftershocks. The result shows that the proposed method stably estimates the parameters of the distribution simultaneously their credible intervals even within three hours after the main shock.
△ Less
Submitted 3 February, 2021; v1 submitted 14 June, 2020;
originally announced June 2020.
-
Study of Deep Generative Models for Inorganic Chemical Compositions
Authors:
Yoshihide Sawada,
Koji Morikawa,
Mikiya Fujii
Abstract:
Generative models based on generative adversarial networks (GANs) and variational autoencoders (VAEs) have been widely studied in the fields of image generation, speech generation, and drug discovery, but, only a few studies have focused on the generation of inorganic materials. Such studies use the crystal structures of materials, but material researchers rarely store this information. Thus, we g…
▽ More
Generative models based on generative adversarial networks (GANs) and variational autoencoders (VAEs) have been widely studied in the fields of image generation, speech generation, and drug discovery, but, only a few studies have focused on the generation of inorganic materials. Such studies use the crystal structures of materials, but material researchers rarely store this information. Thus, we generate chemical compositions without using crystal information. We use a conditional VAE (CondVAE) and a conditional GAN (CondGAN) and show that CondGAN using the bag-of-atom representation with physical descriptors generates better compositions than other generative models. Also, we evaluate the effectiveness of the Metropolis-Hastings-based atomic valency modification and the extrapolation performance, which is important to material discovery.
△ Less
Submitted 24 October, 2019;
originally announced October 2019.
-
Bayesian Semiparametric Modeling of Response Mechanism for Nonignorable Missing Data
Authors:
Shonosuke Sugasawa,
Kosuke Morikawa,
Keisuke Takahata
Abstract:
Statistical inference with nonresponse is quite challenging, especially when the response mechanism is nonignorable. In this case, the validity of statistical inference depends on untestable correct specification of the response model. To avoid the misspecification, we propose semiparametric Bayesian estimation in which an outcome model is parametric, but the response model is semiparametric in th…
▽ More
Statistical inference with nonresponse is quite challenging, especially when the response mechanism is nonignorable. In this case, the validity of statistical inference depends on untestable correct specification of the response model. To avoid the misspecification, we propose semiparametric Bayesian estimation in which an outcome model is parametric, but the response model is semiparametric in that we do not assume any parametric form for the nonresponse variable. We adopt penalized spline methods to estimate the unknown function. We also consider a fully nonparametric approach to modeling the response mechanism by using radial basis function methods. Using Polya-gamma data augmentation, we developed an efficient posterior computation algorithm via Gibbs sampling in which most full conditional distributions can be obtained in familiar forms. The performance of the proposed method is demonstrated in simulation studies and an application to longitudinal data.
△ Less
Submitted 14 January, 2021; v1 submitted 6 September, 2019;
originally announced September 2019.
-
A Profile Likelihood Approach to Semiparametric Estimation with Nonignorable Nonresponse
Authors:
Hejian Sang,
Kosuke Morikawa
Abstract:
Statistical inference with nonresponse is quite challenging, especially when the response mechanism is nonignorable. The existing methods often require correct model specifications for both outcome and response models. However, due to nonresponse, both models cannot be verified from data directly and model misspecification can lead to a seriously biased inference. To overcome this limitation, we d…
▽ More
Statistical inference with nonresponse is quite challenging, especially when the response mechanism is nonignorable. The existing methods often require correct model specifications for both outcome and response models. However, due to nonresponse, both models cannot be verified from data directly and model misspecification can lead to a seriously biased inference. To overcome this limitation, we develop a robust and efficient semiparametric method based on the profile likelihood. The proposed method uses the robust semiparametric response model, in which fully unspecified function of study variable is assumed. An efficient computation algorithm using fractional imputation is developed. A quasi-likelihood approach for testing ignorability is also developed. The consistency and asymptotic normality of the proposed method are established. The finite-sample performance is examined in the extensive simulation studies and an application to the Korean Labor and Income Panel Study dataset is also presented.
△ Less
Submitted 10 September, 2018;
originally announced September 2018.
-
Semiparametric Optimal Estimation With Nonignorable Nonresponse Data
Authors:
Kosuke Morikawa,
Jae Kwang Kim
Abstract:
When the response mechanism is believed to be not missing at random (NMAR), a valid analysis requires stronger assumptions on the response mechanism than standard statistical methods would otherwise require. Semiparametric estimators have been developed under the model assumptions on the response mechanism. In this paper, a new statistical test is proposed to guarantee model identifiability withou…
▽ More
When the response mechanism is believed to be not missing at random (NMAR), a valid analysis requires stronger assumptions on the response mechanism than standard statistical methods would otherwise require. Semiparametric estimators have been developed under the model assumptions on the response mechanism. In this paper, a new statistical test is proposed to guarantee model identifiability without using any instrumental variable. Furthermore, we develop optimal semiparametric estimation for parameters such as the population mean. Specifically, we propose two semiparametric optimal estimators that do not require any model assumptions other than the response mechanism. Asymptotic properties of the proposed estimators are discussed. An extensive simulation study is presented to compare with some existing methods. We present an application of our method using Korean Labor and Income Panel Survey data.
△ Less
Submitted 7 May, 2020; v1 submitted 29 December, 2016;
originally announced December 2016.
-
Statistical Inference with Different Missing-data Mechanisms
Authors:
Kosuke Morikawa,
Yutaka Kano
Abstract:
When data are missing due to at most one cause from some time to next time, we can make sampling distribution inferences about the parameter of the data by modeling the missing-data mechanism correctly. Proverbially, in case its mechanism is missing at random (MAR), it can be ignored, but in case not missing at random (NMAR), it can not be. There are no methods, however, to analyze when missing of…
▽ More
When data are missing due to at most one cause from some time to next time, we can make sampling distribution inferences about the parameter of the data by modeling the missing-data mechanism correctly. Proverbially, in case its mechanism is missing at random (MAR), it can be ignored, but in case not missing at random (NMAR), it can not be. There are no methods, however, to analyze when missing of the data can occur because of several causes despite of there being many such data in practice. Hence the aim of this paper is to propose how to inference on such data. Concretely, we extend the missing-data indicator from usual binary random vectors to discrete random vectors, define missing-data mechanism for every causes and research ignorability of a mixture of missing-data mechanisms such as "MAR & MAR" and "MAR & NMAR". In particular, when the combination of mechanisms is "MAR & NMAR", generally the component of MAR can not be ignored, but in special case, it can be.
△ Less
Submitted 18 July, 2014;
originally announced July 2014.
-
Identification Problem for The Analysis of Binary Data with Non-ignorable Missing
Authors:
Kosuke Morikawa,
Yutaka Kano
Abstract:
When a missing-data mechanism is NMAR or non-ignorable, missingness is itself vital information and it must be taken into the likelihood, which, however, needs to introduce additional parameters to be estimated. The incompleteness of the data and introduction of more parameters can cause the identification problem. When a response variable is binary, it becomes a more serious problem because of le…
▽ More
When a missing-data mechanism is NMAR or non-ignorable, missingness is itself vital information and it must be taken into the likelihood, which, however, needs to introduce additional parameters to be estimated. The incompleteness of the data and introduction of more parameters can cause the identification problem. When a response variable is binary, it becomes a more serious problem because of less information of bi- nary data, however, there are no methods to briefly verify whether a mode is identified or not. Therefore, we provide a new necessary and sufficient condition to easily check model identifiability when analyzing binary data with non-ignorable missing by condi- tional models. This condition can give us what condition is needed for a model to have identifiability as well as make easily check the identifiability of a model.
△ Less
Submitted 14 May, 2014;
originally announced May 2014.
-
Long-term monitoring of the short period SU UMa-type dwarf nova, V844 Herculis
Authors:
Shota Oizumi,
Toshihiro Omodaka,
Hiroyuki Yamamoto,
Shunsuke Tanada,
Tatsuki Yasuda,
Yoshihiro Arao,
Kie Kodama,
Miho Suzuki,
Takeshi Matsuo,
Hiroyuki Maehara,
Kazuhiro Nakajima,
Pavol A. Dubovsky,
Taichi Kato,
Akira Imada,
Kaori Kubota,
Kei Sugiyasu,
Koichi Morikawa,
Ken'ichi Torii,
Makoto Uemura,
Ryoko Ishioka,
Kenji Tanabe,
Daisaku Nogami
Abstract:
We report on time-resolved CCD photometry of four outbursts of a short-period SU UMa-type dwarf nova, V844 Herculis. We successfully determined the mean superhump periods to be 0.05584(64) days, and 0.055883(3) for the 2002 May superoutburst, and the 2006 April-May superoutburst, respectively. During the 2002 October observations, we confirmed that the outburst is a normal outburst, which is the…
▽ More
We report on time-resolved CCD photometry of four outbursts of a short-period SU UMa-type dwarf nova, V844 Herculis. We successfully determined the mean superhump periods to be 0.05584(64) days, and 0.055883(3) for the 2002 May superoutburst, and the 2006 April-May superoutburst, respectively. During the 2002 October observations, we confirmed that the outburst is a normal outburst, which is the first recorded normal outburst in V844 Her. We also examined superhump period changes during 2002 May and 2006 April-May superoutbursts, both of which showed increasing superhump period over the course of the plateau stage. In order to examine the long-term behavior of V844 Her, we analyzed archival data over the past ten years since the discovery of this binary. Although photometry is not satisfactory in some superoutbursts, we found that V844 Her showed no precursors and rebrightenings. Based on the long-term light curve, we further confirmed V844 Her has shown almost no normal outbursts despite the fact that the supercycle of the system is estimated to be about 300 days. In order to explain the long-term light curves of V844 Her, evaporation in the accretion disk may play a role in the avoidance of several normal outbursts, which does not contradict with the relatively large X-ray luminosity of V844 Her.
△ Less
Submitted 28 February, 2007;
originally announced February 2007.
-
A New SU UMa-Type Dwarf Nova, QW Serpentis (= TmzV46)
Authors:
D. Nogami,
M. Uemura,
R. Ishioka,
H. Iwamatsu,
T. Kato,
E. P. Pavlenko,
A. Baklanov,
R. Novak,
S. Kiyota,
K. Tanabe,
G. Masi,
L. M. Cook,
K. Morikawa,
P. Schmeer
Abstract:
We report on the results of the QW Ser campaign which has been continued from 2000 to 2003 by the VSNET collaboration team. Four long outbursts and many short ones were caught during this period. Our intensive photometric observations revealed superhumps with a period of 0.07700(4) d during all four superoutbursts, proving the SU UMa nature of this star. The recurrence cycles of the normal outbu…
▽ More
We report on the results of the QW Ser campaign which has been continued from 2000 to 2003 by the VSNET collaboration team. Four long outbursts and many short ones were caught during this period. Our intensive photometric observations revealed superhumps with a period of 0.07700(4) d during all four superoutbursts, proving the SU UMa nature of this star. The recurrence cycles of the normal outbursts and the superoutbursts were measured to be $\sim$50 days and 240(30) days, respectively. The change rate of the superhump period was -5.8x10^{-5}. The distance and the X-ray luminosity in the range of 0.5-2.4 keV are estimated to be 380(60) pc and log L_x = 31.0 \pm 0.1 erg s^{-1}. These properties have typical values for an SU UMa-type dwarf nova with this superhump period.
△ Less
Submitted 9 October, 2003;
originally announced October 2003.
-
Period change of Superhumps in the WZ Sge-Type Dwarf Nova, HV Virginis
Authors:
R. Ishioka,
T. Kato,
M. Uemura,
J. Pietz,
T. Vanmunster,
T. Krajci,
K. Torii,
K. Tanabe,
S. Kiyota,
K. Kinugasa,
G. Masi,
K. Morikawa,
L. M. Cook,
P. Schmeer,
H. Yamaoka
Abstract:
After 10 years of quiescence, HV Vir underwent a superoutburst in January 2002. We report time-series observations clearly revealing the period change of ordinary superhumps during the superoutburst. We derived a mean superhump period of 0.058260 d and a positive period derivative of $7 \times 10^{-5}$. These results are in good agreement with the value obtained from the 1992 superoutburst. We a…
▽ More
After 10 years of quiescence, HV Vir underwent a superoutburst in January 2002. We report time-series observations clearly revealing the period change of ordinary superhumps during the superoutburst. We derived a mean superhump period of 0.058260 d and a positive period derivative of $7 \times 10^{-5}$. These results are in good agreement with the value obtained from the 1992 superoutburst. We also detected early superhumps, which were not clearly recognized in the past outburst, and a possible rebrightening. Both of them are the common characteristics of WZ Sge-type stars.
△ Less
Submitted 11 March, 2003;
originally announced March 2003.
-
Dwarf Nova-like Outburst of Short Period Intermediate Polar HT Camelopardalis
Authors:
Ryoko Ishioka,
Taichi Kato,
Makoto Uemura,
Gary W. Billings,
Koichi Morikawa,
Ken'ichi Torii,
Kenji Tanabe,
Arto Oksanen,
Harri Hyvonen,
Hitoshi Itoh
Abstract:
We report the first time-series observations of the short outburst of the proposed intermediate polar HT Cam (=RX J0757.0+6306). On 2001 December 29, we detected the object was undergoing a bright outburst at the magnitude of $m_{vis}=12.2$. Following this detection, we started international joint observations through VSNET. The light curve showed a gradual decline for the first 0.5 d. Following…
▽ More
We report the first time-series observations of the short outburst of the proposed intermediate polar HT Cam (=RX J0757.0+6306). On 2001 December 29, we detected the object was undergoing a bright outburst at the magnitude of $m_{vis}=12.2$. Following this detection, we started international joint observations through VSNET. The light curve showed a gradual decline for the first 0.5 d. Following this short plateau phase, the rate of decline dramatically increased to more than 4 mag d$^{-1}$. Within 1.5 d from the outburst detection, the object almost declined to the quiescent level. During the rapidly declining phase, long-term modulations with a period of 86 min and strong pulses with a period of 8.6 min were observed. We concluded that 86 min and 8.6 min are the orbital period and the spin period of HT Cam, respectively. By the detection of the spin period, we confirmed the IP classification of HT Cam. However, its outburst behavior rather resembles that of dwarf novae.
The discrepancy between the declining rates of the total flux and the pulse flux strongly suggests that the disk instabilities were taking place during the outburst.
△ Less
Submitted 12 June, 2002;
originally announced June 2002.
-
1998 Superoutburst of the Large-Amplitude SU UMa-Type Dwarf Nova WX Ceti
Authors:
Taichi Kato,
Katsura Matsumoto,
Daisaku Nogami,
Koichi Morikawa,
Seiichiro Kiyota
Abstract:
We observed the 1998 November superoutburst of WX Cet, a dwarf nova originally proposed as a WZ Sge-like system. The observation established that WX Cet is an SU UMa-type dwarf nova with a mean superhump period of 0.05949(1) d, which is 2.1% longer than the reported orbital period. The lack of early superhumps at the earliest stage of the superoutburst, the rapid development of usual superhumps,…
▽ More
We observed the 1998 November superoutburst of WX Cet, a dwarf nova originally proposed as a WZ Sge-like system. The observation established that WX Cet is an SU UMa-type dwarf nova with a mean superhump period of 0.05949(1) d, which is 2.1% longer than the reported orbital period. The lack of early superhumps at the earliest stage of the superoutburst, the rapid development of usual superhumps, and the possible rapid decay of late superhumps seem to support that WX Cet is a fairly normal large-amplitude SU UMa-type dwarf nova, rather than a WZ Sge-type dwarf nova with a number of peculiarities. However, a period increase of superhumps at a rate dot(P)/P = (+8.5+/-1.0) x 10^-5 was observed, which is one of the largest dot(P)/P ever observed in SU UMa-type dwarf novae. A linear decline of light, with a rate of 0.10 mag/d, was observed in the post-superoutburst stage. This may be an exemplification of the decay of the viscosity in the accretion disk after the termination of a superoutburst, mechanism of which is proposed to explain a variety of post-superoutburst phenomena in some SU UMa-type dwarf novae.
△ Less
Submitted 10 October, 2001;
originally announced October 2001.
-
Large-N_c meson theory
Authors:
Shin-Ichiro Kuroki,
Keiichi Morikawa,
Takuya Morozumi
Abstract:
We derive an effective Lagrangian for meson fields. This is done in the light-cone gauge for two-dimensional large-N_c QCD by using the bilocal auxiliary field method. The auxiliary fields are bilocal on light-cone space and their Fourier transformation determines the parton momentum distribution. As the first test of our method, the 't Hooft equation is derived from the effective Lagrangian.
We derive an effective Lagrangian for meson fields. This is done in the light-cone gauge for two-dimensional large-N_c QCD by using the bilocal auxiliary field method. The auxiliary fields are bilocal on light-cone space and their Fourier transformation determines the parton momentum distribution. As the first test of our method, the 't Hooft equation is derived from the effective Lagrangian.
△ Less
Submitted 15 March, 2001;
originally announced March 2001.
-
Electronic structure of the strongly hybridized ferromagnet CeFe2
Authors:
T. Konishi,
K. Morikawa,
K. Kobayashi,
T. Mizokawa,
A. Fujimori,
K. Mamiya,
F. Iga,
H. Kawanaka,
Y. Nishihara,
A. Delin,
O. Eriksson
Abstract:
We report on results from high-energy spectroscopic measurements on CeFe2, a system of particular interest due to its anomalous ferromagnetism with an unusually low Curie temperature and small magnetization compared to the other rare earth-iron Laves phase compounds. Our experimental results indicate very strong hybridization of the Ce 4f states with the delocalized band states, mainly the Fe 3d…
▽ More
We report on results from high-energy spectroscopic measurements on CeFe2, a system of particular interest due to its anomalous ferromagnetism with an unusually low Curie temperature and small magnetization compared to the other rare earth-iron Laves phase compounds. Our experimental results indicate very strong hybridization of the Ce 4f states with the delocalized band states, mainly the Fe 3d states. In the interpretation and analysis of our measured spectra, we have made use of two different theoretical approaches: The first one is based on the Anderson impurity model, with surface contributions explicitly taken into account. The second method consists of band-structure calculations for bulk CeFe2. The analysis based on the Anderson impurity model gives calculated spectra in good agreement with the whole range of measured spectra, and reveals that the Ce 4f -- Fe 3d hybridization is considerably reduced at the surface, resulting in even stronger hybridization in the bulk than previously thought. The band-structure calculations are ab initio full-potential linear muffin-tin orbital calculations within the local-spin-density approximation of the density functional. The Ce 4f electrons were treated as itinerant band electrons. Interestingly, the Ce 4f partial density of states obtained from the band-structure calculations also agree well with the experimental spectra concerning both the 4f peak position and the 4f bandwidth, if the surface effects are properly taken into account. In addition, results, notably the partial spin magnetic moments, from the band-structure calculations are discussed in some detail and compared to experimental findings and earlier calculations.
△ Less
Submitted 13 September, 2000;
originally announced September 2000.
-
Scheme Dependence of the Wilsonian Effective Action and Sharp Cutoff Limit of the Flow Equation
Authors:
Jun-Ichi Sumi,
Wataru Souma,
Ken-Ichi Aoki,
Haruhiko Terao,
Keiichi Morikawa
Abstract:
The cutoff scheme dependence in the several formulations of the Exact Renormalization Group (ERG) is investigated. It is shown that the cutoff scheme dependence of the Wilsonian effective action is regarded as a certain coordinate transformation on the theory space. From this observation the Wilsonian effective actions are found to suffer from strong dependence on the schemes even in the infra-r…
▽ More
The cutoff scheme dependence in the several formulations of the Exact Renormalization Group (ERG) is investigated. It is shown that the cutoff scheme dependence of the Wilsonian effective action is regarded as a certain coordinate transformation on the theory space. From this observation the Wilsonian effective actions are found to suffer from strong dependence on the schemes even in the infra-red asymptotic region for massive theories. However there is no such scheme dependence in the one particle irreducible parts of them, which is called the effective average actions. We also derive the explicit form of the Polchinski RG equation in the sharp cutoff limit. Finally this equation is shown to be identical with the Wegner-Houghton RG equation.
△ Less
Submitted 28 February, 2000;
originally announced February 2000.
-
Analysis of the Wilsonian Effective Potentials in Dynamical Chiral Symmetry Breaking
Authors:
Ken-Ichi Aoki,
Keiichi Morikawa,
Jun-Ichi Sumi,
Haruhiko Terao,
Masashi Tomoyose
Abstract:
The non-perturbative renormalization group equation for the Wilsonian effective potential is given in a certain simple approximation scheme in order to study chiral symmetry breaking phenomena dynamically induced by strong gauge interactions. The evolving effective potential is found to be non-analytic in infrared, which indicates spontaneous generation of the fermion mass. It is also shown that…
▽ More
The non-perturbative renormalization group equation for the Wilsonian effective potential is given in a certain simple approximation scheme in order to study chiral symmetry breaking phenomena dynamically induced by strong gauge interactions. The evolving effective potential is found to be non-analytic in infrared, which indicates spontaneous generation of the fermion mass. It is also shown that the renormalization group equation gives the identical effective fermion mass with that obtained by solving the Schwinger-Dyson equation in the (improved) ladder approximation. Moreover introduction of the collective field corresponding to the fermion composite into the theory space is found to offer an efficient method to evaluate the order parameters; the dynamical mass and the chiral condensate. The relation between the renormalization group equation incorporating the collective field and the Schwinger-Dyson equation is also clarified.
△ Less
Submitted 4 August, 1999;
originally announced August 1999.
-
Wilson Renormalization Group Equations for the Critical Dynamics of Chiral Symmetry
Authors:
Ken-Ichi Aoki,
Keiichi Morikawa,
Jun-Ichi Sumi,
Haruhiko Terao,
Masashi Tomoyose
Abstract:
The critical dynamics of the chiral symmetry breaking induced by gauge interaction is examined in the Wilson renormalization group framework in comparison with the Schwinger-Dyson approach. We derive the beta functions for the four-fermi couplings in the sharp cutoff renormalzation group scheme, from which the critical couplings and the anomalous dimensions of the fermion composite operators nea…
▽ More
The critical dynamics of the chiral symmetry breaking induced by gauge interaction is examined in the Wilson renormalization group framework in comparison with the Schwinger-Dyson approach. We derive the beta functions for the four-fermi couplings in the sharp cutoff renormalzation group scheme, from which the critical couplings and the anomalous dimensions of the fermion composite operators near criticality are immediately obtained. It is also shown that the beta functions lead to the same critical behavior found by solving the so-called ladder Schwinger-Dyson equation, if we restrict the radiative corrections to a certain limited type.
△ Less
Submitted 4 August, 1999;
originally announced August 1999.
-
Rapidly Converging Truncation Scheme of the Exact Renormalization Group
Authors:
Ken-Ichi Aoki,
Keiichi Morikawa,
Wataru Souma,
Jun-Ichi Sumi,
Haruhiko Terao
Abstract:
The truncation scheme dependence of the exact renormalization group equations is investigated for scalar field theories in three dimensions. The exponents are numerically estimated to the next-to-leading order of the derivative expansion. It is found that the convergence property in various truncations in the number of powers of the fields is remarkably improved if the expansion is made around t…
▽ More
The truncation scheme dependence of the exact renormalization group equations is investigated for scalar field theories in three dimensions. The exponents are numerically estimated to the next-to-leading order of the derivative expansion. It is found that the convergence property in various truncations in the number of powers of the fields is remarkably improved if the expansion is made around the minimum of the effective potential. It is also shown that this truncation scheme is suitable for evaluation of infrared effective potentials. The physical interpretation of this improvement is discussed by considering O(N) symmetric scalar theories in the large-N limit.
△ Less
Submitted 23 March, 1998; v1 submitted 6 March, 1998;
originally announced March 1998.
-
Discovery of X-ray emission rom the distant lensing cluster of galaxies CL2236-04 at z = 0.552
Authors:
M. Hattori,
H. Matuzawa,
K. Morikawa,
J. -P. Kneib,
K. Yamashita,
K. Watanabe,
H. Boehringer,
T. G. Tsuru
Abstract:
X-ray emission from the distant lensing cluster CL2236-04 at $z$ = 0.552 was discovered by ASCA and ROSAT/HRI observations. If the spherical symmetric mass distribution model of the cluster is assumed, the lensing estimate of the cluster mass is a factor of two higher than that obtained from X-ray observations as reported for many distant clusters. However, the elliptical and clumpy lens model p…
▽ More
X-ray emission from the distant lensing cluster CL2236-04 at $z$ = 0.552 was discovered by ASCA and ROSAT/HRI observations. If the spherical symmetric mass distribution model of the cluster is assumed, the lensing estimate of the cluster mass is a factor of two higher than that obtained from X-ray observations as reported for many distant clusters. However, the elliptical and clumpy lens model proposed by Kneib et al.(1993) is surprisingly consistent with the X-ray observations assuming that the X-ray emitting hot gas is isothermal and in a hydrostatic equilibrium state. The existence of the cooling flow in the central region of the cluster is indicated by the short central cooling time and the excess flux detected by ROSAT/HRI compared to the ASCA flux. However, it is shown that even if the AXJ2239-0429 has a cooling flow in the central region, the temperature measured by ASCA which is the mean emission-weighted cluster temperature in this case, should not be cooler than and different from the virial temperature of the cluster. Therefore, we conclude that the effect of the clumpiness and non-zero ellipticity in the mass distribution of the cluster are essential to explain the observed feature of the giant luminous arc, and there is no discrepancy between strong lensing and X-ray estimation of the mass of the cluster in this cluster.
△ Less
Submitted 9 March, 1998;
originally announced March 1998.
-
Non-Perturbative Renormalization Group Analysis of the Chiral Critical Behaviors in QED
Authors:
Ken-Ichi Aoki,
Keiichi Morikawa,
Jun-Ichi Sumi,
Haruhiko Terao,
Masashi Tomoyose
Abstract:
We study the chiral critical behaviors of QED by using Non-Perturbative Renormalization Group (NPRG). Taking account of the non-ladder contributions, our flow equations are free from the gauge parameter ($α$) dependence. We clarify the chiral phase structure, and calculate the anomalous dimension of $\barψψ$, which is enhanced compared to the ladder approximation. We find that the cutoff scheme…
▽ More
We study the chiral critical behaviors of QED by using Non-Perturbative Renormalization Group (NPRG). Taking account of the non-ladder contributions, our flow equations are free from the gauge parameter ($α$) dependence. We clarify the chiral phase structure, and calculate the anomalous dimension of $\barψψ$, which is enhanced compared to the ladder approximation. We find that the cutoff scheme dependence of the physical results is very small.
△ Less
Submitted 29 January, 1997; v1 submitted 24 December, 1996;
originally announced December 1996.
-
The effectiveness of the local potential approximation in the Wegner-Houghton renormalization group
Authors:
Ken-Ichi Aoki,
Keiichi Morikawa,
Wataru Souma,
Jun-ichi Sumi,
Haruhiko Terao
Abstract:
The non-perturbative Wegner-Houghton renormalization group is analyzed by the local potential approximation in O(N) scalar theories in d-dimensions $(3\leq d\leq 4)$. The leading critical exponents νare calculated in order to investigate the effectiveness of the local potential approximation by comparing them with the other non-perturbative methods. We show analytically that the local potential…
▽ More
The non-perturbative Wegner-Houghton renormalization group is analyzed by the local potential approximation in O(N) scalar theories in d-dimensions $(3\leq d\leq 4)$. The leading critical exponents νare calculated in order to investigate the effectiveness of the local potential approximation by comparing them with the other non-perturbative methods. We show analytically that the local potential approximation gives the exact exponents up to $O(ε)$ in ε-expansion and the leading in 1/N-expansion. We claim that this approximation offers fairly accurate results in the whole range of the parameter space of N and d. It is a great advantage of our method that no diverging expansions appear in the procedure.
△ Less
Submitted 24 December, 1996;
originally announced December 1996.