### ZAJĘCIA 1.12.2010 ###

zakupy=read.csv2("zakupy.csv")
attach(zakupy)
zakupy
names(zakupy)
zakupy

p1=subset(Wydatek, PLEC=="K");
mean(p1)
which.max(WIEK)
min(WIEK)
max(WIEK)
quantile(WIEK)
summary(WIEK)

table(WYKSZTAŁCENIE)
table(WIEK)

table(cut(WIEK,b=6))


####szereg rozdzielczy na wykresie słupkowym
k=20
dl=(max(WYDATEK)-min(WYDATEK))/(k-1)
podzial=min(WYDATEK)+dl*0:k-dl/2;
podzial
podzial=seq(min(WYDATEK)-dl/2,max(WYDATEK)+dl/2,,k+1);
podzial=seq(min(WYDATEK)-dl/2,max(WYDATEK)+dl/2,,k+1);
podzial
(szereg1=table(cut(WYDATEK,podzial)))
barplot(szereg1) #wykres slupkowy

hist(WYDATEK,40) #histogram
hist(WYDATEK,breaks=c(0,10,20,55,100,200,300))

plot(density(WYDATEK))

### TEST SHAPIRO-WILKA

shapiro.test(WYDATEK)
mean(WYDATEK)

### TEST T-STUDENTA ###
t.test(WYDATEK,mu=70)

##############################################################################
#ZADANIE 5(pkt1)
Wzrost1=read.csv2("Wzrost1.csv")
attach(Wzrost1)
Wzrost1
names(Wzrost1)
Wzrost1

mi=mean(Wzrost1)
sigma=sd(Wzrost1)
ile=length(Wzrost1)
k=20
w=0:k
podzial=qnorm(w/k,mi,sigma)
podzial
podzial[1]=min(Wzrost1)-3
podzial[length(Wzrost1)]=max(Wzrost1)+3
szereg=table(cut(X212,podzial))
szereg
(t=chisq.test(szereg))
names(t)