HW c3

library(vioplot)

## Warning: пакет 'vioplot' был собран под R версии 4.3.3

## Загрузка требуемого пакета: sm

## Warning: пакет 'sm' был собран под R версии 4.3.3

## Package 'sm', version 2.2-6.0: type help(sm) for summary information

## Загрузка требуемого пакета: zoo

## Warning: пакет 'zoo' был собран под R версии 4.3.3

## 
## Присоединяю пакет: 'zoo'

## Следующие объекты скрыты от 'package:base':
## 
##     as.Date, as.Date.numeric

data("WWWusage")

min_val <- min(WWWusage)
max_val <- max(WWWusage)
mean_val <- mean(WWWusage)
median_val <- median(WWWusage)
sd_val <- sd(WWWusage)
var_val <- var(WWWusage)
summary_stats <- summary(WWWusage)

cat("Minimum:", min_val, "\n")

## Minimum: 83

cat("Maximum:", max_val, "\n")

## Maximum: 228

cat("Mean:", mean_val, "\n")

## Mean: 137.08

cat("Median:", median_val, "\n")

## Median: 138.5

cat("Standard Deviation:", sd_val, "\n")

## Standard Deviation: 39.99941

cat("Variance:", var_val, "\n")

## Variance: 1599.953

cat("Summary:\n")

## Summary:

print(summary_stats)

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

hist(WWWusage, main="Histogram", xlab="Používanie", ylab="Frekvencia", col="lightblue", border="black")

boxplot(WWWusage, main="Boxplot", ylab="Používanie")

plot(WWWusage, type="l", col="blue", lwd=2, main="V priebehu času", xlab="Čas", ylab="Používanie")

q25 <- quantile(WWWusage, probs=0.25)
q50 <- quantile(WWWusage, probs=0.50)
q75 <- quantile(WWWusage, probs=0.75)
iqr_val <- IQR(WWWusage)
iqr_val_2 <- q75-q25

cat("25th Percentile (Q1):", q25, "\n")

## 25th Percentile (Q1): 99

cat("50th Percentile (Median):", q50, "\n")

## 50th Percentile (Median): 138.5

cat("75th Percentile (Q3):", q75, "\n")

## 75th Percentile (Q3): 167.5

cat("Interquartile Range (IQR):", iqr_val, "\n")

## Interquartile Range (IQR): 68.5

cat("Interquartile Range (IQR) 2:", iqr_val_2, "\n")

## Interquartile Range (IQR) 2: 68.5

k <- 0.5
upper_bound <- q75 + k * iqr_val
lower_bound <- q25 - k * iqr_val

outliers <- WWWusage[WWWusage > upper_bound | WWWusage < lower_bound]
cat("Outliers:\n")

## Outliers:

print(outliers)

## [1] 204 208 210 215 222 228 226 222 220

boxplot(WWWusage, main="Boxplot", ylab="Používanie")

plot(WWWusage, type="b", col="darkgreen", pch=19, main="V priebehu času", xlab="Čas", ylab="Používanie")