mnTransform {Transform} | R Documentation |
Manly Transformation for Normality
Description
mnTransform
performs Manly transformation for normality of a variable and provides graphical analysis.
Usage
mnTransform(data, lambda = seq(-3,3,0.01), plot = TRUE, alpha = 0.05,
verbose = TRUE)
Arguments
data |
a numeric vector of data values. |
lambda |
a vector which includes the sequence of candidate lambda values. Default is set to (-3,3) with increment 0.01. |
plot |
a logical to plot histogram with its density line and qqplot of raw and transformed data. Defaults plot = TRUE. |
alpha |
the level of significance to check the normality after transformation. Default is set to alpha = 0.05. |
verbose |
a logical for printing output to R console. |
Details
Denote y
the variable at the original scale and y'
the transformed variable. The Manly power transformation is defined by:
y' = \left\{ \begin{array}{ll}
\frac{e^{\lambda y} - 1}{\lambda} \mbox{ , if $\lambda \neq 0$} \cr
y \mbox{ , if $\lambda = 0$}
\end{array} \right.
Value
A list with class "mn" containing the following elements:
method |
method to estimate Manly transformation parameter |
lambda.hat |
estimate of Manly transformation parameter |
statistic |
Shapiro-Wilk test statistic for transformed data |
p.value |
Shapiro-Wilk test p.value for transformed data |
alpha |
level of significance to assess normality |
tf.data |
transformed data set |
var.name |
variable name |
Author(s)
Muge Coskun Yildirim, Osman Dag
References
Asar, O., Ilk, O., Dag, O. (2017). Estimating Box-Cox Power Transformation Parameter via Goodness of Fit Tests. Communications in Statistics - Simulation and Computation, 46:1, 91–105.
Manly, B.F. (1976). Exponential Data Transformations. Journal of the Royal Statistical Society: Series D (The Statistician), 25:1, 37–42.
Examples
data <- cars$dist
library(Transform)
out <- mnTransform(data)
out$lambda.hat # the estimate of Manly parameter based on Shapiro-Wilk test statistic
out$p.value # p.value of Shapiro-Wilk test for transformed data
out$tf.data # transformed data set