Instruktioner - Matematiska institutionen - Stockholms universitet
Instruktioner - Matematiska institutionen - Stockholms universitet
Instruktioner - Matematiska institutionen - Stockholms universitet
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
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