# 6. cvičení # spustit R v režimu SDI library(Rcmdr) # nechat zavést RcmdrPlugins.TeachingDemos # # modelujme si situaci, kdy testuje hypotétu o střední hodnotě # proti jednostranné alternativě # Demos | Power of the test # Sample size 1 # Alpha level 0.05 # Standard deviation 1 # to znamená, že rozhodujeme pomocí Y~N(0,1) # # při jednostranné alternativě s kritickou hodnotou 1,645 # True difference 0.5 # prohlédni obrázek, reprodukuj výpočet hodnoty1.644854 # qnorm(.95) (nesouvisí s alternativou!) # jakou má test sílu, je-li stř. hodnota o 0.5 vyšší? # jak se projeví, když změníme hladinu na alpha = 10 %? # jak se projeví větší počet pozorování? # jak se projeví menší počet pozorování? # # nastav sogma = 4 # nastav TRUE difference na 1 a znovu porovnej sílu pro n=48, 12 a 100 # kolik potřebuješ pozorování, aby nyní byla síla aspoň 60 %? (58) # ověření téhož vzorečkem z přednášky: ((qnorm(1-0.05)+qnorm(1-0.40))/1)^2*4^2 # # nastavit jako data Kojeni # bud poklepáním na tlačítko vedle Data set: # nebo pomocí Data | Load data set ... # nebo pomocí Data | Import ... # připomenout data např. pomocí View data set # uložit data Data | Active data set | Save active ... # # připomenout popisné statistiky pro delka podle pohlaví: # průměr, medián, směrodatná odhchylka # # dále budeme pracovat pouze s hochy: # Data | Active data set | Subset active data set # nech zaškrtnuto Include all .. (nebo si vyber) # Subset selection např. HochL nebo hoch==1 nebo Hoch=="hoch" # Name for ... např. KojeniN (nepřepsat původní soubor!) # # POZNÁMKA: # občas se mi stalo, že při přípravě takové podmnožiny dat nereagovalo # tlačítko OK; pak stačilo klepnout kamkoliv mimo toto tlačítko # a jeho funkce se obnovila (patří to ke známým problémům knihovny # tcltk, kterou Commander používá) # # testujme hypotézu, že populační průměr proměnné delka je roven 68 # (známá konstanta z minulých let) # proti alternativě, že je větší # Statistics | Means | Single sample t-test # nastav delka, zvol alternativu # Null hypothesis: mu= 68, nech Confidence level .95 # zamítli jsme na 5% hladině? NEZAMÍTLI # na 10% hladině bychom již zamítli # # # ale předpokládali jsme normální rozdělení # ověření normality: # Graphs | Quantile-comparison plot # zvolíme proměnnou (delka), necháme nastavení Distribution # Statistics | Summaries | Shapiro-Wilk test of normality # # rozhodněte o tvrzení, že hmotnost hochů v 24. týdnu života je # rovna 7650 g (rozumí se její populační průměr) # řešením je t-test s oboustrannou alternativou! # # samostatně (možnost od jisté chvíle ukládat do souboru, # který si odnesu a zhodnotím) # popisné statistiky a interval spolehlivosti proměnných: # por.hmotnost, por.delka, hmotnost, vyska.m, vyska.o