# Materialy do zajec ze statystyk opisowych
# http://www.biecek.pl/semestr/stats

#
# ---------------------------------------------------------------
# zaladujmy zbior danych
library(coin)
data(glioma)

# podsumowanie
summary(glioma$sex)
summary(glioma$age)

# tablica wielodzielcza, lub jak kto woli, tabica kontyngencji
table(glioma$sex, glioma$histology)

# rozstep zmiennej
range(glioma$age)

# histogram, czyli przyklad sitowego estymatora gestosci ;-)
hist(glioma$age, 10)

# jadrowy estymator gestosci, czyli wygladzony histogram
plot(density(glioma$age))

# boxplot, czyli wykres pudelko wasy
boxplot(time~histology, data = glioma)
boxplot(glioma$time, glioma$age)

# vioplot, nazywany tez wykresem skrzypcowym
vioplot(glioma$time[glioma$histology=="GBM"],
  glioma$time[glioma$histology=="Grade3"],glioma$time)

# wykres rozrzutu
scatterplot(glioma$age,glioma$time,groups=glioma$group,smooth=F

# wykres mozaikowy
mosaicplot(~sex+histology+group,glioma)