BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Tatsuya Kubokawa

Tatsuya Kubokawa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Methodologymath.STStatistics Theory

Trust snapshot

Quick read

Trust 21 - Emerging
10works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2025arXiv

Priors for second-order unbiased Bayes estimators

Asymptotically unbiased priors, introduced by Hartigan (1965), are designed to achieve second-order unbiasedness of Bayes estimators. This paper extends Hartigan's framework to non-i.i.d. models by deriving a system of partial differential equations that characterizes asymptotically unbiased priors. Furthermore, we establish a necessary and sufficient condition for the existence of such priors and propose a simple procedure for constructing them. The proposed method is applied to the linear regression model and the nested error regression model (also known as the random effects model). Simulation studies evaluate the frequentist properties of the Bayes estimator under the asymptotically unbiased prior for the nested error regression model, highlighting its effectiveness in small-sample settings.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ridge-type Linear Shrinkage Estimation of the Matrix Mean of High-dimensional Normal Distribution

The estimation of the mean matrix of the multivariate normal distribution is addressed in the high dimensional setting. Efron-Morris-type linear shrinkage estimators based on ridge estimators for the precision matrix instead of the Moore-Penrose generalized inverse are considered, and the weights in the ridge-type linear shrinkage estimators are estimated in terms of minimizing the Stein unbiased risk estimators under the quadratic loss. It is shown that the ridge-type linear shrinkage estimators with the estimated weights are minimax, and that the estimated weights converge to the optimal weights in the Bayesian model with high dimension by using the random matrix theory. The performance of the ridge-type linear shrinkage estimators is numerically compared with the existing estimators including the Efron-Morris and James-Stein estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Consistent Estimation in Box-Cox Transformed Linear Mixed Models

The Box-Cox transformation is applied to the linear mixed models for analyzing positive and grouped data. The problem in using Box Cox transformation is that the maximum likelihood estimator of the transformation parameter is generally inconsistent. To fix it, we suggest a simple and consistent estimator for the transformation parameter based on the moment method. The consistent estimator is used to construct consistent estimators of the parameters involved in the model and under some conditions, the estimators of model parameters are shown to be consistent under $m\to\infty$, where $m$ is the number of groups. Moreover, in estimation of the expectation of the (future) observations, it is shown that the resulting estimators also hold consistency owing to the consistent estimator of the transformation parameter. The proposed estimating method is compared with the maximum likelihood method via simulation and it is shown that the proposed estimator works better than the maximum likelihood estimator. Finally, the proposed method is applied to a real data set.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Empirical Uncertain Bayes Methods in Area-level Models

Random effects model can account for the lack of fitting a regression model and increase precision of estimating area-level means. However, in case that the synthetic mean provides accurate estimates, the prior distribution may inflate an estimation error. Thus it is desirable to consider the uncertain prior distribution, which is expressed as the mixture of a one-point distribution and a proper prior distribution. In this paper, we develop an empirical Bayes approach for estimating area-level means, using the uncertain prior distribution in the context of a natural exponential family, which we call the empirical uncertain Bayes (EUB) method. The regression model considered in this paper includes the Poisson-gamma and the binomial-beta, and the normal-normal (Fay-Herriot) model, which are typically used in small area estimation. We obtain the estimators of hyperparameters based on the marginal likelihood by using a well-known Expectation-Maximization algorithm, and propose the EUB estimators of area means. For risk evaluation of the EUB estimator, we derive a second-order unbiased estimator of a conditional mean squared error by using some techniques of numerical calculation. Through simulation studies and real data applications, we evaluate a performance of the EUB estimator and compare it with the usual empirical Bayes estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heteroscedastic Nested Error Regression Models with Variance Functions

The nested error regression model is a useful tool for analyzing clustered (grouped) data, and is especially used in small area estimation. The classical nested error regression model assumes normality of random effects and error terms, and homoscedastic variances. However, these assumptions are often violated in real applications and more flexible models are required. This article proposes a nested error regression model with heteroscedastic variances, where the normality for the underlying distributions is not assumed. We propose the structure of heteroscedastic variances by using some specified variance functions and some covariates with unknown parameters. Under the setting, we construct the moment-type estimators of model parameters and some asymptotic properties including asymptotic biases and variances are derived. For predicting linear quantities including random effects, we suggest the empirical best linear unbiased predictors and the second-order unbiased estimators of mean squared errors are derived in the closed form. We investigate the proposed method with simulation and empirical studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On Conditional Prediction Errors in Mixed Models with Application to Small Area Estimation

