# Beremo podatke z ukazom read.table ter jih preimenujemo v Naloga2 z ukazom <-
Naloga2 <- read.table("VAJE 1/SA_Vaje1_Naloga2.csv", header=TRUE, sep=",")

# Izpiše nam prvih 6 vrstic podatkov, skupaj z imeni
head(Naloga2)
##   Donosnost
## 1       3.4
## 2       4.1
## 3       6.6
## 4      10.0
## 5      10.3
## 6      10.4
#Če izberemo ukaz print namesto head, nam izpiše celotne podatke
#print(Naloga2)
#Za opisno statistiko uporabimo funkcijo stat.desc (predhodno zahteva ukaz library(pastecs)). Če zapišemo stat.desc(Naloga2) nam izpiše opisno statistiko za vse spremenljivke, če zapišemo stat.desc(Naloga2$Donosnost) pa samo za določeno spremenljivko. Ker je privzeta nastavitev te funkcije 5 decimalnih mest, preimenujemo ukaz v Opisna, nato pa z ukazom round dosežemo, da nam izpiše samo dve decimalki.
Opisna <- stat.desc(Naloga2$Donosnost)
round(Opisna, 2)
##      nbr.val     nbr.null       nbr.na          min          max        range 
##        46.00         0.00         0.00         3.40        41.30        37.90 
##          sum       median         mean      SE.mean CI.mean.0.95          var 
##       837.50        15.35        18.21         1.34         2.70        82.86 
##      std.dev     coef.var 
##         9.10         0.50
#Z ukazom boxplot izrišemo škatlasti diagram, s podukazom ylab pa poimenujemo y-os.
boxplot(Naloga2, 
        ylab="Donosnost")

#Z ukazom hist narišemo histogram, s podukazom breaks določimo prikazani interval ter kako podrobno razbijemo stolpce, s podukazom main poimenujemo histogram, s podukazi xlab in ylab pa poimenujemo osi.
hist (Naloga2$Donosnost, 
      breaks = seq(from = 0, to = 50, by = 10), 
      main = "Histogram donosnosti", 
      xlab = "Donosnost", 
      ylab = "Frekvenca")