[OPGELOST] Hypothesetoetsen:oef

2de semestervak op 3 studiepunten
Gebruikersavatar
Ruben
Doctor in de forumwetenschappen
Doctor in de forumwetenschappen
Berichten: 4848
Lid geworden op: 20 dec 2007, 21:15
Locatie: Steenhuffel

Re: [back from opgelost]Hypothesetoetsen:oef

Berichtdoor Ruben » 21 jun 2009, 14:32

Bedankt Pieter!
Ruben - Delivering awesomeness since 1989
Gebruikersavatar
hanssimi
Heeft dit forum graag
Heeft dit forum graag
Berichten: 188
Lid geworden op: 15 dec 2007, 16:59
Contacteer:

Re: [back from opgelost]Hypothesetoetsen:oef

Berichtdoor hanssimi » 21 jun 2009, 15:05

Pieter schreef:

Code: Selecteer alles

H0=false;
periode=1;

while H0==false      periode=periode+1;
   
    try z=(reshape(u,periode,length(u)/periode));
    catch continue
    end
   
Z=mean(z,2)
T=sqrt(periode)*mean(Z)/std(Z);
   
    if abs(T)<5 % kies een waarde uit de tabel
        H0=true;
    end
end
disp(['de periode bedraagt: ',num2str(periode)]);


Dit geeft je de gezochte periode.




ik heb enkele vraagjes



waarom neem je sqrt(periode)?
in de T - test
en wat is die 5
en heb je dat zelf geschreven:p?
Gebruikersavatar
yoachim
Prof in de forumwetenschappen
Prof in de forumwetenschappen
Berichten: 7725
Lid geworden op: 29 nov 2007, 20:44
Locatie: Halle

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor yoachim » 21 jun 2009, 16:19

ook een vraagje:
ge gebruikt ne try-catch, das wrs om te voorkomen dat ge periode neemt dat geen deler is van 25000 (lengte signaal), ma stel dat de periode toch geen deler is (bvb 21)
Afbeelding
like a bouws
Afbeelding
Ur mom is so fat she doesn't need the internet, she's already world wide
Pieter
Regelmatig forumgebruiker
Regelmatig forumgebruiker
Berichten: 88
Lid geworden op: 17 dec 2007, 22:20
Contacteer:

Re: [back from opgelost]Hypothesetoetsen:oef

Berichtdoor Pieter » 21 jun 2009, 19:22



ik heb enkele vraagjes



waarom neem je sqrt(periode)?
in de T - test
en wat is die 5
en heb je dat zelf geschreven:p?


1. omdat de periode gelijk is aan n^1/2
2. die vijf is een getalleke uit de tabel. Eigenlijk moet je elke keer ifv de vrijheidsgraden en de nauwkeurigheid een nieuw getalleke nemen, maar daar had ik geen zin in, ik heb dus maar iet ingevuld
3. ik heb dat helemaal zelf geschreven!
Pieter
Regelmatig forumgebruiker
Regelmatig forumgebruiker
Berichten: 88
Lid geworden op: 17 dec 2007, 22:20
Contacteer:

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor Pieter » 21 jun 2009, 19:25

yoachim schreef:ook een vraagje:
ge gebruikt ne try-catch, das wrs om te voorkomen dat ge periode neemt dat geen deler is van 25000 (lengte signaal), ma stel dat de periode toch geen deler is (bvb 21)


Dan hebt ge inderdaad een probleem.
Je kan altijd een aantal punten schrappen totdat de data wel deelbaar is door je periode.
Maw: je verwijderd mod(lengte,periode) aantal punten.
Dit kan je dan doen voor elke mogelijke periode.
Gebruikersavatar
ideglier
Prof in de forumwetenschappen
Prof in de forumwetenschappen
Berichten: 5400
Lid geworden op: 11 okt 2008, 18:27
Locatie: Asse

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor ideglier » 24 jun 2009, 11:29

Voor de mensen die er niet naar willen zoeken:

A = u (met u matrix 25000x10)
A((length(u(:,1))+1 - mod(length(u(:,1)),periode)):length(u(:,1)),:) = []; (de tweede index in u nl 1 mag je uiteraard vrij kiezen zolang die kleiner is dan 11)

Zo kap je een stuk van de matrix af
I love the smell of petrol in the morning. Smells like ... SPEED!
Afbeelding
Gebruikersavatar
yoachim
Prof in de forumwetenschappen
Prof in de forumwetenschappen
Berichten: 7725
Lid geworden op: 29 nov 2007, 20:44
Locatie: Halle

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor yoachim » 24 jun 2009, 12:35

is dit niet beetje makkelijker?

Xorigineel: 1=> 25000
Xnieuw: 1 => 24653 bvb

Xnieuw(1:24653) = Xorigineel(1:24653)
Afbeelding
like a bouws
Afbeelding
Ur mom is so fat she doesn't need the internet, she's already world wide
Gebruikersavatar
ideglier
Prof in de forumwetenschappen
Prof in de forumwetenschappen
Berichten: 5400
Lid geworden op: 11 okt 2008, 18:27
Locatie: Asse

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor ideglier » 24 jun 2009, 13:18

zo kan het idd ook...
I love the smell of petrol in the morning. Smells like ... SPEED!
Afbeelding
Gebruikersavatar
Tom V
Master in de forumwetenschappen
Master in de forumwetenschappen
Berichten: 2996
Lid geworden op: 28 nov 2007, 20:09
Contacteer:

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor Tom V » 24 jun 2009, 19:02

Een algemenere mogelijkheid is gebruik maken van floor:

Code: Selecteer alles

function bool=periodejuist(u, kolom, periode)   
Y=reshape(u(1:floor(length(u(:,kolom))/periode)*periode,kolom),periode,floor(length(u(:,kolom))/periode));
sinus=mean(Y,2);
plot(sinus)
T=sqrt(periode)*mean(sinus)/std(sinus);
if T<5
    bool=true;
else
    bool=false;
end


Deze kan alle periodes aan. All right!

EDIT: En wie echt even cool wil zijn als ik, kan daarnaast gebruik maken van deze functie die de periode bepaalt:

Code: Selecteer alles

function ZoekPeriode(u,kolom)
periode=25;
bool=false;
while bool==false
    periode=periode+1;
    bool=periodejuist(u, kolom, periode);
end
disp(['De periode is ',num2str(periode)]);
Dit bericht kreeg een Chuck Norris quality label:

Afbeelding
Jorne
Beginnend forumgebruiker
Beginnend forumgebruiker
Berichten: 25
Lid geworden op: 02 dec 2008, 21:13
Locatie: Brugge
Contacteer:

Re: [opgelost]Hypothesetoetsen:oef

Berichtdoor Jorne » 24 jun 2009, 20:00

Goed gezien. Dat komt eigenlijk gewoon neer op afronden op een periode nauwkeurig.

Kleine opmerking: als je de goeie periode (25) wilt vinden, raad ik je aan de beginwaarde kleiner dan 25 te kiezen in jouw programma. Hij zal nu namelijk beginnen bij 25+1, wat niet juist is en stoppen bij 50, de eerstvolgende juiste periode.
What is mind? No matter.
What is matter? Never mind.

Terug naar “Statistiek voor ingenieurs”

Wie is er online

Gebruikers op dit forum: Geen geregistreerde gebruikers en 1 gast

cron