rm(list=ls())
graphics.off()
library(lattice)
library(moments)
setwd('I:\\mpis\\lab4') # prosz� odkomentowa� i wskaza� katalog, gdzie zapisano skrypt

tab <-read.table('ankieta.txt', header = TRUE)

dane <- tab$Wzrost[tab$P�e�=='M']
qqnorm(dane)
qqline(dane)
cat('Zmienne sa skupione w poblizu lini ')
cat ('\nH0: proba pochodzi z rozkladu normalnego\nH1 proba nie pochodzi z rozkladu normalnego')
shapiroTest=shapiro.test(dane)
print(shapiroTest)
if(shapiroTest$p.value<0.01)
{
cat ('Odrzucamy H0\n')
}else
{
cat ('Brak podstaw do odrzucenia H0\n')
}
wzrostTTest<-t.test(dane,mu=179,alternative='greater',conf.level=0.95)

cat('H0: m=179cm\nH1 m>179cm\n')

if(wzrostTTest$p.value<0.05)
{
cat ('Odrzucamy H0\n')
}else
{
cat ('Brak podstaw do odrzucenia H0\n')
}

ci_mean <- function(x,conf.level)
{
n <- length(x)
a=1-conf.level
ta <- qt(1-a/2,n-1)

c(mean(x)-ta*(sd(x)/sqrt(n-1)),mean(x)+ta*(sd(x)/sqrt(n-1)))
}
cat('Przedzial ufnosci: 0.95')
print(ci_mean(dane,0.95))
cat('Przedzial ufnosci: 0.98')
print(ci_mean(dane,0.98))

cat('Im wiekszy poziom ufnosci, tym wieksza szerokosc\n')
ci_sd<-function(x,conf.level){
n=length(x)
a=1-conf.level
s=sd(x)
c1=qchisq(a/2,n-1)
c2=qchisq(1-(a/2),n-1)
sqrt(c((n*(s^2))/c2,(n*(s^2))/c1))
}
cat('Poziom ufnosci: 0.97')
print(ci_sd(dane,0.97))