Introduction

In this section we will:
* identify which Basic Metabolic Panel procedures were excluded and why * identify which different Basic Metabolic Panel subgroups were created and why

If you have questions or concerns about this data please contact Alexander Nielson ()

Load Libraries

Load Libraries

library(data.table)
library(tidyverse)
## Warning: replacing previous import 'vctrs::data_frame' by 'tibble::data_frame'
## when loading 'dplyr'
library(stringi)
library(ggridges)
library(broom)
library(disk.frame)
library(RecordLinkage)
library(googlesheets4)
library(bigrquery)
library(DBI)
devtools::install_github("utah-osa/hcctools2", upgrade="never" )
library(hcctools2)

Establish color palettes

cust_color_pal1 <- c(
        "Anesthesia" = "#f3476f",
        "Facility" = "#e86a33",
        "Medicine Proc" = "#e0a426",
        "Pathology" = "#77bf45",
        "Radiology" = "#617ff7",
        "Surgery" = "#a974ee"
    )

cust_color_pal2 <- c(
        "TRUE" = "#617ff7",
        "FALSE" = "#e0a426"
    )

cust_color_pal3 <- c(
        "above avg" = "#f3476f",
        "avg" = "#77bf45",
        "below avg" = "#a974ee"
    )



  fac_ref_regex <- "(UTAH)|(IHC)|(HOSP)|(HOSPITAL)|(CLINIC)|(ANESTH)|(SCOPY)|(SURG)|(LLC)|(ASSOC)|(MEDIC)|(CENTER)|(ASSOCIATES)|(U OF U)|(HEALTH)|(OLOGY)|(OSCOPY)|(FAMILY)|(VAMC)|(SLC)|(SALT LAKE)|(CITY)|(PROVO)|(OGDEN)|(ENDO )|( VALLEY)|( REGIONAL)|( CTR)|(GRANITE)|( INSTITUTE)|(INSTACARE)|(WASATCH)|(COUNTY)|(PEDIATRIC)|(CORP)|(CENTRAL)|(URGENT)|(CARE)|(UNIV)|(ODYSSEY)|(MOUNTAINSTAR)|( ORTHOPEDIC)|(INSTITUT)|(PARTNERSHIP)|(PHYSICIAN)|(CASTLEVIEW)|(CONSULTING)|(MAGEMENT)|(PRACTICE)|(EMERGENCY)|(SPECIALISTS)|(DIVISION)|(GUT WHISPERER)|(INTERMOUNTAIN)|(OBGYN)"

Connect to GCP database

bigrquery::bq_auth(path = 'D:/gcp_keys/healthcare-compare-prod-95b3b7349c32.json')

# set my project ID and dataset name
project_id <- 'healthcare-compare-prod'
dataset_name <- 'healthcare_compare'

con <- dbConnect(
  bigrquery::bigquery(),
  project = project_id,
  dataset = dataset_name,
  billing = project_id
)

Get NPI table

