Zadanie

Zadanie: Načítajte dáta WWWusage z dát z R. Vykonajte prvotnú štatistickú analýzu + grafickú analýzu.

myData <- datasets::WWWusage
myData
## Time Series:
## Start = 1 
## End = 100 
## Frequency = 1 
##   [1]  88  84  85  85  84  85  83  85  88  89  91  99 104 112 126 138 146 151
##  [19] 150 148 147 149 143 132 131 139 147 150 148 145 140 134 131 131 129 126
##  [37] 126 132 137 140 142 150 159 167 170 171 172 172 174 175 172 172 174 174
##  [55] 169 165 156 142 131 121 112 104 102  99  99  95  88  84  84  87  89  88
##  [73]  85  86  89  91  91  94 101 110 121 135 145 149 156 165 171 175 177 182
##  [91] 193 204 208 210 215 222 228 226 222 220

Sumarna tabulka

summary(myData)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    83.0    99.0   138.5   137.1   167.5   228.0

Aritmetický priemer

mean(myData)
## [1] 137.08

Geometrický priemer

Sposob 1

library('psych')
## 
## Присоединяю пакет: 'psych'
## Следующие объекты скрыты от 'package:ggplot2':
## 
##     %+%, alpha
geometric.mean(myData)
## [1] 131.3419

Sposob 2

suc <- 1

for(i in 1:length(myData)){
  suc <- suc * myData[i] 
}

geom_mean <- (suc)^(1/length(myData))
geom_mean
## [1] 131.3419

Sposob 3

exp(mean(log(myData))) 
## [1] 131.3419

Harmonický priemer

Sposob 1

harmonic.mean(myData)
## [1] 125.7698

Sposob 2

harm_mean <- length(myData)/(sum(1/myData))
harm_mean
## [1] 125.7698

Median

median(myData)
## [1] 138.5

Kvantily(pre 0.25 a 0.75)

kv1 <- quantile(myData,0.25)
kv2 <- quantile(myData,0.75)

cat("kv1: ", kv1)
## kv1:  99
cat("\nkv2: ", kv2)
## 
## kv2:  167.5

Variacn

var(myData)
## [1] 1599.953

SD

sd(myData)
## [1] 39.99941

Dolný a hornú hranicu odĺahlých udajov

k <- 3
horna_hranica <- quantile(myData, probs=0.75) + k*IQR(myData)
cat("Horná hranica: ", horna_hranica) 
## Horná hranica:  373
dolna_hranica <- quantile(myData, probs=0.25) - k*IQR(myData)
cat("\nDolná hranica:", dolna_hranica) 
## 
## Dolná hranica: -106.5

Plot

plot(myData, type = "l", col = "blue")

Histogram

hist(myData, breaks = seq(0, max(myData)+50, 10), 
     main="Histogram", 
     xlab = "Počet návštevníkov", 
     ylab = "Pocet",
     col = "lightblue")

Boxplot

boxplot(myData, main = "Boxplot", 
        col = "lightgreen")

## Vioplot

library(vioplot)
## Загрузка требуемого пакета: sm
## Package 'sm', version 2.2-6.0: type help(sm) for summary information
## Загрузка требуемого пакета: zoo
## 
## Присоединяю пакет: 'zoo'
## Следующие объекты скрыты от 'package:base':
## 
##     as.Date, as.Date.numeric
vioplot(myData)