下記のようなデータフレームで、男女別にBPの四分位変数(BPq)を作成し、それぞれの下位25%を新変数(lowBP)= 1(その他は0)として作成したいです。
df <- data.frame(
ID = 1:4, age = c(43, 62, 54, 55), sex = c(0, 1, 1, 0), group = c(1, 2, 3, 3)
BP = c(120, 130, 132, 110), BMI = c(21, 26, 23, 19))
以下のコードで試みていますが、正確に分類(男女別に)できているかが不明です。変数の確認方法までご教示いただけると助かります。
# 男女別のBP四分位変数を作成
df <- df |> group_by(sex) |> mutate(BPq = ntile (BP,4))
# 男女それぞれの下位25%をlowBP=1と定義
df <- df |> mutate(lowBP = case_when(BPq == 1 ~ 1, TRUE ~ 0 ))