data(iris)
dat <- iris
dat[, "dĺžka"] <- dat[, "Sepal.Length"]
dat[, "šírka"] <- dat[, "Sepal.Width"]
dat <- dat[, !(names(dat) %in% c("Sepal.Length", "Sepal.Width"))]
str(dat)
## 'data.frame': 150 obs. of 5 variables:
## $ Petal.Length: num 1.4 1.4 1.3 1.5 1.4 1.7 1.4 1.5 1.4 1.5 ...
## $ Petal.Width : num 0.2 0.2 0.2 0.2 0.2 0.4 0.3 0.2 0.2 0.1 ...
## $ Species : Factor w/ 3 levels "setosa","versicolor",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ dĺžka : num 5.1 4.9 4.7 4.6 5 5.4 4.6 5 4.4 4.9 ...
## $ šírka : num 3.5 3 3.2 3.1 3.6 3.9 3.4 3.4 2.9 3.1 ...
summary(dat)
## Petal.Length Petal.Width Species dĺžka
## Min. :1.000 Min. :0.100 setosa :50 Min. :4.300
## 1st Qu.:1.600 1st Qu.:0.300 versicolor:50 1st Qu.:5.100
## Median :4.350 Median :1.300 virginica :50 Median :5.800
## Mean :3.758 Mean :1.199 Mean :5.843
## 3rd Qu.:5.100 3rd Qu.:1.800 3rd Qu.:6.400
## Max. :6.900 Max. :2.500 Max. :7.900
## šírka
## Min. :2.000
## 1st Qu.:2.800
## Median :3.000
## Mean :3.057
## 3rd Qu.:3.300
## Max. :4.400
hist(dat$šírka, main = "Rozdelenie šírky kališných lístkov", xlab = "šírka", col = "skyblue")

hist(dat$dĺžka, main = "Rozdelenie dĺžky kališných lístkov", xlab = "dĺžka", col = "lightgreen")

means <- aggregate(. ~ Species, data = dat[, c("šírka", "dĺžka", "Species")], FUN = mean)
means
## Species šírka dĺžka
## 1 setosa 3.428 5.006
## 2 versicolor 2.770 5.936
## 3 virginica 2.974 6.588
boxplot(šírka ~ Species, data = dat, main = "Rozdelenie šírky podľa druhov", xlab = "Typy", ylab = "šírka", col = c("skyblue", "lightgreen", "lightcoral"))

boxplot(dĺžka ~ Species, data = dat, main = "Rozdelenie dĺžky podľa druhov", xlab = "Typy", ylab = "dĺžka", col = c("skyblue", "lightgreen", "lightcoral"))

plot(dat$dĺžka, dat$dĺžka, main = "Vzťah medzi šírkou a dĺžkou kališného lístka", xlab = "šírka", ylab = "dĺžka", col = "blue")

plot(dĺžka ~ šírka, data = dat, col = as.numeric(dat$Species), pch = as.numeric(dat$Species), main = "Vzťah medzi šírkou a dĺžkou kališného lístka podľa druhu", xlab = "šírka", ylab = "dĺžka")
legend("topright", legend = levels(dat$Species), col = 1:length(levels(dat$Species)), pch = 1:length(levels(dat$Species)))
