Tema 4. T de Student

pagos <- c(107,92,97,95,105,101,91,99,95,104)

t.test(pagos,y=NULL,alternative="two.sided",mu=100,paired=FALSE,var.equal=FALSE,conf.level=0.90)
## 
##  One Sample t-test
## 
## data:  pagos
## t = -0.79888, df = 9, p-value = 0.4449
## alternative hypothesis: true mean is not equal to 100
## 90 percent confidence interval:
##   95.38755 101.81245
## sample estimates:
## mean of x 
##      98.6
# Si es de una cola, alternative puede ser: "less" o "greater".

# Si p value es mayor a 0.10 (por significancia del 10%) no se rechaza H0.

# Conclusión: No se rechaza H0.
# Intervalo de Confianza: Entre 95 y 101 dólares. 

# 8-70
# Drive-a-Lemon, renta automóviles en buenas condiciones mecánicas, pero más antiguos que
# aquellos que rentan las grandes cadenas nacionales de renta de coches. Como resultado, anuncia
# que sus tarifas son considerablemente más bajas que las de sus grandes competidores. Una
# encuesta en la industria estableció que el cargo total promedio por renta en una de las mayores
# compañías es de $77.38 dólares. Una muestra aleatoria de 18 transacciones realizadas por Drive-a-
# Lemon mostró un cargo total promedio de $87.61, con una desviación estándar de la muestra de
# $19.48. Verifique que, con un nivel de significancia de 0.025, el cargo total promedio de Drive-a-
# Lemon es más alto que el de las grandes compañías. ¿Indica este resultado que las tarifas de Drive-
# a-Lemon, no son más bajas que las de las grandes cadenas nacionales? Justifique su respuesta.

# Respuesta

# Paso 1. Plantear Hipótesis
# H0: μ = xbar
# H1: μ > xbar

# Paso 2. Nivel de Confianza y Grados de Libertad
# α = 0.025
# GL = 17

# Paso 3. Zona de Aceptación/Rechazo
# t de Tabla =  2.110

# Paso 4. Fórmula
t <- (87.61 - 77.38)/(19.48/sqrt(18))
t
## [1] 2.22804
# Paso 5. Conclusión
# Se rechaza H0

# t=2.228 (x=87.61), tU=2.110 (xbarU= 87.07), de modo que se rechaza la hipótesis nula. Sin embargo,
# si Drive-a-Lemon no tiene la misma presencia a nivel nacional que las principales cadenas del país, 
# entonces una comparación de sus tasas promedio con el promedio nacional de las grandes cadenas puede
# conducir a una conclusión errónea.

Tema 5. ANOVA

# install.packages("stats")
library(stats)

# Ejercicio 1
# file.choose()
resistencia <- read.csv("C:\\Users\\raulc\\Downloads\\Resistencia (1).csv")
resistencia$Caja <- as.factor(resistencia$Caja)
qf(.95,df1=3,df2=20)
## [1] 3.098391
anova1 <- aov(Resistencia ~ Caja, data=resistencia)
summary(anova1)
##             Df Sum Sq Mean Sq F value   Pr(>F)    
## Caja         3 127375   42458   25.09 5.53e-07 ***
## Residuals   20  33839    1692                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
# Ejercicio 2
# file.choose()
manchado <- read.csv("C:\\Users\\raulc\\Downloads\\Grado.csv")
manchado$Mezcla <- as.factor(manchado$Mezcla)
qf(.99,df1=2,df2=12)
## [1] 6.926608
anova2 <- aov(Grado.de.manchado ~ Mezcla, data=manchado)
summary(anova2)
##             Df Sum Sq Mean Sq F value Pr(>F)
## Mezcla       2 0.0609 0.03043   0.987  0.401
## Residuals   12 0.3701 0.03084

Tema 6. Regresión Lineal Simple

# Ejercicio 1
x <- c(0.2,0.5,1,2,3)
y <- c(8,10,18,35,60)
regresion <- lm(y~x)
summary(regresion)
## 
## Call:
## lm(formula = y ~ x)
## 
## Residuals:
##       1       2       3       4       5 
##  3.0591 -0.5354 -1.8596 -3.5079  2.8437 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)   
## (Intercept)    1.211      2.451   0.494  0.65510   
## x             18.648      1.450  12.863  0.00101 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.341 on 3 degrees of freedom
## Multiple R-squared:  0.9822, Adjusted R-squared:  0.9763 
## F-statistic: 165.5 on 1 and 3 DF,  p-value: 0.001014
# Ejercicio 2
a <- c(-6,-3,0,3,6,9,12,15,20,25)
b <- c(2,2.8,3.9,4.2,5.8,6.2,7.5,8.2,9.3,10.9)
regresion2 <- lm(b~a)
summary(regresion2)
## 
## Call:
## lm(formula = b ~ a)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.41088 -0.12507 -0.03329  0.14807  0.32493 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 3.746694   0.101162   37.04 3.10e-10 ***
## a           0.288062   0.008087   35.62 4.22e-10 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.2438 on 8 degrees of freedom
## Multiple R-squared:  0.9937, Adjusted R-squared:  0.993 
## F-statistic:  1269 on 1 and 8 DF,  p-value: 4.223e-10