The empirical Bayes estimators in mixed models are useful for small area estimation in the sense of increasing precision of prediction for small area means, and one wants to know the prediction errors of the empirical Bayes estimators based on the data. This paper is concerned with conditional prediction errors in the mixed models instead of conventional unconditional prediction errors. In the mixed models based on natural exponential families with quadratic variance functions, it is shown that the difference between the conditional and unconditional prediction errors is significant under distributions far from normality. Especially for the binomial-beta mixed and the Poisson-gamma mixed models, the leading terms in the conditional prediction errors are, respectively, a quadratic concave function and an increasing function of the direct estimate in the small area, while the corresponding leading terms in the unconditional prediction errors are constants. Second-order unbiased estimators of the conditional prediction errors are also derived and their performances are examined through simulation and empirical studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Small Area Predictors with Dual Shrinkage of Means and Variances

The paper concerns small-area estimation in the Fay-Herriot type area-level model with random dispersions, which models the case that the sampling errors change from area to area. The resulting Bayes estimator shrinks both means and variances, but needs numerical computation to provide the estimates. In this paper, an approximated empirical Bayes (AEB) estimator with a closed form is suggested. The model parameters are estimated via the moment method, and the mean squared error of the AEB is estimated via the single parametric bootstrap method. The benchmarked estimator and a second-order unbiased estimator of the mean squared error are also derived.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Predictive Density Estimation for Location Families under Integrated $L_2$ and $L_1$ Losses

Our investigation concerns the estimation of predictive densities and a study of efficiency as measured by the frequentist risk of such predictive densities with integrated $L_2$ and $L_1$ losses. Our findings relate to a $p-$variate spherically symmetric observable $X \sim p_X(\|x-μ\|^2)$ and the objective of estimating the density of $Y \sim q_Y(\|y-μ\|^2)$ based on $X$. For $L_2$ loss, we describe Bayes estimation, minimum risk equivariant estimation (MRE), and minimax estimation. We focus on the risk performance of the benchmark minimum risk equivariant estimator, plug-in estimators, and plug-in type estimators with expanded scale. For the multivariate normal case, we make use of a duality result with a point estimation problem bringing into play reflected normal loss. In three of more dimensions (i.e., $p \geq 3$), we show that the MRE estimator is inadmissible under $L_2$ loss and provide dominating estimators. This brings into play Stein-type results for estimating a multivariate normal mean with a loss which is a concave and increasing function of $\|\hatμ-μ\|^2$. We also study the phenomenon of improvement on the plug-in density estimator of the form $q_Y(\|y-aX\|^2)\,, 0<a \leq 1\,,$ by a subclass of scale expansions $\frac{1}{c^p} \, q_Y(\|(y -aX)/c \|^2)$ with $c>1$, showing in some cases, inevitably for large enough $p$, that all choices $c>1$ are dominating estimators. Extensions are obtained for scale mixture of normals including a general inadmissibility result of the MRE estimator for $p \geq 3$. Finally, we describe and expand on analogous plug-in dominance results for spherically symmetric distributions with $p \geq 4$ under $L_1$ loss.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parametric Transformed Fay-Herriot Model for Small Area Estimation

In this paper, we consider parametric transformed Fay-Herriot models, and clarify conditions on transformations under which the estimator of the transformation is consistent. It is shown that the dual power transformation satisfies the conditions. Based on asymptotic properties for estimators of parameters, we derive a second-order approximation of the prediction error of the empirical best linear unbiased predictors (EBLUP) and obtain a second-order unbiased estimator of the prediction error. Finally, performances of the proposed procedures are investigated through simulation and empirical studies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dominance properties of constrained Bayes and empirical Bayes estimators

This paper studies decision theoretic properties of benchmarked estimators which are of some importance in small area estimation problems. Benchmarking is intended to improve certain aggregate properties (such as study-wide averages) when model based estimates have been applied to individual small areas. We study decision-theoretic properties of such estimators by reducing the problem to one of studying these problems in a related derived problem. For certain such problems, we show that unconstrained solutions in the original (unbenchmarked) problem give unconstrained Bayes and improved estimators which automatically satisfy the benchmark constraint. Also, dominance properties of constrained empirical Bayes estimators are shown in the Fay-Herriot model, a frequently used model in small area estimation.

0 citations0 reviewsNo code linkedOpen access