Блог

Я получил оценки и соответствующие 95% — ные доверительные интервалы из регрессионной модели, которую я пытаюсь построить.

Данные выглядят следующим образом:

 library(tidyverse)

mydata <- structure(list(term = structure(c(1L, 7L, 18L, 19L, 20L), .Label = c("Age (years)", 
"Sex (male)", "Never smoker (reference)", "Current smoker", "Former smoker", 
"Obesity", "BMI (kg/m^2)", "Diabetes", "Glucose (mmol/L)", "Glucose lowering medication use", 
"Hypertension", "Systolic blood pressure (mmHg)", "Diastolic blood pressure (mmHg)", 
"Antihypertensive medication use", "Hypercholesterolemia", "LDL cholesterol (mmol/L)", 
"Lipid lowering medication use", ">90 (reference)", "60-89", 
"<60"), class = c("ordered", "factor")), estimate = c(0.4, 0.9, 
1, 1.5, 1.9), conf_low = c(0.2, 1.4, 0.95, 1, 1.7), conf_high = c(0.6, 
0.6, 1.05, 2, 2.1)), row.names = c(NA, -5L), class = c("tbl_df", 
"tbl", "data.frame"))

mydata
# A tibble: 5 x 4
  term            estimate conf_low conf_high
  <ord>              <dbl>    <dbl>     <dbl>
1 Age (years)          0.4     0.2       0.6 
2 BMI (kg/m^2)         0.9     1.4       0.6 
3 >90 (reference)      1       0.95      1.05
4 60-89                1.5     1         2   
5 <60                  1.9     1.7       2.1
 

Я построил эти цифры следующим образом:

 ggplot(data=mydata,
  aes(x=estimate, y=fct_rev(term)))  
  geom_point()  
  geom_errorbarh(aes(xmin=conf_low, xmax=conf_high, height=0.15))  
  annotation_custom(grob=grid::textGrob(label="Chronic kidney disease", 
                                        gp=grid::gpar(fontface="bold", fontsize=11), 
                                        hjust=1.0), 
                    xmin=-Inf, xmax=-Inf, ymin=3.2, ymax=3.2)   
  coord_cartesian(clip="off")   
  #scale_x_continuous(expand=expansion(mult=c(.05, .3)))  
  theme(axis.text.y=element_text(margin=margin(t=0, r=2.2, b=0, l=40, "pt")))  
  scale_y_discrete(name="") 
 

Я хотел бы добиться следующих вещей, которые я не могу заставить работать:

1. В метке BMI на оси Y ^2 должен быть надстрочным индексом. Я пробовал expression() , но, похоже, это не работает.
2. В обозначении переменной Хронического заболевания почек на оси y >90 должно быть знаком «больше или равно», а > Здесь я попробовал выражение и однократные коды, но оба, похоже, не работают.

Отредактируйте добавленную информацию о сеансе, так как это может быть полезно:

 library(utils)
sessionInfo(package=NULL)

R version 4.1.2 (2021-11-01)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

Matrix products: default

locale:
[1] LC_COLLATE=Dutch_Netherlands.1252  LC_CTYPE=Dutch_Netherlands.1252    LC_MONETARY=Dutch_Netherlands.1252
[4] LC_NUMERIC=C                       LC_TIME=Dutch_Netherlands.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] forcats_0.5.1   stringr_1.4.0   dplyr_1.0.7     purrr_0.3.4     readr_2.1.0     tidyr_1.1.4     tibble_3.1.6   
[8] ggplot2_3.3.5   tidyverse_1.3.1

loaded via a namespace (and not attached):
 [1] Rcpp_1.0.7       cellranger_1.1.0 pillar_1.6.4     compiler_4.1.2   dbplyr_2.1.1     tools_4.1.2     
 [7] digest_0.6.28    jsonlite_1.7.2   lubridate_1.8.0  lifecycle_1.0.1  gtable_0.3.0     pkgconfig_2.0.3 
[13] rlang_0.4.12     reprex_2.0.1     cli_3.1.0        rstudioapi_0.13  DBI_1.1.1        haven_2.4.3     
[19] xml2_1.3.2       withr_2.4.2      httr_1.4.2       fs_1.5.0         generics_0.1.1   vctrs_0.3.8     
[25] hms_1.1.1        grid_4.1.2       tidyselect_1.1.1 glue_1.5.0       R6_2.5.1         fansi_0.5.0     
[31] readxl_1.3.1     farver_2.1.0     tzdb_0.2.0       modelr_0.1.8     magrittr_2.0.1   backports_1.3.0 
[37] scales_1.1.1     ellipsis_0.3.2   rvest_1.0.2      assertthat_0.2.1 colorspace_2.0-2 labeling_0.4.2  
[43] utf8_1.2.2       stringi_1.7.5    munsell_0.5.0    broom_0.7.10     crayon_1.4.2 

 

Как вы добавили юникод?

 "BMI (kg/mu00b2)"
"u2265 90 (reference)" or "u226590 (reference)"
"u2264 60" or "u226460"
 

Похоже, это работает.