Statystyka, prognozowanie, ekonometria, data mining Strona Główna
Reklama pqstat.pl
Statystyka, prognozowanie, ekonometria, data mining
Forum miłośników statystyki - Portal Statystyczny

FAQFAQ  SzukajSzukaj  UżytkownicyUżytkownicy  GrupyGrupy  StatystykiStatystyki
RejestracjaRejestracja  ZalogujZaloguj  Chat   Regulamin  Kadra forum
PORTAL STATYSTYCZNY
 Ogłoszenie 
FORUM STATYSTYCZNE MA JUŻ 10 LAT

Znasz statystykę lub ekonometrię, metody prognozowania, data mining i chcesz pomóc w rozwoju forum statystycznego ?
Pisz na: administrator(małpa)statystycy.pl

Rozpoczął swoją działalność portal statystyczny - masz pomysł na jego rozwój ?

Drogi forumowiczu! Zanim napiszesz posta zapoznaj się z regulaminem forum i przedstaw się
The International Year of Statistics (Statistics2013) Free statistics help forum. Discuss statistical research, statistical consulting Smarter Poland Portal statystyczny

Poprzedni temat «» Następny temat

Tagi tematu: Brak tagów.

Przesunięty przez: mathkit
2015-03-23, 16:03
[R] t.test()
Autor Wiadomość
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-01, 19:12   [R] t.test()

Witam.

Czy ja mógłbym prosić o pomoc w rozwiązaniu pewnego zadania? Oto jego treść:


Zadanie mam do rozwiązania z programie R. Tzn to nie jest zadanie które muszę zrobić ale przygotowuje się do egzaminu poprawkowego i rozwiązuje sobie zadania i właśnie z tym zadaniem mam problem a nie robiliśmy nic podobnego na zajęciach. Nie wiem nawet jak mam wprowadzić tą tabelkę do programu R i to uporządkować aby to było sensownie ułożone czyli tak jak jest w zadaniu. Z góry dziękuje za pomoc.

chi-square.jpg
chi-square
Plik ściągnięto 242 raz(y) 69,44 KB

Ostatnio zmieniony przez mathkit 2011-11-19, 00:15, w całości zmieniany 2 razy  
 
 
     
Google

Wysłany:    Reklama google.

 
 
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-01, 19:54   

> m=cbind(c(86,17,54),c(31,64,39),c(132,43,132))
> m
[,1] [,2] [,3]
[1,] 86 31 132
[2,] 17 64 43
[3,] 54 39 132

> chisq.test(m)

Pearson's Chi-squared test

data: m
X-squared = 84.3745, df = 4, p-value < 2.2e-16
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-01, 20:15   

I to jest naprawdę tylko tyle? Kurcze, aż nie wiedziałem że to takie proste. A jak jest z hipotezami? Tzn ja hipotezy postawię a Ty ja możesz to powiedz czy dobre :).

H -> nie ma zależności poglądów
K -> istnieje zależność poglądów.

pv<alfy zatem odrzucamy hipotezę H.

Dobrze napisałem czy popełniłem błąd?
 
 
     
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-01, 20:21   

OK

zmienne są zależne
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-01, 21:16   

Czyli hipotezy postawiłem poprawne, tak?

Jeszcze mam pytanie, może mi pomożesz, nie chce już zakładac nowego tematu żeby nie zaśmiecać forum. Chodzi o funkcje "R". Bo np żeby policzyć np 95% przedział ufności dla średniej to używam funkcji "t.test(dane, conf.level=0.95)". A czy sa jakieś gotowe funkcję gdy chcę policzyć np 95% przedział dla wariancji i 95% przedział dla odchylenia standardowego?
 
 
     
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-02, 07:43   

Hipotezy postawiłeś poprawne
--------------------
library("TeachingDemos")

sigma.test(x,sigma=2) # dla odchylenia standardowego

sigma.test(x,sigmasq=2) # dla wariancji
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-02, 20:32   

Jak mam użyć tych testów dla odchylenia i wariancji? U mnie wywala komunikat "Error: could not find function "sigma.test" ". Podejrzewam że trzeba coś doinstalować ale czy możesz powiedzieć jak to zrobić?
 
 
     
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-02, 20:48   

> install.packages("TeachingDemos")
> library("TeachingDemos")
> sigma.test(......)
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-02, 20:50   

A możesz jeszcze powiedzieć jak mam w tym teście uwzględnić ten 95% przedział ufności? Bo przy średniej miałem parametr conf.level a w tym teście nic takiego nie mam.

Ok, on mi standardowo testuje dla poziomu ufności 95%. A jeśli mam poziom ufności 98%?

Ok, chyba wiem. Podałem parametr conf.level=0.98 ale czy to się tak właśnie robi?
 
 
     
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-02, 20:57   

Jak chcesz się dowiedzieć o danym np. teście skorzystaj z dokumentacji np:
> ?sigma.test
------
> sigma.test(...... , conf.level=0.98)
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-02, 21:08   