query <-  paste0("SELECT npi, clean_name, osa_group, osa_class, osa_specialization
                 FROM `healthcare-compare-prod.healthcare_compare.npi_full`")
                       
#bq_project_query(billing, query) # uncomment to determine billing price for above query.

npi_full <- dbGetQuery(con, query) %>%
  data.table() 

get a subset of the NPI providers based upon taxonomy groups

gs4_auth(email="alexnielson@utah.gov")
  npi_prov_pair <-  npi_full %>% 
    .[,.(npi=npi,
         clean_name = clean_name
         )
      ] 

Load Data

bun_proc <-  disk.frame("full_apcd.df")
bmp <- bun_proc[proc_code_str_sorted %>% stri_detect_regex("80048") & 
                  tp_med > 0 & 
                  tp_med < 2000 &
                  surg_bun_sum_med == 0 &
                  anes_bun_sum_med == 0 &
                  radi_bun_sum_med == 0 &
                # stri_length(proc_code_str_sorted) < 27 &
                duration_mean==0]
bmp %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 338 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 faci_bun_t_dept            0.705
##  2 faci_bun_t_emergency       0.705
##  3 medi_bun_t_inject          0.413
##  4 medi_bun_t_add-on          0.380
##  5 medi_bun_t_hydrate         0.336
##  6 faci_bun_t_visit           0.240
##  7 path_bun_t_complete        0.237
##  8 medi_bun_t_hydration       0.234
##  9 medi_bun_t_normal          0.200
## 10 medi_bun_t_saline          0.200
## # ... with 328 more rows
cols_to_check <-  bmp %>% select_if(is.logical)
tag_frequency <- apply(cols_to_check, 2, function(i) sum(i > 0)) %>% enframe() %>% as.data.table() %>% .[,.(tag=name, freq=value)]
outlier_tag_list <- tag_frequency[freq < 100] %>% .[["tag"]]
bmp_clean <- remove_bun_proc_with_rare_tags(outlier_tag_list=outlier_tag_list, df=bmp)
## [1] "Start of rare tag removal"
## [1] "nrow:34263 , ncol:1404"
## [1] "nrow:34263"
## [1] "removing: surg_bun_t_a-grft"
## [1] 34263
## [1] "removing: surg_bun_t_abd"
## [1] 34263
## [1] "removing: surg_bun_t_abdom"
## [1] 34263
## [1] "removing: surg_bun_t_abdomen"
## [1] 34263
## [1] "removing: surg_bun_t_ablate"
## [1] 34263
## [1] "removing: surg_bun_t_abortion"
## [1] 34263
## [1] "removing: surg_bun_t_abscess"
## [1] 34263
## [1] "removing: surg_bun_t_acell"
## [1] 34263
## [1] "removing: surg_bun_t_acellular"
## [1] 34263
## [1] "removing: surg_bun_t_acne"
## [1] 34263
## [1] "removing: surg_bun_t_acute"
## [1] 34263
## [1] "removing: surg_bun_t_add-on"
## [1] 34263
## [1] "removing: surg_bun_t_adenoid"
## [1] 34263
## [1] "removing: surg_bun_t_adjust"
## [1] 34263
## [1] "removing: surg_bun_t_airway"
## [1] 34263
## [1] "removing: surg_bun_t_amput"
## [1] 34263
## [1] "removing: surg_bun_t_anal"
## [1] 34263
## [1] "removing: surg_bun_t_analysis"
## [1] 34263
## [1] "removing: surg_bun_t_aneurysm"
## [1] 34263
## [1] "removing: surg_bun_t_ankle"
## [1] 34263
## [1] "removing: surg_bun_t_anoscopy"
## [1] 34263
## [1] "removing: surg_bun_t_ant"
## [1] 34263
## [1] "removing: surg_bun_t_antepartum"
## [1] 34263
## [1] "removing: surg_bun_t_aorta"
## [1] 34263
## [1] "removing: surg_bun_t_aortic"
## [1] 34263
## [1] "removing: surg_bun_t_aply"
## [1] 34263
## [1] "removing: surg_bun_t_appendectomy"
## [1] 34263
## [1] "removing: surg_bun_t_application"
## [1] 34263
## [1] "removing: surg_bun_t_apply"
## [1] 34263
## [1] "removing: surg_bun_t_approach"
## [1] 34263
## [1] "removing: surg_bun_t_arm"
## [1] 34263
## [1] "removing: surg_bun_t_art"
## [1] 34263
## [1] "removing: surg_bun_t_arterial"
## [1] 34263
## [1] "removing: surg_bun_t_artery"
## [1] 34263
## [1] "removing: surg_bun_t_arthro"
## [1] 34263
## [1] "removing: surg_bun_t_arthroplasty"
## [1] 34263
## [1] "removing: surg_bun_t_arthroscopy"
## [1] 34263
## [1] "removing: surg_bun_t_artificial"
## [1] 34263
## [1] "removing: surg_bun_t_augmentation"
## [1] 34263
## [1] "removing: surg_bun_t_autograft"
## [1] 34263
## [1] "removing: surg_bun_t_b9"
## [1] 34263
## [1] "removing: surg_bun_t_back"
## [1] 34263
## [1] "removing: surg_bun_t_balloon"
## [1] 34263
## [1] "removing: surg_bun_t_bed"
## [1] 34263
## [1] "removing: surg_bun_t_biceps"
## [1] 34263
## [1] "removing: surg_bun_t_bile"
## [1] 34263
## [1] "removing: surg_bun_t_bilobectomy"
## [1] 34263
## [1] "removing: surg_bun_t_biopsy"
## [1] 34263
## [1] "removing: surg_bun_t_bladder"
## [1] 34263
## [1] "removing: surg_bun_t_bleeding"
## [1] 34263
## [1] "removing: surg_bun_t_block"
## [1] 34263
## [1] "removing: surg_bun_t_blockage"
## [1] 34263
## [1] "removing: surg_bun_t_blood"
## [1] 34262
## [1] "removing: surg_bun_t_body"
## [1] 34257
## [1] "removing: surg_bun_t_bone"
## [1] 34257
## [1] "removing: surg_bun_t_bowel"
## [1] 34257
## [1] "removing: surg_bun_t_brace"
## [1] 34257
## [1] "removing: surg_bun_t_brain"
## [1] 34257
## [1] "removing: surg_bun_t_breast"
## [1] 34257
## [1] "removing: surg_bun_t_bronchial"
## [1] 34257
## [1] "removing: surg_bun_t_bronchoscopy"
## [1] 34257
## [1] "removing: surg_bun_t_brow"
## [1] 34257
## [1] "removing: surg_bun_t_bulge"
## [1] 34257
## [1] "removing: surg_bun_t_burn"
## [1] 34257
## [1] "removing: surg_bun_t_burr"
## [1] 34257
## [1] "removing: surg_bun_t_bursa"
## [1] 34257
## [1] "removing: surg_bun_t_buttck"
## [1] 34257
## [1] "removing: surg_bun_t_buttock"
## [1] 34257
## [1] "removing: surg_bun_t_byp"
## [1] 34257
## [1] "removing: surg_bun_t_bypass"
## [1] 34257
## [1] "removing: surg_bun_t_cabg"
## [1] 34257
## [1] "removing: surg_bun_t_canal"
## [1] 34257
## [1] "removing: surg_bun_t_cannula"
## [1] 34257
## [1] "removing: surg_bun_t_capillary"
## [1] 34257
## [1] "removing: surg_bun_t_capsule"
## [1] 34257
## [1] "removing: surg_bun_t_care"
## [1] 34257
## [1] "removing: surg_bun_t_carotid"
## [1] 34257
## [1] "removing: surg_bun_t_cartilage"
## [1] 34257
## [1] "removing: surg_bun_t_cast"
## [1] 34257
## [1] "removing: surg_bun_t_cath"
## [1] 34257
## [1] "removing: surg_bun_t_cavity"
## [1] 34257
## [1] "removing: surg_bun_t_cell"
## [1] 34257
## [1] "removing: surg_bun_t_cervix"
## [1] 34257
## [1] "removing: surg_bun_t_cesarean"
## [1] 34257
## [1] "removing: surg_bun_t_chamber"
## [1] 34257
## [1] "removing: surg_bun_t_change"
## [1] 34257
## [1] "removing: surg_bun_t_chemical"
## [1] 34257
## [1] "removing: surg_bun_t_chemodenerv"
## [1] 34257
## [1] "removing: surg_bun_t_chest"
## [1] 34257
## [1] "removing: surg_bun_t_chin"
## [1] 34257
## [1] "removing: surg_bun_t_cholangiogram"
## [1] 34257
## [1] "removing: surg_bun_t_cholangiopancreatograph"
## [1] 34257
## [1] "removing: surg_bun_t_cholecystoenterostomy"
## [1] 34257
## [1] "removing: surg_bun_t_chronic"
## [1] 34257
## [1] "removing: surg_bun_t_circulation"
## [1] 34257
## [1] "removing: surg_bun_t_circumcision"
## [1] 34257
## [1] "removing: surg_bun_t_clamp"
## [1] 34257
## [1] "removing: surg_bun_t_clavicle"
## [1] 34257
## [1] "removing: surg_bun_t_clearance"
## [1] 34257
## [1] "removing: surg_bun_t_cleft"
## [1] 34257
## [1] "removing: surg_bun_t_clitoris"
## [1] 34257
## [1] "removing: surg_bun_t_cloacal"
## [1] 34257
## [1] "removing: surg_bun_t_closed"
## [1] 34257
## [1] "removing: surg_bun_t_closure"
## [1] 34257
## [1] "removing: surg_bun_t_clot"
## [1] 34257
## [1] "removing: surg_bun_t_collar"
## [1] 34257
## [1] "removing: surg_bun_t_colon"
## [1] 34256
## [1] "removing: surg_bun_t_colonoscopy"
## [1] 34256
## [1] "removing: surg_bun_t_coloproctostomy"
## [1] 34256
## [1] "removing: surg_bun_t_colostomy"
## [1] 34256
## [1] "removing: surg_bun_t_complex"
## [1] 34256
## [1] "removing: surg_bun_t_compression"
## [1] 34256
## [1] "removing: surg_bun_t_contriction"
## [1] 34256
## [1] "removing: surg_bun_t_control"
## [1] 34256
## [1] "removing: surg_bun_t_cord"
## [1] 34256
## [1] "removing: surg_bun_t_cornea"
## [1] 34256
## [1] "removing: surg_bun_t_coronary"
## [1] 34256
## [1] "removing: surg_bun_t_coronoid"
## [1] 34256
## [1] "removing: surg_bun_t_correct"
## [1] 34256
## [1] "removing: surg_bun_t_cranial"
## [1] 34256
## [1] "removing: surg_bun_t_craniofacial"
## [1] 34256
## [1] "removing: surg_bun_t_cryotherapy"
## [1] 34256
## [1] "removing: surg_bun_t_cuff"
## [1] 34256
## [1] "removing: surg_bun_t_cult"
## [1] 34256
## [1] "removing: surg_bun_t_cyst"
## [1] 34256
## [1] "removing: surg_bun_t_cysto"
## [1] 34256
## [1] "removing: surg_bun_t_cystometrogram"
## [1] 34256
## [1] "removing: surg_bun_t_cystoscopy"
## [1] 34256
## [1] "removing: surg_bun_t_cystouretero"
## [1] 34256
## [1] "removing: surg_bun_t_debrid"
## [1] 34256
## [1] "removing: surg_bun_t_debride"
## [1] 34256
## [1] "removing: surg_bun_t_decompression"
## [1] 34256
## [1] "removing: surg_bun_t_deep"
## [1] 34256
## [1] "removing: surg_bun_t_defect"
## [1] 34256
## [1] "removing: surg_bun_t_defib"
## [1] 34256
## [1] "removing: surg_bun_t_deformity"
## [1] 34256
## [1] "removing: surg_bun_t_degenerated"
## [1] 34256
## [1] "removing: surg_bun_t_delivery"
## [1] 34256
## [1] "removing: surg_bun_t_deplet"
## [1] 34256
## [1] "removing: surg_bun_t_dermabrasion"
## [1] 34256
## [1] "removing: surg_bun_t_design"
## [1] 34256
## [1] "removing: surg_bun_t_destroy"
## [1] 34256
## [1] "removing: surg_bun_t_destruction"
## [1] 34256
## [1] "removing: surg_bun_t_dev"
## [1] 34256
## [1] "removing: surg_bun_t_device"
## [1] 34256
## [1] "removing: surg_bun_t_diagnostic"
## [1] 34256
## [1] "removing: surg_bun_t_dialysis"
## [1] 34256
## [1] "removing: surg_bun_t_diaphragm"
## [1] 34256
## [1] "removing: surg_bun_t_digit"
## [1] 34256
## [1] "removing: surg_bun_t_dilat"
## [1] 34256
## [1] "removing: surg_bun_t_dilation"
## [1] 34256
## [1] "removing: surg_bun_t_disc"
## [1] 34256
## [1] "removing: surg_bun_t_diskectomy"
## [1] 34256
## [1] "removing: surg_bun_t_dislocat"
## [1] 34256
## [1] "removing: surg_bun_t_dissect"
## [1] 34256
## [1] "removing: surg_bun_t_donar"
## [1] 34256
## [1] "removing: surg_bun_t_donor"
## [1] 34256
## [1] "removing: surg_bun_t_drain"
## [1] 34256
## [1] "removing: surg_bun_t_draw"
## [1] 34256
## [1] "removing: surg_bun_t_dressing"
## [1] 34256
## [1] "removing: surg_bun_t_drill"
## [1] 34256
## [1] "removing: surg_bun_t_drug"
## [1] 34256
## [1] "removing: surg_bun_t_duct"
## [1] 34256
## [1] "removing: surg_bun_t_duoden"
## [1] 34256
## [1] "removing: surg_bun_t_dura"
## [1] 34256
## [1] "removing: surg_bun_t_ear"
## [1] 34256
## [1] "removing: surg_bun_t_eardrum"
## [1] 34256
## [1] "removing: surg_bun_t_ectopic"
## [1] 34256
## [1] "removing: surg_bun_t_egd"
## [1] 34256
## [1] "removing: surg_bun_t_elbow"
## [1] 34256
## [1] "removing: surg_bun_t_electrical"
## [1] 34256
## [1] "removing: surg_bun_t_electrode"
## [1] 34256
## [1] "removing: surg_bun_t_electroejaculation"
## [1] 34256
## [1] "removing: surg_bun_t_electrolysis"
## [1] 34256
## [1] "removing: surg_bun_t_embolize"
## [1] 34256
## [1] "removing: surg_bun_t_embryo"
## [1] 34256
## [1] "removing: surg_bun_t_endo"
## [1] 34256
## [1] "removing: surg_bun_t_endoscopy"
## [1] 34256
## [1] "removing: surg_bun_t_endovas"
## [1] 34256
## [1] "removing: surg_bun_t_enlarge"
## [1] 34256
## [1] "removing: surg_bun_t_enterectomy"
## [1] 34256
## [1] "removing: surg_bun_t_epiderm"
## [1] 34256
## [1] "removing: surg_bun_t_epididymis"
## [1] 34256
## [1] "removing: surg_bun_t_epiphysis"
## [1] 34256
## [1] "removing: surg_bun_t_episiotomy"
## [1] 34256
## [1] "removing: surg_bun_t_esoph"
## [1] 34256
## [1] "removing: surg_bun_t_esophagoscopy"
## [1] 34256
## [1] "removing: surg_bun_t_esophagus"
## [1] 34256
## [1] "removing: surg_bun_t_exam"
## [1] 34255
## [1] "removing: surg_bun_t_exc"
## [1] 34255
## [1] "removing: surg_bun_t_excessive"
## [1] 34255
## [1] "removing: surg_bun_t_excise"
## [1] 34255
## [1] "removing: surg_bun_t_excision"
## [1] 34255
## [1] "removing: surg_bun_t_exostosis"
## [1] 34255
## [1] "removing: surg_bun_t_expander"
## [1] 34255
## [1] "removing: surg_bun_t_explora"
## [1] 34255
## [1] "removing: surg_bun_t_explore"
## [1] 34255
## [1] "removing: surg_bun_t_extra"
## [1] 34255
## [1] "removing: surg_bun_t_extremity"
## [1] 34255
## [1] "removing: surg_bun_t_eye"
## [1] 34255
## [1] "removing: surg_bun_t_eyelashes"
## [1] 34255
## [1] "removing: surg_bun_t_eyelid"
## [1] 34255
## [1] "removing: surg_bun_t_face"
## [1] 34255
## [1] "removing: surg_bun_t_fallopian"
## [1] 34255
## [1] "removing: surg_bun_t_fascia"
## [1] 34255
## [1] "removing: surg_bun_t_fat"
## [1] 34255
## [1] "removing: surg_bun_t_femur"
## [1] 34255
## [1] "removing: surg_bun_t_fibula"
## [1] 34255
## [1] "removing: surg_bun_t_finger"
## [1] 34255
## [1] "removing: surg_bun_t_fissure"
## [1] 34255
## [1] "removing: surg_bun_t_fistula"
## [1] 34255
## [1] "removing: surg_bun_t_fixation"
## [1] 34255
## [1] "removing: surg_bun_t_flap"
## [1] 34255
## [1] "removing: surg_bun_t_flow"
## [1] 34255
## [1] "removing: surg_bun_t_fluid"
## [1] 34255
## [1] "removing: surg_bun_t_fna"
## [1] 34255
## [1] "removing: surg_bun_t_follow-up"
## [1] 34255
## [1] "removing: surg_bun_t_foot"
## [1] 34255
## [1] "removing: surg_bun_t_forearm"
## [1] 34255
## [1] "removing: surg_bun_t_forehead"
## [1] 34255
## [1] "removing: surg_bun_t_foreign"
## [1] 34255
## [1] "removing: surg_bun_t_fracture"
## [1] 34255
## [1] "removing: surg_bun_t_free"
## [1] 34255
## [1] "removing: surg_bun_t_frontal"
## [1] 34255
## [1] "removing: surg_bun_t_fundoplasty"
## [1] 34255
## [1] "removing: surg_bun_t_fusion"
## [1] 34255
## [1] "removing: surg_bun_t_gallbladder"
## [1] 34255
## [1] "removing: surg_bun_t_gastrointestinal"
## [1] 34255
## [1] "removing: surg_bun_t_gastrostomy"
## [1] 34255
## [1] "removing: surg_bun_t_gentialia"
## [1] 34255
## [1] "removing: surg_bun_t_gi"
## [1] 34255
## [1] "removing: surg_bun_t_gland"
## [1] 34255
## [1] "removing: surg_bun_t_glaucoma"
## [1] 34255
## [1] "removing: surg_bun_t_graft"
## [1] 34255
## [1] "removing: surg_bun_t_growth"
## [1] 34255
## [1] "removing: surg_bun_t_gum"
## [1] 34255
## [1] "removing: surg_bun_t_hair"
## [1] 34255
## [1] "removing: surg_bun_t_hallux"
## [1] 34255
## [1] "removing: surg_bun_t_halo"
## [1] 34255
## [1] "removing: surg_bun_t_hand"
## [1] 34255
## [1] "removing: surg_bun_t_harvest"
## [1] 34255
## [1] "removing: surg_bun_t_head"
## [1] 34255
## [1] "removing: surg_bun_t_heart"
## [1] 34255
## [1] "removing: surg_bun_t_hematoma"
## [1] 34255
## [1] "removing: surg_bun_t_hemorrhoid"
## [1] 34255
## [1] "removing: surg_bun_t_hemorrhoidopexy"
## [1] 34255
## [1] "removing: surg_bun_t_hernia"
## [1] 34255
## [1] "removing: surg_bun_t_hip"
## [1] 34255
## [1] "removing: surg_bun_t_hole"
## [1] 34255
## [1] "removing: surg_bun_t_hormone"
## [1] 34255
## [1] "removing: surg_bun_t_humer"
## [1] 34255
## [1] "removing: surg_bun_t_hydrocele"
## [1] 34255
## [1] "removing: surg_bun_t_hyst"
## [1] 34255
## [1] "removing: surg_bun_t_hysterectomy"
## [1] 34255
## [1] "removing: surg_bun_t_iliopsoas"
## [1] 34255
## [1] "removing: surg_bun_t_illiac"
## [1] 34255
## [1] "removing: surg_bun_t_imag"
## [1] 34255
## [1] "removing: surg_bun_t_implant"
## [1] 34255
## [1] "removing: surg_bun_t_incision"
## [1] 34255
## [1] "removing: surg_bun_t_infect"
## [1] 34255
## [1] "removing: surg_bun_t_injection"
## [1] 34254
## [1] "removing: surg_bun_t_inner"
## [1] 34254
## [1] "removing: surg_bun_t_insert"
## [1] 34254
## [1] "removing: surg_bun_t_interlaminar"
## [1] 34254
## [1] "removing: surg_bun_t_internal"
## [1] 34254
## [1] "removing: surg_bun_t_intestine"
## [1] 34254
## [1] "removing: surg_bun_t_intracorporeal"
## [1] 34254
## [1] "removing: surg_bun_t_intracranial"
## [1] 34254
## [1] "removing: surg_bun_t_intranasal"
## [1] 34254
## [1] "removing: surg_bun_t_inverted"
## [1] 34254
## [1] "removing: surg_bun_t_iris"
## [1] 34254
## [1] "removing: surg_bun_t_jaw"
## [1] 34254
## [1] "removing: surg_bun_t_joint"
## [1] 34254
## [1] "removing: surg_bun_t_kidney"
## [1] 34254
## [1] "removing: surg_bun_t_knee"
## [1] 34254
## [1] "removing: surg_bun_t_kneecap"
## [1] 34254
## [1] "removing: surg_bun_t_knuckle"
## [1] 34254
## [1] "removing: surg_bun_t_kyphectomy"
## [1] 34254
## [1] "removing: surg_bun_t_laceration"
## [1] 34254
## [1] "removing: surg_bun_t_lamina"
## [1] 34254
## [1] "removing: surg_bun_t_laminotomy"
## [1] 34254
## [1] "removing: surg_bun_t_lap"
## [1] 34254
## [1] "removing: surg_bun_t_laparoscop"
## [1] 34254
## [1] "removing: surg_bun_t_large"
## [1] 34254
## [1] "removing: surg_bun_t_laryngoscopy"
## [1] 34254
## [1] "removing: surg_bun_t_larynx"
## [1] 34254
## [1] "removing: surg_bun_t_laser"
## [1] 34254
## [1] "removing: surg_bun_t_lavage"
## [1] 34254
## [1] "removing: surg_bun_t_lefort"
## [1] 34254
## [1] "removing: surg_bun_t_leg"
## [1] 34254
## [1] "removing: surg_bun_t_lengthen"
## [1] 34254
## [1] "removing: surg_bun_t_lesion"
## [1] 34254
## [1] "removing: surg_bun_t_ligament"
## [1] 34254
## [1] "removing: surg_bun_t_ligation"
## [1] 34254
## [1] "removing: surg_bun_t_lining"
## [1] 34254
## [1] "removing: surg_bun_t_lip"
## [1] 34254
## [1] "removing: surg_bun_t_lipectomy"
## [1] 34254
## [1] "removing: surg_bun_t_liver"
## [1] 34254
## [1] "removing: surg_bun_t_lobectomy"
## [1] 34254
## [1] "removing: surg_bun_t_long"
## [1] 34254
## [1] "removing: surg_bun_t_lower"
## [1] 34254
## [1] "removing: surg_bun_t_lumbar"
## [1] 34254
## [1] "removing: surg_bun_t_lung"
## [1] 34254
## [1] "removing: surg_bun_t_lymph"
## [1] 34254
## [1] "removing: surg_bun_t_major"
## [1] 34254
## [1] "removing: surg_bun_t_mandible"
## [1] 34254
## [1] "removing: surg_bun_t_marrow"
## [1] 34254
## [1] "removing: surg_bun_t_mast"
## [1] 34254
## [1] "removing: surg_bun_t_mastectomy"
## [1] 34254
## [1] "removing: surg_bun_t_mastoid"
## [1] 34254
## [1] "removing: surg_bun_t_material"
## [1] 34254
## [1] "removing: surg_bun_t_maxilla"
## [1] 34254
## [1] "removing: surg_bun_t_mcp"
## [1] 34254
## [1] "removing: surg_bun_t_mesh"
## [1] 34254
## [1] "removing: surg_bun_t_metacarpal"
## [1] 34254
## [1] "removing: surg_bun_t_microvasc"
## [1] 34254
## [1] "removing: surg_bun_t_middle"
## [1] 34254
## [1] "removing: surg_bun_t_midfoot"
## [1] 34254
## [1] "removing: surg_bun_t_mitral"
## [1] 34254
## [1] "removing: surg_bun_t_mohs"
## [1] 34254
## [1] "removing: surg_bun_t_mouth"
## [1] 34254
## [1] "removing: surg_bun_t_musc"
## [1] 34253
## [1] "removing: surg_bun_t_myelography"
## [1] 34253
## [1] "removing: surg_bun_t_myomectomy"
## [1] 34253
## [1] "removing: surg_bun_t_nail"
## [1] 34253
## [1] "removing: surg_bun_t_nasal"
## [1] 34253
## [1] "removing: surg_bun_t_nasolacrimal"
## [1] 34253
## [1] "removing: surg_bun_t_neck"
## [1] 34253
## [1] "removing: surg_bun_t_needle"
## [1] 34253
## [1] "removing: surg_bun_t_nephrectomy"
## [1] 34253
## [1] "removing: surg_bun_t_nerve"
## [1] 34253
## [1] "removing: surg_bun_t_neuroelectrode"
## [1] 34253
## [1] "removing: surg_bun_t_neuroendoscopy"
## [1] 34253
## [1] "removing: surg_bun_t_nipple"
## [1] 34253
## [1] "removing: surg_bun_t_njx"
## [1] 34253
## [1] "removing: surg_bun_t_node"
## [1] 34253
## [1] "removing: surg_bun_t_nodule"
## [1] 34253
## [1] "removing: surg_bun_t_nonfacial"
## [1] 34253
## [1] "removing: surg_bun_t_nose"
## [1] 34253
## [1] "removing: surg_bun_t_nosebleed"
## [1] 34253
## [1] "removing: surg_bun_t_nuero"
## [1] 34253
## [1] "removing: surg_bun_t_obstetrical"
## [1] 34253
## [1] "removing: surg_bun_t_open"
## [1] 34253
## [1] "removing: surg_bun_t_opening"
## [1] 34253
## [1] "removing: surg_bun_t_oral"
## [1] 34253
## [1] "removing: surg_bun_t_orbitocranial"
## [1] 34253
## [1] "removing: surg_bun_t_orchiopexy"
## [1] 34253
## [1] "removing: surg_bun_t_other"
## [1] 34253
## [1] "removing: surg_bun_t_ovarian"
## [1] 34253
## [1] "removing: surg_bun_t_ovary"
## [1] 34253
## [1] "removing: surg_bun_t_pacemaker"
## [1] 34253
## [1] "removing: surg_bun_t_pad"
## [1] 34253
## [1] "removing: surg_bun_t_palate"
## [1] 34253
## [1] "removing: surg_bun_t_palm"
## [1] 34253
## [1] "removing: surg_bun_t_palsy"
## [1] 34253
## [1] "removing: surg_bun_t_pancreas"
## [1] 34253
## [1] "removing: surg_bun_t_pancreatectomy"
## [1] 34253
## [1] "removing: surg_bun_t_paravert"
## [1] 34253
## [1] "removing: surg_bun_t_part"
## [1] 34253
## [1] "removing: surg_bun_t_partial"
## [1] 34253
## [1] "removing: surg_bun_t_pedcle"
## [1] 34253
## [1] "removing: surg_bun_t_pedicle"
## [1] 34253
## [1] "removing: surg_bun_t_peel"
## [1] 34253
## [1] "removing: surg_bun_t_pelvic"
## [1] 34253
## [1] "removing: surg_bun_t_pelvis"
## [1] 34253
## [1] "removing: surg_bun_t_penile"
## [1] 34253
## [1] "removing: surg_bun_t_penis"
## [1] 34253
## [1] "removing: surg_bun_t_percut"
## [1] 34253
## [1] "removing: surg_bun_t_percutaneous"
## [1] 34253
## [1] "removing: surg_bun_t_perineum"
## [1] 34253
## [1] "removing: surg_bun_t_pernineum"
## [1] 34253
## [1] "removing: surg_bun_t_pharynx"
## [1] 34253
## [1] "removing: surg_bun_t_pierce"
## [1] 34253
## [1] "removing: surg_bun_t_pilonidal"
## [1] 34253
## [1] "removing: surg_bun_t_place"
## [1] 34253
## [1] "removing: surg_bun_t_plastic"
## [1] 34253
## [1] "removing: surg_bun_t_plate"
## [1] 34253
## [1] "removing: surg_bun_t_pneumonectomy"
## [1] 34253
## [1] "removing: surg_bun_t_polyp"
## [1] 34253
## [1] "removing: surg_bun_t_pouch"
## [1] 34253
## [1] "removing: surg_bun_t_pregnancy"
## [1] 34253
## [1] "removing: surg_bun_t_prep"
## [1] 34253
## [1] "removing: surg_bun_t_prepare"
## [1] 34253
## [1] "removing: surg_bun_t_pressure"
## [1] 34253
## [1] "removing: surg_bun_t_process"
## [1] 34253
## [1] "removing: surg_bun_t_proctosigmoidoscopy"
## [1] 34253
## [1] "removing: surg_bun_t_prolapse"
## [1] 34253
## [1] "removing: surg_bun_t_prostat"
## [1] 34253
## [1] "removing: surg_bun_t_prostate"
## [1] 34253
## [1] "removing: surg_bun_t_prosthesis"
## [1] 34253
## [1] "removing: surg_bun_t_prosthetic"
## [1] 34253
## [1] "removing: surg_bun_t_pulmonary"
## [1] 34253
## [1] "removing: surg_bun_t_punch"
## [1] 34253
## [1] "removing: surg_bun_t_puncture"
## [1] 34253
## [1] "removing: surg_bun_t_pyeloplasty"
## [1] 34253
## [1] "removing: surg_bun_t_radius"
## [1] 34253
## [1] "removing: surg_bun_t_realignment"
## [1] 34253
## [1] "removing: surg_bun_t_rechannel"
## [1] 34253
## [1] "removing: surg_bun_t_reconstruct"
## [1] 34253
## [1] "removing: surg_bun_t_reconstruction"
## [1] 34253
## [1] "removing: surg_bun_t_rectal"
## [1] 34253
## [1] "removing: surg_bun_t_rectum"
## [1] 34253
## [1] "removing: surg_bun_t_reduction"
## [1] 34253
## [1] "removing: surg_bun_t_reinforce"
## [1] 34253
## [1] "removing: surg_bun_t_release"
## [1] 34253
## [1] "removing: surg_bun_t_remov"
## [1] 34253
## [1] "removing: surg_bun_t_removal"
## [1] 34253
## [1] "removing: surg_bun_t_remove"
## [1] 34253
## [1] "removing: surg_bun_t_renal"
## [1] 34253
## [1] "removing: surg_bun_t_repair"
## [1] 34253
## [1] "removing: surg_bun_t_replantation"
## [1] 34253
## [1] "removing: surg_bun_t_resect"
## [1] 34253
## [1] "removing: surg_bun_t_reset"
## [1] 34253
## [1] "removing: surg_bun_t_retina"
## [1] 34253
## [1] "removing: surg_bun_t_retinopathy"
## [1] 34253
## [1] "removing: surg_bun_t_revise"
## [1] 34253
## [1] "removing: surg_bun_t_revision"
## [1] 34253
## [1] "removing: surg_bun_t_rib"
## [1] 34253
## [1] "removing: surg_bun_t_ridge"
## [1] 34253
## [1] "removing: surg_bun_t_ring"
## [1] 34253
## [1] "removing: surg_bun_t_roof"
## [1] 34253
## [1] "removing: surg_bun_t_rotator"
## [1] 34253
## [1] "removing: surg_bun_t_routine"
## [1] 34253
## [1] "removing: surg_bun_t_sacrum"
## [1] 34253
## [1] "removing: surg_bun_t_salivary"
## [1] 34253
## [1] "removing: surg_bun_t_scapula"
## [1] 34253
## [1] "removing: surg_bun_t_scar"
## [1] 34253
## [1] "removing: surg_bun_t_scope"
## [1] 34253
## [1] "removing: surg_bun_t_scrotum"
## [1] 34253
## [1] "removing: surg_bun_t_segment"
## [1] 34253
## [1] "removing: surg_bun_t_segmental"
## [1] 34253
## [1] "removing: surg_bun_t_segmentectomy"
## [1] 34253
## [1] "removing: surg_bun_t_septum"
## [1] 34253
## [1] "removing: surg_bun_t_sex"
## [1] 34253
## [1] "removing: surg_bun_t_shave"
## [1] 34253
## [1] "removing: surg_bun_t_sheath"
## [1] 34253
## [1] "removing: surg_bun_t_shorten"
## [1] 34253
## [1] "removing: surg_bun_t_shoulder"
## [1] 34253
## [1] "removing: surg_bun_t_shunt"
## [1] 34253
## [1] "removing: surg_bun_t_sigmoidoscopy"
## [1] 34253
## [1] "removing: surg_bun_t_simple"
## [1] 34253
## [1] "removing: surg_bun_t_sinus"
## [1] 34253
## [1] "removing: surg_bun_t_skeletal"
## [1] 34253
## [1] "removing: surg_bun_t_skin"
## [1] 34253
## [1] "removing: surg_bun_t_skull"
## [1] 34253
## [1] "removing: surg_bun_t_sleeve"
## [1] 34253
## [1] "removing: surg_bun_t_slipped"
## [1] 34253
## [1] "removing: surg_bun_t_small"
## [1] 34253
## [1] "removing: surg_bun_t_socket"
## [1] 34253
## [1] "removing: surg_bun_t_soft"
## [1] 34253
## [1] "removing: surg_bun_t_sore"
## [1] 34253
## [1] "removing: surg_bun_t_sperm"
## [1] 34253
## [1] "removing: surg_bun_t_sphenoid"
## [1] 34253
## [1] "removing: surg_bun_t_sphincter"
## [1] 34253
## [1] "removing: surg_bun_t_spine"
## [1] 34253
## [1] "removing: surg_bun_t_spleen"
## [1] 34253
## [1] "removing: surg_bun_t_splint"
## [1] 34253
## [1] "removing: surg_bun_t_split"
## [1] 34253
## [1] "removing: surg_bun_t_stem"
## [1] 34253
## [1] "removing: surg_bun_t_stenosis"
## [1] 34253
## [1] "removing: surg_bun_t_stent"
## [1] 34253
## [1] "removing: surg_bun_t_sternal"
## [1] 34253
## [1] "removing: surg_bun_t_sternum"
## [1] 34253
## [1] "removing: surg_bun_t_stimulation"
## [1] 34253
## [1] "removing: surg_bun_t_stomach"
## [1] 34253
## [1] "removing: surg_bun_t_stone"
## [1] 34253
## [1] "removing: surg_bun_t_stop"
## [1] 34253
## [1] "removing: surg_bun_t_strapping"
## [1] 34253
## [1] "removing: surg_bun_t_stricture"
## [1] 34253
## [1] "removing: surg_bun_t_styloid"
## [1] 34253
## [1] "removing: surg_bun_t_sub"
## [1] 34253
## [1] "removing: surg_bun_t_suction"
## [1] 34253
## [1] "removing: surg_bun_t_superficial"
## [1] 34253
## [1] "removing: surg_bun_t_suspension"
## [1] 34253
## [1] "removing: surg_bun_t_sutures"
## [1] 34253
## [1] "removing: surg_bun_t_sweat"
## [1] 34253
## [1] "removing: surg_bun_t_tags"
## [1] 34253
## [1] "removing: surg_bun_t_tail"
## [1] 34253
## [1] "removing: surg_bun_t_tear"
## [1] 34253
## [1] "removing: surg_bun_t_temporal"
## [1] 34253
## [1] "removing: surg_bun_t_tendon"
## [1] 34253
## [1] "removing: surg_bun_t_testis"
## [1] 34253
## [1] "removing: surg_bun_t_therapy"
## [1] 34253
## [1] "removing: surg_bun_t_thigh"
## [1] 34253
## [1] "removing: surg_bun_t_thoracoscopy"
## [1] 34253
## [1] "removing: surg_bun_t_thoracotomy"
## [1] 34253
## [1] "removing: surg_bun_t_thorax"
## [1] 34253
## [1] "removing: surg_bun_t_throat"
## [1] 34253
## [1] "removing: surg_bun_t_thrombin"
## [1] 34253
## [1] "removing: surg_bun_t_thrombo"
## [1] 34253
## [1] "removing: surg_bun_t_thumb"
## [1] 34253
## [1] "removing: surg_bun_t_thymus"
## [1] 34253
## [1] "removing: surg_bun_t_thyroid"
## [1] 34253
## [1] "removing: surg_bun_t_tip"
## [1] 34253
## [1] "removing: surg_bun_t_tiss"
## [1] 34253
## [1] "removing: surg_bun_t_tkr"
## [1] 34253
## [1] "removing: surg_bun_t_toe"
## [1] 34253
## [1] "removing: surg_bun_t_tongue"
## [1] 34253
## [1] "removing: surg_bun_t_tonsil"
## [1] 34253
## [1] "removing: surg_bun_t_torsion"
## [1] 34253
## [1] "removing: surg_bun_t_total"
## [1] 34253
## [1] "removing: surg_bun_t_tract"
## [1] 34253
## [1] "removing: surg_bun_t_transcath"
## [1] 34253
## [1] "removing: surg_bun_t_transfer"
## [1] 34253
## [1] "removing: surg_bun_t_transform"
## [1] 34253
## [1] "removing: surg_bun_t_transplant"
## [1] 34253
## [1] "removing: surg_bun_t_treat"
## [1] 34253
## [1] "removing: surg_bun_t_triceps"
## [1] 34253
## [1] "removing: surg_bun_t_tricuspid"
## [1] 34253
## [1] "removing: surg_bun_t_trigger"
## [1] 34253
## [1] "removing: surg_bun_t_trim"
## [1] 34253
## [1] "removing: surg_bun_t_tube"
## [1] 34253
## [1] "removing: surg_bun_t_tum"
## [1] 34253
## [1] "removing: surg_bun_t_tumor"
## [1] 34253
## [1] "removing: surg_bun_t_tunneled"
## [1] 34253
## [1] "removing: surg_bun_t_ulcer"
## [1] 34253
## [1] "removing: surg_bun_t_ulna"
## [1] 34253
## [1] "removing: surg_bun_t_umbil"
## [1] 34253
## [1] "removing: surg_bun_t_unit"
## [1] 34253
## [1] "removing: surg_bun_t_unstable"
## [1] 34253
## [1] "removing: surg_bun_t_upper"
## [1] 34253
## [1] "removing: surg_bun_t_ureter"
## [1] 34253
## [1] "removing: surg_bun_t_urethra"
## [1] 34253
## [1] "removing: surg_bun_t_urethrovaginal"
## [1] 34253
## [1] "removing: surg_bun_t_uterus"
## [1] 34253
## [1] "removing: surg_bun_t_uvula"
## [1] 34253
## [1] "removing: surg_bun_t_vag"
## [1] 34253
## [1] "removing: surg_bun_t_vagotomy"
## [1] 34253
## [1] "removing: surg_bun_t_valgus"
## [1] 34253
## [1] "removing: surg_bun_t_valve"
## [1] 34253
## [1] "removing: surg_bun_t_vascular"
## [1] 34253
## [1] "removing: surg_bun_t_vein"
## [1] 34253
## [1] "removing: surg_bun_t_venipuncture"
## [1] 34253
## [1] "removing: surg_bun_t_ventricular"
## [1] 34253
## [1] "removing: surg_bun_t_vertebra"
## [1] 34253
## [1] "removing: surg_bun_t_vessel"
## [1] 34253
## [1] "removing: surg_bun_t_vulva"
## [1] 34253
## [1] "removing: surg_bun_t_wall"
## [1] 34253
## [1] "removing: surg_bun_t_wax"
## [1] 34253
## [1] "removing: surg_bun_t_web"
## [1] 34253
## [1] "removing: surg_bun_t_wedge"
## [1] 34253
## [1] "removing: surg_bun_t_window"
## [1] 34253
## [1] "removing: surg_bun_t_windpipe"
## [1] 34253
## [1] "removing: surg_bun_t_wound"
## [1] 34253
## [1] "removing: surg_bun_t_wrinkle"
## [1] 34253
## [1] "removing: surg_bun_t_wrist"
## [1] 34253
## [1] "removing: surg_bun_t_x-ray"
## [1] 34253
## [1] "removing: surg_bun_t_xenograft"
## [1] 34253
## [1] "removing: surg_bun_t_xenogrft"
## [1] 34253
## [1] "removing: surg_bun_t_xgraft"
## [1] 34253
## [1] "removing: medi_bun_t_applicat"
## [1] 34253
## [1] "removing: medi_bun_t_apply"
## [1] 34253
## [1] "removing: medi_bun_t_asses"
## [1] 34168
## [1] "removing: medi_bun_t_benzathine"
## [1] 34161
## [1] "removing: medi_bun_t_betamethasone"
## [1] 34159
## [1] "removing: medi_bun_t_boot"
## [1] 34147
## [1] "removing: medi_bun_t_candida"
## [1] 34147
## [1] "removing: medi_bun_t_cast"
## [1] 34147
## [1] "removing: medi_bun_t_cathet"
## [1] 34147
## [1] "removing: medi_bun_t_catheterisation"
## [1] 34147
## [1] "removing: medi_bun_t_chiropract"
## [1] 34081
## [1] "removing: medi_bun_t_choline"
## [1] 34081
## [1] "removing: medi_bun_t_clostridium"
## [1] 34081
## [1] "removing: medi_bun_t_cochlear"
## [1] 34081
## [1] "removing: medi_bun_t_cold"
## [1] 34075
## [1] "removing: medi_bun_t_consult"
## [1] 34075
## [1] "removing: medi_bun_t_consultation"
## [1] 34075
## [1] "removing: medi_bun_t_crutch"
## [1] 34075
## [1] "removing: medi_bun_t_cyst"
## [1] 34075
## [1] "removing: medi_bun_t_derma"
## [1] 34075
## [1] "removing: medi_bun_t_dermal"
## [1] 34075
## [1] "removing: medi_bun_t_dialys"
## [1] 34048
## [1] "removing: medi_bun_t_discharg"
## [1] 34048
## [1] "removing: medi_bun_t_drain"
## [1] 34048
## [1] "removing: medi_bun_t_draw"
## [1] 33971
## [1] "removing: medi_bun_t_dtap"
## [1] 33959
## [1] "removing: medi_bun_t_dtap-ipv"
## [1] 33959
## [1] "removing: medi_bun_t_endocrine"
## [1] 33959
## [1] "removing: medi_bun_t_epstein-barr"
## [1] 33959
## [1] "removing: medi_bun_t_equipment"
## [1] 33959
## [1] "removing: medi_bun_t_establish"
## [1] 33929
## [1] "removing: medi_bun_t_estradiol"
## [1] 33928
## [1] "removing: medi_bun_t_evaluat"
## [1] 33865
## [1] "removing: medi_bun_t_follow-up"
## [1] 33865
## [1] "removing: medi_bun_t_gastroenterology"
## [1] 33865
## [1] "removing: medi_bun_t_gonorrhoeae"
## [1] 33865
## [1] "removing: medi_bun_t_haloperidol"
## [1] 33858
## [1] "removing: medi_bun_t_hearing"
## [1] 33794
## [1] "removing: medi_bun_t_hemorrhoid"
## [1] 33794
## [1] "removing: medi_bun_t_hepatitis"
## [1] 33787
## [1] "removing: medi_bun_t_hiv - 2"
## [1] 33787
## [1] "removing: medi_bun_t_hiv - 3"
## [1] 33787
## [1] "removing: medi_bun_t_home"
## [1] 33786
## [1] "removing: medi_bun_t_hormone"
## [1] 33786
## [1] "removing: medi_bun_t_hpv"
## [1] 33700
## [1] "removing: medi_bun_t_immune"
## [1] 33698
## [1] "removing: medi_bun_t_immunoassay"
## [1] 33657
## [1] "removing: medi_bun_t_immunohisto"
## [1] 33657
## [1] "removing: medi_bun_t_immunotherapy"
## [1] 33576
## [1] "removing: medi_bun_t_insulin"
## [1] 33553
## [1] "removing: medi_bun_t_intrauterine"
## [1] 33553
## [1] "removing: medi_bun_t_iodine"
## [1] 33553
## [1] "removing: medi_bun_t_laiv4"
## [1] 33552
## [1] "removing: medi_bun_t_medroxyprogesterone"
## [1] 33533
## [1] "removing: medi_bun_t_mental"
## [1] 33529
## [1] "removing: medi_bun_t_metabolic"
## [1] 33529
## [1] "removing: medi_bun_t_methylprednisolone"
## [1] 33497
## [1] "removing: medi_bun_t_mgmt"
## [1] 33475
## [1] "removing: medi_bun_t_mmr"
## [1] 33416
## [1] "removing: medi_bun_t_muscle"
## [1] 33416
## [1] "removing: medi_bun_t_needle"
## [1] 33409
## [1] "removing: medi_bun_t_non-invasive"
## [1] 33409
## [1] "removing: medi_bun_t_non-stress"
## [1] 33409
## [1] "removing: medi_bun_t_nursing"
## [1] 33409
## [1] "removing: medi_bun_t_nutrition"
## [1] 33375
## [1] "removing: medi_bun_t_observ"
## [1] 33367
## [1] "removing: medi_bun_t_oral"
## [1] 33323
## [1] "removing: medi_bun_t_outpatient"
## [1] 33323
## [1] "removing: medi_bun_t_pap"
## [1] 33314
## [1] "removing: medi_bun_t_physical"
## [1] 33312
## [1] "removing: medi_bun_t_portable"
## [1] 33312
## [1] "removing: medi_bun_t_prednisolone"
## [1] 33312
## [1] "removing: medi_bun_t_pressure"
## [1] 33311
## [1] "removing: medi_bun_t_progesterone"
## [1] 33311
## [1] "removing: medi_bun_t_psych"
## [1] 33296
## [1] "removing: medi_bun_t_psychiat"
## [1] 33296
## [1] "removing: medi_bun_t_psycho"
## [1] 33296
## [1] "removing: medi_bun_t_pulmonary"
## [1] 33280
## [1] "removing: medi_bun_t_rabies"
## [1] 33279
## [1] "removing: medi_bun_t_radiation"
## [1] 33278
## [1] "removing: medi_bun_t_radio"
## [1] 33278
## [1] "removing: medi_bun_t_reeducation"
## [1] 33222
## [1] "removing: medi_bun_t_reeval"
## [1] 33222
## [1] "removing: medi_bun_t_remove"
## [1] 33222
## [1] "removing: medi_bun_t_rheumatoid"
## [1] 33222
## [1] "removing: medi_bun_t_riv4"
## [1] 33127
## [1] "removing: medi_bun_t_rubella"
## [1] 33127
## [1] "removing: medi_bun_t_sciatic"
## [1] 33127
## [1] "removing: medi_bun_t_serum"
## [1] 33127
## [1] "removing: medi_bun_t_services"
## [1] 33125
## [1] "removing: medi_bun_t_shave"
## [1] 33125
## [1] "removing: medi_bun_t_simple"
## [1] 33124
## [1] "removing: medi_bun_t_sleep"
## [1] 33113
## [1] "removing: medi_bun_t_smear"
## [1] 33113
## [1] "removing: medi_bun_t_soft"
## [1] 33113
## [1] "removing: medi_bun_t_sound"
## [1] 33097
## [1] "removing: medi_bun_t_speech"
## [1] 33096
## [1] "removing: medi_bun_t_splint"
## [1] 33096
## [1] "removing: medi_bun_t_stent"
## [1] 33096
## [1] "removing: medi_bun_t_sterile"
## [1] 33096
## [1] "removing: medi_bun_t_stimulation"
## [1] 33053
## [1] "removing: medi_bun_t_stocking"
## [1] 33053
## [1] "removing: medi_bun_t_stone"
## [1] 33053
## [1] "removing: medi_bun_t_stool"
## [1] 33053
## [1] "removing: medi_bun_t_strep"
## [1] 33053
## [1] "removing: medi_bun_t_stress"
## [1] 33026
## [1] "removing: medi_bun_t_subseq"
## [1] 33013
## [1] "removing: medi_bun_t_substitute"
## [1] 33013
## [1] "removing: medi_bun_t_supplies"
## [1] 33013
## [1] "removing: medi_bun_t_suppls"
## [1] 33013
## [1] "removing: medi_bun_t_support"
## [1] 33013
## [1] "removing: medi_bun_t_susceptible"
## [1] 33013
## [1] "removing: medi_bun_t_swallow"
## [1] 33013
## [1] "removing: medi_bun_t_synth"
## [1] 33013
## [1] "removing: medi_bun_t_syphilis"
## [1] 33013
## [1] "removing: medi_bun_t_syringe"
## [1] 33013
## [1] "removing: medi_bun_t_system"
## [1] 33013
## [1] "removing: medi_bun_t_testing"
## [1] 33013
## [1] "removing: medi_bun_t_testosterone"
## [1] 32972
## [1] "removing: medi_bun_t_tests"
## [1] 32963
## [1] "removing: medi_bun_t_tetanus"
## [1] 32963
## [1] "removing: medi_bun_t_thyroxine"
## [1] 32963
## [1] "removing: medi_bun_t_tobacco"
## [1] 32963
## [1] "removing: medi_bun_t_toe"
## [1] 32963
## [1] "removing: medi_bun_t_tolerance"
## [1] 32963
## [1] "removing: medi_bun_t_topical"
## [1] 32963
## [1] "removing: medi_bun_t_toxin"
## [1] 32962
## [1] "removing: medi_bun_t_trach"
## [1] 32962
## [1] "removing: medi_bun_t_traction"
## [1] 32954
## [1] "removing: medi_bun_t_training"
## [1] 32944
## [1] "removing: medi_bun_t_transfusion"
## [1] 32935
## [1] "removing: medi_bun_t_treat"
## [1] 32915
## [1] "removing: medi_bun_t_typhoid"
## [1] 32914
## [1] "removing: medi_bun_t_ultrasound"
## [1] 32914
## [1] "removing: medi_bun_t_urinalysis"
## [1] 32914
## [1] "removing: medi_bun_t_urinary"
## [1] 32914
## [1] "removing: medi_bun_t_urine"
## [1] 32914
## [1] "removing: medi_bun_t_uterine"
## [1] 32914
## [1] "removing: medi_bun_t_valve"
## [1] 32914
## [1] "removing: medi_bun_t_vein"
## [1] 32914
## [1] "removing: medi_bun_t_visual"
## [1] 32898
## [1] "removing: medi_bun_t_vitamin"
## [1] 32842
## [1] "removing: medi_bun_t_wbc"
## [1] 32842
## [1] "removing: medi_bun_t_wheelchair"
## [1] 32842
## [1] "removing: medi_bun_t_wheezing"
## [1] 32842
## [1] "removing: medi_bun_t_whirlpool"
## [1] 32841
## [1] "removing: medi_bun_t_wire"
## [1] 32841
## [1] "removing: medi_bun_t_with"
## [1] 32799
## [1] "removing: medi_bun_t_without"
## [1] 32799
## [1] "removing: medi_bun_t_wkend"
## [1] 32799
## [1] "removing: medi_bun_t_wound"
## [1] 32791
## [1] "removing: medi_bun_t_menb"
## [1] 32783
## [1] "removing: medi_bun_t_xgraft"
## [1] 32783
## [1] "removing: medi_bun_t_pcv13"
## [1] 32719
## [1] "removing: medi_bun_t_var_vaccine_live_subq"
## [1] 32713
## [1] "removing: medi_bun_t_menacwy"
## [1] 32680
## [1] "removing: radi_bun_t_abdom"
## [1] 32677
## [1] "removing: radi_bun_t_adrenal"
## [1] 32677
## [1] "removing: radi_bun_t_angiography"
## [1] 32677
## [1] "removing: radi_bun_t_ankle"
## [1] 32677
## [1] "removing: radi_bun_t_aorta"
## [1] 32677
## [1] "removing: radi_bun_t_area"
## [1] 32677
## [1] "removing: radi_bun_t_arm"
## [1] 32677
## [1] "removing: radi_bun_t_artery"
## [1] 32677
## [1] "removing: radi_bun_t_bile"
## [1] 32677
## [1] "removing: radi_bun_t_bladder"
## [1] 32677
## [1] "removing: radi_bun_t_blade"
## [1] 32677
## [1] "removing: radi_bun_t_blockage"
## [1] 32677
## [1] "removing: radi_bun_t_body"
## [1] 32677
## [1] "removing: radi_bun_t_bone"
## [1] 32677
## [1] "removing: radi_bun_t_brain"
## [1] 32676
## [1] "removing: radi_bun_t_breast"
## [1] 32676
## [1] "removing: radi_bun_t_cardiac"
## [1] 32676
## [1] "removing: radi_bun_t_cast"
## [1] 32676
## [1] "removing: radi_bun_t_chest"
## [1] 32676
## [1] "removing: radi_bun_t_complex"
## [1] 32676
## [1] "removing: radi_bun_t_contrast"
## [1] 32676
## [1] "removing: radi_bun_t_cranial"
## [1] 32676
## [1] "removing: radi_bun_t_ct"
## [1] 32673
## [1] "removing: radi_bun_t_ct_scan"
## [1] 32673
## [1] "removing: radi_bun_t_delivery"
## [1] 32673
## [1] "removing: radi_bun_t_diagnostic"
## [1] 32673
## [1] "removing: radi_bun_t_dislocation"
## [1] 32673
## [1] "removing: radi_bun_t_duct"
## [1] 32673
## [1] "removing: radi_bun_t_echo"
## [1] 32673
## [1] "removing: radi_bun_t_elbow"
## [1] 32673
## [1] "removing: radi_bun_t_epidurograpy"
## [1] 32673
## [1] "removing: radi_bun_t_evaluation"
## [1] 32673
## [1] "removing: radi_bun_t_exam"
## [1] 32672
## [1] "removing: radi_bun_t_eye"
## [1] 32672
## [1] "removing: radi_bun_t_femur"
## [1] 32672
## [1] "removing: radi_bun_t_fetal"
## [1] 32672
## [1] "removing: radi_bun_t_finger"
## [1] 32672
## [1] "removing: radi_bun_t_fistula"
## [1] 32672
## [1] "removing: radi_bun_t_fluoroscopy"
## [1] 32672
## [1] "removing: radi_bun_t_fmri"
## [1] 32672
## [1] "removing: radi_bun_t_foot"
## [1] 32672
## [1] "removing: radi_bun_t_forearm"
## [1] 32672
## [1] "removing: radi_bun_t_fracture"
## [1] 32672
## [1] "removing: radi_bun_t_gadobutrol"
## [1] 32672
## [1] "removing: radi_bun_t_gastric"
## [1] 32672
## [1] "removing: radi_bun_t_general"
## [1] 32672
## [1] "removing: radi_bun_t_gential"
## [1] 32672
## [1] "removing: radi_bun_t_gestation"
## [1] 32672
## [1] "removing: radi_bun_t_gi_tract"
## [1] 32672
## [1] "removing: radi_bun_t_glands"
## [1] 32672
## [1] "removing: radi_bun_t_guide"
## [1] 32672
## [1] "removing: radi_bun_t_hand"
## [1] 32672
## [1] "removing: radi_bun_t_head"
## [1] 32672
## [1] "removing: radi_bun_t_heart"
## [1] 32672
## [1] "removing: radi_bun_t_heel"
## [1] 32672
## [1] "removing: radi_bun_t_humerus"
## [1] 32672
## [1] "removing: radi_bun_t_hyperthermia"
## [1] 32672
## [1] "removing: radi_bun_t_imag"
## [1] 32672
## [1] "removing: radi_bun_t_image"
## [1] 32672
## [1] "removing: radi_bun_t_imaging"
## [1] 32672
## [1] "removing: radi_bun_t_infant"
## [1] 32672
## [1] "removing: radi_bun_t_interpret"
## [1] 32672
## [1] "removing: radi_bun_t_jaw"
## [1] 32672
## [1] "removing: radi_bun_t_joint"
## [1] 32672
## [1] "removing: radi_bun_t_kidney"
## [1] 32672
## [1] "removing: radi_bun_t_leg"
## [1] 32672
## [1] "removing: radi_bun_t_lesion"
## [1] 32672
## [1] "removing: radi_bun_t_liver"
## [1] 32672
## [1] "removing: radi_bun_t_lower"
## [1] 32672
## [1] "removing: radi_bun_t_lumbar"
## [1] 32671
## [1] "removing: radi_bun_t_lung"
## [1] 32671
## [1] "removing: radi_bun_t_mammo"
## [1] 32671
## [1] "removing: radi_bun_t_mammogra"
## [1] 32671
## [1] "removing: radi_bun_t_management"
## [1] 32671
## [1] "removing: radi_bun_t_marrow"
## [1] 32671
## [1] "removing: radi_bun_t_mastoids"
## [1] 32671
## [1] "removing: radi_bun_t_mra"
## [1] 32671
## [1] "removing: radi_bun_t_mri"
## [1] 32671
## [1] "removing: radi_bun_t_myelography"
## [1] 32671
## [1] "removing: radi_bun_t_neck"
## [1] 32671
## [1] "removing: radi_bun_t_nuclear"
## [1] 32671
## [1] "removing: radi_bun_t_oncology"
## [1] 32671
## [1] "removing: radi_bun_t_pancreas"
## [1] 32671
## [1] "removing: radi_bun_t_pelvi"
## [1] 32671
## [1] "removing: radi_bun_t_penis"
## [1] 32671
## [1] "removing: radi_bun_t_pet"
## [1] 32671
## [1] "removing: radi_bun_t_pituitary"
## [1] 32671
## [1] "removing: radi_bun_t_plan"
## [1] 32671
## [1] "removing: radi_bun_t_plasma"
## [1] 32671
## [1] "removing: radi_bun_t_radiation"
## [1] 32671
## [1] "removing: radi_bun_t_radio"
## [1] 32671
## [1] "removing: radi_bun_t_repair"
## [1] 32671
## [1] "removing: radi_bun_t_rib"
## [1] 32671
## [1] "removing: radi_bun_t_scan"
## [1] 32671
## [1] "removing: radi_bun_t_shoulder"
## [1] 32671
## [1] "removing: radi_bun_t_sinuses"
## [1] 32671
## [1] "removing: radi_bun_t_skull"
## [1] 32671
## [1] "removing: radi_bun_t_socket"
## [1] 32671
## [1] "removing: radi_bun_t_speech"
## [1] 32671
## [1] "removing: radi_bun_t_spine"
## [1] 32671
## [1] "removing: radi_bun_t_spleen"
## [1] 32671
## [1] "removing: radi_bun_t_tear"
## [1] 32671
## [1] "removing: radi_bun_t_teeth"
## [1] 32671
## [1] "removing: radi_bun_t_therapy"
## [1] 32671
## [1] "removing: radi_bun_t_thigh"
## [1] 32671
## [1] "removing: radi_bun_t_thorac"
## [1] 32671
## [1] "removing: radi_bun_t_thorax"
## [1] 32671
## [1] "removing: radi_bun_t_throat"
## [1] 32671
## [1] "removing: radi_bun_t_thrombosis"
## [1] 32671
## [1] "removing: radi_bun_t_thyroid"
## [1] 32671
## [1] "removing: radi_bun_t_toe"
## [1] 32671
## [1] "removing: radi_bun_t_tract"
## [1] 32671
## [1] "removing: radi_bun_t_treatment"
## [1] 32671
## [1] "removing: radi_bun_t_trunk"
## [1] 32671
## [1] "removing: radi_bun_t_tumor"
## [1] 32671
## [1] "removing: radi_bun_t_upper"
## [1] 32671
## [1] "removing: radi_bun_t_urethra"
## [1] 32671
## [1] "removing: radi_bun_t_urinary"
## [1] 32671
## [1] "removing: radi_bun_t_vein"
## [1] 32671
## [1] "removing: radi_bun_t_ventilation"
## [1] 32671
## [1] "removing: radi_bun_t_veous"
## [1] 32671
## [1] "removing: radi_bun_t_vessel"
## [1] 32671
## [1] "removing: radi_bun_t_view"
## [1] 32671
## [1] "removing: radi_bun_t_wrist"
## [1] 32671
## [1] "removing: radi_bun_t_with_dye"
## [1] 32671
## [1] "removing: radi_bun_t_without_dye"
## [1] 32671
## [1] "removing: radi_bun_t_x-ray"
## [1] 32671
## [1] "removing: radi_bun_t_xray"
## [1] 32671
## [1] "removing: radi_bun_t_xtremity"
## [1] 32671
## [1] "removing: path_bun_t_acet"
## [1] 32633
## [1] "removing: path_bun_t_acute"
## [1] 32586
## [1] "removing: path_bun_t_allerg"
## [1] 32557
## [1] "removing: path_bun_t_amikacin"
## [1] 32555
## [1] "removing: path_bun_t_amitripytline"
## [1] 32555
## [1] "removing: path_bun_t_ammonia"
## [1] 32542
## [1] "removing: path_bun_t_amniotic"
## [1] 32542
## [1] "removing: path_bun_t_anabolic"
## [1] 32542
## [1] "removing: path_bun_t_analgesics"
## [1] 32538
## [1] "removing: path_bun_t_analysis"
## [1] 32536
## [1] "removing: path_bun_t_antidepressant"
## [1] 32534
## [1] "removing: path_bun_t_antiepileptics"
## [1] 32534
## [1] "removing: path_bun_t_antipsychotic"
## [1] 32534
## [1] "removing: path_bun_t_aortic"
## [1] 32534
## [1] "removing: path_bun_t_autopsy"
## [1] 32534
## [1] "removing: path_bun_t_barbiturate"
## [1] 32534
## [1] "removing: path_bun_t_benzo"
## [1] 32534
## [1] "removing: path_bun_t_bile"
## [1] 32531
## [1] "removing: path_bun_t_biopsy"
## [1] 32531
## [1] "removing: path_bun_t_bleeding"
## [1] 32531
## [1] "removing: path_bun_t_cadmium"
## [1] 32531
## [1] "removing: path_bun_t_cannabinoids"
## [1] 32530
## [1] "removing: path_bun_t_carbamazepine"
## [1] 32524
## [1] "removing: path_bun_t_carbon"
## [1] 32524
## [1] "removing: path_bun_t_catecholamine"
## [1] 32514
## [1] "removing: path_bun_t_chlamydia"
## [1] 32514
## [1] "removing: path_bun_t_chloride"
## [1] 32487
## [1] "removing: path_bun_t_chromosome"
## [1] 32487
## [1] "removing: path_bun_t_chromotograph"
## [1] 32481
## [1] "removing: path_bun_t_class"
## [1] 32472
## [1] "removing: path_bun_t_clot"
## [1] 32463
## [1] "removing: path_bun_t_clozapine"
## [1] 32463
## [1] "removing: path_bun_t_coagulation"
## [1] 32463
## [1] "removing: path_bun_t_cocaine"
## [1] 32463
## [1] "removing: path_bun_t_cryopreservd"
## [1] 32463
## [1] "removing: path_bun_t_cyanide"
## [1] 32463
## [1] "removing: path_bun_t_cyclosporine"
## [1] 32406
## [1] "removing: path_bun_t_cystines"
## [1] 32406
## [1] "removing: path_bun_t_cytopathol"
## [1] 32406
## [1] "removing: path_bun_t_desipramine"
## [1] 32406
## [1] "removing: path_bun_t_digoxin"
## [1] 32390
## [1] "removing: path_bun_t_dihydro"
## [1] 32382
## [1] "removing: path_bun_t_dilation"
## [1] 32382
## [1] "removing: path_bun_t_dioxide"
## [1] 32382
## [1] "removing: path_bun_t_dipropylacetic"
## [1] 32382
## [1] "removing: path_bun_t_disease"
## [1] 32373
## [1] "removing: path_bun_t_draw"
## [1] 32373
## [1] "removing: path_bun_t_dysfunction"
## [1] 32373
## [1] "removing: path_bun_t_embryo"
## [1] 32361
## [1] "removing: path_bun_t_encephalitis"
## [1] 32361
## [1] "removing: path_bun_t_endocrine"
## [1] 32361
## [1] "removing: path_bun_t_epstein-barr"
## [1] 32299
## [1] "removing: path_bun_t_estrogen"
## [1] 32256
## [1] "removing: path_bun_t_ethanol"
## [1] 32256
## [1] "removing: path_bun_t_evaluation"
## [1] 32256
## [1] "removing: path_bun_t_fecal"
## [1] 32180
## [1] "removing: path_bun_t_fentanyl"
## [1] 32180
## [1] "removing: path_bun_t_fetal"
## [1] 32180
## [1] "removing: path_bun_t_fibrinogen"
## [1] 32175
## [1] "removing: path_bun_t_fibronlytic"
## [1] 32175
## [1] "removing: path_bun_t_flouride"
## [1] 32175
## [1] "removing: path_bun_t_galacto"
## [1] 32175
## [1] "removing: path_bun_t_gases"
## [1] 32149
## [1] "removing: path_bun_t_gastrin"
## [1] 32145
## [1] "removing: path_bun_t_gene"
## [1] 32067
## [1] "removing: path_bun_t_general"
## [1] 32067
## [1] "removing: path_bun_t_genome"
## [1] 32067
## [1] "removing: path_bun_t_genotpying"
## [1] 32067
## [1] "removing: path_bun_t_giardia"
## [1] 32057
## [1] "removing: path_bun_t_glucagon"
## [1] 32057
## [1] "removing: path_bun_t_growth"
## [1] 32032
## [1] "removing: path_bun_t_gynecologic"
## [1] 32032
## [1] "removing: path_bun_t_haloperidol"
## [1] 32032
## [1] "removing: path_bun_t_health"
## [1] 32032
## [1] "removing: path_bun_t_hearing"
## [1] 32032
## [1] "removing: path_bun_t_heavy"
## [1] 32027
## [1] "removing: path_bun_t_heinz"
## [1] 32027
## [1] "removing: path_bun_t_heridatary"
## [1] 32027
## [1] "removing: path_bun_t_heroin"
## [1] 32027
## [1] "removing: path_bun_t_herpes"
## [1] 32012
## [1] "removing: path_bun_t_hold"
## [1] 32012
## [1] "removing: path_bun_t_hpa"
## [1] 32012
## [1] "removing: path_bun_t_hsv"
## [1] 32009
## [1] "removing: path_bun_t_imipramine"
## [1] 32009
## [1] "removing: path_bun_t_immunohisto"
## [1] 32007
## [1] "removing: path_bun_t_iodine"
## [1] 32007
## [1] "removing: path_bun_t_lab"
## [1] 32007
## [1] "removing: path_bun_t_lamotrigine"
## [1] 32000
## [1] "removing: path_bun_t_levetiracetam"
## [1] 32000
## [1] "removing: path_bun_t_lidocaine"
## [1] 32000
## [1] "removing: path_bun_t_loss"
## [1] 32000
## [1] "removing: path_bun_t_low"
## [1] 31998
## [1] "removing: path_bun_t_lung"
## [1] 31998
## [1] "removing: path_bun_t_lyme"
## [1] 31998
## [1] "removing: path_bun_t_lympg"
## [1] 31998
## [1] "removing: path_bun_t_malaria"
## [1] 31998
## [1] "removing: path_bun_t_mercury"
## [1] 31996
## [1] "removing: path_bun_t_metal"
## [1] 31996
## [1] "removing: path_bun_t_microscop"
## [1] 31930
## [1] "removing: path_bun_t_microscopic"
## [1] 31930
## [1] "removing: path_bun_t_microscopy"
## [1] 31930
## [1] "removing: path_bun_t_mitocondrial"
## [1] 31930
## [1] "removing: path_bun_t_mopath"
## [1] 31930
## [1] "removing: path_bun_t_mumps"
## [1] 31901
## [1] "removing: path_bun_t_mycophenolate"
## [1] 31901
## [1] "removing: path_bun_t_nicotine"
## [1] 31901
## [1] "removing: path_bun_t_nile"
## [1] 31898
## [1] "removing: path_bun_t_non-opioid"
## [1] 31898
## [1] "removing: path_bun_t_nucleic"
## [1] 31898
## [1] "removing: path_bun_t_obstretric"
## [1] 31898
## [1] "removing: path_bun_t_occult"
## [1] 31868
## [1] "removing: path_bun_t_onco"
## [1] 31860
## [1] "removing: path_bun_t_oocyte"
## [1] 31860
## [1] "removing: path_bun_t_opioid"
## [1] 31860
## [1] "removing: path_bun_t_ovulation"
## [1] 31860
## [1] "removing: path_bun_t_oxycodone"
## [1] 31859
## [1] "removing: path_bun_t_paper"
## [1] 31859
## [1] "removing: path_bun_t_patho"
## [1] 31851
## [1] "removing: path_bun_t_pathologist"
## [1] 31851
## [1] "removing: path_bun_t_period"
## [1] 31851
## [1] "removing: path_bun_t_phenobarbital"
## [1] 31849
## [1] "removing: path_bun_t_phenytoin"
## [1] 31848
## [1] "removing: path_bun_t_pituitary"
## [1] 31848
## [1] "removing: path_bun_t_plasminogen"
## [1] 31848
## [1] "removing: path_bun_t_poliovirus"
## [1] 31848
## [1] "removing: path_bun_t_potassium"
## [1] 31830
## [1] "removing: path_bun_t_pregabalin"
## [1] 31830
## [1] "removing: path_bun_t_primidone"
## [1] 31830
## [1] "removing: path_bun_t_procainamide"
## [1] 31830
## [1] "removing: path_bun_t_procedure"
## [1] 31830
## [1] "removing: path_bun_t_prophobulinogen"
## [1] 31830
## [1] "removing: path_bun_t_prophyrins"
## [1] 31830
## [1] "removing: path_bun_t_propoxyphene"
## [1] 31830
## [1] "removing: path_bun_t_prostate"
## [1] 31830
## [1] "removing: path_bun_t_quinide"
## [1] 31830
## [1] "removing: path_bun_t_quinine"
## [1] 31830
## [1] "removing: path_bun_t_reagent"
## [1] 31778
## [1] "removing: path_bun_t_renal"
## [1] 31763
## [1] "removing: path_bun_t_retinal"
## [1] 31763
## [1] "removing: path_bun_t_rickettsia"
## [1] 31763
## [1] "removing: path_bun_t_rotovirus"
## [1] 31763
## [1] "removing: path_bun_t_rubella"
## [1] 31758
## [1] "removing: path_bun_t_salicylate"
## [1] 31758
## [1] "removing: path_bun_t_salmonella"
## [1] 31758
## [1] "removing: path_bun_t_scan"
## [1] 31758
## [1] "removing: path_bun_t_sedative"
## [1] 31758
## [1] "removing: path_bun_t_semen"
## [1] 31755
## [1] "removing: path_bun_t_sequence"
## [1] 31755
## [1] "removing: path_bun_t_serotonin"
## [1] 31752
## [1] "removing: path_bun_t_silica"
## [1] 31752
## [1] "removing: path_bun_t_sirolimus"
## [1] 31680
## [1] "removing: path_bun_t_skin"
## [1] 31680
## [1] "removing: path_bun_t_slenium"
## [1] 31680
## [1] "removing: path_bun_t_specimen"
## [1] 31671
## [1] "removing: path_bun_t_sperm"
## [1] 31670
## [1] "removing: path_bun_t_steroid"
## [1] 31670
## [1] "removing: path_bun_t_stimulation"
## [1] 31668
## [1] "removing: path_bun_t_syringe"
## [1] 31668
## [1] "removing: path_bun_t_tacrolims"
## [1] 31668
## [1] "removing: path_bun_t_thawing"
## [1] 31668
## [1] "removing: path_bun_t_theophylline"
## [1] 31668
## [1] "removing: path_bun_t_tissue"
## [1] 31664
## [1] "removing: path_bun_t_tobramycin"
## [1] 31664
## [1] "removing: path_bun_t_tolerance"
## [1] 31662
## [1] "removing: path_bun_t_toxoplasma"
## [1] 31661
## [1] "removing: path_bun_t_tpiramate"
## [1] 31661
## [1] "removing: path_bun_t_transfusion"
## [1] 31661
## [1] "removing: path_bun_t_tumor"
## [1] 31631
## [1] "removing: path_bun_t_vancomycin"
## [1] 31565
## [1] "removing: path_bun_t_vein"
## [1] 31565
## [1] "removing: path_bun_t_venipuncture"
## [1] 31565
## [1] "removing: path_bun_t_virus"
## [1] 31512
## [1] "removing: path_bun_t_west"
## [1] 31512
## [1] "removing: path_bun_t_with_scope"
## [1] 31512
## [1] "removing: path_bun_t_zonisamide"
## [1] 31511
## [1] "removing: anes_bun_t_abdom"
## [1] 31511
## [1] "removing: anes_bun_t_add - on"
## [1] 31511
## [1] "removing: anes_bun_t_adrenal"
## [1] 31511
## [1] "removing: anes_bun_t_amputation"
## [1] 31511
## [1] "removing: anes_bun_t_analg"
## [1] 31511
## [1] "removing: anes_bun_t_anes"
## [1] 31511
## [1] "removing: anes_bun_t_anesth"
## [1] 31511
## [1] "removing: anes_bun_t_area"
## [1] 31511
## [1] "removing: anes_bun_t_arest"
## [1] 31511
## [1] "removing: anes_bun_t_arteriography"
## [1] 31511
## [1] "removing: anes_bun_t_artery"
## [1] 31511
## [1] "removing: anes_bun_t_biopsy"
## [1] 31511
## [1] "removing: anes_bun_t_bladder"
## [1] 31511
## [1] "removing: anes_bun_t_blood"
## [1] 31511
## [1] "removing: anes_bun_t_bone"
## [1] 31511
## [1] "removing: anes_bun_t_breast"
## [1] 31511
## [1] "removing: anes_bun_t_bronchi"
## [1] 31511
## [1] "removing: anes_bun_t_burn"
## [1] 31511
## [1] "removing: anes_bun_t_casting"
## [1] 31511
## [1] "removing: anes_bun_t_cervix"
## [1] 31511
## [1] "removing: anes_bun_t_chemonucleolysis"
## [1] 31511
## [1] "removing: anes_bun_t_chest"
## [1] 31511
## [1] "removing: anes_bun_t_cleft"
## [1] 31511
## [1] "removing: anes_bun_t_cord"
## [1] 31511
## [1] "removing: anes_bun_t_corneal"
## [1] 31511
## [1] "removing: anes_bun_t_cs"
## [1] 31511
## [1] "removing: anes_bun_t_debride"
## [1] 31511
## [1] "removing: anes_bun_t_deliver"
## [1] 31511
## [1] "removing: anes_bun_t_drainage"
## [1] 31511
## [1] "removing: anes_bun_t_duct"
## [1] 31511
## [1] "removing: anes_bun_t_ear"
## [1] 31511
## [1] "removing: anes_bun_t_elbow"
## [1] 31511
## [1] "removing: anes_bun_t_electroshock"
## [1] 31511
## [1] "removing: anes_bun_t_embolectomy"
## [1] 31511
## [1] "removing: anes_bun_t_endoscopy"
## [1] 31511
## [1] "removing: anes_bun_t_epidural"
## [1] 31511
## [1] "removing: anes_bun_t_esophag"
## [1] 31511
## [1] "removing: anes_bun_t_eyelid"
## [1] 31511
## [1] "removing: anes_bun_t_facial"
## [1] 31511
## [1] "removing: anes_bun_t_fat"
## [1] 31511
## [1] "removing: anes_bun_t_femoral"
## [1] 31511
## [1] "removing: anes_bun_t_fentanyl"
## [1] 31511
## [1] "removing: anes_bun_t_genitalia"
## [1] 31511
## [1] "removing: anes_bun_t_glands"
## [1] 31511
## [1] "removing: anes_bun_t_head"
## [1] 31511
## [1] "removing: anes_bun_t_hernia"
## [1] 31511
## [1] "removing: anes_bun_t_humer"
## [1] 31511
## [1] "removing: anes_bun_t_hydromorphone"
## [1] 31511
## [1] "removing: anes_bun_t_hysterectomy"
## [1] 31511
## [1] "removing: anes_bun_t_image"
## [1] 31511
## [1] "removing: anes_bun_t_insertion"
## [1] 31511
## [1] "removing: anes_bun_t_intracranial"
## [1] 31511
## [1] "removing: anes_bun_t_iredectomy"
## [1] 31511
## [1] "removing: anes_bun_t_joint"
## [1] 31511
## [1] "removing: anes_bun_t_ketorolac"
## [1] 31511
## [1] "removing: anes_bun_t_kidney"
## [1] 31511
## [1] "removing: anes_bun_t_knee"
## [1] 31511
## [1] "removing: anes_bun_t_larynx"
## [1] 31511
## [1] "removing: anes_bun_t_layer"
## [1] 31511
## [1] "removing: anes_bun_t_leg"
## [1] 31511
## [1] "removing: anes_bun_t_lip"
## [1] 31511
## [1] "removing: anes_bun_t_liver"
## [1] 31511
## [1] "removing: anes_bun_t_local"
## [1] 31511
## [1] "removing: anes_bun_t_lower"
## [1] 31511
## [1] "removing: anes_bun_t_lumbar"
## [1] 31511
## [1] "removing: anes_bun_t_lung"
## [1] 31511
## [1] "removing: anes_bun_t_major"
## [1] 31511
## [1] "removing: anes_bun_t_manipulation"
## [1] 31511
## [1] "removing: anes_bun_t_midazolam"
## [1] 31511
## [1] "removing: anes_bun_t_mod_sed"
## [1] 31511
## [1] "removing: anes_bun_t_neck"
## [1] 31511
## [1] "removing: anes_bun_t_nerve"
## [1] 31511
## [1] "removing: anes_bun_t_nose"
## [1] 31511
## [1] "removing: anes_bun_t_open"
## [1] 31511
## [1] "removing: anes_bun_t_organ"
## [1] 31511
## [1] "removing: anes_bun_t_pacemaker"
## [1] 31511
## [1] "removing: anes_bun_t_palate"
## [1] 31511
## [1] "removing: anes_bun_t_pelvic"
## [1] 31511
## [1] "removing: anes_bun_t_pelvis"
## [1] 31511
## [1] "removing: anes_bun_t_penis"
## [1] 31511
## [1] "removing: anes_bun_t_percutaneous"
## [1] 31511
## [1] "removing: anes_bun_t_plastic"
## [1] 31511
## [1] "removing: anes_bun_t_propofol"
## [1] 31511
## [1] "removing: anes_bun_t_prostate"
## [1] 31511
## [1] "removing: anes_bun_t_pump"
## [1] 31511
## [1] "removing: anes_bun_t_puncture"
## [1] 31511
## [1] "removing: anes_bun_t_radical"
## [1] 31511
## [1] "removing: anes_bun_t_release"
## [1] 31511
## [1] "removing: anes_bun_t_removal"
## [1] 31511
## [1] "removing: anes_bun_t_repair"
## [1] 31511
## [1] "removing: anes_bun_t_replace"
## [1] 31511
## [1] "removing: anes_bun_t_rib"
## [1] 31511
## [1] "removing: anes_bun_t_sedation"
## [1] 31511
## [1] "removing: anes_bun_t_shoulder"
## [1] 31511
## [1] "removing: anes_bun_t_sinus"
## [1] 31511
## [1] "removing: anes_bun_t_skull"
## [1] 31511
## [1] "removing: anes_bun_t_sperm"
## [1] 31511
## [1] "removing: anes_bun_t_spine"
## [1] 31511
## [1] "removing: anes_bun_t_spit"
## [1] 31511
## [1] "removing: anes_bun_t_sternal"
## [1] 31511
## [1] "removing: anes_bun_t_stone"
## [1] 31511
## [1] "removing: anes_bun_t_surg"
## [1] 31511
## [1] "removing: anes_bun_t_tendon"
## [1] 31511
## [1] "removing: anes_bun_t_testis"
## [1] 31511
## [1] "removing: anes_bun_t_thigh"
## [1] 31511
## [1] "removing: anes_bun_t_thyriod"
## [1] 31511
## [1] "removing: anes_bun_t_trachea"
## [1] 31511
## [1] "removing: anes_bun_t_transplant"
## [1] 31511
## [1] "removing: anes_bun_t_tromethamine"
## [1] 31511
## [1] "removing: anes_bun_t_tubal"
## [1] 31511
## [1] "removing: anes_bun_t_tympanotomy"
## [1] 31511
## [1] "removing: anes_bun_t_upper"
## [1] 31511
## [1] "removing: anes_bun_t_ureter"
## [1] 31511
## [1] "removing: anes_bun_t_urethra"
## [1] 31511
## [1] "removing: anes_bun_t_vag"
## [1] 31511
## [1] "removing: anes_bun_t_vein"
## [1] 31511
## [1] "removing: anes_bun_t_ventilation"
## [1] 31511
## [1] "removing: anes_bun_t_vesectomy"
## [1] 31511
## [1] "removing: anes_bun_t_vessel"
## [1] 31511
## [1] "removing: anes_bun_t_viteroretinal"
## [1] 31511
## [1] "removing: anes_bun_t_wall"
## [1] 31511
## [1] "removing: anes_bun_t_wrist"
## [1] 31511
## [1] "removing: faci_bun_t_advanced"
## [1] 31511
## [1] "removing: faci_bun_t_analysis"
## [1] 31511
## [1] "removing: faci_bun_t_annual"
## [1] 31511
## [1] "removing: faci_bun_t_assessmnt"
## [1] 31511
## [1] "removing: faci_bun_t_audit"
## [1] 31511
## [1] "removing: faci_bun_t_basic"
## [1] 31511
## [1] "removing: faci_bun_t_behav"
## [1] 31501
## [1] "removing: faci_bun_t_bus"
## [1] 31501
## [1] "removing: faci_bun_t_chng"
## [1] 31501
## [1] "removing: faci_bun_t_computer"
## [1] 31501
## [1] "removing: faci_bun_t_conf"
## [1] 31501
## [1] "removing: faci_bun_t_consultation"
## [1] 31489
## [1] "removing: faci_bun_t_counseling"
## [1] 31483
## [1] "removing: faci_bun_t_critical"
## [1] 31482
## [1] "removing: faci_bun_t_data"
## [1] 31482
## [1] "removing: faci_bun_t_direct"
## [1] 31482
## [1] "removing: faci_bun_t_disability"
## [1] 31482
## [1] "removing: faci_bun_t_discharge"
## [1] 31468
## [1] "removing: faci_bun_t_dommicl"
## [1] 31468
## [1] "removing: faci_bun_t_exam"
## [1] 31468
## [1] "removing: faci_bun_t_fac"
## [1] 31433
## [1] "removing: faci_bun_t_facility"
## [1] 31433
## [1] "removing: faci_bun_t_follow-up"
## [1] 31433
## [1] "removing: faci_bun_t_handling"
## [1] 31398
## [1] "removing: faci_bun_t_home"
## [1] 31371
## [1] "removing: faci_bun_t_hospice"
## [1] 31371
## [1] "removing: faci_bun_t_hospital"
## [1] 31365
## [1] "removing: faci_bun_t_initial"
## [1] 31362
## [1] "removing: faci_bun_t_inpatient"
## [1] 31362
## [1] "removing: faci_bun_t_interfacility"
## [1] 31362
## [1] "removing: faci_bun_t_interpof"
## [1] 31362
## [1] "removing: faci_bun_t_lab"
## [1] 31362
## [1] "removing: faci_bun_t_life"
## [1] 31362
## [1] "removing: faci_bun_t_nonemergency"
## [1] 31362
## [1] "removing: faci_bun_t_noner"
## [1] 31362
## [1] "removing: faci_bun_t_nursing"
## [1] 31362
## [1] "removing: faci_bun_t_observation"
## [1] 31362
## [1] "removing: faci_bun_t_on_call"
## [1] 31362
## [1] "removing: faci_bun_t_online"
## [1] 31361
## [1] "removing: faci_bun_t_phlebotomy"
## [1] 31355
## [1] "removing: faci_bun_t_phone"
## [1] 31355
## [1] "removing: faci_bun_t_physician"
## [1] 31355
## [1] "removing: faci_bun_t_postop"
## [1] 31355
## [1] "removing: faci_bun_t_preventine"
## [1] 31355
## [1] "removing: faci_bun_t_prolong"
## [1] 31342
## [1] "removing: faci_bun_t_service"
## [1] 31337
## [1] "removing: faci_bun_t_smoking"
## [1] 31337
## [1] "removing: faci_bun_t_specimen"
## [1] 31337
## [1] "removing: faci_bun_t_standby"
## [1] 31337
## [1] "removing: faci_bun_t_subsequent"
## [1] 31337
## [1] "removing: faci_bun_t_support"
## [1] 31337
## [1] "removing: faci_bun_t_taxi"
## [1] 31337
## [1] "removing: faci_bun_t_team"
## [1] 31337
## [1] "removing: faci_bun_t_transort"
## [1] 31337
## [1] "removing: faci_bun_t_travel"
## [1] 31337
## [1] "removing: faci_bun_t_van"
## [1] 31337
bmp_clean %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 120 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 faci_bun_t_dept            0.729
##  2 faci_bun_t_emergency       0.729
##  3 medi_bun_t_inject          0.456
##  4 medi_bun_t_add-on          0.393
##  5 medi_bun_t_hydrate         0.356
##  6 faci_bun_t_visit           0.248
##  7 medi_bun_t_hydration       0.230
##  8 path_bun_t_complete        0.229
##  9 medi_bun_t_normal          0.207
## 10 medi_bun_t_saline          0.207
## # ... with 110 more rows
bmp_clean %>% get_tag_density_information("faci_bun_t_emergency") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ faci_bun_t_emergency"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 53.7
## $dist_plots

## 
## $stat_tables

bmp_clean <- bmp_clean[faci_bun_t_emergency==F]
bmp_clean %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 118 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 medi_bun_t_add-on             0.242
##  2 faci_bun_t_visit              0.240
##  3 path_bun_t_assay              0.177
##  4 medi_bun_t_inject             0.17 
##  5 medi_bun_t_vacc               0.158
##  6 medi_bun_t_immun              0.156
##  7 medi_bun_t_immunization       0.156
##  8 faci_bun_t_est                0.152
##  9 path_bun_t_antibod            0.149
## 10 medi_bun_t_sulfate            0.146
## # ... with 108 more rows
bmp_clean %>% get_tag_density_information("medi_bun_t_inject") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_inject"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 61.5
## $dist_plots

## 
## $stat_tables

bmp_clean <- bmp_clean[medi_bun_t_inject==F]
bmp_clean %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 117 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_visit              0.264
##  2 medi_bun_t_add-on             0.250
##  3 path_bun_t_assay              0.191
##  4 faci_bun_t_est                0.169
##  5 medi_bun_t_vacc               0.168
##  6 medi_bun_t_immun              0.166
##  7 medi_bun_t_immunization       0.166
##  8 path_bun_t_antibod            0.156
##  9 faci_bun_t_pat                0.156
## 10 faci_bun_t_office             0.156
## # ... with 107 more rows
bmp_clean[`medi_bun_t_add-on`==T] %>% nrow()
## [1] 74
bmp_clean <- bmp_clean[`medi_bun_t_add-on`==F]
bmp_clean %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 115 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_visit              0.272
##  2 path_bun_t_assay              0.192
##  3 medi_bun_t_vacc               0.180
##  4 medi_bun_t_immun              0.178
##  5 medi_bun_t_immunization       0.178
##  6 faci_bun_t_est                0.171
##  7 path_bun_t_antibod            0.166
##  8 faci_bun_t_pat                0.155
##  9 faci_bun_t_office             0.155
## 10 path_bun_t_vitamin            0.155
## # ... with 105 more rows
bmp_clean %>% get_tag_density_information("faci_bun_t_visit") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ faci_bun_t_visit"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 14
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean[faci_bun_t_visit==T]
bmp_clean_n_visit <- bmp_clean[faci_bun_t_visit==F]

##bmp_clean_w_visit

bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 114 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 path_bun_t_assay              0.198
##  2 medi_bun_t_vacc               0.181
##  3 medi_bun_t_immun              0.179
##  4 medi_bun_t_immunization       0.179
##  5 path_bun_t_antibod            0.173
##  6 path_bun_t_immuno             0.159
##  7 path_bun_t_vitamin            0.157
##  8 path_bun_t_immunoassay        0.155
##  9 path_bun_t_complete           0.153
## 10 path_bun_t_dna                0.152
## # ... with 104 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_vacc") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_vacc"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 16.6
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_vacc==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 108 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay             0.200
##  2 path_bun_t_antibod           0.180
##  3 path_bun_t_complete          0.176
##  4 path_bun_t_immuno            0.164
##  5 path_bun_t_vitamin           0.163
##  6 path_bun_t_immunoassay       0.158
##  7 faci_bun_t_est               0.151
##  8 medi_bun_t_tracing           0.148
##  9 path_bun_t_dna               0.147
## 10 path_bun_t_agent             0.144
## # ... with 98 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_antibod") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_antibod"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 30.8
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_antibod==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 106 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 path_bun_t_assay           0.195
##  2 path_bun_t_complete        0.173
##  3 medi_bun_t_tracing         0.160
##  4 path_bun_t_vitamin         0.152
##  5 medi_bun_t_stud            0.142
##  6 path_bun_t_magnesium       0.142
##  7 path_bun_t_hormone         0.137
##  8 path_bun_t_dna             0.135
##  9 faci_bun_t_est             0.132
## 10 path_bun_t_free            0.130
## # ... with 96 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_vitamin") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_vitamin"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 16.4
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_vitamin==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 105 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 medi_bun_t_tracing         0.175
##  2 path_bun_t_assay           0.167
##  3 path_bun_t_complete        0.157
##  4 medi_bun_t_stud            0.151
##  5 path_bun_t_time            0.142
##  6 path_bun_t_dna             0.142
##  7 medi_bun_t_monitor         0.138
##  8 path_bun_t_magnesium       0.137
##  9 medi_bun_t_visit           0.132
## 10 medi_bun_t_clinic          0.129
## # ... with 95 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_stud") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_stud"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 63.9
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_stud==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 104 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 medi_bun_t_tracing         0.180
##  2 path_bun_t_assay           0.170
##  3 path_bun_t_complete        0.157
##  4 path_bun_t_dna             0.144
##  5 path_bun_t_time            0.144
##  6 medi_bun_t_monitor         0.139
##  7 path_bun_t_magnesium       0.132
##  8 medi_bun_t_visit           0.132
##  9 medi_bun_t_sulfate         0.130
## 10 medi_bun_t_clinic          0.129
## # ... with 94 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_dna") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_dna"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 28.7
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_dna==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 102 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 path_bun_t_assay           0.181
##  2 medi_bun_t_tracing         0.173
##  3 path_bun_t_complete        0.153
##  4 medi_bun_t_monitor         0.145
##  5 path_bun_t_magnesium       0.136
##  6 medi_bun_t_sulfate         0.134
##  7 path_bun_t_time            0.132
##  8 medi_bun_t_visit           0.132
##  9 medi_bun_t_clinic          0.129
## 10 medi_bun_t_outpt           0.129
## # ... with 92 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_monitor") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_monitor"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 40.7
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_monitor==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 101 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 path_bun_t_assay           0.180
##  2 medi_bun_t_tracing         0.165
##  3 path_bun_t_complete        0.155
##  4 medi_bun_t_sulfate         0.137
##  5 path_bun_t_time            0.133
##  6 faci_bun_t_est             0.127
##  7 path_bun_t_magnesium       0.127
##  8 medi_bun_t_visit           0.122
##  9 medi_bun_t_clinic          0.119
## 10 medi_bun_t_outpt           0.119
## # ... with 91 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_sulfate") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_sulfate"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 102
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_sulfate==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 99 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 path_bun_t_assay              0.183
##  2 medi_bun_t_tracing            0.167
##  3 path_bun_t_complete           0.153
##  4 path_bun_t_time               0.132
##  5 faci_bun_t_est                0.13 
##  6 path_bun_t_magnesium          0.124
##  7 medi_bun_t_visit              0.121
##  8 path_bun_t_testosterone       0.118
##  9 medi_bun_t_clinic             0.117
## 10 medi_bun_t_outpt              0.117
## # ... with 89 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_time") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_time"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 25.9
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_time==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 98 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 path_bun_t_assay              0.201
##  2 path_bun_t_complete           0.136
##  3 faci_bun_t_est                0.131
##  4 path_bun_t_magnesium          0.128
##  5 medi_bun_t_visit              0.128
##  6 path_bun_t_hormone            0.127
##  7 medi_bun_t_clinic             0.126
##  8 medi_bun_t_outpt              0.126
##  9 path_bun_t_testosterone       0.126
## 10 medi_bun_t_tracing            0.122
## # ... with 88 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_magnesium") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_magnesium"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.1
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_magnesium==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 94 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 path_bun_t_assay              0.190
##  2 faci_bun_t_est                0.139
##  3 path_bun_t_lipid              0.132
##  4 path_bun_t_complete           0.130
##  5 medi_bun_t_tracing            0.129
##  6 path_bun_t_hormone            0.129
##  7 path_bun_t_testosterone       0.121
##  8 path_bun_t_thyroid            0.119
##  9 faci_bun_t_office             0.112
## 10 faci_bun_t_pat                0.112
## # ... with 84 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_lipid") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_lipid"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 11
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_lipid==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 91 x 2
##    name                              correlation
##    <chr>                                   <dbl>
##  1 path_bun_t_assay                        0.166
##  2 path_bun_t_complete                     0.159
##  3 medi_bun_t_tracing                      0.153
##  4 faci_bun_t_est                          0.148
##  5 path_bun_t_testosterone                 0.120
##  6 path_bun_t_free                         0.120
##  7 path_bun_t_hormone                      0.114
##  8 path_bun_t_dehydroepiandrosterone       0.107
##  9 medi_bun_t_exercise                     0.105
## 10 path_bun_t_protein                      0.104
## # ... with 81 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_tracing") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_tracing"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 23.5
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_tracing==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 90 x 2
##    name                              correlation
##    <chr>                                   <dbl>
##  1 path_bun_t_assay                        0.177
##  2 path_bun_t_complete                     0.146
##  3 faci_bun_t_est                          0.146
##  4 path_bun_t_free                         0.130
##  5 path_bun_t_testosterone                 0.129
##  6 path_bun_t_hormone                      0.12 
##  7 path_bun_t_dehydroepiandrosterone       0.112
##  8 path_bun_t_iron                         0.107
##  9 path_bun_t_thyroid                      0.105
## 10 path_bun_t_protein                      0.103
## # ... with 80 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_testosterone") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_testosterone"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 27.9
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_testosterone==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 89 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay             0.161
##  2 path_bun_t_complete          0.153
##  3 faci_bun_t_est               0.139
##  4 path_bun_t_iron              0.108
##  5 path_bun_t_protein           0.107
##  6 path_bun_t_hormone           0.107
##  7 path_bun_t_free              0.106
##  8 medi_bun_t_therapeutic       0.104
##  9 medi_bun_t_exercise          0.103
## 10 medi_bun_t_clinic            0.100
## # ... with 79 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_hormone") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_hormone"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 12.5
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_hormone==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 86 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_complete         0.153 
##  2 path_bun_t_assay            0.133 
##  3 faci_bun_t_est              0.121 
##  4 path_bun_t_protein          0.115 
##  5 medi_bun_t_therapeutic      0.110 
##  6 medi_bun_t_exercise         0.109 
##  7 medi_bun_t_clinic           0.104 
##  8 medi_bun_t_outpt            0.104 
##  9 medi_bun_t_visit            0.103 
## 10 medi_bun_t_test             0.0984
## # ... with 76 more rows
bmp_clean_w_visit %>% get_tag_density_information("faci_bun_t_office") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ faci_bun_t_office"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 13.9
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[faci_bun_t_office==T]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 84 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_complete          0.160
##  2 path_bun_t_assay             0.132
##  3 faci_bun_t_est               0.121
##  4 path_bun_t_protein           0.118
##  5 medi_bun_t_therapeutic       0.112
##  6 medi_bun_t_exercise          0.111
##  7 medi_bun_t_clinic            0.108
##  8 medi_bun_t_outpt             0.108
##  9 medi_bun_t_visit             0.107
## 10 medi_bun_t_test              0.102
## # ... with 74 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_clinic") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_clinic"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 44.8
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_clinic==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 82 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_complete         0.159 
##  2 path_bun_t_assay            0.123 
##  3 faci_bun_t_est              0.122 
##  4 path_bun_t_protein          0.117 
##  5 medi_bun_t_therapeutic      0.117 
##  6 medi_bun_t_exercise         0.116 
##  7 medi_bun_t_test             0.101 
##  8 path_bun_t_count            0.0976
##  9 path_bun_t_smear            0.0964
## 10 medi_bun_t_manual           0.0893
## # ... with 72 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_therapeutic") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_therapeutic"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 29.3
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_therapeutic==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 80 x 2
##    name                  correlation
##    <chr>                       <dbl>
##  1 path_bun_t_complete        0.160 
##  2 path_bun_t_assay           0.128 
##  3 faci_bun_t_est             0.124 
##  4 path_bun_t_protein         0.119 
##  5 medi_bun_t_test            0.103 
##  6 path_bun_t_count           0.100 
##  7 path_bun_t_smear           0.0984
##  8 path_bun_t_urinalysis      0.0844
##  9 medi_bun_t_psytx           0.0837
## 10 medi_bun_t_hydration       0.0797
## # ... with 70 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_protein") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_protein"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.4
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_protein==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 76 x 2
##    name                  correlation
##    <chr>                       <dbl>
##  1 path_bun_t_complete        0.136 
##  2 path_bun_t_assay           0.127 
##  3 faci_bun_t_est             0.124 
##  4 medi_bun_t_test            0.108 
##  5 path_bun_t_urinalysis      0.0848
##  6 path_bun_t_iron            0.0844
##  7 medi_bun_t_hydration       0.0843
##  8 path_bun_t_count           0.084 
##  9 path_bun_t_ferritin        0.0795
## 10 medi_bun_t_psytx           0.0781
## # ... with 66 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_psa") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_psa"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 17.9
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[path_bun_t_psa==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 75 x 2
##    name                  correlation
##    <chr>                       <dbl>
##  1 path_bun_t_complete        0.137 
##  2 path_bun_t_assay           0.122 
##  3 faci_bun_t_est             0.120 
##  4 medi_bun_t_test            0.110 
##  5 path_bun_t_urinalysis      0.0862
##  6 path_bun_t_iron            0.0858
##  7 medi_bun_t_hydration       0.0853
##  8 path_bun_t_count           0.0822
##  9 path_bun_t_ferritin        0.0799
## 10 medi_bun_t_psytx           0.0792
## # ... with 65 more rows
bmp_clean_w_visit %>% get_tag_density_information("faci_bun_t_est") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ faci_bun_t_est"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 14.8
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[faci_bun_t_est==T]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 74 x 2
##    name                  correlation
##    <chr>                       <dbl>
##  1 path_bun_t_complete        0.139 
##  2 medi_bun_t_test            0.124 
##  3 path_bun_t_assay           0.119 
##  4 path_bun_t_urinalysis      0.0879
##  5 medi_bun_t_psytx           0.0867
##  6 path_bun_t_iron            0.083 
##  7 path_bun_t_count           0.0824
##  8 path_bun_t_smear           0.0744
##  9 medi_bun_t_hydration       0.0743
## 10 path_bun_t_ferritin        0.0742
## # ... with 64 more rows
bmp_clean_w_visit %>% get_tag_density_information("medi_bun_t_test") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_test"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 105
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_test==F]
bmp_clean_w_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_w_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 73 x 2
##    name                  correlation
##    <chr>                       <dbl>
##  1 path_bun_t_complete        0.144 
##  2 path_bun_t_assay           0.123 
##  3 path_bun_t_urinalysis      0.0915
##  4 medi_bun_t_psytx           0.089 
##  5 path_bun_t_iron            0.0854
##  6 path_bun_t_count           0.0851
##  7 path_bun_t_ferritin        0.0764
##  8 path_bun_t_smear           0.0763
##  9 medi_bun_t_hydration       0.0762
## 10 medi_bun_t_normal          0.0751
## # ... with 63 more rows
bmp_clean_w_visit %>% get_tag_density_information("path_bun_t_complete") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_complete"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 11.4
## $dist_plots

## 
## $stat_tables

bmp_clean_w_visit <- bmp_clean_w_visit[medi_bun_t_test==T]
bmp_clean_w_visit <- bmp_clean_w_visit[,`:=`(
  surg_sp_name_clean = surg_sp_npi %>% map_chr(get_npi_standard_name),
  surg_bp_name_clean = surg_bp_npi %>% map_chr(get_npi_standard_name),
  
  medi_sp_name_clean = medi_sp_npi %>% map_chr(get_npi_standard_name),
  medi_bp_name_clean = medi_bp_npi %>% map_chr(get_npi_standard_name),
  
  radi_sp_name_clean = radi_sp_npi %>% map_chr(get_npi_standard_name),
  radi_bp_name_clean = radi_bp_npi %>% map_chr(get_npi_standard_name),
  
  path_sp_name_clean = path_sp_npi %>% map_chr(get_npi_standard_name),
  path_bp_name_clean = path_bp_npi %>% map_chr(get_npi_standard_name),
  
  anes_sp_name_clean = anes_sp_npi %>% map_chr(get_npi_standard_name),
  anes_bp_name_clean = anes_bp_npi %>% map_chr(get_npi_standard_name),
  
  faci_sp_name_clean = faci_sp_npi %>% map_chr(get_npi_standard_name),
  faci_bp_name_clean = faci_bp_npi %>% map_chr(get_npi_standard_name)
                           )]
bmp_clean_w_visit_btbv5 <- bmp_clean_w_visit %>% btbv5_medi_fac_only(
                                  osa_group_p=NA,
                                  osa_class_p = NA,
                                  osa_specialization_p = NA
                                  ) %>% 
  .[,`:=`(procedure_type=26, 
          procedure_modifier="Basic Metabolic Panel + Visit", 
          ingest_date =Sys.Date())]
## Warning in max(most_important_fac_npi): no non-missing arguments to max;
## returning -Inf
bmp_clean_w_visit_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## 
spine_q1 <- bmp_clean_w_visit_btbv5$tp_med_med %>% quantile(.25)
spine_q3 <- bmp_clean_w_visit_btbv5$tp_med_med %>% quantile(.75)
spine_iqr <- bmp_clean_w_visit_btbv5$tp_med_med %>% IQR()
bmp_clean_w_visit_btbv5_final <- bmp_clean_w_visit_btbv5[tp_med_med >spine_q1 -1.5*spine_iqr & tp_med_med < spine_q3 + 1.5*spine_iqr ]
bmp_clean_w_visit_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## 
bmp_clean_w_visit_btbv5_final$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## 
bq_table_upload(x=procedure_dev_table, values= bmp_clean_w_visit_btbv5_final, create_disposition='CREATE_IF_NEEDED', write_disposition='WRITE_APPEND')

##bmp_clean_n_visit

bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 113 x 2
##    name                              correlation
##    <chr>                                   <dbl>
##  1 medi_bun_t_sulfate                      0.226
##  2 medi_bun_t_sodium                       0.206
##  3 path_bun_t_dehydroepiandrosterone       0.188
##  4 path_bun_t_progesterone                 0.184
##  5 path_bun_t_antibod                      0.175
##  6 medi_bun_t_hepa                         0.165
##  7 path_bun_t_magnesium                    0.158
##  8 path_bun_t_vitamin                      0.156
##  9 path_bun_t_assay                        0.154
## 10 path_bun_t_free                         0.148
## # ... with 103 more rows
bmp_clean_n_visit %>% get_tag_density_information("medi_bun_t_sulfate") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_sulfate"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 161
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[medi_bun_t_sulfate==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 112 x 2
##    name                              correlation
##    <chr>                                   <dbl>
##  1 path_bun_t_dehydroepiandrosterone       0.198
##  2 path_bun_t_progesterone                 0.194
##  3 path_bun_t_antibod                      0.185
##  4 path_bun_t_vitamin                      0.166
##  5 path_bun_t_free                         0.159
##  6 medi_bun_t_tracing                      0.156
##  7 path_bun_t_insulin                      0.153
##  8 path_bun_t_complete                     0.152
##  9 path_bun_t_testosterone                 0.152
## 10 path_bun_t_assay                        0.15 
## # ... with 102 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_dehydroepiandrosterone") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_dehydroepiandrosterone"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 60.8
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_dehydroepiandrosterone==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 111 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 medi_bun_t_tracing         0.166
##  2 path_bun_t_complete        0.143
##  3 path_bun_t_antibod         0.140
##  4 path_bun_t_vitamin         0.135
##  5 path_bun_t_assay           0.134
##  6 medi_bun_t_stud            0.114
##  7 path_bun_t_free            0.110
##  8 path_bun_t_magnesium       0.109
##  9 path_bun_t_screen          0.104
## 10 path_bun_t_hormone         0.102
## # ... with 101 more rows
bmp_clean_n_visit %>% get_tag_density_information("medi_bun_t_tracing") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_tracing"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 27.2
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[medi_bun_t_tracing==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 110 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 path_bun_t_assay           0.167
##  2 path_bun_t_vitamin         0.156
##  3 path_bun_t_antibod         0.149
##  4 path_bun_t_complete        0.128
##  5 path_bun_t_hormone         0.127
##  6 path_bun_t_free            0.126
##  7 medi_bun_t_stud            0.122
##  8 path_bun_t_enzyme          0.11 
##  9 path_bun_t_immuno          0.109
## 10 path_bun_t_magnesium       0.107
## # ... with 100 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_vitamin") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_vitamin"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 24
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_vitamin==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 109 x 2
##    name                correlation
##    <chr>                     <dbl>
##  1 path_bun_t_antibod       0.142 
##  2 path_bun_t_assay         0.133 
##  3 medi_bun_t_stud          0.122 
##  4 path_bun_t_complete      0.119 
##  5 path_bun_t_screen        0.115 
##  6 path_bun_t_typing        0.114 
##  7 path_bun_t_enzyme        0.105 
##  8 path_bun_t_dna           0.0938
##  9 path_bun_t_immuno        0.0884
## 10 medi_bun_t_test          0.0836
## # ... with 99 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_antibod") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_antibod"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 35.7
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_antibod==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 107 x 2
##    name                 correlation
##    <chr>                      <dbl>
##  1 path_bun_t_assay          0.130 
##  2 medi_bun_t_stud           0.121 
##  3 path_bun_t_complete       0.105 
##  4 path_bun_t_enzyme         0.104 
##  5 medi_bun_t_test           0.0875
##  6 path_bun_t_dna            0.0835
##  7 path_bun_t_magnesium      0.079 
##  8 faci_bun_t_care           0.0773
##  9 path_bun_t_free           0.0767
## 10 medi_bun_t_monitor        0.0757
## # ... with 97 more rows
bmp_clean_n_visit %>% get_tag_density_information("medi_bun_t_stud") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_stud"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 96.8
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[medi_bun_t_stud==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 106 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 path_bun_t_assay             0.133 
##  2 path_bun_t_enzyme            0.106 
##  3 path_bun_t_complete          0.104 
##  4 path_bun_t_dna               0.0856
##  5 path_bun_t_magnesium         0.081 
##  6 path_bun_t_free              0.0792
##  7 faci_bun_t_care              0.0788
##  8 medi_bun_t_monitor           0.0772
##  9 medi_bun_t_psytx             0.0765
## 10 path_bun_t_testosterone      0.0763
## # ... with 96 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_enzyme") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_enzyme"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 39.1
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_enzyme==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 105 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 path_bun_t_assay             0.124 
##  2 path_bun_t_complete          0.0915
##  3 path_bun_t_free              0.084 
##  4 faci_bun_t_care              0.0813
##  5 path_bun_t_hormone           0.0808
##  6 path_bun_t_testosterone      0.0801
##  7 medi_bun_t_monitor           0.0797
##  8 medi_bun_t_psytx             0.0793
##  9 medi_bun_t_rehab             0.0783
## 10 path_bun_t_dna               0.0776
## # ... with 95 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_free") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_free"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 21.7
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_free==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 104 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay            0.108 
##  2 path_bun_t_complete         0.0929
##  3 faci_bun_t_care             0.086 
##  4 medi_bun_t_psytx            0.0844
##  5 medi_bun_t_monitor          0.0842
##  6 medi_bun_t_rehab            0.0828
##  7 medi_bun_t_hydration        0.0815
##  8 path_bun_t_phosphorus       0.0804
##  9 path_bun_t_dna              0.0726
## 10 medi_bun_t_therapeutic      0.0714
## # ... with 94 more rows
bmp_clean_n_visit %>% get_tag_density_information("medi_bun_t_psytx") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_psytx"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.3
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[medi_bun_t_psytx==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 103 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay            0.108 
##  2 path_bun_t_complete         0.0957
##  3 faci_bun_t_care             0.0872
##  4 medi_bun_t_monitor          0.0864
##  5 medi_bun_t_rehab            0.085 
##  6 medi_bun_t_hydration        0.0834
##  7 path_bun_t_phosphorus       0.0832
##  8 medi_bun_t_sodium           0.0727
##  9 path_bun_t_dna              0.071 
## 10 medi_bun_t_therapeutic      0.0702
## # ... with 93 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_testosterone") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_testosterone"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 23.9
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_testosterone==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 102 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay            0.100 
##  2 path_bun_t_complete         0.0958
##  3 faci_bun_t_care             0.0892
##  4 medi_bun_t_monitor          0.0884
##  5 medi_bun_t_rehab            0.0869
##  6 path_bun_t_phosphorus       0.0858
##  7 medi_bun_t_hydration        0.0851
##  8 medi_bun_t_sodium           0.074 
##  9 path_bun_t_dna              0.0731
## 10 medi_bun_t_therapeutic      0.072 
## # ... with 92 more rows
bmp_clean_n_visit %>% get_tag_density_information("medi_bun_t_monitor") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ medi_bun_t_monitor"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 27.3
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[medi_bun_t_monitor==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 100 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay            0.102 
##  2 path_bun_t_complete         0.0984
##  3 faci_bun_t_care             0.0895
##  4 path_bun_t_phosphorus       0.0878
##  5 medi_bun_t_hydration        0.0864
##  6 path_bun_t_dna              0.075 
##  7 medi_bun_t_sodium           0.075 
##  8 medi_bun_t_therapeutic      0.0733
##  9 path_bun_t_hiv              0.071 
## 10 path_bun_t_gonorr           0.0669
## # ... with 90 more rows
bmp_clean_n_visit %>% get_tag_density_information("faci_bun_t_care") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ faci_bun_t_care"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 24.3
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[faci_bun_t_care==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 99 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay            0.105 
##  2 path_bun_t_complete         0.0969
##  3 path_bun_t_phosphorus       0.0895
##  4 medi_bun_t_hydration        0.0873
##  5 path_bun_t_dna              0.0767
##  6 medi_bun_t_sodium           0.0756
##  7 medi_bun_t_therapeutic      0.0744
##  8 path_bun_t_hiv              0.0722
##  9 path_bun_t_gonorr           0.0679
## 10 medi_bun_t_normal           0.0641
## # ... with 89 more rows
bmp_clean_n_visit %>% get_tag_density_information("path_bun_t_dna") %>% print()
## Scale for 'y' is already present. Adding another scale for 'y', which will
## replace the existing scale.
## Scale for 'x' is already present. Adding another scale for 'x', which will
## replace the existing scale.
## [1] "tp_med ~ path_bun_t_dna"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 20.3
## $dist_plots

## 
## $stat_tables

bmp_clean_n_visit <- bmp_clean_n_visit[path_bun_t_dna==F]
bmp_clean_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

bmp_clean_n_visit%>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 96 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 path_bun_t_assay            0.105 
##  2 path_bun_t_phosphorus       0.0922
##  3 medi_bun_t_hydration        0.0902
##  4 path_bun_t_complete         0.09  
##  5 medi_bun_t_sodium           0.0778
##  6 medi_bun_t_therapeutic      0.0774
##  7 medi_bun_t_normal           0.0662
##  8 medi_bun_t_saline           0.0662
##  9 medi_bun_t_solution         0.0662
## 10 path_bun_t_protein          0.0643
## # ... with 86 more rows
bmp_clean_n_visit <- bmp_clean_n_visit[,`:=`(
  surg_sp_name_clean = surg_sp_npi %>% map_chr(get_npi_standard_name),
  surg_bp_name_clean = surg_bp_npi %>% map_chr(get_npi_standard_name),
  
  medi_sp_name_clean = medi_sp_npi %>% map_chr(get_npi_standard_name),
  medi_bp_name_clean = medi_bp_npi %>% map_chr(get_npi_standard_name),
  
  radi_sp_name_clean = radi_sp_npi %>% map_chr(get_npi_standard_name),
  radi_bp_name_clean = radi_bp_npi %>% map_chr(get_npi_standard_name),
  
  path_sp_name_clean = path_sp_npi %>% map_chr(get_npi_standard_name),
  path_bp_name_clean = path_bp_npi %>% map_chr(get_npi_standard_name),
  
  anes_sp_name_clean = anes_sp_npi %>% map_chr(get_npi_standard_name),
  anes_bp_name_clean = anes_bp_npi %>% map_chr(get_npi_standard_name),
  
  faci_sp_name_clean = faci_sp_npi %>% map_chr(get_npi_standard_name),
  faci_bp_name_clean = faci_bp_npi %>% map_chr(get_npi_standard_name)
                           )]
bmp_clean_n_visit_btbv5 <- bmp_clean_n_visit %>% btbv5_medi_fac_only(
                                  osa_group_p=NA,
                                  osa_class_p = NA,
                                  osa_specialization_p = NA
                                  ) %>% 
  .[,`:=`(procedure_type=26, 
          procedure_modifier="Standard", 
          ingest_date =Sys.Date())]
bmp_clean_n_visit_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    2.40   36.72  122.50  168.55  196.82 1640.89
spine_q1 <- bmp_clean_n_visit_btbv5$tp_med_med %>% quantile(.25)
spine_q3 <- bmp_clean_n_visit_btbv5$tp_med_med %>% quantile(.75)
spine_iqr <- bmp_clean_n_visit_btbv5$tp_med_med %>% IQR()
bmp_clean_n_visit_btbv5_final <- bmp_clean_n_visit_btbv5[tp_med_med >spine_q1 -1.5*spine_iqr & tp_med_med < spine_q3 + 1.5*spine_iqr ]
bmp_clean_n_visit_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    2.40   36.72  122.50  168.55  196.82 1640.89
bmp_clean_n_visit_btbv5_final$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    2.40   33.07  111.53  127.94  185.02  433.39
bq_table_upload(x=procedure_dev_table, values= bmp_clean_n_visit_btbv5_final, create_disposition='CREATE_IF_NEEDED', write_disposition='WRITE_APPEND')