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