2

R3.3.2とパッケージtidyverseを用いて前処理中です。
以下のようなデータなのですが、同じidについて複数回調査(chousaが調査回数)しています。
同じ調査項目(今回はQ1)が0から1に変わってから1回目、2回目、…4回目までのデータがほしいです。
どのようにmutateすればよいでしょうか?
よろしくお願いします。

ID Q1 Q6 chousa …
1    0   3    1
1    1  23     2 
1    1
2    0  27     3 
2    0  28     4 
2    1  29     5

これの場合は

ID change1 change2 … 
 1     0        0
 1   1        0
 1      0        1
 2      0        0
 2      0        0
 2      1        0

となるようにしたいです。

追記
同一IDについてQ1が0,0,1,1,1,1,0と来たら

change1:0,0,1,0,0,0,0
change2:0,0,0,1,0,0,0
change3:0,0,0,0,1,0,0
change4:0,0,0,0,0,1,0
change5:0,0,0,0,0,0,0

となるようにしたいです。

mutate(change1=ifelse(Q1==1&Q1-lag(Q1)==1,1,0),
change2=ifelse(Q1==1&lag(change1)==1,1,0),
change3=ifelse(Q1==1&lag(change2)==1,1,0),
change4=ifelse(Q1==1&lag(change3)==1,1,0))

で行けるような気がしますが

id chousa change1 change2 change3
 5      2       1       0       0
 5      3       0       1       0
 5      4       0       0       1

となってほしいものが

id chousa change1 change2 change3
 5      2       0       0       1
 5      3       0       1       0
 5      4       1       0       0

と反対になってしまいます。group_byした時にIDの昇順降順が逆になっているようなのですがどうしたら良いのでしょうか…

→group_by(ID)だけでなくその後にarrange(ID,chousa) で解決しました。

3
  • 求めてるデータについてですが,change1はQ1がそのままだと思うのですが,change2がどのようなデータなのでしょうか。change2はchange1(Q1)が1に切り替わった際のchousaを0とし,以降の同一IDをそこから4つまで数え上げる,という感じでしょうか。
    – kazutan
    Commented 2017年1月25日 11:16
  • すいません。間違いなので修正します。
    – user87562
    Commented 2017年1月25日 11:17
  • change2は1に切り替わって2回連続で1になったときに1になります。そうでないと0です。(1回めや3回目以降は0)。 change3は3回連続で1になったとき1,そうでないときは0です。
    – user87562
    Commented 2017年1月25日 11:54

1 件の回答 1

0

質問がきちんと理解できていませんが、次のようにすれば希望の形が出力されます。
(change1、change2の判定がわかりにくい)

IDごとに処理するのがミソになるのだと思います。

library(tidyverse)
data_frame(
     ID    = c(rep(1, 3), rep(2, 3)),
     Q1    = c(0, 1, 1, 0, 0, 1),
     chosa = c(1, 2, NA, 3, 4, 5)
) %>% 
  # IDごとに集計
  group_by(ID) %>% 
    # change1: 対象行のQ1と対象行から次の行のQ1の値を引く
    #   このときIDの先頭行はNAとなり、連続値の時には0となる
    # change2: change1の値とQ1の値の差分をとり、Q1に同じ値が続く場合に1となるようにする
  mutate(change1 = Q1 - lag(Q1),
         change2 = Q1 - change1) %>% 
  ungroup() %>% 
  # NAを0として処理する
  mutate_at(num_range("change", 1:2), funs(ifelse(test = is.na(.), yes = 0, no = .))) %>% 
  select(ID, num_range("change", 1:2))

ロジックが理解できていないので、データ数が増えるとダメになるかもしれません。ご参考までに。

3
  • ありがとうございます。試してみましたが、やはりロジックに漏れがあり現在なんとかして改良しようと試みています。出来たら書き込ませていただきます。
    – user87562
    Commented 2017年1月25日 13:37
  • funsのところがchangeだけでなくデータすべてのNAを0にしてしまっている気がします…あまり使ったことがないので自信ないですが…
    – user87562
    Commented 2017年1月25日 14:43
  • やはりmutate_atのところで全てのNAを0にしてしまっていました。
    – user87562
    Commented 2017年1月26日 9:14

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

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