如何將 ggplot 與 Haven 數據集一起使用

Question

編碼和 R 的新手，但有一個 STATA 數據集，我想使用 ggplot 來查看我的數據，但是，我遇到了多個錯誤，例如

沒有適用於“重新縮放”的方法應用於 class 的 object "c('haven_labelled', 'vctrs_vctr', 'double')"

我不知道如何轉換它們，所以我可以 plot 進行可視化，

代碼行如下：

Data <- read_dta("longitudinal_td.dta")


Data <- Data %>%
  select(pidp,wave,age_dv,sex_dv,ethn_dv,sf1_dv,bmi_dv,sf12pcs_dv,fihhmnnet1_dv,sf12mcs_dv) %>%
  filter(wave == "1", age_dv<=50)%>%
  mutate(pipd = row_number(),age=age_dv, sex=sex_dv, ethnicity = ethn_dv, general_health=sf1_dv,
        bmi=bmi_dv, physical_component_score=sf12pcs_dv, mental_component_score=sf12mcs_dv, household_income=fihhmnnet1_dv)%>%
select(-pipd,-age_dv,-sex_dv,-ethn_dv,-sf1_dv,-bmi_dv,-sf12pcs_dv,-sf12mcs_dv,-fihhmnnet1_dv)

我希望這是正確的，這是 dput：

本質上我只是想探索 BMI，但我不知道我是否只能 plot 這些或必須將數字分配給 label 就像它已經在避風港標簽中完成一樣

dput(head(Data))
structure(list(pidp = structure(c(68001367, 68006127, 68008167, 
68009527, 68010207, 68010887), label = "cross-wave person identifier (public release)", format.stata = "%12.0g"), 
    wave = structure(c(1, 1, 1, 1, 1, 1), label = "interview wave", format.stata = "%8.0g"), 
    age = structure(c(39, 39, 38, 31, 24, 45), label = "Age, derived from dob_dv and intdat_dv", format.stata = "%8.0g"), 
    sex = structure(c(1, 2, 2, 1, 2, 2), label = "Sex, derived", format.stata = "%8.0g", labels = c(Male = 1, 
    Female = 2), class = c("haven_labelled", "vctrs_vctr", "double"
    )), ethnicity = structure(c(1, 1, 1, 1, 1, 1), label = "Ethnic group (derived from multiple sources)", format.stata = "%8.0g", labels = c(`white uk` = 1, 
    irish = 2, `gypsy or irish traveller` = 3, `any other white background` = 4, 
    `white and black caribbean` = 5, `white and black african` = 6, 
    `white and asian` = 7, `any other mixed background` = 8, 
    indian = 9, pakistani = 10, bangladeshi = 11, chinese = 12, 
    `any other asian background` = 13, caribbean = 14, african = 15, 
    `any other black background` = 16, arab = 17, `any other ethnic group` = 97
    ), class = c("haven_labelled", "vctrs_vctr", "double")), 
    general_health = structure(c(2, 4, 5, 3, 1, 1), label = "General health", format.stata = "%8.0g", labels = c(excellent = 1, 
    `very good` = 2, good = 3, fair = 4, `or Poor?` = 5), class = c("haven_labelled", 
    "vctrs_vctr", "double")), bmi = structure(c(29.6, 38.8, 21.5, 
    24.2, 25, 25.5), label = "Body Mass Index", format.stata = "%12.0g")

Answer 1

感謝您使用dput()發布數據示例。 您發布的數據格式表明它以某種方式變成了列表而不是數據框。 您需要將其轉換為數據框 - 當您使用tidyverse haven使用as_tibble()來完成。

同樣，您需要標簽而不是基礎整數。 您可以簡單地將as_factor應用於整個數據框來執行此操作。

然后，您的數據就可以通過管道傳輸到ggplot2 。 例如：

library(dplyr)
library(ggplot2)
library(haven)


Data  |> 
    as_tibble()   |>
    as_factor()  |>
    ggplot() +
    geom_boxplot(aes(x=sex, y=bmi))