1 Estandarización

1.1 Score-Z

\[ Z=\frac{x-\mu}{\sigma} \]

Nota X= variable de interes
\(\mu\)=media de los datos
\(\sigma\)=desviación estandar

1.2 Usando R para estandarizar una variable

set.seed(123)

x= rnorm(60,3,0.3)
head(x)
## [1] 2.831857 2.930947 3.467612 3.021153 3.038786 3.514519
hist(x, ylab = "Frecuencia", col="skyblue")
abline(v = mean(x),
       col="red",
       lty=3,
       lwd=3)

boxplot(x, col="darkgreen")
points(mean(x), col="red", cex= 1.5, pch=16)

1.3 Score-Z de la variable X

ZX = scale(x)
head(ZX)
##              [,1]
## [1,] -0.687729949
## [2,] -0.324914908
## [3,]  1.640081452
## [4,]  0.005372616
## [5,]  0.069938613
## [6,]  1.811830980
par(mfrow=c(1,2))
hist(x)
hist(ZX, breaks=6)

plot(x, ZX, col="blue", main= "Gráfico de dispersión")

1.4 Correlacción de pearson

cor(x, ZX, method="pearson")
##      [,1]
## [1,]    1

1.5 Transformaciones lineales

1.6 Demostración

1.7 Demostración Matemática

La transformación \(Z\) se define mediante la siguiente expresión:

\[Z = \frac{x - \mu}{\sigma}\]

Podemos descomponer esta fracción para observar su estructura algebraica:

\[Z(x) = \underbrace{\left( \frac{1}{\sigma} \right)}_{a} x + \underbrace{\left( -\frac{\mu}{\sigma} \right)}_{b}\]

Al compararla con la forma general de una línea recta \(y = ax + b\), observamos que:

  • Pendiente (\(a\)): \(\frac{1}{\sigma}\), que representa un cambio de escala.
  • Intersección (\(b\)): \(-\frac{\mu}{\sigma}\), que representa una traslación.

Por lo tanto, al ser una función de primer grado respecto a \(x\), el Z-score es una transformación lineal (técnicamente afín).

1.8 Demostración: Incumplimiento de la Propiedad Aditiva

Para que el \(Z\)-score sea una transformación lineal estricta, debe cumplir que \(Z(x + y) = Z(x) + Z(y)\).

Sea la función del \(Z\)-score: \[Z(x) = \frac{x - \mu}{\sigma}\]

1. Evaluamos el lado izquierdo \(Z(x + y)\): \[Z(x + y) = \frac{(x + y) - \mu}{\sigma}\]

2. Evaluamos el lado derecho \(Z(x) + Z(y)\): \[Z(x) + Z(y) = \frac{x - \mu}{\sigma} + \frac{y - \mu}{\sigma} = \frac{x + y - 2\mu}{\sigma}\]

3. Comparación: Comparamos ambos resultados: \[\frac{x + y - \mu}{\sigma} \neq \frac{x + y - 2\mu}{\sigma}\]

La igualdad solo se cumpliría si \(\mu = 0\). Como en la mayoría de los conjuntos de datos la media \(\mu \neq 0\), entonces: \[Z(x + y) \neq Z(x) + Z(y)\]

Conclusión: El \(Z\)-score no es una transformación lineal estricta bajo la propiedad aditiva. Matemáticamente, es una transformación afín, ya que incluye una traslación (el término \(-\mu\)) que “rompe” la linealidad en el sentido del álgebra lineal pura.

