setwd(' ')

# 1)
data <- read.table('newborns.txt', header=TRUE)

x<-na.omit(data$weight.C)
x<- as.numeric(x) # technicka uprava dat pro odstraneni NA hodnot

# zopakujme si z minula, kdy jsme porodni delku modelovali pomoci normalniho rozdeleni

mu<-mean(x) # stredni hodnotu odhadneme pomoci prumeru
sigma_squared<-var(x) # rozptyl odhadneme pomoci vyberoveho rozptylu
sigma<-sd(x) # a tedy smerodatnou odchylku  pomoci vyberove smerodatne odchylky

hist(x, prob = T, density = 30,
col = 'blue', xlab = 'porodni hmotnost',ylab = 'hustota', main = '')
curve(dnorm(x,mu,sigma),add=T,col='red',lwd=2)

# 1a)


# 1cde)

?t.test


# 2)
IQ<-read.table("IQ.txt")
x<-IQ$V1 #opet technikalie pro odstaneni data.framu

# priste si ukazeme poradne, ale vetsina testu predpoklada, ze data jsou normalni
# tak to aspon lehce overime

mu<-mean(x) # stredni hodnotu odhadneme pomoci prumeru
sigma<-sd(x) # a smerodatnou odchylku  pomoci vyberove smerodatne odchylky

hist(x, prob = T, density = 30,
col = 'blue', xlab = 'IQ',ylab = 'hustota', main = '')
curve(dnorm(x,mu,sigma),add=T,col='red',lwd=2)

# 2a)
library(EnvStats)
?varTest


# 3) meli bychom overit, ze data pochazi z normalniho rozdeleneni se sm.odchylkou 15
# opet zacneme s vizualizaci dat

mu<-mean(x) # stredni hodnotu odhadneme pomoci prumeru
sigma<-15 # a smerodatna odchylka je presne 15

hist(x, prob = T, density = 30,
col = 'blue', xlab = 'IQ',ylab = 'hustota', main = '')
curve(dnorm(x,mu,sigma),add=T,col='red',lwd=2)

# 3a)
library(BSDA)
?z.test


# 4) nacteme data a vybereme jen muze
data <- read.table('klicnikosti.txt', header=TRUE)

x<-data$length.L[data$sex=='m'] # delka leve klicni kosti u muzu
y<-data$length.R[data$sex=='m'] # delka prave klicni kosti u muzu

plot(x,y,asp=1,xlab='delka leve klicni kosti',ylab='delka prave klicni kosti')
# provedeme parovy t-test
# ale nejprve overime predpoklad normality rozdilu

z<-x-y
mu<-mean(z) # stredni hodnotu odhadneme pomoci prumeru
sigma<-sd(z) # a tedy smerodatnou odchylku  pomoci vyberove smerodatne odchylky

hist(z, prob = T, density = 30,
col = 'blue', xlab = 'rozdil delek',ylab = 'hustota', main = '')
curve(dnorm(x,mu,sigma),add=T,col='red',lwd=2)

?t.test