Instruktioner - Matematiska institutionen - Stockholms universitet

More documents

Recommendations

Info

när man startar det nästa gång.Uppgifterna 1–5 ska besvaras skriftligen. Bifoga den R-kod du skrivit, samt några grafer somvisar vad du har gjort och som illustrerar ditt resonemang.UPPGIFT 1.Plotta temperaturen mot tid (=observationsnummer här). Beskriv grafen i ord!UPPGIFT 2.Skapa en ny variabel innehållande årsvisa medelvärden (se R-tips nedan för hur man gördetta), och gör en graf över dem. Verkar temperaturen ha ökat eller minskat under perioden?Spara grafen!UPPGIFT 3.a) Ange ett viktsystem för centrerade glidande medelvärden. Viktsystemet ska ha egenskapenatt jämna ut den årsvisa säsongvariationen. Hur många värden innehåller viktsystemet?b) Hur många värden kommer serien av centrerade glidande medelvärden att ha?c) Beräkna serien av centrerade glidande medelvärden. Gör en graf över denna serie, ochjämför den med grafen över årsvisa medelvärden du gjorde i Uppgift 2. Vilka skillnader ochlikheter ser du?UPPGIFT 4.Beräkna de justerade säsongkomponenterna. Vilken månad är varmast? Vilken är kallast?UPPGIFT 5.Genomför en säsongrensning av data. Gör en graf över den rensade serien. Lägg in rensadserie, trend och ursprunglig serie i en och samma graf. (Se R-tips I nedan.)- 2 -
R-TIPS.I. Grafisk framställning av tidsserier.I vårt fall har vi ingen variabel i datasetet som specifikt anger löpande tid, så vad vi vill göra äratt plotta olika variabler mot radnummer. Om variabeln heter Y, så kan detta göras genomkommandot> plot(Y,type="l")som ritar serien som en sammanhängande kurva. Vill man ha färg så kan man skrivaexempelvis> plot(Y,type="l", col=”blue”)Bokstaven l (”ell”, inte ”ett”) står för ”line”.Om man vill lägga in två tidsserier Y och Z i samma graf, och färga grafen för Z röd, så kanman skriva> plot(Y,type="l")> lines(Z,col="red")Om de båda tidsserierna har olika längd och/eller ligger ur fas , så måste man komplettera denena serien eller båda med NA (som är R:s beteckning för saknade värden) i början och/ellerslutet. Med hjälp av rep (se nedan under IV) kan man åstadkomma så många NA som behövs,och sedan kombinera ihop dem med tidsserien med hjälp av c(…,…,…). (Se Lab 1, Statistiskanalys GN, ”Enkla numeriska operationer”.)II. Att spara figurer.Om man vill spara en figur — exempelvis den som skapats av det första plot-kommandotovan — måste man först ange kommandona> jpeg("figur.jpg")> plot(Y,type="l")> dev.off()som skapar filen figur.jpg i underbiblioteket tsa1. För att skriva ut filen, öppna först ettnytt terminalfönster, byt katalog till tsa1, och skriv$ lpr figur.jpgPå liknande vis görs en bild i ps-format.III. Att sammanfatta data.Kommandot tapply kan användas för att beräkna medelvärden (och även andra statistikor)över undergrupper av datasetets observationer. Om Y är den variabel vi vill beräknamedelvärden av (exempelvis temperatur), och X för varje observation anger vilkenundergrupp den tillhör (exempelvis månad eller år), så kan gruppvisa medelvärden beräknasgenom> MedelY
Page 1: MATEMATISKA INSTITUTIONENTillämpad
Page 5 and 6: På samma vis kan vi förvandla M t

Instruktioner - Matematiska institutionen - Stockholms universitet

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?