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