setwd('C:/Disk D/ND-Skola/02-Vyuka/01-Aplikovaná statistika/Data_cviceni_txt')
data<-read.delim('ocel.txt',sep='',dec='.')
head(data)
plasticita<-data$mez_plasticity
pevnost<-data$mez_pevnosti

n<-dim(data)[1]
n<-length(plasticita)
round(Sturges<-1+3.3*log10(n),2) #7 tridicich intervalu


#pevnost
range(pevnost)
dh_pev<-c(50,70,90,110,130,150,170)
hh_pev<-c(70,90,110,130,150,170,190)

#plasticita
range(plasticita)
dh_pl<-c(30,50,70,90,110,130,150)
hh_pl<-c(50,70,90,110,130,150,170)

m<-length(dh_pl)
nj<-rep(0,m)
for(j in 1:m){
  nj[j]<-sum(pevnost > dh_pev[j] & pevnost <= hh_pev[j])
}

nj                        #absolutni cetnost
pj<-nj/n                  #relativni cetnost
Nj<-cumsum(nj)            #absolutni kumulativni cetnost
Fj<-cumsum(pj)            #relativni kumulativni cetnost
xj<-(dh_pev+hh_pev)/2     #stred intervalu
dj<-hh_pev-dh_pev         #delka intervalu

var.rada<-data.frame(dh=dh_pev, hh=hh_pev, stred=xj, delka=dj, 
                     nj=nj, pj=pj, Nj=Nj,Fj=Fj)
round(var.rada,2)
hist(pevnost,xlab='pevnost oceli', ylab='',main='Histogram - pevnost oceli', col='bisque')

#Dvourozmerna intervalova data
KT<-matrix(NA,m,m)
for(i in 1:m){
  for(j in 1:m){
    KT[i,j]<-sum((plasticita > dh_pl[i] & plasticita <=hh_pl[i]) &
                   pevnost > dh_pev[j] & pevnost <= hh_pev[j])
    }
}
round(KT/n,2)
round(prop.table(KT,1),2)
round(prop.table(KT,2),2)
#====================================================================================
#Charakteristiky
#Nominalni data
library(lsr)
data<-matrix(c(13,15,14,11,7,2,19,9,10),3,3,byrow=T)
cramersV(data)
#INTERPRETACE: Mezi barvou vlasu a barvou oci existuje slaby stupen zavislosti