0

Rで生成したデータの前処理
で以前質問したのですが、dplyrを用いた場合の最後のmutateのところがうまく行きません。

params <- expand.grid(p = 1:2, q = 1:1, r = 1:2, s = 1:1)
ans <- params %>% rowwise() %>% mutate(aic = f(p, q, r, s))

関数fの部分は単独で変数を指定(コンソールにf(1,1,1,1)などと打ち込む)するとしっかり

返ってきて問題はないようですが、返り値がmatrixになっているためas.data.frameで直すと変なくっつき方をしました(具体的にうまく説明できない)。
そのため、しっかり横につなげようとmatrixを転地したあとくっつけようとすると

Error: incompatible size (4), expecting 1 (the group size) or 1

のエラーが出ました。そのため、matrixの中身だけ取り出そうと

ans <- params %>% rowwise() %>% mutate(aic = f(p, q, r, s)$[1,1]) %>% mutate(bic=f(p, q, r, s)$[1,2])

とやるとこの場合は$[]を用いてはいけないようでエラーになりました。

詳しい方、ヘルプお願い致します。

1 件の回答 1

1

そのエラーはおそらく、dplyrは関係なく文法上のミスです。$[]ではなく[]が正しいです。

fはこういう結果を返す関数ということであっていますか?

f <- function(x,y,z,v) {
  m <- matrix(c(4.108074, 4.112013, 4.108074, 4.109363))
  rownames(m) <- c("Akaike", "Bayes", "Shibata", "Hannan-Quinn")
  m
}
#> f(1,1,1,1)
#>                  [,1]
#> Akaike       4.108074
#> Bayes        4.112013
#> Shibata      4.108074
#> Hannan-Quinn 4.109363

もしそうなら、↓のように[]を使えば私の手元ではエラーなく動きましたが、いかがでしょうか。

ans <- params %>%
  rowwise() %>%
  mutate(aic = f(p, q, r, s)[1,1]) %>%
  mutate(bic = f(p, q, r, s)[1,2])

ちなみに、この書き方だとf(p, q, r, s)を複数回計算しなおすことになるので、あまり効率よくありません。以下のようにすればいちど計算結果を格納して、そこから値を取り出すことができます。
あんまりスマートな書き方じゃないですが、今のところこれしか思いつかないのでひとまずこれを回答とさせてください。すみません…

ans <- params %>%
  rowwise() %>%
  mutate(result = list(f(p, q, r, s)),
         aic = result[1,1],
         bic = result[2,1])
ans
#> Source: local data frame [4 x 7]
#> Groups: <by row>
#> 
#>       p     q     r     s     result      aic      bic
#>   (int) (int) (int) (int)     (list)    (dbl)    (dbl)
#> 1     1     1     1     1 <dbl[4,1]> 4.108074 4.112013
#> 2     2     1     1     1 <dbl[4,1]> 4.108074 4.112013
#> 3     1     1     2     1 <dbl[4,1]> 4.108074 4.112013
#> 4     2     1     2     1 <dbl[4,1]> 4.108074 4.112013
2
  • ありがとうございます。無事動作しました。 また、最後のmutateのところですがbic = result[1,2]ですね。
    – user87562
    Commented 2016年3月7日 7:29
  • あ、すみません。。指摘ありがとうございます。 Commented 2016年3月8日 9:22

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。