Ooo, nawet nie wiedziałem że taka dokumentacja istnieje :). Dzięki wielkie. A czy ten parametr: sigma=2 to jest jakas stała czy ona wynika z zadania?

A jest jakiś test dzięki któremu mogę na poziomie istotności 0.05 sprawdzić czy średnia waga przekracza 110? Dane z wagami oczywiście mam.
 
 
     
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-03, 07:51   

sigma = 2: test dla odchylenia standardowego = 2.
sigmasq=2: test dla wariancji = 2.
Wartości te wynikają z treści zadania.
-------------
Test dla jednej średniej możesz przeprowadzić za pomocą testu t-Studenta. Przeczytaj dokumentację i wszystko będzie jasne:
?t.test
------------
t.test(waga, alternative="less", mu=110)

[ Dodano: 2010-09-03, 07:53 ]
t.test(waga, alternative="greater", mu=110) # przekracza 110.
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-03, 14:45   

Czyli gdyby było bez opcji "greater" wówczas było by że nie przekracza, tak? A jeśli wstawię ta opcję to wtedy liczę ile przekracza?
 
 
     
mikos 
Podporucznik


Pomógł: 45 razy
Posty: 161
Skąd: Kutno
Wysłany: 2010-09-03, 14:59   

Powtarzam. W dokumentacji jest wszystko o opcjach testów napisane. ;-)
----------
mu=110 , alternative = "less" (pol. "mniej") - mniej niż 110
mu=110, alternative = "greater" (pol. "więcej") - więcej niż 110
mu=110, alternative = "two.sided" (różny) - różny od 110
 
     
lukpio3 
Starszy Szeregowy


Posty: 20
Skąd: Warszawa
Wysłany: 2010-09-03, 15:10   

Ok, teraz właśnie zauważyłem Tylko wiesz co, coś mi chyba nie wychodzi. Bo mam zadanie w którym pytają czy na poziomie istotności można przyjąć że średnia waga kurczaków przekracza 110. Ty też dane na pewno też masz. Otwierasz je komendą ChickWeight i tam w kolumnie weight jest własnie ta średnia waga kurczaków. Ja to robię tak:

> w=ChickWeight$weight
> t.test(w, alternative="greater", mu=110)

One Sample t-test

data: w
t = 3.9978, df = 577, p-value = 3.612e-05
alternative hypothesis: true mean is greater than 110
95 percent confidence interval:
116.948 Inf
sample estimates:
mean of x
121.8183

H0 -> średnia waga wynosi 110
K0 -> średnia waga przekracza 110

No i wyszło pv<alfy bo wyszło jakieś strasznie małe zatem odrzucamy H0 i przyjmujemy K ale to jest chyba raczej źle zrobione. Tak mi się przynajmniej wydaje i nie wiem gdzie popełniłem błąd.
 
 
     
Wyświetl posty z ostatnich:   
Odpowiedz do tematu
Nie możesz pisać nowych tematów
Możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Możesz ściągać załączniki na tym forum
Dodaj temat do Ulubionych zakładek(IE)
Wersja do druku

Skocz do:  

Podobne Tematy
Temat Autor Forum Odpowiedzi Ostatni post
Brak nowych postów Przyklejony: Chi2 test (vs test dla dwóch wskaźników struktury)
Ania Testowanie hipotez statystycznych 103 2013-05-06, 17:21
prosie23
Brak nowych postów Przyklejony: Dokładny test Fishera (Fisher's Exact Test)
Badanie zmiennej binarnej w przekroju klas
wariancja Testowanie hipotez statystycznych 59 2015-03-16, 22:08
mathkit
Brak nowych postów Przyklejony: Test istotności wahań sezonowych-test Kendalla
tomsew Testowanie hipotez statystycznych 3 2009-10-21, 18:08
Crunchy
Brak nowych postów Przyklejony: Test ADF
[Gretl]
barrel Modelowanie ekonometryczne 16 2011-05-24, 16:57
Madrich
Brak nowych postów Przyklejony: Test F
m|chał Testowanie hipotez statystycznych 32 2010-05-03, 11:00
mareczekfm

Ideą przyświecającą istnieniu forum statystycznego jest stworzenie możliwości wymiany informacji, poglądów i doświadczeń osób związanych ze statystyką, mierzenie się z różnego rodzaju problemami statystycznymi i aktuarialnymi. Poruszane problemy: Statystyka w badaniach sondażowych rynku, metody reprezentacyjne, Teoria i rachunek prawdopodobieństwa, statystyka opisowa, teoria estymacji, testowanie hipotez statystycznych, ekonometria, prognozowanie, metody data mining.
Copyright (C) 2006-2015 Statystycy.pl
Powered by phpBB modified by Przemo © 2003 phpBB Group
Strona wygenerowana w 0,06 sekundy. Zapytań do SQL: 30