Statystyka, prognozowanie, ekonometria, data mining

haver · Plutonowy Posty: 64 Skąd: Lubaczów

Witam,
mam problem z funkcją ts(). Dokładniej chodzi mi o dopasowanie dat do moich danych podczas zmieniania danych z wektora na klasę ts

ts(data = NA, start = 1, end = numeric(0), frequency = 1,
deltat = 1, ts.eps = getOption("ts.eps"), class = , names = )

Mam problem z danymi giełdowymi (wartosci akcj itd.), ponieważ notowania dotyczą jedynie 5 dni w tygodniu....Co mam wpisac we frequency?

Dzieki z góry za pomoc.
B.

Google

piotrek · Plutonowy **Pomógł:** 1 raz Posty: 60

Sugeruję wykorzystanie bibliotek: timeDate oraz timeSeries (dawniej fSeries) oraz innych z fmetrics, są tam m.in. specjalne klasy do obsługi (nie)regularnych serii czasowych (np. giełdowych).

haver · Plutonowy Posty: 64 Skąd: Lubaczów

Szukalem i nie moge znaleźć...tam jest tego za dużo :(. Wczytuje dane giełdowe i daty. Dane są numeric a daty szybko i sprawnie można zamienić na format Date. Tylko jak teraz narysować wykres np. wartości akcji, żeby na osi x był numer roku (np co kilka lat), ale tak zeby wszystko bylo zsynchronizowane...

Pomóżcie proszę!

Pozdrawiam,
B.

piotrek · Plutonowy **Pomógł:** 1 raz Posty: 60

zobacz:

Kod:

library(timeSeries)
?timeSeries::plot
?timeSeries::as
data(MSFT)
head(MSFT)

haver · Plutonowy Posty: 64 Skąd: Lubaczów

Znowu nie wychodzi....mam szesc kolumn...wybieram kolumny z data i kursem zamknięcia. Kurs zamknięcia zamieniam na numeric...date na date i potem robie z tego timeSeries za pomocą funkcji as.timeSeries. Wszytko sie chyba udaje bo jak dam dane[,1] to wyswietla mi kurs zamkniecia (jakby nie bylo timeSeriec to wyswietliłaby sie kolumna z datami)...daje plot(dane) i na osi x mam numery obserwacji a nie daty...yhhh. R potrafi zdenerwowac jak sie nie umie czegos zrobic :) w gretlu wczytuje dane ustawiam tydzien 5 dniowy i zaznaczam punkt startowy i wszystko smiga...

Crunchy · Wysłany: 2009-10-18, 19:57

piotrek · Plutonowy **Pomógł:** 1 raz Posty: 60

haver · Plutonowy Posty: 64 Skąd: Lubaczów

http://www.parkiet.com/dane/danesesji/bazytxt/akcje/vistula.zip

Kod:

dane<-read.table("........",sep=";",h=T)
dane<-dane[,c[2,6])
dane[,1]<-as.Date(as.character(dane[,1]),"%Y%m%d")
dane[,2]<-as.numeric(dane[,2])

Dzieki za pomoc i pozdrawiam : )

Crunchy · Wysłany: 2009-10-20, 00:28

Kod:

library(zoo)
dane<-read.table("vistula.txt",sep=";",dec=",",h=T)
dane.z.datami<- zoo(as.matrix(dane[3:7]), as.Date((as.character(dane[[2]])),format="%Y%m%d"))
plot(dane.z.datami)

haver · Plutonowy Posty: 64 Skąd: Lubaczów

Dziala...nie wiedziałem że istnieje zoo : ) Dzieki wielkie!

piotrek · Plutonowy **Pomógł:** 1 raz Posty: 60

Trochę inaczej:

Kod:

temp <- read.table(file="vistula.txt", header = T, sep = ";", dec = ",")
colnames(temp)<-c("TICKER", "DATA", "OPEN", "HIGH", "LOW", "CLOSE", "VOL")
library(xts)
a <- xts(zoo(temp[, 3:7], as.Date(as.character(temp[, 2]), format = "%Y %m %d")))
plot(a) # na osi X daty
plot(a$HIGH) # na osi X daty

haver · Plutonowy Posty: 64 Skąd: Lubaczów

Hmmm coś nie działa. Gdy korzystam z funkcji plot wyskakuje komunikat:

Warning messages:
1: In plot.xts(a) : only the univariate series will be plotted
2: In as.double.xts(y) : NAs introduced by coercion

Co mam zrobic?

Dzieki za pomoc.
Pozdrawiam,
Bartek

bstq · Chorąży **Pomógł:** 9 razy Posty: 108 Skąd: Warszawa

univariate adj. having one variable (in Statistics - about a distribution) - szereg musi byc jedowymiarowy zeby go narysowal

a w tym drugim:
znalazlem w: http://www.mail-archive.c...h/msg15767.html

"My guess would be you're trying to convert from character to numeric and are unable to do so. As in,

> as.numeric("A")
[1] NA
Warning message:
NAs introduced by coercion
> as.numeric("1")
[1] 1
>

But without more information from you it's impossible to tell.

See the posting guide at

http://www.R-project.org/posting-guide.html

Regards,
Sundar"

google się kłania :-D

pozdrawiam