dúvida function fts to create functional time series objects
Caros colegas, Estou a tentar fazer uns graficos de series temporais com a libraria rainbow, a qual precisa de ter os dados arranjados numa functional time series (FTS). Eu ja sei transformar os dados numa ts normal, mas agora gostaria de mostrar todo num só gráfico com esta libraria rainbow: 1.- o comportamento da serie temporal por meses (no eixo x) dos diferentes barcos a longo dos anos, 2.- o comportamento dos desembarques por portos de desembarque, 3.- simplesmente os desembarques mensais por ano para todos os barcos e portos juntos (diferentes cores para cada ano). Os dados que tenho, ficheiro OCTOPUS.csv, são desembarques em kilogramos (QVENDA) de polvo, em preço (QVALOR, €), porto de desembarque (IPORTO), codigo de cada embarcação (CodEmbarc), DAY, MONTH, YEAR. Realmente não consigo perceber a ficheiro de ajuda quando escrevo ?fts, alguma ajuda disponivel? data= http://www.datafilehost.com/download-5bbb581a.html obrigado pelo tempo disponivel, Carlos -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605
Carlos, Eu nao verifiquei detalhes do teu conjunto de dados, entretanto, imagino que necessites formata-lo apropriadamente (uma formatacao distinta para cada questao). Para a primeira, assumindo que teu interesse seja 'QVENDA', eu comecaria com algo como: dadosOriginais = read.csv('OCTOPUS.csv') library(reshape2) porBarco = dcast(dadosOriginais, DAY + MONTH + YEAR ~ CodEmbarc, value.var='QVENDA', fun.aggregate=mean) x0 = 1:nrow(porBarco) y0 = as.matrix(porBarco[, -(1:3)]) library(rainbow) obj0 = fts(x0, y0) plot(obj0) E usaria a mesma estrategia para os demais. benilton
ok benilton, obrigado, mas a intenção de usar a libraria rainbow e para fazer um grafico como esse daqui que conseguí num paper do R Journal: #os dados ElNinosmooth vem incluidos no pacote rainbow plot(ElNinosmooth, plot.type = 'functions', plotlegend = TRUE) Então eu quero fazer um grafico similar, onde as diferentes series sejam os anos da serie ( de cores diferentes), o eixo* x =* meses do ano; e o eixo *y*=QVENDA. #eu tentei transformar meus dados numa fts assim: OCTOPUS.fts <- fts(OCTOPUS$MONTH, OCTOPUS$QVENDA, frequency = 12, start = 2003, xname= 'Month', yname = 'Tonnes') # e logo tentei fazer o grafico como o do exemplo assim: plot(OCTOPUS.fts, plot.type = 'functions', plotlegend = TRUE) # mas deu-me um error Error in is.language(x) : 'x' is missing mesmo assim saiu um grafico, nem com cores nem nada de jeito! Alguma dica, ajuda? obrigado Carlos El 4 de abril de 2012 17:37, Benilton Carvalho <beniltoncarvalho@gmail.com>escribió:
Carlos,
Eu nao verifiquei detalhes do teu conjunto de dados, entretanto, imagino que necessites formata-lo apropriadamente (uma formatacao distinta para cada questao).
Para a primeira, assumindo que teu interesse seja 'QVENDA', eu comecaria com algo como:
dadosOriginais = read.csv('OCTOPUS.csv') library(reshape2) porBarco = dcast(dadosOriginais, DAY + MONTH + YEAR ~ CodEmbarc, value.var='QVENDA', fun.aggregate=mean) x0 = 1:nrow(porBarco) y0 = as.matrix(porBarco[, -(1:3)]) library(rainbow) obj0 = fts(x0, y0) plot(obj0)
E usaria a mesma estrategia para os demais.
benilton _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605
A dica e' transformar os dados apropriadamente como mostrei na mensagem anterior. Ao inves de transformar como funcao de embarcacao, transforme como funcao de ano. b 2012/4/5 carlos pombo sonderblohm <c.sonderblohm@gmail.com>:
ok benilton, obrigado, mas a intenção de usar a libraria rainbow e para fazer um grafico como esse daqui que conseguí num paper do R Journal:
#os dados ElNinosmooth vem incluidos no pacote rainbow
plot(ElNinosmooth, plot.type = 'functions', plotlegend = TRUE)
Então eu quero fazer um grafico similar, onde as diferentes series sejam os anos da serie ( de cores diferentes), o eixo x = meses do ano; e o eixo y=QVENDA.
#eu tentei transformar meus dados numa fts assim:
OCTOPUS.fts <- fts(OCTOPUS$MONTH, OCTOPUS$QVENDA, frequency = 12, start = 2003, xname= 'Month', yname = 'Tonnes')
# e logo tentei fazer o grafico como o do exemplo assim:
plot(OCTOPUS.fts, plot.type = 'functions', plotlegend = TRUE)
# mas deu-me um error
Error in is.language(x) : 'x' is missing
mesmo assim saiu um grafico, nem com cores nem nada de jeito!
Alguma dica, ajuda?
obrigado
Carlos
El 4 de abril de 2012 17:37, Benilton Carvalho <beniltoncarvalho@gmail.com> escribió:
Carlos,
Eu nao verifiquei detalhes do teu conjunto de dados, entretanto, imagino que necessites formata-lo apropriadamente (uma formatacao distinta para cada questao).
Para a primeira, assumindo que teu interesse seja 'QVENDA', eu comecaria com algo como:
dadosOriginais = read.csv('OCTOPUS.csv') library(reshape2) porBarco = dcast(dadosOriginais, DAY + MONTH + YEAR ~ CodEmbarc, value.var='QVENDA', fun.aggregate=mean) x0 = 1:nrow(porBarco) y0 = as.matrix(porBarco[, -(1:3)]) library(rainbow) obj0 = fts(x0, y0) plot(obj0)
E usaria a mesma estrategia para os demais.
benilton _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Novamente, um outro exemplo: dadosOriginais = read.csv('OCTOPUS.csv') library(reshape2) porAno = dcast(dadosOriginais, MONTH ~ YEAR, value.var='QVENDA', fun.aggregate=mean) x0 = 1:nrow(porAno) y0 = as.matrix(porAno[, -(1:3)]) library(rainbow) obj0 = fts(x0, y0) plot(obj0, plot.type='functions', plotlegend=TRUE) altere o codigo acima para os seus casos especificos. benilton
Este arquivo OCTOPUS está disponível em algum lugar? Novamente, um outro exemplo: dadosOriginais = read.csv('OCTOPUS.csv') library(reshape2) porAno = dcast(dadosOriginais, MONTH ~ YEAR, value.var='QVENDA', fun.aggregate=mean) x0 = 1:nrow(porAno) y0 = as.matrix(porAno[, -(1:3)]) library(rainbow) obj0 = fts(x0, y0) plot(obj0, plot.type='functions', plotlegend=TRUE) altere o codigo acima para os seus casos especificos. benilton
participantes (3)
-
Benilton Carvalho -
carlos pombo sonderblohm -
Mauro Sznelwar