# Ejercicio 1: prueba de hipotesis de dos medias independientes
planT=c(75,72,73,76,78,72,80,74,76,75)
planN=c(80,83,84,78,79,81,84,85,78,56)
# parte 1-varianzas
var(planT)= 6.544444
var(planN)= 70.84444
var.test(planN,planT)
F= 10.825
p-value = 0.001534
# parte 2
t.test(planN,planT,mu=0,alternative = "g",var.equal = F,paired = F)
t= 1.33
p.valor=0.1056
# ejercicio 2 Prueba de hipotesis de dos medias dependientes
compañia1=c(10,12.2,15.3,9.6,8.6,9.4,12.5,7.3,9.4,8.7,9.1)
compañia2=c(9.2,10,9.2,10.5,9.5,8.4,7.2,8.4,10.5,6.2,8.1)
t.test(compañia1, compañia2,mu=0,alternative="g",paired=T)
t=1.7915
p.valor= 0.05174
# ejercicio 3: Muestras independientes con varianza poblacional desconocida
accion1=c(5.6,7.2,6.3,6.3,7.1,8.2,7.9,5.3,6.2,6.2)
# Parte 1
accion2=c(7.5,7.3,6.2,8.3,8.2,8.0,8.1,7.3,5.9,5.3)
var(accion1)=0.8934444
var(accion2)=1.118778
var.test(accion2,accion1)
F = 1.2522
p-value = 0.7431
# Parte 2
t.test(accion1,accion2,mu=0,alternative = "l",var.equal = T,paired = F)
t= -1.293
p.valor= 0.1062
# Ejercicio 4 : Prueba de hipotesis de dos medias pareadas o dependientes
antes=c(25,25,27,44,30,67,53,53,52)
despues=c(27,29,37,56,46,82,57,80,61)
t.test(antes,despues,mu=0,alternative = "l",paired = T)
t = -4.2514
p-value = 0.001397
# ejercicio de examen
despues=c(35.44,35.17,31.07,31.57,26.48,23.11,25.18,32.39)
antes=c(33.76,34.30,29.55,30.90,24.92,21.78,24.30,31.25)
t.test(antes,despues,mu=0,alternative = "l",paired = T)
#Ejercicio 2 del examen
# Ejercicio 1: prueba de hipotesis de dos medias independientes
planT=c(75,72,73,76,78,72,80,74,76,75)
planN=c(80,83,84,78,79,81,84,85,78,56)
# parte 1-varianzas
var(planT)= 6.544444
var(planN)= 70.84444
var.test(planN,planT)
F= 10.825
# parte 2
t.test(planN,planT,mu=0,alternative = "g",var.equal = F,paired = F)
p-value = 0.001534
Rating1=c(78,86,86,80,92,89,87,91,85,78,79,90)
Rating2=c(83,80,78,82,84,76,77,81,81,78)
var(Rating1)=26.44697
var(Rating2)=7.111111
t.test(Rating1,Rating2,mu=0,alternative = "g",var.equal = F,paired = F)
var.test(Rating1,Rating2)
library(readxl)
datos <- read_excel("C:/Users/ALEXIA/Desktop/datosanova.xlsx")
View(datos)
attach(datos)
# primero analisis exploratorio
boxplot(Venta~Tratamientos)
#hallar las medias por cada tratamiento
tapply(Venta,Tratamientos,mean)
# la desviacion estandar por cada tratamiento
tapply(Venta,Tratamientos,sd )
# Anova de una Via
anva=aov(Venta~Tratamientos)
summary(anva)
# Prueba de post hoc(tukey)
TukeyHSD(anva)
plot(TukeyHSD(anva))
# distribucion multinomial
# calculo del pvalor
pchisq(23.156,df=3,lower.tail = FALSE)
# PRUEBA DE INDEPENDENCIA- CHICUADRADO
# pvalor
pchisq(1.41959496744289,df=2,lower.tail = FALSE)
# EJEMPLO CON SOLAMENTE R
matriz1 = matrix(c(12,18,6,14,23,14), nrow = 3)
matriz1
resultado = chisq.test(matriz1)
resultado$expected
# Ejemplo 2 con R.
matriz2=matrix(c(40,60,55,45,45,55,50,50),nrow=2)
matriz2
chisq.test(matriz2)
# ejemplo 3
matriz3 = matrix(c(20,40,15,11,8,6), nrow = 3)
matriz3
chisq.test(matriz3)
# ejercicio 3 de la practica
matriz4 = matrix(c(83,52,63,62,71,58,37,49,63), nrow = 3)
matriz4
chisq.test(matriz4)
# kruskal wallis
library(readxl)
datos_KW <- read_excel("C:/Users/ALEXIA/Desktop/datos K.xlsx")
View(datos_KW)
library(agricolae)
attach(datos_KW)
salida=kruskal(disenos,tratamientos,group = F)
salida
salida2=kruskal(disenos,tratamientos,group =T)
salida2
# Ejemplo 12 de la ficha
library(readxl)
datos_KW2 <- read_excel("C:/Users/ALEXIA/Desktop/datos K.xlsx",
sheet = "pre 12 ")
View(datos_KW2)
attach(datos_KW2)
salida=kruskal(y,sector,group = F)
salida
salida2=kruskal(y,sector,group =T)
salida2
