Paquetes de uso frecuente en el curso
library(ggplot2)
library(plotly)
library(tidyverse)
Tenemos dos muestras, n1 y n2. Queremos
saber si las dos pertenecen a la misma población (N). Supongamos que
nuestro marco teórico nos da suficiente evidencia para predecir (según
la hipótesis de trabajo o investigación) que las muestras son de
diferentes poblaciones (N1 y N2), y
específicamente, que mu1<mu2. Entonces,
nuestra hipótesis estadística (hipótesis Nula), es que ambas, n1 y n2,
provienen de la misma N.
n1
Promedio = 49.99
desv.est = 6.39
error.est = 1.01
n2
Promedio = 55.39
desv.est = 7.52
error.est = 1.16
Muestro los límites de confianza para el promedio de la población,
para cada muestra (usando los grados de libertad gl
corrspondientes).
IC95 = ȳ +/- (qt() * std.error)
En ROJO
n1=40, gl= 39, q= 2.023
50 +/- (2.023* 1.01)
En AZUL
n1=42, gl= 41, q= 2.0195
55.4 +/- (2.0195* 1.16)

T TEST
Ya desde la gráfica anterior, podemos deducir que la prueba de T nos
va a decir que las muestras vienen de dos poblaciones difernetes,
N1y N2. ¿Cómo lo sabemos?
En el gráfico de arriba, podemos ver que los IC.95% para los promedios
poblacionales de cada muestra, no se traslapan. ENtonces lo que nos va a
decir la prueba de T, es que la cantidad de errores estándar que separa
a las muestras, es mayor que lo que hubiéramos esperado, si las dos
muestras hubieran salido de la misma población.
t <- t.test(n1, n2, var.equal = TRUE, )
t$statistic ## T
t
-3.494702
t$p.value ## P Value
[1] 0.0007767692
Atención a la P y a
t.
El estadístico de t= -3.495 , se ve en rojo en el gráfico de abajo.
La P = 8^{-4} es el área a la izquierda de la linea roja de
t. El área sombreada es P= 0.95, y el error alfa (1-P=0.05)
es el área sin sombrear

