Y, ¿cuál trasformación uso? (en R)

Reemplazando Y à Y^l para “normalizar” variables. La potencia l proporciona posibles transformaciones para normalidad (Box-Cox).

Si l =	-2	-1	-0.5	0	0.5	1	2
Transformación	1/Y2	1/Y	1/Y0.5	log(Y)	Y0.5	Ninguna	Y2

# Este programa realiza la transformacion Box-Cox de la Tabla 3.9 Niveles de Plasma en Neter, Kutner, Nachtsheim and Wasserman, 4th ed. pp. 134-7

# Cargando datos de 'orig.txt'
orig = read.table("C:/Users/Administrator/Desktop/orig.txt", header = T)
orig

   age plasma
1    0  13.44
2    0  12.84
3    0  11.91
4    0  20.09
5    0  15.60
6    1  10.11
7    1  11.38
8    1  10.28
9    1   8.96
10   1   8.59
11   2   9.83
12   2   9.00
13   2   8.65
14   2   7.85
15   2   8.88
16   3   7.94
17   3   6.01
18   3   5.14
19   3   6.90
20   3   6.77
21   4   4.86
22   4   5.10
23   4   5.67
24   4   5.75
25   4   6.23

# Cargando 'car' a la biblioteca
library(car)

Loading required package: MASS

Loading required package: nnet

summary(powerTransform(orig$plasma))

bcPower Transformation to Normality 

            Est.Power Std.Err. Wald Lower Bound Wald Upper Bound
orig$plasma    -0.455   0.5376           -1.509           0.5987

Likelihood ratio tests about transformation parameters
                         LRT df     pval
LR test, lambda = (0) 0.7434  1 0.388585
LR test, lambda = (1) 8.0795  1 0.004477

Lamda = -0.455 sugiere la trasnformación de Y' =  1/Y^0.5 Ademas, existe evidencia que la transformación es necesaria (p = 0.004477).