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
library(datos)
library(modeest)
library(descriptr)
library(knitr)
library(fdth)
## 
## Attaching package: 'fdth'
## The following object is masked from 'package:modeest':
## 
##     mfv
## The following objects are masked from 'package:stats':
## 
##     sd, var

Selección de la muestra de datos

set.seed(6094)
datos.d <- sample_n(diamantes,size = 120,replace = FALSE)

Medidas de tendencia central para la variable quilate

media_a <- mean(datos.d$quilate)
mediana <-  median(datos.d$quilate)
moda <- mfv(datos.d$quilate)


df <- data.frame (media_a, mediana, moda)

df
##     media_a mediana moda
## 1 0.7814167   0.655  0.3

Medidas de variabilidad para la variable profundidad

*Rango

rango_diam <- max(datos.d$profundidad) -min(datos.d$profundidad)
cat("El rango de la profundidad de los diamantes es: ",rango_diam)
## El rango de la profundidad de los diamantes es:  11.1

*Varianza

var_diam <-  var(datos.d$profundidad)
cat("la varianza de la  profundidad de los diamantes es: ",var_diam)
## la varianza de la  profundidad de los diamantes es:  1.906856

*La desviación estándar

ds_diam <- sqrt(var(datos.d$profundidad))
cat("La desviación estandar de la profundidad de los diamantes es: ",ds_diam)
## La desviación estandar de la profundidad de los diamantes es:  1.38089

*Coeficiente de variación:

CV <-  (sd(datos.d$profundidad) / mean(datos.d$profundidad))*100
cat("El coeficiente de variación de la profundidad de los diamantes es: ",CV)
## El coeficiente de variación de la profundidad de los diamantes es:  2.237678

*Medidas de Posición

percentiles <- quantile(datos.d$tabla, probs = c(0.25, 0.3, 0.75, 0.78, 0.9), type= 6)
percentiles
## 25% 30% 75% 78% 90% 
##  56  56  59  59  60

*Medidas de forma

asim_diam <- ds_skewness(datos.d$x)
curt_diam <- ds_kurtosis(datos.d$x)

cat("El coeficiente de asimetría de la variable x de los diamantes es: ", asim_diam) 
## El coeficiente de asimetría de la variable x de los diamantes es:  0.2836175
cat("El coeficiente de curtosis de la variable x de los diamanates es: ",curt_diam)
## El coeficiente de curtosis de la variable x de los diamanates es:  -0.9827952

*Resumen

resumen_precio.datos.d <- ds_tidy_stats(datos.d, profundidad)
kable(resumen_precio.datos.d)
vars min max mean t_mean median mode range variance stdev skew kurtosis coeff_var q1 q3 iqrange
profundidad 56.5 67.6 61.71083 61.70556 61.75 61.6 11.1 1.906856 1.38089 0.0955529 3.532202 2.237678 60.9 62.425 1.525

*Tabla de frecuencias e histograma para la variable y

tabla1 <- fdt(x =datos.d$y,k =10)

tabla1
##   Class limits  f   rf rf(%)  cf  cf(%)
##   [3.95,4.407) 19 0.16 15.83  19  15.83
##  [4.407,4.863) 22 0.18 18.33  41  34.17
##  [4.863,5.319) 12 0.10 10.00  53  44.17
##  [5.319,5.776)  9 0.07  7.50  62  51.67
##  [5.776,6.232)  6 0.05  5.00  68  56.67
##  [6.232,6.689) 31 0.26 25.83  99  82.50
##  [6.689,7.145) 10 0.08  8.33 109  90.83
##  [7.145,7.601)  5 0.04  4.17 114  95.00
##  [7.601,8.058)  3 0.03  2.50 117  97.50
##  [8.058,8.514)  3 0.03  2.50 120 100.00
plot(tabla1, col=c("blue4","blue3","blue2","blue", "deepskyblue2","deepskyblue4","dodgerblue", "dodgerblue4", "dodgerblue3", "royalblue4"),type = "fh", main = "Histograma para la variable y")

*Tabla de frecuencias, diagrama de barras, y diagrama de torta para la variable corte

tabla2 <- table(datos.d$corte)

tabla2
## 
##   Regular     Bueno Muy bueno   Premium     Ideal 
##         4         7        29        37        43
barplot(tabla2 , col = c("dodgerblue3", "royalblue4","dodgerblue","blue4", "dodgerblue4" ),
        horiz=FALSE, density= NULL,ylab= "Frecuencia absoluta" , border = TRUE,
        ylim=c(0,75),xlab="Tipos de corte",main=" Gráfico de barras para la variable corte")

pie(tabla2,labels = names(tabla2),
    density=NULL,edges = 500, radius = 0.8,
    col=c("dodgerblue3", "royalblue4","dodgerblue","blue4", "dodgerblue4"),
    clockwise = TRUE, main = "Gráfico de torta para la variable corte ",
    border = "blue")