LS0tDQp0aXRsZTogIlRfVGVzdCINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCioqKlBhcXVldGVzIGRlIHVzbyBmcmVjdWVudGUgZW4gZWwgY3Vyc28qKiogIA0KDQpgbGlicmFyeShnZ3Bsb3QyKWAgIA0KYGxpYnJhcnkocGxvdGx5KWAgIA0KYGxpYnJhcnkodGlkeXZlcnNlKWAgIA0KDQpgYGB7ciBlY2hvPUZBTFNFfQ0KbXUxID0gICAgNTAuMg0Kc2lnbWExID0gNi44DQoNCm11MiA9ICAgIDU1LjQNCnNpZ21hMiA9ICA3LjENCmBgYCAgDQoNCmBgYHtyIGVjaG89RkFMU0V9DQpzZXQuc2VlZCgxMzQ0KQ0KbjEgPC0gcm5vcm0oIG49IDQwLCBtZWFuPSBtdTEsIHNkPSBzaWdtYTEpICAgIyMgbXVlc3RyYSBuMSANCiMjIyMjIyMjIyMjIyMNCnNldC5zZWVkKDEyMzUpDQpuMiA8LSBybm9ybSggbj0gNDIsIG1lYW49IG11Miwgc2Q9IHNpZ21hMikgICAjIyBtdWVzdHJhIG4yDQojIyMjIyMjIyMjIyMjDQpuPW4xDQplZTEgPC0gc2Qobikvc3FydChsZW5ndGgobikpDQpuPSBuMg0KZWUyIDwtIHNkKG4pL3NxcnQobGVuZ3RoKG4pKQ0KYGBgDQoNClRlbmVtb3MgZG9zIG11ZXN0cmFzLCBgbjFgIHkgYG4yYC4gUXVlcmVtb3Mgc2FiZXIgc2kgbGFzIGRvcyBwZXJ0ZW5lY2VuIGEgbGEgbWlzbWEgcG9ibGFjacOzbiAoTikuIFN1cG9uZ2Ftb3MgcXVlIG51ZXN0cm8gbWFyY28gdGXDs3JpY28gbm9zIGRhIHN1ZmljaWVudGUgZXZpZGVuY2lhIHBhcmEgcHJlZGVjaXIgKHNlZ8O6biBsYSBoaXDDs3Rlc2lzIGRlIHRyYWJham8gbyBpbnZlc3RpZ2FjacOzbikgcXVlIGxhcyBtdWVzdHJhcyBzb24gZGUgZGlmZXJlbnRlcyBwb2JsYWNpb25lcyAoYE4xYCB5IGBOMmApLCB5IGVzcGVjw61maWNhbWVudGUsIHF1ZSBgbXUxYDxgbXUyYC4gRW50b25jZXMsIG51ZXN0cmEgaGlww7N0ZXNpcyBlc3RhZMOtc3RpY2EgKGhpcMOzdGVzaXMgTnVsYSksIGVzIHF1ZSBhbWJhcywgbjEgeSBuMiwgcHJvdmllbmVuIGRlIGxhIG1pc21hIGBOYC4NCg0KKioqbjEqKiogIA0KUHJvbWVkaW8gPSBgciByb3VuZChtZWFuKG4xKSwyKWAgIA0KZGVzdi5lc3QgPSBgciByb3VuZChzZChuMSksMilgICANCmVycm9yLmVzdCA9IGByIHJvdW5kKGVlMSwyKWAgIA0KDQoqKipuMioqKiAgDQpQcm9tZWRpbyA9IGByIHJvdW5kKG1lYW4objIpLDIpYCAgDQpkZXN2LmVzdCA9IGByIHJvdW5kKHNkKG4yKSwgMilgICANCmVycm9yLmVzdCAgPSBgciByb3VuZChlZTIsMilgICANCiAgDQoNCmBgYHtyIGVjaG89RkFMU0V9DQpkYXQgPC0gZGF0YS5mcmFtZShjYmluZCh5ID0gYXMubnVtZXJpYyhjKG4xLG4yKSkgKSApDQpkYXQkeCA8LSBhcy5mYWN0b3IoYyggcmVwKCJuMSIsIDQwKSwgcmVwKCJuMiIsIDQyKSApKQ0KYGBgICANCg0KYGBge3IgZWNobz0gRkFMU0UsIHJlc3VsdHM9RkFMU0V9DQpsaWJyYXJ5KGdncGxvdDIpDQpgYGANCmBgYHtyIGVjaG89RkFMU0UscmVzdWx0cz1GQUxTRSxldmFsPUZBTFNFfQ0KZ2dwbG90KGFlcyh4PXkpLGRhdGE9ZGF0KSArIGdlb21faGlzdG9ncmFtKGFlcyh5PS4uZGVuc2l0eS4uKSwgYmlud2lkdGg9LjUsY29sb3VyPSJibGFjayIsIGZpbGw9IndoaXRlIikgKyBnZW9tX2RlbnNpdHkoYWxwaGE9LjIsIGZpbGw9IiNGRjY2NjYiKQ0KYGBgICANCg0KTXVlc3RybyBsb3MgbMOtbWl0ZXMgZGUgY29uZmlhbnphIHBhcmEgZWwgcHJvbWVkaW8gZGUgbGEgcG9ibGFjacOzbiwgcGFyYSBjYWRhIG11ZXN0cmEgKHVzYW5kbyBsb3MgZ3JhZG9zIGRlIGxpYmVydGFkIGBnbGAgY29ycnNwb25kaWVudGVzKS4NCg0KYElDOTVgID0gYMizICsvLSAocXQoKSAqIHN0ZC5lcnJvcilgDQoNCkVuIFJPSk8gIA0KYG4xYD00MCwgYGdsYD0gMzksIGBxYD0gYHIgcm91bmQocXQoIDAuOTc1LCBkZj0gMzkpLDMpYCAgDQoNCmByIHJvdW5kKG1lYW4objEpLDEpYCArLy0gIChgciByb3VuZChxdCggMC45NzUsIGRmPSAzOSksMylgKiBgciByb3VuZChlZTEsMilgKQ0KDQpFbiBBWlVMICANCmBuMWA9NDIsIGBnbGA9IDQxLCBgcWA9IGByIHJvdW5kKHF0KCAwLjk3NSwgZGY9IDQxKSw0KWAgIA0KDQpgciByb3VuZChtZWFuKG4yKSwxKWAgKy8tICAoYHIgcm91bmQocXQoIDAuOTc1LCBkZj0gNDEpLDQpYCogYHIgcm91bmQoZWUyLDIpYCkNCg0KYGBge3IgZWNobz1GQUxTRX0NCnNpemU9IDEuMg0KZ2dwbG90KGRhdGE9ZGF0LCBhZXMoeD0geCwgeT15LCBjb2w9eCkpICsgDQogIGdlb21fYm94cGxvdChtaWRkbGU9RkFMU0UsIGFscGhhPTAuNSkgKyANCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gbWVhbihuMSkgKyAocXQoIDAuOTc1LGRmPWxlbmd0aChuMSktMSkqZWUxICogKGMoLTEsMSkpKSAsICBjb2w9ICJyZWQiKSArIA0KICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSBtZWFuKG4yKSsgKHF0KCAwLjk3NSxkZj1sZW5ndGgobjIpLTEpKmVlMiAqIChjKC0xLDEpKSksIGNvbD0gImJsdWUiICkgKyANCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gbWVhbihuMSksIGNvbD0gInJlZCIsIGxpbmV0eXBlPSJkYXNoZWQiLCBzaXplPSBzaXplKSArDQogIGdlb21faGxpbmUoeWludGVyY2VwdD0gbWVhbihuMiksIGNvbD0gImJsdWUiLCBsaW5ldHlwZT0iZGFzaGVkIiwgc2l6ZT0gc2l6ZSkgKw0KICBnZW9tX3BvaW50KHBvc2l0aW9uID0gcG9zaXRpb25faml0dGVyKHdpZHRoPSAwLjIpKQ0KYGBgICANCg0KIyMjIFQgVEVTVA0KDQpZYSBkZXNkZSBsYSBncsOhZmljYSBhbnRlcmlvciwgcG9kZW1vcyBkZWR1Y2lyIHF1ZSBsYSBwcnVlYmEgZGUgVCBub3MgdmEgYSBkZWNpciBxdWUgbGFzIG11ZXN0cmFzIHZpZW5lbiBkZSBkb3MgcG9ibGFjaW9uZXMgZGlmZXJuZXRlcywgYE4xYHkgYE4yYC4gwr9Dw7NtbyBsbyBzYWJlbW9zPyAgDQpFbiBlbCAgZ3LDoWZpY28gZGUgYXJyaWJhLCBwb2RlbW9zIHZlciBxdWUgbG9zIElDLjk1JSBwYXJhIGxvcyBwcm9tZWRpb3MgcG9ibGFjaW9uYWxlcyBkZSBjYWRhIG11ZXN0cmEsIG5vIHNlIHRyYXNsYXBhbi4gRU50b25jZXMgbG8gcXVlIG5vcyB2YSBhIGRlY2lyIGxhIHBydWViYSBkZSBULCBlcyBxdWUgbGEgY2FudGlkYWQgZGUgZXJyb3JlcyBlc3TDoW5kYXIgcXVlIHNlcGFyYSBhIGxhcyBtdWVzdHJhcywgZXMgbWF5b3IgcXVlIGxvIHF1ZSBodWJpw6lyYW1vcyBlc3BlcmFkbywgc2kgbGFzIGRvcyBtdWVzdHJhcyBodWJpZXJhbiBzYWxpZG8gZGUgbGEgbWlzbWEgcG9ibGFjacOzbi4NCg0KYGBge3J9DQp0IDwtIHQudGVzdChuMSwgbjIsIHZhci5lcXVhbCA9IFRSVUUsICkNCnQkc3RhdGlzdGljICAgICAgICAgICAgICAgICAgICAgICAgICAgIyMgVA0KdCRwLnZhbHVlICAgICAgICAgICAgICAgICAgICAgICAgICAgICAjIyBQIFZhbHVlDQpgYGAgIA0KDQoqKkF0ZW5jacOzbiBhIGxhIF9gUGBfIHkgYSBfYHRgXy4qKiAgDQoNCkVsIGVzdGFkw61zdGljbyBkZSB0PSBgciByb3VuZCh0JHN0YXRpc3RpYywzKWAgLCBzZSB2ZSBlbiByb2pvIGVuIGVsIGdyw6FmaWNvIGRlIGFiYWpvLiBMYSBQID0gYHIgcm91bmQodCRwLnZhbHVlLDQpYCBlcyBlbCDDoXJlYSBhIGxhIGl6cXVpZXJkYSBkZSBsYSBsaW5lYSByb2phIGRlIGB0YC4gRWwgw6FyZWEgc29tYnJlYWRhIGVzIFA9IDAuOTUsIHkgZWwgZXJyb3IgYWxmYSAoMS1QPTAuMDUpIGVzIGVsIMOhcmVhIHNpbiBzb21icmVhcg0KIA0KYGBge3IgZWNobz1GQUxTRX0NCiMjIyMgVCBkaXN0cmlidXRpb24gUGxvdA0Kbj0gODINCm1lYW49MDsgc2Q9MQ0KZ2w9IG4tMg0KDQpsYj0gcXQocD0wLjA1LCBkZj1nbCApOyB1Yj0gcXQocD0xLCBkZj1nbCApIA0KDQp4IDwtIHNlcSgtNSw1LGxlbmd0aD0xMDApDQpoeCA8LSBkdCh4LCBkZj1nbCkNCg0KcGxvdCh4LCBoeCwgdHlwZT0gImwiLCB4bGFiPSJ0IHZhbHVlcyIsIHlsYWI9ImRlbnNpdHkiLA0KICBtYWluPSwgYXhlcz0gVFJVRSkNCg0KaSA8LSB4ID49IGxiICYgeCA8PSB1Yg0KYWJsaW5lKGg9MCkNCmxpbmVzKHgsIGh4KQ0KDQpwb2x5Z29uKGMobGIseFtpXSx1YiksIGMoMCxoeFtpXSwwKSwgY29sPSJncmF5IikNCg0KDQphcmVhIDwtIHB0KHViLCBkZj1nbCkgLSBwdChsYiwgZGY9Z2wpDQpyZXN1bHQgPC0gcGFzdGUoIlAoIixyb3VuZChsYiwzKSwiPCBJUSA8Iixyb3VuZCh1YiwzKSwiKSA9IiwNCiAgIHNpZ25pZihhcmVhLCBkaWdpdHM9MykpDQptdGV4dChyZXN1bHQsMykNCmF4aXMoMSwgYXQ9c2VxKDQwLCAxNjAsIDIwKSwgcG9zPTApDQoNCmFibGluZSh2PXQkc3RhdGlzdGljLCBjb2w9ICJyZWQiICkNCmBgYA0KDQo=