# Small Area Estimation using Hierarchical Bayesian under Student-t (non central) Distribution This function is implemented to variable of interest $(y)$ that assumed to be a Student-t (non central) Distribution. The range of data is $(-\infty < y < \infty)$ ```r Student_tnc( formula, iter.update = 3, iter.mcmc = 10000, coef, var.coef, thin = 2, burn.in = 2000, tau.u = 1, data ) ``` ## Arguments - `formula`: Formula that describe the fitted model - `iter.update`: Number of updates with default `3` - `iter.mcmc`: Number of total iterations per chain with default `10000` - `coef`: a vector contains prior initial value of Coefficient of Regression Model for fixed effect with default vector of `0` with the length of the number of regression coefficients - `var.coef`: a vector contains prior initial value of variance of Coefficient of Regression Model with default vector of `1` with the length of the number of regression coefficients - `thin`: Thinning rate, must be a positive integer with default `2` - `burn.in`: Number of iterations to discard at the beginning with default `2000` - `tau.u`: Prior initial value of inverse of Variance of area random effect with default `1` - `data`: The data frame ## Returns This function returns a list of the following objects: - **Est**: A vector with the values of Small Area mean Estimates using Hierarchical bayesian method - **refVar**: Estimated random effect variances - **coefficient**: A dataframe with the estimated model coefficient - **plot**: Trace, Dencity, Autocorrelation Function Plot of MCMC samples ## Examples ```r ## Data Generation set.seed(123) m <- 30 x1 <- runif(m, 1, 100) x2 <- runif(m, 10, 15) b0 <- b1 <- b2 <- 0.5 u <- rnorm(m, 0, 1) MU <- b0 + b1 * x1 + b2 * x2 + u k <- rgamma(1, 10, 1) y <- rt(m, k, MU) vardir <- k / (k - 1) vardir <- sd(y)^2 datatnc <- as.data.frame(cbind(y, x1, x2, vardir)) datatncNs <- datatnc datatncNs$y[c(3, 14, 22, 29, 30)] <- NA datatncNs$vardir[c(3, 14, 22, 29, 30)] <- NA ## Compute Fitted Model ## y ~ x1 +x2 ## For data without any nonsampled area formula <- y ~ x1+x2 v <- c(1, 1, 1) c <- c(0, 0, 0) dat <- datatnc ## Using parameter coef and var.coef saeHBtnc <- Student_tnc(formula, coef = c, var.coef = v, iter.update = 10, data = dat) saeHBtnc$Est # Small Area mean Estimates saeHBtnc$refVar # Random effect variance saeHBtnc$coefficient # coefficient # Load Library 'coda' to execute the plot # autocorr.plot(saeHBtnc$plot[[3]]) is used to generate ACF Plot # plot(saeHBtnc$plot[[3]]) is used to generate Density and trace plot ## Do not using parameter coef and var.coef saeHBtnc <- Student_tnc(formula, data = datatnc) ## For data with nonsampled area use datatncNs ```

Student_tnc function