Ejercicios del Mundo Real

11-26

En un estudio se compararon los efectos de cuatro promociones mensuales sobre las ventas. A continuación, presentamos las ventas unitarias de 5 tiendas que utilizaron las 4 promociones en meses diferentes:

Tipo de Promoción Ventas Tienda 1 Ventas Tienda 2 Ventas Tienda 3 Ventas Tienda 4 Ventas Tienda 5
Muestra gratis 78 87 81 89 85
Producto Extra 94 91 87 90 88
Descuento 73 78 69 83 76
Reembolso 79 83 78 69 81

Al nivel de significancia de 0.01 ¿las promociones producen diferentes efectos sobre las ventas?

LS0tDQp0aXRsZTogIldvcmtzaG9wIDIiDQphdXRob3I6ICJSYXVsIENhbnR1IEEwMTA4NzY4MyINCmRhdGU6ICIyMDIzLTA1LTI0Ig0Kb3V0cHV0OiANCiAgaHRtbF9kb2N1bWVudDoNCiAgICB0aGVtZTogdW5pdGVkDQogICAgaGlnaGxpZ2h0OiB0YW5nbw0KICAgIHRvYzogVFJVRQ0KICAgIHRvY19mbG9hdDogVFJVRQ0KICAgIGNvZGVfZG93bmxvYWQ6IFRSVUUNCi0tLQ0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogZ29sZDsiPlRlbWEgNC4gVCBkZSBTdHVkZW50PC9zcGFuPg0KDQohW10oQzpcXFVzZXJzXFxyYXVsY1xcT25lRHJpdmVcXEVzY3JpdG9yaW9cXGF1dG8ucG5nKQ0KDQpgYGB7cn0NCnBhZ29zIDwtIGMoMTA3LDkyLDk3LDk1LDEwNSwxMDEsOTEsOTksOTUsMTA0KQ0KDQp0LnRlc3QocGFnb3MseT1OVUxMLGFsdGVybmF0aXZlPSJ0d28uc2lkZWQiLG11PTEwMCxwYWlyZWQ9RkFMU0UsdmFyLmVxdWFsPUZBTFNFLGNvbmYubGV2ZWw9MC45MCkNCg0KIyBTaSBlcyBkZSB1bmEgY29sYSwgYWx0ZXJuYXRpdmUgcHVlZGUgc2VyOiAibGVzcyIgbyAiZ3JlYXRlciIuDQoNCiMgU2kgcCB2YWx1ZSBlcyBtYXlvciBhIDAuMTAgKHBvciBzaWduaWZpY2FuY2lhIGRlbCAxMCUpIG5vIHNlIHJlY2hhemEgSDAuDQoNCiMgQ29uY2x1c2nDs246IE5vIHNlIHJlY2hhemEgSDAuDQojIEludGVydmFsbyBkZSBDb25maWFuemE6IEVudHJlIDk1IHkgMTAxIGTDs2xhcmVzLiANCg0KIyA4LTcwDQojIERyaXZlLWEtTGVtb24sIHJlbnRhIGF1dG9tw7N2aWxlcyBlbiBidWVuYXMgY29uZGljaW9uZXMgbWVjw6FuaWNhcywgcGVybyBtw6FzIGFudGlndW9zIHF1ZQ0KIyBhcXVlbGxvcyBxdWUgcmVudGFuIGxhcyBncmFuZGVzIGNhZGVuYXMgbmFjaW9uYWxlcyBkZSByZW50YSBkZSBjb2NoZXMuIENvbW8gcmVzdWx0YWRvLCBhbnVuY2lhDQojIHF1ZSBzdXMgdGFyaWZhcyBzb24gY29uc2lkZXJhYmxlbWVudGUgbcOhcyBiYWphcyBxdWUgbGFzIGRlIHN1cyBncmFuZGVzIGNvbXBldGlkb3Jlcy4gVW5hDQojIGVuY3Vlc3RhIGVuIGxhIGluZHVzdHJpYSBlc3RhYmxlY2nDsyBxdWUgZWwgY2FyZ28gdG90YWwgcHJvbWVkaW8gcG9yIHJlbnRhIGVuIHVuYSBkZSBsYXMgbWF5b3Jlcw0KIyBjb21wYcOxw61hcyBlcyBkZSAkNzcuMzggZMOzbGFyZXMuIFVuYSBtdWVzdHJhIGFsZWF0b3JpYSBkZSAxOCB0cmFuc2FjY2lvbmVzIHJlYWxpemFkYXMgcG9yIERyaXZlLWEtDQojIExlbW9uIG1vc3Ryw7MgdW4gY2FyZ28gdG90YWwgcHJvbWVkaW8gZGUgJDg3LjYxLCBjb24gdW5hIGRlc3ZpYWNpw7NuIGVzdMOhbmRhciBkZSBsYSBtdWVzdHJhIGRlDQojICQxOS40OC4gVmVyaWZpcXVlIHF1ZSwgY29uIHVuIG5pdmVsIGRlIHNpZ25pZmljYW5jaWEgZGUgMC4wMjUsIGVsIGNhcmdvIHRvdGFsIHByb21lZGlvIGRlIERyaXZlLWEtDQojIExlbW9uIGVzIG3DoXMgYWx0byBxdWUgZWwgZGUgbGFzIGdyYW5kZXMgY29tcGHDscOtYXMuIMK/SW5kaWNhIGVzdGUgcmVzdWx0YWRvIHF1ZSBsYXMgdGFyaWZhcyBkZSBEcml2ZS0NCiMgYS1MZW1vbiwgbm8gc29uIG3DoXMgYmFqYXMgcXVlIGxhcyBkZSBsYXMgZ3JhbmRlcyBjYWRlbmFzIG5hY2lvbmFsZXM/IEp1c3RpZmlxdWUgc3UgcmVzcHVlc3RhLg0KDQojIFJlc3B1ZXN0YQ0KDQojIFBhc28gMS4gUGxhbnRlYXIgSGlww7N0ZXNpcw0KIyBIMDogzrwgPSB4YmFyDQojIEgxOiDOvCA+IHhiYXINCg0KIyBQYXNvIDIuIE5pdmVsIGRlIENvbmZpYW56YSB5IEdyYWRvcyBkZSBMaWJlcnRhZA0KIyDOsSA9IDAuMDI1DQojIEdMID0gMTcNCg0KIyBQYXNvIDMuIFpvbmEgZGUgQWNlcHRhY2nDs24vUmVjaGF6bw0KIyB0IGRlIFRhYmxhID0gIDIuMTEwDQoNCiMgUGFzbyA0LiBGw7NybXVsYQ0KdCA8LSAoODcuNjEgLSA3Ny4zOCkvKDE5LjQ4L3NxcnQoMTgpKQ0KdA0KDQojIFBhc28gNS4gQ29uY2x1c2nDs24NCiMgU2UgcmVjaGF6YSBIMA0KDQojIHQ9Mi4yMjggKHg9ODcuNjEpLCB0VT0yLjExMCAoeGJhclU9IDg3LjA3KSwgZGUgbW9kbyBxdWUgc2UgcmVjaGF6YSBsYSBoaXDDs3Rlc2lzIG51bGEuIFNpbiBlbWJhcmdvLA0KIyBzaSBEcml2ZS1hLUxlbW9uIG5vIHRpZW5lIGxhIG1pc21hIHByZXNlbmNpYSBhIG5pdmVsIG5hY2lvbmFsIHF1ZSBsYXMgcHJpbmNpcGFsZXMgY2FkZW5hcyBkZWwgcGHDrXMsIA0KIyBlbnRvbmNlcyB1bmEgY29tcGFyYWNpw7NuIGRlIHN1cyB0YXNhcyBwcm9tZWRpbyBjb24gZWwgcHJvbWVkaW8gbmFjaW9uYWwgZGUgbGFzIGdyYW5kZXMgY2FkZW5hcyBwdWVkZQ0KIyBjb25kdWNpciBhIHVuYSBjb25jbHVzacOzbiBlcnLDs25lYS4NCmBgYA0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogYnJvd247Ij5UZW1hIDUuIEFOT1ZBPC9zcGFuPg0KDQohW10oQzpcXFVzZXJzXFxyYXVsY1xcT25lRHJpdmVcXEVzY3JpdG9yaW9cXGNhamFzLnBuZykNCmBgYHtyfQ0KIyBpbnN0YWxsLnBhY2thZ2VzKCJzdGF0cyIpDQpsaWJyYXJ5KHN0YXRzKQ0KDQojIEVqZXJjaWNpbyAxDQojIGZpbGUuY2hvb3NlKCkNCnJlc2lzdGVuY2lhIDwtIHJlYWQuY3N2KCJDOlxcVXNlcnNcXHJhdWxjXFxEb3dubG9hZHNcXFJlc2lzdGVuY2lhICgxKS5jc3YiKQ0KcmVzaXN0ZW5jaWEkQ2FqYSA8LSBhcy5mYWN0b3IocmVzaXN0ZW5jaWEkQ2FqYSkNCnFmKC45NSxkZjE9MyxkZjI9MjApDQphbm92YTEgPC0gYW92KFJlc2lzdGVuY2lhIH4gQ2FqYSwgZGF0YT1yZXNpc3RlbmNpYSkNCnN1bW1hcnkoYW5vdmExKQ0KDQojIEVqZXJjaWNpbyAyDQojIGZpbGUuY2hvb3NlKCkNCm1hbmNoYWRvIDwtIHJlYWQuY3N2KCJDOlxcVXNlcnNcXHJhdWxjXFxEb3dubG9hZHNcXEdyYWRvLmNzdiIpDQptYW5jaGFkbyRNZXpjbGEgPC0gYXMuZmFjdG9yKG1hbmNoYWRvJE1lemNsYSkNCnFmKC45OSxkZjE9MixkZjI9MTIpDQphbm92YTIgPC0gYW92KEdyYWRvLmRlLm1hbmNoYWRvIH4gTWV6Y2xhLCBkYXRhPW1hbmNoYWRvKQ0Kc3VtbWFyeShhbm92YTIpDQpgYGANCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6IHJlZDsiPlRlbWEgNi4gUmVncmVzacOzbiBMaW5lYWwgU2ltcGxlPC9zcGFuPg0KDQohW10oQzpcXFVzZXJzXFxyYXVsY1xcT25lRHJpdmVcXEVzY3JpdG9yaW9cXGp1cmFzc2ljLmdpZikNCmBgYHtyfQ0KIyBFamVyY2ljaW8gMQ0KeCA8LSBjKDAuMiwwLjUsMSwyLDMpDQp5IDwtIGMoOCwxMCwxOCwzNSw2MCkNCnJlZ3Jlc2lvbiA8LSBsbSh5fngpDQpzdW1tYXJ5KHJlZ3Jlc2lvbikNCg0KIyBFamVyY2ljaW8gMg0KYSA8LSBjKC02LC0zLDAsMyw2LDksMTIsMTUsMjAsMjUpDQpiIDwtIGMoMiwyLjgsMy45LDQuMiw1LjgsNi4yLDcuNSw4LjIsOS4zLDEwLjkpDQpyZWdyZXNpb24yIDwtIGxtKGJ+YSkNCnN1bW1hcnkocmVncmVzaW9uMikNCmBgYA0KDQojIEVqZXJjaWNpb3MgZGVsIE11bmRvIFJlYWwNCg0KIyMgKioxMS0yNioqDQpFbiB1biBlc3R1ZGlvIHNlIGNvbXBhcmFyb24gbG9zIGVmZWN0b3MgZGUgY3VhdHJvIHByb21vY2lvbmVzIG1lbnN1YWxlcyBzb2JyZSBsYXMgdmVudGFzLiBBIGNvbnRpbnVhY2nDs24sIHByZXNlbnRhbW9zIGxhcyB2ZW50YXMgdW5pdGFyaWFzIGRlIDUgdGllbmRhcyBxdWUgdXRpbGl6YXJvbiBsYXMgNCBwcm9tb2Npb25lcyBlbiBtZXNlcyBkaWZlcmVudGVzOg0KDQpUaXBvIGRlIFByb21vY2nDs24gfCBWZW50YXMgVGllbmRhIDEgfCBWZW50YXMgVGllbmRhIDIgfCBWZW50YXMgVGllbmRhIDMgfCBWZW50YXMgVGllbmRhIDQgfCBWZW50YXMgVGllbmRhIDUNCi0tLS0tLS0tLS0tLS0gfCAtLS0tLS0tLS0tLS0tIHwgLS0tLS0tLS0tLS0tLSB8IC0tLS0tLS0tLS0tLS0gfCAtLS0tLS0tLS0tLS0tIHwgLS0tLS0tLS0tLS0tLQ0KTXVlc3RyYSBncmF0aXMgfCA3OCB8IDg3IHwgODEgfCA4OSB8IDg1DQpQcm9kdWN0byBFeHRyYSB8IDk0IHwgOTEgfCA4NyB8IDkwIHwgODgNCkRlc2N1ZW50byB8IDczIHwgNzggfCA2OSB8IDgzIHwgNzYNClJlZW1ib2xzbyB8IDc5IHwgODMgfCA3OCB8IDY5IHwgODENCkFsIG5pdmVsIGRlIHNpZ25pZmljYW5jaWEgZGUgMC4wMSDCv2xhcyBwcm9tb2Npb25lcyBwcm9kdWNlbiBkaWZlcmVudGVzIGVmZWN0b3Mgc29icmUgbGFzIHZlbnRhcz8NCg0K