setwd('C:/Disk D/ND-Skola/02-Vyuka/01-Aplikovaná statistika/Data_cviceni_txt')
data<-read.delim('znamky.txt',sep='',dec='.')
head(data)
math<-data$math
english<-data$english

#variacni rada pro znamky z matematiky
n<-length(math)
nj<-rep(0,4)
for(j in 1:4){
  nj[j]<-sum(math==j)
}
(pj<-nj/n)
(Nj<-cumsum(nj))
(Fj<-cumsum(pj))
(Var.rada.M<-data.frame(nj=nj,pj=pj,Nj=Nj,Fj=Fj, 
                        row.names=c('Vyborne','Velmi_dobre',
                                    'Prospel','Neprospel')))

#-------------------------------------------------------------------------------------------
barplot(nj,col='olivedrab1', names.arg=c('vyborne','velmi dobre','prospel','neprospel'),
        main='Abs. cetnosti - Matematika',ylim=c(0,8),border='green4',
        xlab='znamka',ylab='absolutni cetnost',space=0)
#---------------------------------------------------------------------------------------
plot(1:4,nj,type='l',lwd=3,col='green4',xlab='znamka',ylab='abs. cetnost',
     main='Polygon abs.cetn. - Matematika',axes=F)
box()
points(1:4,nj,col='green4',pch=20)
axis(1, at=1:5, lab=c('vyborne', 'velmi dobre', 'dostatecne', 'nedostatecne',''))
axis(2, at=1:10)
#=====================================================================================
#Dvourozmetna bodova data
#Kontingencni tabulka simultannich absolutnich cetnosti
Ktab<-matrix(NA,4,4)
for(i in 1:4){
  for(j in 1:4){
    Ktab[i,j]<-sum(data$math==i & data$english==j)
  }
}
Ktab
K.tabulka<-data.frame(round(Ktab,digits=2),row.names=c('M_Vyborne','M_Velmi_dobre','M_Prospel',
                                                            'M_Neprospel'))
names(K.tabulka)<-c('E_Vyborne','E_Velmi_dobre','E_Prospel',
                    'E_Neprospel')
K.tabulka #tabulka typu data.frame
Ktab      #matice
#-------------------------------------------------------------------------------
#Radkove podminene relativni cetnosti (znamka z Aj podminena znamkou z matematiky)
KtabR<-prop.table(Ktab,1)         #do prikazu prop.table musime vlozit matici, tedy Ktab
KtabR

#Pojmenovani radku a sloupcu - nepovinne - nedelalo se na hodine
KtabR<-data.frame(round(KtabR,digits=2),row.names=c('M_Vyborne','M_Velmi_dobre','M_Prospel',
                                                    'M_Neprospel'))
names(KtabR)<-c('E_Vyborne','E_Velmi_dobre','E_Prospel',
                'E_Neprospel')
KtabR
#-------------------------------------------------------------------------------
#Sloupcove podminene relativni cetnosti (znamka z M podminena znamkou z Aj)
KtabS<-prop.table(Ktab,2)  #do prikazu prop.table musime vlozit matici, tedy Ktab
KtabS

#Pojmenovani radku a sloupcu - nepovinne - nedelalo se na hodine
KtabS<-data.frame(round(KtabS,digits=2),row.names=c('M_Vyborne','M_Velmi_dobre','M_Prospel',
                                                    'M_Neprospel'))
names(KtabS)<-c('E_Vyborne','E_Velmi_dobre','E_Prospel',
                'E_Neprospel')
KtabS