2

...スプライン-（最良）または線形補間（OK）、または四半期全体で繰り返される値（細かい）のいずれか。問題は、返されるデータ型を使用可能なものgetFin()に変換する方法がわからないことです。これが私のコードです：viewFin()timeSeries

getFin('F')
x <- viewFin(F.f, "BS", period="Q")["Total Common Shares Outstanding",]*1000

私の希望する出力は

> x
GMT     x.ts
2011-09-01  3816000
2011-08-01  3816000
2011-07-01  3816000
2011-06-01  3815000
2011-05-01  3815000
2011-04-01  3815000
2011-03-01  3813000
2011-02-01  3813000
2011-01-01  3813000
2010-12-01  3778000
2010-11-01  3778000
2010-10-01  3778000
2010-09-01  3484000

ただし、実際の出力は次のとおりです。

> x
2011-09-30 2011-06-30 2011-03-31 2010-12-31 2010-09-30 
   3816000    3815000    3813000    3778000    3484000 
> str(x)
 Named num [1:5] 3816000 3815000 3813000 3778000 3484000
 - attr(*, "names")= chr [1:5] "2011-09-30" "2011-06-30" "2011-03-31" "2010-12-31" ...

xオブジェクトは奇妙な逆形式であるように見えます。ここで、キーは数値であり、値は日付の文字列です。日付または数値コンポーネントを抽出しようとすると、数値部分を分離して時系列オブジェクトを生成できません。

理想的には、私の希望する出力に到達するために、私は言うことができるでしょう

mydates <- timeSequence(from = "2011-01-01", to=Sys.Date(), by = "month")
series <- timeSeries(x$data, mydates)

しかし、数値データ部分を抽出できないようです。

アップデート

こことここから、私は次のコードを適応させました：

getFin('F')
x <- viewFin(F.f, "BS", period="Q")["Total Common Shares Outstanding",]*1000
zoox = zoo(x, order.by=as.Date(names(x)))
x2 <- na.spline(merge(zoox, foo=zoo(NA, order.by=seq(start(zoox), end(zoox), "month")))[, 1])

ただし、私の出力は日付を少し混乱させ、補間を台無しにします。

>x2
2010-09-30 2010-10-30 2010-11-30 2010-12-30 2010-12-31 2011-01-30 2011-03-02 
   3484000    3623591    3720509    3776671    3778000    3804738    3813071 
2011-03-30 2011-03-31 2011-04-30 2011-05-30 2011-06-30 2011-07-30 2011-08-30 
   3813025    3813000    3813100    3814116    3815000    3814976    3814884 
2011-09-30 
   3816000

ご覧のとおり、2011年3月には12-30と12-31の両方の値がありますが、2月などはありません。これを解決するにはどうすればよいですか。

4

2 に答える 2

1

私はちょうどこれと同じ問題を抱えていました-私はそれを実現するために2つの関数を書きました。

まず、日付関数で日付を移動する必要があります。おそらく月末の日付が必要ですが、Rは最初の日付をはるかに簡単に処理します。

日付関数

toLastDay <- function(dateObj, monAdv=0, ToFirst = FALSE)
# takes date object, transforms and returns a dat object
{
  tt <- as.POSIXlt(dateObj)
  tt$mon <- tt$mon + monAdv # moves the month
  tt$mday <- 1L             # make date the first
  if(ToFirst) {             
      tt <- as.Date(tt)
  } else { 
      tt$mon <- tt$mon + 1L # go to the first of the next month
      tt <- as.Date(tt) - 1L # subtract one day, yielding the last of prior month
  }
  return(tt)
}

ここで、データを取得して線形interpを実行する関数が必要です。私はとでna.approx関数を使用しzooますxts。ここでは、おそらくqtrの中間の日付が必要です。私の経験では、ほとんどのqtrlyデータはqtrの中間と考えるのが最適です。

# make Qtrly monthly -- with obsv in the mid qtr
qtr2Mon <- function(QD) 
{
    fromD <- toLastDay(index(first(QD)), -2L, ToFirst = TRUE)
    toD <- toLastDay(index(last(QD)), ToFirst = TRUE)
    q2m_dates <- toLastDay(seq(fromD, toD, by = 'mon'))
    emptyX <- xts(, q2m_dates)
    QD_adj <- QD
    index(QD_adj) <- toLastDay(index(QD), monAdv = -1L)
    mm <- merge(emptyX, QD_adj)
    mm_filled <- na.approx(mm)
    return(mm_filled)
}

まず、qtrlyxtsオブジェクトを作成します-qtrlyGDPなどと考えてください

Sys.setenv(TZ = 'GMT')
require(xts)

qdates <- seq(as.Date("2000-03-01"), as.Date("2013-06-01"), by = "3 mon")
qdates <- toLastDay(qdates)
qdata <- rnorm(length(qdates), mean = 1)
qtr_XTS <- xts(qdata, order.by = qdates)

これで、線形成長などを想定して、上記のライブラリ関数を使用して月次に変換できます。

mon_fromQtr <- qtr2Mon(qtr_XTS)

終わり！

于 2013-05-26T10:33:54.340 に答える
0

更新2：

より良い答えを投稿してください！これは本当に醜いですが、これが私が受け入れられるものを得た方法です：

getFin('F')
x <- viewFin(F.f, "BS", period="Q")["Total Common Shares Outstanding",]*1000
zoox = zoo(x, order.by=as.Date(names(x)))
foo=zoo(NA, order.by=seq(as.Date(as.character(timeFirstDayInMonth(start(zoox)))), as.Date(as.character(timeFirstDayInMonth(end(zoox)))), "month"))
foo2 <- na.approx(merge(zoox, foo)[, 1])
fx <- merge(foo2,foo, all=FALSE)[,1]

変換する

> x
2011-09-30 2011-06-30 2011-03-31 2010-12-31 2010-09-30 
   3816000    3815000    3813000    3778000    3484000

の中へ

> fx
2010-10-01 2010-11-01 2010-12-01 2011-01-01 2011-02-01 2011-03-01 2011-04-01 
   3487196    3586261    3682130    3778389    3790444    3801333    3813022 
2011-05-01 2011-06-01 2011-07-01 2011-08-01 2011-09-01 
   3813681    3814363    3815011    3815348    3815685

私はこれがあまりにも醜くて真実ではないと思うので、私はこの答えを投稿しますが、他の誰かにもっとハンサムなものを投稿してもらいたいです。

于 2011-11-18T17:13:30.133 に答える