cat( "Media de X:", mean(x));
## Media de X: 3.019685
cat("Media de ZX:", round(mean(ZX),1))
## Media de ZX: 0
var(x)
## [1] 0.07459062
var(ZX)
##      [,1]
## [1,]    1
LS0tDQp0aXRsZTogIkNsYXNlIDMgLSBDRSINCnN1YnRpdGxlOiAiQW7DoWxpc2lzIEVzdGFkw61zdGljbyBBdmFuemFkbyINCmF1dGhvcjogIlR1IE5vbWJyZSINCmRhdGU6ICIyMDI2LTAyLTA5Ig0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50Og0KICAgIHRoZW1lOiBmbGF0bHkgICAgICMgQ2FtYmlhIGVsIGVzdGlsbyAob3Ryb3M6IGZsYXRseSwgbHVtZW4sIHlldGkpDQogICAgaGlnaGxpZ2h0OiB0YW5nbyAgICAjIEVzdGlsbyBkZSBjb2xvcmVzIGRlbCBjw7NkaWdvDQogICAgdG9jOiB0cnVlICAgICAgICAgICAjIEluY2x1eWUgdGFibGEgZGUgY29udGVuaWRvcw0KICAgIHRvY19mbG9hdDogdHJ1ZSAgICAgIyBFbCDDrW5kaWNlIHRlIHNpZ3VlIGFsIGhhY2VyIHNjcm9sbA0KICAgIGNvZGVfZm9sZGluZzogc2hvdyAgIyBQZXJtaXRlIG9jdWx0YXIvbW9zdHJhciBsb3MgYmxvcXVlcyBkZSBjw7NkaWdvDQogICAgZGZfcHJpbnQ6IHBhZ2VkICAgICAjIExhcyB0YWJsYXMgZ3JhbmRlcyBzZSB2ZW4gY29uIHBhZ2luYWNpw7NuIGludGVyYWN0aXZhDQogICAgY29kZV9kb3dubG9hZDogdHJ1ZSAjIEFncmVnYSB1biBib3TDs24gcGFyYSBiYWphciBlbCBhcmNoaXZvIC5SbWQNCiAgICBudW1iZXJfc2VjdGlvbnM6IHllcw0KLS0tDQoNCmBgYHs9aHRtbH0NCjxzdHlsZT4NCi5ub3RhLWF6dWwgew0KICBiYWNrZ3JvdW5kLWNvbG9yOiAjZjBmOGZmOw0KICBib3JkZXI6IDFweCBzb2xpZCAjYWRkOGU2Ow0KICBwYWRkaW5nOiAxMHB4Ow0KICBib3JkZXItcmFkaXVzOiA4cHg7DQogIGNvbG9yOiAjMmMzZTUwOw0KfQ0KPC9zdHlsZT4NCmBgYA0KDQojIyBFc3RhbmRhcml6YWNpw7NuDQoNCiMjIyBTY29yZS1aDQoNCiQkDQpaPVxmcmFje3gtXG11fXtcc2lnbWF9DQokJA0KDQo6Ojoge3N0eWxlPSJiYWNrZ3JvdW5kLWNvbG9yOiAjZTZmN2ZmOyBib3JkZXItbGVmdDogNnB4IHNvbGlkICMwMDdiZmY7IHBhZGRpbmc6IDE1cHg7IG1hcmdpbjogMTBweCAwOyBib3JkZXItcmFkaXVzOiA1cHg7In0NCioqTm90YSoqIFg9IHZhcmlhYmxlIGRlIGludGVyZXM8YnI+ICRcbXUkPW1lZGlhIGRlIGxvcyBkYXRvcyA8YnI+ICRcc2lnbWEkPWRlc3ZpYWNpw7NuIGVzdGFuZGFyDQo6OjoNCg0KIyMjIFVzYW5kbyBSIHBhcmEgZXN0YW5kYXJpemFyIHVuYSB2YXJpYWJsZQ0KDQpgYGB7cn0NCnNldC5zZWVkKDEyMykNCg0KeD0gcm5vcm0oNjAsMywwLjMpDQpoZWFkKHgpDQoNCg0KYGBgDQoNCmBgYHtyfQ0KaGlzdCh4LCB5bGFiID0gIkZyZWN1ZW5jaWEiLCBjb2w9InNreWJsdWUiKQ0KYWJsaW5lKHYgPSBtZWFuKHgpLA0KICAgICAgIGNvbD0icmVkIiwNCiAgICAgICBsdHk9MywNCiAgICAgICBsd2Q9MykNCg0KYGBgDQoNCmBgYHtyfQ0KDQpib3hwbG90KHgsIGNvbD0iZGFya2dyZWVuIikNCnBvaW50cyhtZWFuKHgpLCBjb2w9InJlZCIsIGNleD0gMS41LCBwY2g9MTYpDQpgYGANCg0KIyMjIFNjb3JlLVogZGUgbGEgdmFyaWFibGUgWA0KDQpgYGB7cn0NCg0KWlggPSBzY2FsZSh4KQ0KaGVhZChaWCkNCmBgYA0KDQpgYGB7cn0NCnBhcihtZnJvdz1jKDEsMikpDQpoaXN0KHgpDQpoaXN0KFpYLCBicmVha3M9NikNCmBgYA0KDQpgYGB7cn0NCnBsb3QoeCwgWlgsIGNvbD0iYmx1ZSIsIG1haW49ICJHcsOhZmljbyBkZSBkaXNwZXJzacOzbiIpDQpgYGANCg0KIyMjIENvcnJlbGFjY2nDs24gZGUgcGVhcnNvbg0KDQpgYGB7cn0NCmNvcih4LCBaWCwgbWV0aG9kPSJwZWFyc29uIikNCmBgYA0KDQojIyMgVHJhbnNmb3JtYWNpb25lcyBsaW5lYWxlcyANCg0KIyMjIERlbW9zdHJhY2nDs24NCg0KIyMjIERlbW9zdHJhY2nDs24gTWF0ZW3DoXRpY2ENCg0KTGEgdHJhbnNmb3JtYWNpw7NuICRaJCBzZSBkZWZpbmUgbWVkaWFudGUgbGEgc2lndWllbnRlIGV4cHJlc2nDs246DQoNCiQkWiA9IFxmcmFje3ggLSBcbXV9e1xzaWdtYX0kJA0KDQpQb2RlbW9zIGRlc2NvbXBvbmVyIGVzdGEgZnJhY2Npw7NuIHBhcmEgb2JzZXJ2YXIgc3UgZXN0cnVjdHVyYSBhbGdlYnJhaWNhOg0KDQokJFooeCkgPSBcdW5kZXJicmFjZXtcbGVmdCggXGZyYWN7MX17XHNpZ21hfSBccmlnaHQpfV97YX0geCArIFx1bmRlcmJyYWNle1xsZWZ0KCAtXGZyYWN7XG11fXtcc2lnbWF9IFxyaWdodCl9X3tifSQkDQoNCkFsIGNvbXBhcmFybGEgY29uIGxhIGZvcm1hIGdlbmVyYWwgZGUgdW5hIGzDrW5lYSByZWN0YSAkeSA9IGF4ICsgYiQsIG9ic2VydmFtb3MgcXVlOg0KDQotICAgKipQZW5kaWVudGUgKCoqJGEkKTogJFxmcmFjezF9e1xzaWdtYX0kLCBxdWUgcmVwcmVzZW50YSB1biBjYW1iaW8gZGUgZXNjYWxhLg0KLSAgICoqSW50ZXJzZWNjacOzbiAoKiokYiQpOiAkLVxmcmFje1xtdX17XHNpZ21hfSQsIHF1ZSByZXByZXNlbnRhIHVuYSB0cmFzbGFjacOzbi4NCg0KUG9yIGxvIHRhbnRvLCBhbCBzZXIgdW5hIGZ1bmNpw7NuIGRlIHByaW1lciBncmFkbyByZXNwZWN0byBhICR4JCwgZWwgKipaLXNjb3JlIGVzIHVuYSB0cmFuc2Zvcm1hY2nDs24gbGluZWFsKiogKHTDqWNuaWNhbWVudGUgYWbDrW4pLg0KDQojIyMgRGVtb3N0cmFjacOzbjogSW5jdW1wbGltaWVudG8gZGUgbGEgUHJvcGllZGFkIEFkaXRpdmENCg0KUGFyYSBxdWUgZWwgJFokLXNjb3JlIHNlYSB1bmEgdHJhbnNmb3JtYWNpw7NuIGxpbmVhbCBlc3RyaWN0YSwgZGViZSBjdW1wbGlyIHF1ZSAkWih4ICsgeSkgPSBaKHgpICsgWih5KSQuDQoNClNlYSBsYSBmdW5jacOzbiBkZWwgJFokLXNjb3JlOiAkJFooeCkgPSBcZnJhY3t4IC0gXG11fXtcc2lnbWF9JCQNCg0KKioxLiBFdmFsdWFtb3MgZWwgbGFkbyBpenF1aWVyZG8qKiAkWih4ICsgeSkkOiAkJFooeCArIHkpID0gXGZyYWN7KHggKyB5KSAtIFxtdX17XHNpZ21hfSQkDQoNCioqMi4gRXZhbHVhbW9zIGVsIGxhZG8gZGVyZWNobyoqICRaKHgpICsgWih5KSQ6ICQkWih4KSArIFooeSkgPSBcZnJhY3t4IC0gXG11fXtcc2lnbWF9ICsgXGZyYWN7eSAtIFxtdX17XHNpZ21hfSA9IFxmcmFje3ggKyB5IC0gMlxtdX17XHNpZ21hfSQkDQoNCioqMy4gQ29tcGFyYWNpw7NuOioqIENvbXBhcmFtb3MgYW1ib3MgcmVzdWx0YWRvczogJCRcZnJhY3t4ICsgeSAtIFxtdX17XHNpZ21hfSBcbmVxIFxmcmFje3ggKyB5IC0gMlxtdX17XHNpZ21hfSQkDQoNCkxhIGlndWFsZGFkIHNvbG8gc2UgY3VtcGxpcsOtYSBzaSAkXG11ID0gMCQuIENvbW8gZW4gbGEgbWF5b3LDrWEgZGUgbG9zIGNvbmp1bnRvcyBkZSBkYXRvcyBsYSBtZWRpYSAkXG11IFxuZXEgMCQsIGVudG9uY2VzOiAkJFooeCArIHkpIFxuZXEgWih4KSArIFooeSkkJA0KDQoqKkNvbmNsdXNpw7NuOioqIEVsICRaJC1zY29yZSAqKm5vIGVzIHVuYSB0cmFuc2Zvcm1hY2nDs24gbGluZWFsIGVzdHJpY3RhKiogYmFqbyBsYSBwcm9waWVkYWQgYWRpdGl2YS4gTWF0ZW3DoXRpY2FtZW50ZSwgZXMgdW5hICoqdHJhbnNmb3JtYWNpw7NuIGFmw61uKiosIHlhIHF1ZSBpbmNsdXllIHVuYSB0cmFzbGFjacOzbiAoZWwgdMOpcm1pbm8gJC1cbXUkKSBxdWUgInJvbXBlIiBsYSBsaW5lYWxpZGFkIGVuIGVsIHNlbnRpZG8gZGVsIMOhbGdlYnJhIGxpbmVhbCBwdXJhLg0KDQpgYGB7cn0NCmNhdCggIk1lZGlhIGRlIFg6IiwgbWVhbih4KSk7DQpjYXQoIk1lZGlhIGRlIFpYOiIsIHJvdW5kKG1lYW4oWlgpLDEpKQ0KDQpgYGANCg0KYGBge3J9DQp2YXIoeCkNCnZhcihaWCkNCmBgYA0K