lesson2

enter the data

 x<-c(20,37,29,40,20,26,22,23,40,45,30,20,33,28,23,47,50,46,33,22,28,33,46,24,36,25,50,33,35,28) # Ввод данных
 # просмотр данных

Ввод и просмотр данных

 x

##  [1] 20 37 29 40 20 26 22 23 40 45 30 20 33 28 23 47 50 46 33 22 28 33 46 24 36
## [26] 25 50 33 35 28

 print(x)

##  [1] 20 37 29 40 20 26 22 23 40 45 30 20 33 28 23 47 50 46 33 22 28 33 46 24 36
## [26] 25 50 33 35 28

 head(x)

## [1] 20 37 29 40 20 26

 str(x)

##  num [1:30] 20 37 29 40 20 26 22 23 40 45 ...

 # View(x)
 # fix(x)

Простейшие статистики

 x_mean <- mean(x); x_mean

## [1] 32.4

 x_median <- median(x); x_median

## [1] 31.5

 x_min <- min(x); x_min

## [1] 20

 x_max <- max(x); x_max

## [1] 50

 x_range <- range(x); x_range

## [1] 20 50

 x_rank <- rank(x); x_rank

##  [1]  2.0 22.0 14.0 23.5  2.0 10.0  4.5  6.5 23.5 25.0 15.0  2.0 17.5 12.0  6.5
## [16] 28.0 29.5 26.5 17.5  4.5 12.0 17.5 26.5  8.0 21.0  9.0 29.5 17.5 20.0 12.0

 x_moda <- which.max(table(x)); x_moda

## 33 
## 10

Встроенные команды R для описательных статистик

  var(x)

## [1] 90.17931

  quantile(x)

##    0%   25%   50%   75%  100% 
## 20.00 24.25 31.50 39.25 50.00

  IQR(x)

## [1] 15

R функций для вычисления нескольких описательных характеристик -summary(), fivenum()

  summary(x)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   20.00   24.25   31.50   32.40   39.25   50.00

  fivenum(x)

## [1] 20.0 24.0 31.5 40.0 50.0

Возможности пакета -psych

#install.packages("physch")
#library(physch)
#describe(x)

Визуализация данных

plot(x, co1="red", pch=19, type = "p")

## Warning in plot.window(...): "co1" is not a graphical parameter

## Warning in plot.xy(xy, type, ...): "co1" is not a graphical parameter

## Warning in axis(side = side, at = at, labels = labels, ...): "co1" is not a
## graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "co1" is not a
## graphical parameter

## Warning in box(...): "co1" is not a graphical parameter

## Warning in title(...): "co1" is not a graphical parameter

hist(x, col="lightgreen")

boxplot(x, col="lightblue")