datos<-read.table("D:/Informacion/Desktop/biotinidasa/Biotinidasa Final.txt",header=T)

Analisis descriptivo muestra total

library(knitr)
## Warning: package 'knitr' was built under R version 3.4.4
x <- summary(datos)
kable(x)
Genero Edad Bioitin Amilas
F:52 Min. : 0.08 Min. : 4.920 Min. :162.0
M:50 1st Qu.: 2.00 1st Qu.: 7.293 1st Qu.:177.9
NA Median : 5.00 Median : 8.225 Median :185.9
NA Mean : 4.76 Mean : 8.201 Mean :186.5
NA 3rd Qu.: 7.00 3rd Qu.: 9.030 3rd Qu.:192.1
NA Max. :11.00 Max. :13.460 Max. :227.0
NA NA NA NA’s :5 
par(mfrow = c(1,2))
hist(datos$Bioitin, ylab="Frecuencia", xlab= "Act Biotinidasa", main="Histograma de Frecuencia")
boxplot(datos$Bioitin, ylab="Act Biotinidasa", xlab= "Muestra Total")

Pruebas de Normalidad de los datos

Ho= La población está distribuida normalmente Ha= La poblacion no esta distribuidad normalmente

shapiro.test(datos$Bioitin)
## 
##  Shapiro-Wilk normality test
## 
## data:  datos$Bioitin
## W = 0.9804, p-value = 0.1343
qqnorm(datos$Bioitin)
qqline(datos$Bioitin)

Dado que el valor de p-value es mayor que alfa, no es posible rechazar la hipotesis nula y se concluye que los datos siguen una distribución normal.

Analisis descriptivo por genero

library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
Masculino<- filter(datos,Genero=="M")
## Warning: package 'bindrcpp' was built under R version 3.4.4
Femenino<- filter(datos,Genero=="F")

summary(Masculino$Bioitin)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   4.920   6.968   7.950   7.922   8.685  11.620
summary(Femenino$Bioitin)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   4.990   7.545   8.520   8.469   9.315  13.460
par(mfrow = c(2,2))
hist(Masculino$Bioitin, ylab="Frecuencia", xlab= "Act Biotinidasa Masculino", main="Histograma de Frecuencia")
hist(Femenino$Bioitin, ylab="Frecuencia", xlab= "Act Biotinidasa Femenino", main="Histograma de Frecuencia")
boxplot(Masculino$Bioitin, ylab="Act Biotinidasa", xlab= "Masculino")
boxplot(Femenino$Bioitin, ylab="Act Biotinidasa", xlab= "Femenino")

Analisis Descriptivo por Edad

Jardin<-filter(datos,Edad<=3)

Prescolar<-filter(datos,Edad>3&Edad<=7)

Primaria<-filter(datos,Edad>7)

summary(Jardin$Bioitin)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   5.030   7.317   8.305   8.395   9.537  12.720
summary(Prescolar$Bioitin)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   4.920   7.335   8.060   8.178   9.090  13.460
summary(Primaria$Bioitin)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   4.990   6.875   8.070   7.844   8.755  10.700
par(mfrow = c(1,3))
hist(Jardin$Bioitin, ylab="Frecuencia", xlab= "Act Biotinidasa Jardin", main="Histograma de Frecuencia")
hist(Prescolar$Bioitin, ylab="Frecuencia", xlab= "Act Biotinidasa Preescolar", main="Histograma de Frecuencia")
hist(Primaria$Bioitin, ylab="Frecuencia", xlab= "Act Biotinidasa Primaria", main="Histograma de Frecuencia")

par(mfrow = c(1,3))
boxplot(Jardin$Bioitin, ylab="Act Biotinidasa", xlab= "Jardin")
boxplot(Prescolar$Bioitin, ylab="Act Biotinidasa", xlab= "Preescolar")
boxplot(Primaria$Bioitin, ylab="Act Biotinidasa", xlab= "Primaria")

Pruebas de Hipotesis

Se plantean dos hipotesis

Ho= La actividad enzimatica de los pacientes con deficiencia enzimatica es menos o igual a de 4.4 Ha= La actividad enzimatica de los pacienes sanos es mayor a 4.4

t.test(datos$Bioitin,mu=4.4, alternative = "less")
## 
##  One Sample t-test
## 
## data:  datos$Bioitin
## t = 24.752, df = 101, p-value = 1
## alternative hypothesis: true mean is less than 4.4
## 95 percent confidence interval:
##      -Inf 8.455806
## sample estimates:
## mean of x 
##  8.200882