diff --git a/zajecia5/.Rhistory b/zajecia5/.Rhistory new file mode 100644 index 0000000..e69de29 diff --git a/zajecia5/Estimators-formulae.pdf b/zajecia5/Estimators-formulae.pdf new file mode 100644 index 0000000..8d058ba Binary files /dev/null and b/zajecia5/Estimators-formulae.pdf differ diff --git a/zajecia5/Konstrukcja przedziałów ufności.pdf b/zajecia5/Konstrukcja przedziałów ufności.pdf new file mode 100644 index 0000000..43d0dd6 Binary files /dev/null and b/zajecia5/Konstrukcja przedziałów ufności.pdf differ diff --git a/zajecia5/README.md b/zajecia5/README.md new file mode 100644 index 0000000..7201078 --- /dev/null +++ b/zajecia5/README.md @@ -0,0 +1,27 @@ +# Zajęcia 5 +Estymacja przedziałowa - np jakieś urządzenie moze działać na pewnym przedziale wartości. + + + +## Omówienie +Chcemy "złapać" jakąś wartość w przedział. Jest to lepsze niz próba oszacowania dokładnej wartości. +Jezeli konstruujemy jakiś przedział z poziomem ufności 0,95 to na 100 prób w 95 nasz parametr jest w przedziale. + +Na podstawie funkcji centralnej mozemy stworzyć przedziały ufności. Funkcję centralną bierzemy z tabelki. +Podstawiamy funkcję centralną do prawdopodobieństwa oraz 1-a, wsadzamy parametr pomiędzy funkcję i wyliczamy a i b. +A i b to przedział ufności. Jest przykład w [Konstrukcja przedziałów ufności](Konstrukcja-przedziałów-ufności.pdf) + +**Rozkład t-Studenta** po prostu typ rozkładu. Podobny to rozkładu normalnego. + + + +## 6.1 Przykład +Ten przedział ufności, który mamy podany, tworzy się analogicznie jak w punkcie omówienie. +Konstrukcja jest opisana w pdfie. + + +## Zadania +Zrobiłem z plikiem pomocniczym + + + diff --git a/zajecia5/Zajęcia5.pdf b/zajecia5/Zajęcia5.pdf new file mode 100644 index 0000000..b031be6 Binary files /dev/null and b/zajecia5/Zajęcia5.pdf differ diff --git a/zajecia5/zadania.R b/zajecia5/zadania.R new file mode 100644 index 0000000..9666e6a --- /dev/null +++ b/zajecia5/zadania.R @@ -0,0 +1,38 @@ +# ZAD1 + +load(url("http://ls.home.amu.edu.pl/data_sets/Centrala.RData")) +library(EnvStats) + +epois(Centrala$Liczba, + method = "mle/mme/mvue", + ci = TRUE, ci.type = "two-sided", conf.level = 0.95, + ci.method = "exact")$interval$limits + +epois(Centrala$Liczba, + method = "mle/mme/mvue", + ci = TRUE, ci.type = "two-sided", conf.level = 0.95, + ci.method = "pearson.hartley.approx")$interval$limits + +epois(Centrala$Liczba, + method = "mle/mme/mvue", + ci = TRUE, ci.type = "two-sided", conf.level = 0.95, + ci.method = "normal.approx")$interval$limits + +# ZAD 2 + +awarie <- read.table("http://ls.home.amu.edu.pl/data_sets/awarie.txt") + +lambda_limits <- eexp(awarie$V1, + method = "mle/mme", + ci = TRUE, ci.type = "two-sided", conf.level = 0.95, + ci.method = "exact")$interval$limits +rev(1 / lambda_limits) +rev(1 / lambda_limits^2) + + +# ZAD 3 + + + + + diff --git a/zajecia5/zajecia5.Rproj b/zajecia5/zajecia5.Rproj new file mode 100644 index 0000000..8e3c2eb --- /dev/null +++ b/zajecia5/zajecia5.Rproj @@ -0,0 +1,13 @@ +Version: 1.0 + +RestoreWorkspace: Default +SaveWorkspace: Default +AlwaysSaveHistory: Default + +EnableCodeIndexing: Yes +UseSpacesForTab: Yes +NumSpacesForTab: 2 +Encoding: UTF-8 + +RnwWeave: Sweave +LaTeX: pdfLaTeX