Other
Oefensessie 6
- Course
- Institution
Oefeningen uit oefensessie 6 uitgewerkt
[Show more]
Seller
Followcathoschelkens
Reviews received
Content preview
Oefeningen TijdreeksanalyseP0S37A Thema’s uit de kwantitatieve methodologie
Academiejaar 2022-2023
0. Start R en kopieer het R script en de gegevensbestanden voor de oefeningen vanuit Toledo naar de desktop (of
een andere plaats op de laptop waar deze bestanden gemakkelijk bereikbaar zijn). De gegevens en het R script zijn
terug te vinden op Toledo
Cursusdocumenten / Gegevensbestanden / Oefeningen tijdreeksanalyse
Het R script heeft als naam ‘Joos_script.txt’ en kan in de R Graphic User Interface (RGui) worden geopend via
File / Open script…
Standaard geeft R de extensie .R aan scripts, maar in dit geval is de extensie .txt, dus om de volledige lijst van
beschikbare bestanden op de desktop te zien moet er nog worden aangegeven:
Files of type: All files (*.*)
Ik geef er de voorkeur aan om met de extensie .txt te werken om duidelijk te maken dat het over een ongeformatteerd
tekstbestand gaat dat gemakkelijk met Notepad kan worden bewerkt.
Om het R script stap-voor-stap te doorlopen, selecteert u de lijn (of de lijnen) die u wilt uitvoeren, klik rechts op de
muis en selecteer
Run line or selection
De gegevens voor deze oefening bevinden zich in het csv-bestand ‘Joosdata123.csv’. Vergelijk het resultaat van het
vijfde model vanuit het script (“model5”) met het resultaat op bladzijde 104 van Joos, Onghena, en Pieters (1996). De
resultaten moeten gelijkaardig zijn. Kleine verschillen in de parameterschattingen kunnen voorkomen omwille van
nieuwe software en verbeterde schattingsalgoritmen.
Werken met R
Iedere lijn in R Console (script) een voor een aanduiden en drukken op ‘Run line’
, 2
1. Schat de lag 1 autocorrelatiecoëfficiënt voor de voorbeeldgegevens op dia 15 van het hoorcollege over
tijdreeksanalyse. De gegevens kunnen worden gevonden in het bestand ‘Example slide 15.csv’. Welke waarde komt u
uit? Waarom wijkt deze waarde af van de correlatie die op dia 18 werd berekend?
Tip: Schrijf een nieuw script op basis van de eerste lijnen uit Joos_script.txt. Geef
dat nieuwe script een andere naam, bijvoorbeeld Oef1script.txt, om het later
gemakkelijk terug te vinden. Gebruik een gelijkaardige eerste lijn dan de eerste lijn
van Joos_script.txt om de gegevens in te lezen maar geef een andere willekeurige naam
aan het bestand. Gebruik de acf functie om de lag 1 autocorrelatiecoëfficiënt te
berekenen.
#Stap 1: Script Schrijven in notepad, opslaan als [Oef1script]
#Stap 2: Script openen in R console
> Oef1 <- read.csv(file.choose(new=FALSE),sep=";")
> acf(Oef1$Y,plot=FALSE)
(Gegevensbestand: Example slide 15.cvs gebruiken)
De lag 1 autocorrelatiecoëfficiënt is dus gelijk aan 0.749. Deze waarde wijkt (licht) af van de waarde .7632 die op dia
18 werd berekend omdat op dia 18 een gewone Pearson product-moment correlatiecoëfficiënt werd berekend voor
de verschoven tijdreeks. Deze gewone correlatiecoëfficiënt geeft slechts een ruwe benadering van de door R
gebruikte gecorrigeerde schatting van de populatie-autocorrelatie die, onder andere, rekening houdt met het feit dat
er minder gegevensparen kunnen worden gevormd (slechts 19 gegevensparen voor 20 observaties bij stapwaarde
1). Voor meer uitleg over de problematiek van de zuivere schatting van de populatie-autocorrelatie bij korte
tijdreeksen, zie Huitema, B. E., & McKean, J. W. (1991). Autocorrelation estimation and inference with small samples.
Psychological Bulletin, 110(2), 291–304. https://doi.org/10.1037/0033-2909.110.2.291
, 3
2. Schat de lag 2 partiële autocorrelatiecoëfficiënt voor de voorbeeldgegevens op dia 15 van het hoorcol lege over
tijdreeksanalyse. Welke waarde komt u uit? Waarom wijkt deze waarde af van de partiële correlatie die op dia 18 werd
berekend?
Tip: Gebruik de pacf functie.
Voortgaan op script uit oefening 1
# Stap 3: Schat de lag 2 partiële autocorrelatiecoëfficiënt
> pacf(Oef1$Y,plot=FALSE)
De lag 2 partiële autocorrelatiecoëfficiënt is dus gelijk aan −0.405. Deze waarde wijkt (licht) af van de waarde −.3593
die op dia 18 werd berekend omdat op dia 18 een gewone partiële correlatiecoëfficiënt werd berekend voor de met
twee stappen verschoven tijdreeks, gecorrigeerd voor de met één stap verschoven tijdreeks. Deze gewone partiële
correlatiecoëfficiënt geeft slechts een ruwe benadering van de correcte schatting van de partiële
autocorrelatiecoëfficiënt die rekening houdt met het feit dat er minder gegevensparen kunnen worden gevormd
(slechts 18 gegevensparen voor 20 observaties bij stapwaarde 2).
The benefits of buying summaries with Stuvia:
Guaranteed quality through customer reviews
Stuvia customers have reviewed more than 700,000 summaries. This how you know that you are buying the best documents.
Quick and easy check-out
You can quickly pay through credit card or Stuvia-credit for the summaries. There is no membership needed.
Focus on what matters
Your fellow students write the study notes themselves, which is why the documents are always reliable and up-to-date. This ensures you quickly get to the core!
Frequently asked questions
What do I get when I buy this document?
You get a PDF, available immediately after your purchase. The purchased document is accessible anytime, anywhere and indefinitely through your profile.
Satisfaction guarantee: how does it work?
Our satisfaction guarantee ensures that you always find a study document that suits you well. You fill out a form, and our customer service team takes care of the rest.
Who am I buying these notes from?
Stuvia is a marketplace, so you are not buying this document from us, but from seller cathoschelkens. Stuvia facilitates payment to the seller.
Will I be stuck with a subscription?
No, you only buy these notes for $6.30. You're not tied to anything after your purchase.
Can Stuvia be trusted?
4.6 stars on Google & Trustpilot (+1000 reviews)
77973 documents were sold in the last 30 days
Founded in 2010, the go-to place to buy study notes for 14 years now