library(tidyverse)
## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.2 ──
## ✔ ggplot2 3.3.6      ✔ purrr   0.3.4 
## ✔ tibble  3.1.8      ✔ dplyr   1.0.10
## ✔ tidyr   1.2.1      ✔ stringr 1.4.1 
## ✔ readr   2.1.2      ✔ forcats 0.5.2 
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
library(summarytools)
## 
## Attaching package: 'summarytools'
## 
## The following object is masked from 'package:tibble':
## 
##     view
n=1000
a1=rnorm(n, 10, 1)
a2=rnorm(n, 10, 1)
a3=rnorm(n, 10, 1)
a4=rnorm(n, 10, 1)
a5=rnorm(n, 10, 1)
a6=rnorm(n, 10, 1)
cable=apply(data.frame(a1,a2,a3,a4,a5,a6), 1, min)*6
summarytools::descr(cable)
## Descriptive Statistics  
## cable  
## N: 1000  
## 
##                       cable
## ----------------- ---------
##              Mean     52.56
##           Std.Dev      4.05
##               Min     36.37
##                Q1     50.11
##            Median     52.68
##                Q3     55.38
##               Max     64.46
##               MAD      3.93
##               IQR      5.26
##                CV      0.08
##          Skewness     -0.34
##       SE.Skewness      0.08
##          Kurtosis      0.40
##           N.Valid   1000.00
##         Pct.Valid    100.00
hist(cable)

dnorm(70, 10, 1)
## [1] 0