eRko

Posledná ôsma úloha v rámci predmetu VIKMB42 Nástroje a metódy dátovej analytiky.

# — — — Kod — — — — #

# Nacteme hdp ze souboru
hdp <- read.csv(“hdp.csv”)

# Vytvorime si v tabulce ‘hdp’ dalsi sloupecek ‘koruny’
hdp$koruny <- (hdp$dolary * 22)

# Vykreslime graf s popsanymi osami
plot(
hdp$koruny,
main = “HDP svetovych statu”,
xlab = “Staty”,
ylab = “HDP v milionech korun”,
pch = 1, # Tvar bodu
col = “black” # Barva bodu
)

hdp.prumer <- mean(hdp$koruny)
hdp.median <- median(hdp$koruny)

# Ted do grafu vykreslime dve cary,
# modrou pro prumer a cervenou pro median
abline(a = hdp.prumer, b = 1, col = “blue”)
abline(a = hdp.median, b = 1, col = “red”)

# \n znamena ‘vypsat novy radek’
message(“\n”)

# nrow nam spocita pocet radku v data framu
pocet.statu <- nrow(hdp)

message(“Celkem mame udaje pro “, pocet.statu, “ statu.\n”)

# — — — — — — — — — — -
# Statistiky pro prumer

# floor nam zaokrouhli cislo, ve vystupu to vypada lepe

hdp.prumer <- floor(hdp.prumer)

message(“Prumer svetoveho HDP: “, hdp.prumer, “ milionu Kc.”)

# nrow nam vrati pocet radku
pod.prumerem <- sum(hdp.prumer > hdp$koruny)

nad.prumerem <- sum(hdp.prumer < hdp$koruny)

message(“Nad prumerem svetoveho HPD je “, nad.prumerem, “ statu.”)
message(“Pod prumerem svetoveho HPD je “, pod.prumerem, “ statu.”)

jmeno.statu.prumer <- hdp[nad.prumerem,]$stat

message(“To znamena, ze prumerne se staty maji jako ‘“, jmeno.statu.prumer, “‘\n”)

# — — — — — — — — — — -
# Statistiky pro median
hdp.median <- median(hdp$koruny)

message(“Median svetoveho HDP: “, hdp.median, “ milionu Kc.”)

pod.medianem <- sum(hdp.median > hdp$koruny)

nad.medianem <- sum(hdp.median < hdp$koruny)

message(“Nad medianem svetoveho HPD je “, nad.medianem, “ statu.”)
message(“A pod medianem svetoveho HPD je “, pod.medianem, “ statu.”)

jmeno.statu.median <- hdp$stat[hdp.median == hdp$koruny]

message(“To znamena, ze dle medianu stredni stat urcuje ‘“, jmeno.statu.median, “‘\n”)