Introduction

In this section we will:
* identify which Spine X-ray procedures were excluded and why * identify which different Spine X-ray 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")
xray_spine <- bun_proc[(radi_bun_t_xray==T | `radi_bun_t_x-ray`) & radi_bun_t_spine==T &
                        tp_med >0 &
                        duration_max ==0 &
                        duration_min == 0 & 
                         tp_med < 10000]
xray_spine %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

xray_spine %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 612 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 faci_bun_t_dept              0.534
##  2 faci_bun_t_emergency         0.534
##  3 radi_bun_t_ct                0.480
##  4 radi_bun_t_without_dye       0.464
##  5 radi_bun_t_brain             0.442
##  6 medi_bun_t_inject            0.430
##  7 anes_bun_t_anes              0.430
##  8 radi_bun_t_head              0.429
##  9 anes_bun_t_anesth            0.427
## 10 anes_bun_t_surg              0.412
## # ... with 602 more rows
cols_to_check <-  xray_spine %>% 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 < 1000] %>% .[["tag"]]
xray_spine_c <- remove_bun_proc_with_rare_tags(outlier_tag_list=outlier_tag_list, df=xray_spine)
## [1] "Start of rare tag removal"
## [1] "nrow:35624 , ncol:1404"
## [1] "nrow:35624"
## [1] "removing: surg_bun_t_a-grft"
## [1] 35624
## [1] "removing: surg_bun_t_abd"
## [1] 35623
## [1] "removing: surg_bun_t_abdom"
## [1] 35623
## [1] "removing: surg_bun_t_abdomen"
## [1] 35623
## [1] "removing: surg_bun_t_ablate"
## [1] 35622
## [1] "removing: surg_bun_t_abortion"
## [1] 35622
## [1] "removing: surg_bun_t_abscess"
## [1] 35620
## [1] "removing: surg_bun_t_acell"
## [1] 35620
## [1] "removing: surg_bun_t_acellular"
## [1] 35620
## [1] "removing: surg_bun_t_acne"
## [1] 35620
## [1] "removing: surg_bun_t_acute"
## [1] 35620
## [1] "removing: surg_bun_t_add-on"
## [1] 35497
## [1] "removing: surg_bun_t_adenoid"
## [1] 35497
## [1] "removing: surg_bun_t_adjust"
## [1] 35497
## [1] "removing: surg_bun_t_airway"
## [1] 35493
## [1] "removing: surg_bun_t_amput"
## [1] 35493
## [1] "removing: surg_bun_t_anal"
## [1] 35484
## [1] "removing: surg_bun_t_analysis"
## [1] 35484
## [1] "removing: surg_bun_t_aneurysm"
## [1] 35484
## [1] "removing: surg_bun_t_ankle"
## [1] 35480
## [1] "removing: surg_bun_t_anoscopy"
## [1] 35479
## [1] "removing: surg_bun_t_ant"
## [1] 35466
## [1] "removing: surg_bun_t_antepartum"
## [1] 35466
## [1] "removing: surg_bun_t_aorta"
## [1] 35466
## [1] "removing: surg_bun_t_aortic"
## [1] 35466
## [1] "removing: surg_bun_t_aply"
## [1] 35466
## [1] "removing: surg_bun_t_appendectomy"
## [1] 35466
## [1] "removing: surg_bun_t_application"
## [1] 35453
## [1] "removing: surg_bun_t_apply"
## [1] 35416
## [1] "removing: surg_bun_t_approach"
## [1] 35416
## [1] "removing: surg_bun_t_arm"
## [1] 35416
## [1] "removing: surg_bun_t_art"
## [1] 35394
## [1] "removing: surg_bun_t_arterial"
## [1] 35394
## [1] "removing: surg_bun_t_artery"
## [1] 35394
## [1] "removing: surg_bun_t_arthro"
## [1] 35394
## [1] "removing: surg_bun_t_arthroplasty"
## [1] 35394
## [1] "removing: surg_bun_t_arthroscopy"
## [1] 35394
## [1] "removing: surg_bun_t_artificial"
## [1] 35394
## [1] "removing: surg_bun_t_augmentation"
## [1] 35393
## [1] "removing: surg_bun_t_autograft"
## [1] 35393
## [1] "removing: surg_bun_t_b9"
## [1] 35356
## [1] "removing: surg_bun_t_back"
## [1] 35301
## [1] "removing: surg_bun_t_balloon"
## [1] 35301
## [1] "removing: surg_bun_t_bed"
## [1] 35295
## [1] "removing: surg_bun_t_biceps"
## [1] 35295
## [1] "removing: surg_bun_t_bile"
## [1] 35295
## [1] "removing: surg_bun_t_bilobectomy"
## [1] 35295
## [1] "removing: surg_bun_t_biopsy"
## [1] 35268
## [1] "removing: surg_bun_t_bladder"
## [1] 35262
## [1] "removing: surg_bun_t_bleeding"
## [1] 35262
## [1] "removing: surg_bun_t_block"
## [1] 35240
## [1] "removing: surg_bun_t_blockage"
## [1] 35240
## [1] "removing: surg_bun_t_blood"
## [1] 35240
## [1] "removing: surg_bun_t_body"
## [1] 35211
## [1] "removing: surg_bun_t_bone"
## [1] 35201
## [1] "removing: surg_bun_t_bowel"
## [1] 35201
## [1] "removing: surg_bun_t_brace"
## [1] 35201
## [1] "removing: surg_bun_t_brain"
## [1] 35197
## [1] "removing: surg_bun_t_breast"
## [1] 35195
## [1] "removing: surg_bun_t_bronchial"
## [1] 35195
## [1] "removing: surg_bun_t_bronchoscopy"
## [1] 35195
## [1] "removing: surg_bun_t_brow"
## [1] 35195
## [1] "removing: surg_bun_t_bulge"
## [1] 35195
## [1] "removing: surg_bun_t_burn"
## [1] 35193
## [1] "removing: surg_bun_t_burr"
## [1] 35193
## [1] "removing: surg_bun_t_bursa"
## [1] 34788
## [1] "removing: surg_bun_t_buttck"
## [1] 34788
## [1] "removing: surg_bun_t_buttock"
## [1] 34788
## [1] "removing: surg_bun_t_byp"
## [1] 34788
## [1] "removing: surg_bun_t_bypass"
## [1] 34788
## [1] "removing: surg_bun_t_cabg"
## [1] 34788
## [1] "removing: surg_bun_t_canal"
## [1] 34788
## [1] "removing: surg_bun_t_cannula"
## [1] 34788
## [1] "removing: surg_bun_t_capillary"
## [1] 34788
## [1] "removing: surg_bun_t_capsule"
## [1] 34788
## [1] "removing: surg_bun_t_care"
## [1] 34788
## [1] "removing: surg_bun_t_carotid"
## [1] 34788
## [1] "removing: surg_bun_t_cartilage"
## [1] 34788
## [1] "removing: surg_bun_t_cast"
## [1] 34776
## [1] "removing: surg_bun_t_cath"
## [1] 34775
## [1] "removing: surg_bun_t_cavity"
## [1] 34775
## [1] "removing: surg_bun_t_cell"
## [1] 34775
## [1] "removing: surg_bun_t_cervix"
## [1] 34775
## [1] "removing: surg_bun_t_cesarean"
## [1] 34775
## [1] "removing: surg_bun_t_chamber"
## [1] 34775
## [1] "removing: surg_bun_t_change"
## [1] 34775
## [1] "removing: surg_bun_t_chemical"
## [1] 34775
## [1] "removing: surg_bun_t_chemodenerv"
## [1] 34765
## [1] "removing: surg_bun_t_chest"
## [1] 34764
## [1] "removing: surg_bun_t_chin"
## [1] 34764
## [1] "removing: surg_bun_t_cholangiogram"
## [1] 34764
## [1] "removing: surg_bun_t_cholangiopancreatograph"
## [1] 34764
## [1] "removing: surg_bun_t_cholecystoenterostomy"
## [1] 34764
## [1] "removing: surg_bun_t_chronic"
## [1] 34764
## [1] "removing: surg_bun_t_circulation"
## [1] 34764
## [1] "removing: surg_bun_t_circumcision"
## [1] 34764
## [1] "removing: surg_bun_t_clamp"
## [1] 34764
## [1] "removing: surg_bun_t_clavicle"
## [1] 34747
## [1] "removing: surg_bun_t_clearance"
## [1] 34747
## [1] "removing: surg_bun_t_cleft"
## [1] 34747
## [1] "removing: surg_bun_t_clitoris"
## [1] 34747
## [1] "removing: surg_bun_t_cloacal"
## [1] 34747
## [1] "removing: surg_bun_t_closed"
## [1] 34740
## [1] "removing: surg_bun_t_closure"
## [1] 34740
## [1] "removing: surg_bun_t_clot"
## [1] 34740
## [1] "removing: surg_bun_t_collar"
## [1] 34740
## [1] "removing: surg_bun_t_colon"
## [1] 34733
## [1] "removing: surg_bun_t_colonoscopy"
## [1] 34733
## [1] "removing: surg_bun_t_coloproctostomy"
## [1] 34733
## [1] "removing: surg_bun_t_colostomy"
## [1] 34733
## [1] "removing: surg_bun_t_complex"
## [1] 34733
## [1] "removing: surg_bun_t_compression"
## [1] 34733
## [1] "removing: surg_bun_t_contriction"
## [1] 34733
## [1] "removing: surg_bun_t_control"
## [1] 34732
## [1] "removing: surg_bun_t_cord"
## [1] 34732
## [1] "removing: surg_bun_t_cornea"
## [1] 34732
## [1] "removing: surg_bun_t_coronary"
## [1] 34732
## [1] "removing: surg_bun_t_coronoid"
## [1] 34732
## [1] "removing: surg_bun_t_correct"
## [1] 34732
## [1] "removing: surg_bun_t_cranial"
## [1] 34732
## [1] "removing: surg_bun_t_craniofacial"
## [1] 34732
## [1] "removing: surg_bun_t_cryotherapy"
## [1] 34732
## [1] "removing: surg_bun_t_cuff"
## [1] 34732
## [1] "removing: surg_bun_t_cult"
## [1] 34732
## [1] "removing: surg_bun_t_cyst"
## [1] 34725
## [1] "removing: surg_bun_t_cysto"
## [1] 34725
## [1] "removing: surg_bun_t_cystometrogram"
## [1] 34725
## [1] "removing: surg_bun_t_cystoscopy"
## [1] 34725
## [1] "removing: surg_bun_t_cystouretero"
## [1] 34725
## [1] "removing: surg_bun_t_debrid"
## [1] 34718
## [1] "removing: surg_bun_t_debride"
## [1] 34718
## [1] "removing: surg_bun_t_decompression"
## [1] 34718
## [1] "removing: surg_bun_t_deep"
## [1] 34718
## [1] "removing: surg_bun_t_defect"
## [1] 34718
## [1] "removing: surg_bun_t_defib"
## [1] 34718
## [1] "removing: surg_bun_t_deformity"
## [1] 34718
## [1] "removing: surg_bun_t_degenerated"
## [1] 34718
## [1] "removing: surg_bun_t_delivery"
## [1] 34718
## [1] "removing: surg_bun_t_deplet"
## [1] 34718
## [1] "removing: surg_bun_t_dermabrasion"
## [1] 34718
## [1] "removing: surg_bun_t_design"
## [1] 34718
## [1] "removing: surg_bun_t_destroy"
## [1] 34693
## [1] "removing: surg_bun_t_destruction"
## [1] 34688
## [1] "removing: surg_bun_t_dev"
## [1] 34646
## [1] "removing: surg_bun_t_device"
## [1] 34646
## [1] "removing: surg_bun_t_diagnostic"
## [1] 34602
## [1] "removing: surg_bun_t_dialysis"
## [1] 34602
## [1] "removing: surg_bun_t_diaphragm"
## [1] 34602
## [1] "removing: surg_bun_t_digit"
## [1] 34602
## [1] "removing: surg_bun_t_dilat"
## [1] 34601
## [1] "removing: surg_bun_t_dilation"
## [1] 34601
## [1] "removing: surg_bun_t_disc"
## [1] 34589
## [1] "removing: surg_bun_t_diskectomy"
## [1] 34589
## [1] "removing: surg_bun_t_dislocat"
## [1] 34582
## [1] "removing: surg_bun_t_dissect"
## [1] 34582
## [1] "removing: surg_bun_t_donar"
## [1] 34582
## [1] "removing: surg_bun_t_donor"
## [1] 34582
## [1] "removing: surg_bun_t_drain"
## [1] 34578
## [1] "removing: surg_bun_t_draw"
## [1] 34578
## [1] "removing: surg_bun_t_dressing"
## [1] 34578
## [1] "removing: surg_bun_t_drill"
## [1] 34578
## [1] "removing: surg_bun_t_drug"
## [1] 34575
## [1] "removing: surg_bun_t_duct"
## [1] 34574
## [1] "removing: surg_bun_t_duoden"
## [1] 34574
## [1] "removing: surg_bun_t_dura"
## [1] 34533
## [1] "removing: surg_bun_t_ear"
## [1] 34525
## [1] "removing: surg_bun_t_eardrum"
## [1] 34523
## [1] "removing: surg_bun_t_ectopic"
## [1] 34523
## [1] "removing: surg_bun_t_egd"
## [1] 34523
## [1] "removing: surg_bun_t_elbow"
## [1] 34523
## [1] "removing: surg_bun_t_electrical"
## [1] 34523
## [1] "removing: surg_bun_t_electrode"
## [1] 34523
## [1] "removing: surg_bun_t_electroejaculation"
## [1] 34523
## [1] "removing: surg_bun_t_electrolysis"
## [1] 34523
## [1] "removing: surg_bun_t_embolize"
## [1] 34523
## [1] "removing: surg_bun_t_embryo"
## [1] 34523
## [1] "removing: surg_bun_t_endo"
## [1] 34470
## [1] "removing: surg_bun_t_endoscopy"
## [1] 34470
## [1] "removing: surg_bun_t_endovas"
## [1] 34470
## [1] "removing: surg_bun_t_enlarge"
## [1] 34470
## [1] "removing: surg_bun_t_enterectomy"
## [1] 34470
## [1] "removing: surg_bun_t_epiderm"
## [1] 34470
## [1] "removing: surg_bun_t_epididymis"
## [1] 34470
## [1] "removing: surg_bun_t_epiphysis"
## [1] 34470
## [1] "removing: surg_bun_t_episiotomy"
## [1] 34470
## [1] "removing: surg_bun_t_esoph"
## [1] 34469
## [1] "removing: surg_bun_t_esophagoscopy"
## [1] 34469
## [1] "removing: surg_bun_t_esophagus"
## [1] 34469
## [1] "removing: surg_bun_t_exam"
## [1] 34469
## [1] "removing: surg_bun_t_exc"
## [1] 34466
## [1] "removing: surg_bun_t_excessive"
## [1] 34466
## [1] "removing: surg_bun_t_excise"
## [1] 34466
## [1] "removing: surg_bun_t_excision"
## [1] 34466
## [1] "removing: surg_bun_t_exostosis"
## [1] 34466
## [1] "removing: surg_bun_t_expander"
## [1] 34466
## [1] "removing: surg_bun_t_explora"
## [1] 34466
## [1] "removing: surg_bun_t_explore"
## [1] 34466
## [1] "removing: surg_bun_t_extra"
## [1] 34466
## [1] "removing: surg_bun_t_extremity"
## [1] 34466
## [1] "removing: surg_bun_t_eye"
## [1] 34466
## [1] "removing: surg_bun_t_eyelashes"
## [1] 34466
## [1] "removing: surg_bun_t_eyelid"
## [1] 34466
## [1] "removing: surg_bun_t_face"
## [1] 34448
## [1] "removing: surg_bun_t_fallopian"
## [1] 34448
## [1] "removing: surg_bun_t_fascia"
## [1] 34448
## [1] "removing: surg_bun_t_fat"
## [1] 34448
## [1] "removing: surg_bun_t_femur"
## [1] 34448
## [1] "removing: surg_bun_t_fibula"
## [1] 34447
## [1] "removing: surg_bun_t_finger"
## [1] 34446
## [1] "removing: surg_bun_t_fissure"
## [1] 34446
## [1] "removing: surg_bun_t_fistula"
## [1] 34446
## [1] "removing: surg_bun_t_fixation"
## [1] 34445
## [1] "removing: surg_bun_t_flap"
## [1] 34445
## [1] "removing: surg_bun_t_flow"
## [1] 34445
## [1] "removing: surg_bun_t_fluid"
## [1] 34445
## [1] "removing: surg_bun_t_fna"
## [1] 34445
## [1] "removing: surg_bun_t_follow-up"
## [1] 34445
## [1] "removing: surg_bun_t_foot"
## [1] 34444
## [1] "removing: surg_bun_t_forearm"
## [1] 34444
## [1] "removing: surg_bun_t_forehead"
## [1] 34444
## [1] "removing: surg_bun_t_foreign"
## [1] 34444
## [1] "removing: surg_bun_t_fracture"
## [1] 34429
## [1] "removing: surg_bun_t_free"
## [1] 34429
## [1] "removing: surg_bun_t_frontal"
## [1] 34429
## [1] "removing: surg_bun_t_fundoplasty"
## [1] 34429
## [1] "removing: surg_bun_t_fusion"
## [1] 34427
## [1] "removing: surg_bun_t_gallbladder"
## [1] 34427
## [1] "removing: surg_bun_t_gastrointestinal"
## [1] 34426
## [1] "removing: surg_bun_t_gastrostomy"
## [1] 34426
## [1] "removing: surg_bun_t_gentialia"
## [1] 34426
## [1] "removing: surg_bun_t_gi"
## [1] 34426
## [1] "removing: surg_bun_t_gland"
## [1] 34426
## [1] "removing: surg_bun_t_glaucoma"
## [1] 34426
## [1] "removing: surg_bun_t_graft"
## [1] 34426
## [1] "removing: surg_bun_t_growth"
## [1] 34426
## [1] "removing: surg_bun_t_gum"
## [1] 34426
## [1] "removing: surg_bun_t_hair"
## [1] 34426
## [1] "removing: surg_bun_t_hallux"
## [1] 34426
## [1] "removing: surg_bun_t_halo"
## [1] 34426
## [1] "removing: surg_bun_t_hand"
## [1] 34426
## [1] "removing: surg_bun_t_harvest"
## [1] 34426
## [1] "removing: surg_bun_t_head"
## [1] 34426
## [1] "removing: surg_bun_t_heart"
## [1] 34426
## [1] "removing: surg_bun_t_hematoma"
## [1] 34426
## [1] "removing: surg_bun_t_hemorrhoid"
## [1] 34424
## [1] "removing: surg_bun_t_hemorrhoidopexy"
## [1] 34424
## [1] "removing: surg_bun_t_hernia"
## [1] 34424
## [1] "removing: surg_bun_t_hip"
## [1] 34422
## [1] "removing: surg_bun_t_hole"
## [1] 34422
## [1] "removing: surg_bun_t_hormone"
## [1] 34422
## [1] "removing: surg_bun_t_humer"
## [1] 34422
## [1] "removing: surg_bun_t_hydrocele"
## [1] 34422
## [1] "removing: surg_bun_t_hyst"
## [1] 34422
## [1] "removing: surg_bun_t_hysterectomy"
## [1] 34422
## [1] "removing: surg_bun_t_iliopsoas"
## [1] 34422
## [1] "removing: surg_bun_t_illiac"
## [1] 34422
## [1] "removing: surg_bun_t_imag"
## [1] 34421
## [1] "removing: surg_bun_t_implant"
## [1] 34421
## [1] "removing: surg_bun_t_incision"
## [1] 34421
## [1] "removing: surg_bun_t_infect"
## [1] 34421
## [1] "removing: surg_bun_t_injection"
## [1] 34398
## [1] "removing: surg_bun_t_inner"
## [1] 34398
## [1] "removing: surg_bun_t_insert"
## [1] 34398
## [1] "removing: surg_bun_t_interlaminar"
## [1] 34302
## [1] "removing: surg_bun_t_internal"
## [1] 34302
## [1] "removing: surg_bun_t_intestine"
## [1] 34302
## [1] "removing: surg_bun_t_intracorporeal"
## [1] 34302
## [1] "removing: surg_bun_t_intracranial"
## [1] 34302
## [1] "removing: surg_bun_t_intranasal"
## [1] 34302
## [1] "removing: surg_bun_t_inverted"
## [1] 34302
## [1] "removing: surg_bun_t_iris"
## [1] 34302
## [1] "removing: surg_bun_t_jaw"
## [1] 34302
## [1] "removing: surg_bun_t_joint"
## [1] 34282
## [1] "removing: surg_bun_t_kidney"
## [1] 34282
## [1] "removing: surg_bun_t_knee"
## [1] 34282
## [1] "removing: surg_bun_t_kneecap"
## [1] 34282
## [1] "removing: surg_bun_t_knuckle"
## [1] 34282
## [1] "removing: surg_bun_t_kyphectomy"
## [1] 34282
## [1] "removing: surg_bun_t_laceration"
## [1] 34281
## [1] "removing: surg_bun_t_lamina"
## [1] 34277
## [1] "removing: surg_bun_t_laminotomy"
## [1] 34276
## [1] "removing: surg_bun_t_lap"
## [1] 34276
## [1] "removing: surg_bun_t_laparoscop"
## [1] 34276
## [1] "removing: surg_bun_t_large"
## [1] 34276
## [1] "removing: surg_bun_t_laryngoscopy"
## [1] 34274
## [1] "removing: surg_bun_t_larynx"
## [1] 34274
## [1] "removing: surg_bun_t_laser"
## [1] 34272
## [1] "removing: surg_bun_t_lavage"
## [1] 34272
## [1] "removing: surg_bun_t_lefort"
## [1] 34272
## [1] "removing: surg_bun_t_leg"
## [1] 34272
## [1] "removing: surg_bun_t_lengthen"
## [1] 34272
## [1] "removing: surg_bun_t_lesion"
## [1] 34226
## [1] "removing: surg_bun_t_ligament"
## [1] 34226
## [1] "removing: surg_bun_t_ligation"
## [1] 34226
## [1] "removing: surg_bun_t_lining"
## [1] 34226
## [1] "removing: surg_bun_t_lip"
## [1] 34226
## [1] "removing: surg_bun_t_lipectomy"
## [1] 34226
## [1] "removing: surg_bun_t_liver"
## [1] 34226
## [1] "removing: surg_bun_t_lobectomy"
## [1] 34226
## [1] "removing: surg_bun_t_long"
## [1] 34226
## [1] "removing: surg_bun_t_lower"
## [1] 34226
## [1] "removing: surg_bun_t_lumbar"
## [1] 34225
## [1] "removing: surg_bun_t_lung"
## [1] 34225
## [1] "removing: surg_bun_t_lymph"
## [1] 34225
## [1] "removing: surg_bun_t_major"
## [1] 34225
## [1] "removing: surg_bun_t_mandible"
## [1] 34225
## [1] "removing: surg_bun_t_marrow"
## [1] 34225
## [1] "removing: surg_bun_t_mast"
## [1] 34225
## [1] "removing: surg_bun_t_mastectomy"
## [1] 34225
## [1] "removing: surg_bun_t_mastoid"
## [1] 34225
## [1] "removing: surg_bun_t_material"
## [1] 34225
## [1] "removing: surg_bun_t_maxilla"
## [1] 34225
## [1] "removing: surg_bun_t_mcp"
## [1] 34225
## [1] "removing: surg_bun_t_mesh"
## [1] 34225
## [1] "removing: surg_bun_t_metacarpal"
## [1] 34225
## [1] "removing: surg_bun_t_microvasc"
## [1] 34225
## [1] "removing: surg_bun_t_middle"
## [1] 34225
## [1] "removing: surg_bun_t_midfoot"
## [1] 34225
## [1] "removing: surg_bun_t_mitral"
## [1] 34225
## [1] "removing: surg_bun_t_mohs"
## [1] 34224
## [1] "removing: surg_bun_t_mouth"
## [1] 34224
## [1] "removing: surg_bun_t_musc"
## [1] 34021
## [1] "removing: surg_bun_t_myelography"
## [1] 34021
## [1] "removing: surg_bun_t_myomectomy"
## [1] 34021
## [1] "removing: surg_bun_t_nail"
## [1] 34017
## [1] "removing: surg_bun_t_nasal"
## [1] 34016
## [1] "removing: surg_bun_t_nasolacrimal"
## [1] 34016
## [1] "removing: surg_bun_t_neck"
## [1] 34015
## [1] "removing: surg_bun_t_needle"
## [1] 34015
## [1] "removing: surg_bun_t_nephrectomy"
## [1] 34015
## [1] "removing: surg_bun_t_nerve"
## [1] 34015
## [1] "removing: surg_bun_t_neuroelectrode"
## [1] 34015
## [1] "removing: surg_bun_t_neuroendoscopy"
## [1] 34015
## [1] "removing: surg_bun_t_nipple"
## [1] 34015
## [1] "removing: surg_bun_t_njx"
## [1] 34015
## [1] "removing: surg_bun_t_node"
## [1] 34015
## [1] "removing: surg_bun_t_nodule"
## [1] 34015
## [1] "removing: surg_bun_t_nonfacial"
## [1] 34015
## [1] "removing: surg_bun_t_nose"
## [1] 34015
## [1] "removing: surg_bun_t_nosebleed"
## [1] 34015
## [1] "removing: surg_bun_t_nuero"
## [1] 34015
## [1] "removing: surg_bun_t_obstetrical"
## [1] 34015
## [1] "removing: surg_bun_t_open"
## [1] 34015
## [1] "removing: surg_bun_t_opening"
## [1] 34015
## [1] "removing: surg_bun_t_oral"
## [1] 34015
## [1] "removing: surg_bun_t_orbitocranial"
## [1] 34015
## [1] "removing: surg_bun_t_orchiopexy"
## [1] 34015
## [1] "removing: surg_bun_t_other"
## [1] 34015
## [1] "removing: surg_bun_t_ovarian"
## [1] 34015
## [1] "removing: surg_bun_t_ovary"
## [1] 34015
## [1] "removing: surg_bun_t_pacemaker"
## [1] 34015
## [1] "removing: surg_bun_t_pad"
## [1] 34015
## [1] "removing: surg_bun_t_palate"
## [1] 34015
## [1] "removing: surg_bun_t_palm"
## [1] 34015
## [1] "removing: surg_bun_t_palsy"
## [1] 34015
## [1] "removing: surg_bun_t_pancreas"
## [1] 34015
## [1] "removing: surg_bun_t_pancreatectomy"
## [1] 34015
## [1] "removing: surg_bun_t_paravert"
## [1] 34002
## [1] "removing: surg_bun_t_part"
## [1] 34002
## [1] "removing: surg_bun_t_partial"
## [1] 34002
## [1] "removing: surg_bun_t_pedcle"
## [1] 34002
## [1] "removing: surg_bun_t_pedicle"
## [1] 34002
## [1] "removing: surg_bun_t_peel"
## [1] 34002
## [1] "removing: surg_bun_t_pelvic"
## [1] 34001
## [1] "removing: surg_bun_t_pelvis"
## [1] 34001
## [1] "removing: surg_bun_t_penile"
## [1] 34001
## [1] "removing: surg_bun_t_penis"
## [1] 34001
## [1] "removing: surg_bun_t_percut"
## [1] 34001
## [1] "removing: surg_bun_t_percutaneous"
## [1] 34001
## [1] "removing: surg_bun_t_perineum"
## [1] 34001
## [1] "removing: surg_bun_t_pernineum"
## [1] 34001
## [1] "removing: surg_bun_t_pharynx"
## [1] 34001
## [1] "removing: surg_bun_t_pierce"
## [1] 34001
## [1] "removing: surg_bun_t_pilonidal"
## [1] 34001
## [1] "removing: surg_bun_t_place"
## [1] 34001
## [1] "removing: surg_bun_t_plastic"
## [1] 34001
## [1] "removing: surg_bun_t_plate"
## [1] 34001
## [1] "removing: surg_bun_t_pneumonectomy"
## [1] 34001
## [1] "removing: surg_bun_t_polyp"
## [1] 34001
## [1] "removing: surg_bun_t_pouch"
## [1] 34001
## [1] "removing: surg_bun_t_pregnancy"
## [1] 34001
## [1] "removing: surg_bun_t_prep"
## [1] 34001
## [1] "removing: surg_bun_t_prepare"
## [1] 34001
## [1] "removing: surg_bun_t_pressure"
## [1] 34001
## [1] "removing: surg_bun_t_process"
## [1] 34001
## [1] "removing: surg_bun_t_proctosigmoidoscopy"
## [1] 34001
## [1] "removing: surg_bun_t_prolapse"
## [1] 34001
## [1] "removing: surg_bun_t_prostat"
## [1] 34001
## [1] "removing: surg_bun_t_prostate"
## [1] 34001
## [1] "removing: surg_bun_t_prosthesis"
## [1] 34001
## [1] "removing: surg_bun_t_prosthetic"
## [1] 34001
## [1] "removing: surg_bun_t_pulmonary"
## [1] 34001
## [1] "removing: surg_bun_t_punch"
## [1] 34001
## [1] "removing: surg_bun_t_puncture"
## [1] 34001
## [1] "removing: surg_bun_t_pyeloplasty"
## [1] 34001
## [1] "removing: surg_bun_t_radius"
## [1] 34001
## [1] "removing: surg_bun_t_realignment"
## [1] 34001
## [1] "removing: surg_bun_t_rechannel"
## [1] 34001
## [1] "removing: surg_bun_t_reconstruct"
## [1] 34001
## [1] "removing: surg_bun_t_reconstruction"
## [1] 34001
## [1] "removing: surg_bun_t_rectal"
## [1] 34001
## [1] "removing: surg_bun_t_rectum"
## [1] 34001
## [1] "removing: surg_bun_t_reduction"
## [1] 34001
## [1] "removing: surg_bun_t_reinforce"
## [1] 34001
## [1] "removing: surg_bun_t_release"
## [1] 34001
## [1] "removing: surg_bun_t_remov"
## [1] 33997
## [1] "removing: surg_bun_t_removal"
## [1] 33997
## [1] "removing: surg_bun_t_remove"
## [1] 33997
## [1] "removing: surg_bun_t_renal"
## [1] 33997
## [1] "removing: surg_bun_t_repair"
## [1] 33997
## [1] "removing: surg_bun_t_replantation"
## [1] 33997
## [1] "removing: surg_bun_t_resect"
## [1] 33996
## [1] "removing: surg_bun_t_reset"
## [1] 33996
## [1] "removing: surg_bun_t_retina"
## [1] 33996
## [1] "removing: surg_bun_t_retinopathy"
## [1] 33996
## [1] "removing: surg_bun_t_revise"
## [1] 33996
## [1] "removing: surg_bun_t_revision"
## [1] 33996
## [1] "removing: surg_bun_t_rib"
## [1] 33996
## [1] "removing: surg_bun_t_ridge"
## [1] 33996
## [1] "removing: surg_bun_t_ring"
## [1] 33996
## [1] "removing: surg_bun_t_roof"
## [1] 33996
## [1] "removing: surg_bun_t_rotator"
## [1] 33996
## [1] "removing: surg_bun_t_routine"
## [1] 33996
## [1] "removing: surg_bun_t_sacrum"
## [1] 33996
## [1] "removing: surg_bun_t_salivary"
## [1] 33996
## [1] "removing: surg_bun_t_scapula"
## [1] 33996
## [1] "removing: surg_bun_t_scar"
## [1] 33996
## [1] "removing: surg_bun_t_scope"
## [1] 33995
## [1] "removing: surg_bun_t_scrotum"
## [1] 33995
## [1] "removing: surg_bun_t_segment"
## [1] 33995
## [1] "removing: surg_bun_t_segmental"
## [1] 33995
## [1] "removing: surg_bun_t_segmentectomy"
## [1] 33995
## [1] "removing: surg_bun_t_septum"
## [1] 33995
## [1] "removing: surg_bun_t_sex"
## [1] 33995
## [1] "removing: surg_bun_t_shave"
## [1] 33995
## [1] "removing: surg_bun_t_sheath"
## [1] 33995
## [1] "removing: surg_bun_t_shorten"
## [1] 33995
## [1] "removing: surg_bun_t_shoulder"
## [1] 33994
## [1] "removing: surg_bun_t_shunt"
## [1] 33993
## [1] "removing: surg_bun_t_sigmoidoscopy"
## [1] 33993
## [1] "removing: surg_bun_t_simple"
## [1] 33993
## [1] "removing: surg_bun_t_sinus"
## [1] 33993
## [1] "removing: surg_bun_t_skeletal"
## [1] 33993
## [1] "removing: surg_bun_t_skin"
## [1] 33993
## [1] "removing: surg_bun_t_skull"
## [1] 33993
## [1] "removing: surg_bun_t_sleeve"
## [1] 33993
## [1] "removing: surg_bun_t_slipped"
## [1] 33993
## [1] "removing: surg_bun_t_small"
## [1] 33993
## [1] "removing: surg_bun_t_socket"
## [1] 33993
## [1] "removing: surg_bun_t_soft"
## [1] 33993
## [1] "removing: surg_bun_t_sore"
## [1] 33993
## [1] "removing: surg_bun_t_sperm"
## [1] 33993
## [1] "removing: surg_bun_t_sphenoid"
## [1] 33993
## [1] "removing: surg_bun_t_sphincter"
## [1] 33993
## [1] "removing: surg_bun_t_spine"
## [1] 33992
## [1] "removing: surg_bun_t_spleen"
## [1] 33992
## [1] "removing: surg_bun_t_splint"
## [1] 33992
## [1] "removing: surg_bun_t_split"
## [1] 33992
## [1] "removing: surg_bun_t_stem"
## [1] 33992
## [1] "removing: surg_bun_t_stenosis"
## [1] 33992
## [1] "removing: surg_bun_t_stent"
## [1] 33992
## [1] "removing: surg_bun_t_sternal"
## [1] 33992
## [1] "removing: surg_bun_t_sternum"
## [1] 33992
## [1] "removing: surg_bun_t_stimulation"
## [1] 33992
## [1] "removing: surg_bun_t_stomach"
## [1] 33992
## [1] "removing: surg_bun_t_stone"
## [1] 33992
## [1] "removing: surg_bun_t_stop"
## [1] 33992
## [1] "removing: surg_bun_t_strapping"
## [1] 33992
## [1] "removing: surg_bun_t_stricture"
## [1] 33992
## [1] "removing: surg_bun_t_styloid"
## [1] 33992
## [1] "removing: surg_bun_t_sub"
## [1] 33984
## [1] "removing: surg_bun_t_suction"
## [1] 33984
## [1] "removing: surg_bun_t_superficial"
## [1] 33984
## [1] "removing: surg_bun_t_suspension"
## [1] 33984
## [1] "removing: surg_bun_t_sutures"
## [1] 33984
## [1] "removing: surg_bun_t_sweat"
## [1] 33984
## [1] "removing: surg_bun_t_tags"
## [1] 33984
## [1] "removing: surg_bun_t_tail"
## [1] 33984
## [1] "removing: surg_bun_t_tear"
## [1] 33984
## [1] "removing: surg_bun_t_temporal"
## [1] 33984
## [1] "removing: surg_bun_t_tendon"
## [1] 33984
## [1] "removing: surg_bun_t_testis"
## [1] 33984
## [1] "removing: surg_bun_t_therapy"
## [1] 33984
## [1] "removing: surg_bun_t_thigh"
## [1] 33984
## [1] "removing: surg_bun_t_thoracoscopy"
## [1] 33984
## [1] "removing: surg_bun_t_thoracotomy"
## [1] 33984
## [1] "removing: surg_bun_t_thorax"
## [1] 33984
## [1] "removing: surg_bun_t_throat"
## [1] 33984
## [1] "removing: surg_bun_t_thrombin"
## [1] 33984
## [1] "removing: surg_bun_t_thrombo"
## [1] 33984
## [1] "removing: surg_bun_t_thumb"
## [1] 33984
## [1] "removing: surg_bun_t_thymus"
## [1] 33984
## [1] "removing: surg_bun_t_thyroid"
## [1] 33984
## [1] "removing: surg_bun_t_tip"
## [1] 33984
## [1] "removing: surg_bun_t_tiss"
## [1] 33983
## [1] "removing: surg_bun_t_tkr"
## [1] 33983
## [1] "removing: surg_bun_t_toe"
## [1] 33983
## [1] "removing: surg_bun_t_tongue"
## [1] 33983
## [1] "removing: surg_bun_t_tonsil"
## [1] 33983
## [1] "removing: surg_bun_t_torsion"
## [1] 33983
## [1] "removing: surg_bun_t_total"
## [1] 33983
## [1] "removing: surg_bun_t_tract"
## [1] 33983
## [1] "removing: surg_bun_t_transcath"
## [1] 33983
## [1] "removing: surg_bun_t_transfer"
## [1] 33983
## [1] "removing: surg_bun_t_transform"
## [1] 33983
## [1] "removing: surg_bun_t_transplant"
## [1] 33983
## [1] "removing: surg_bun_t_treat"
## [1] 33983
## [1] "removing: surg_bun_t_triceps"
## [1] 33983
## [1] "removing: surg_bun_t_tricuspid"
## [1] 33983
## [1] "removing: surg_bun_t_trigger"
## [1] 33846
## [1] "removing: surg_bun_t_trim"
## [1] 33846
## [1] "removing: surg_bun_t_tube"
## [1] 33846
## [1] "removing: surg_bun_t_tum"
## [1] 33846
## [1] "removing: surg_bun_t_tumor"
## [1] 33846
## [1] "removing: surg_bun_t_tunneled"
## [1] 33846
## [1] "removing: surg_bun_t_ulcer"
## [1] 33846
## [1] "removing: surg_bun_t_ulna"
## [1] 33846
## [1] "removing: surg_bun_t_umbil"
## [1] 33846
## [1] "removing: surg_bun_t_unit"
## [1] 33846
## [1] "removing: surg_bun_t_unstable"
## [1] 33846
## [1] "removing: surg_bun_t_upper"
## [1] 33846
## [1] "removing: surg_bun_t_ureter"
## [1] 33846
## [1] "removing: surg_bun_t_urethra"
## [1] 33846
## [1] "removing: surg_bun_t_urethrovaginal"
## [1] 33846
## [1] "removing: surg_bun_t_uterus"
## [1] 33846
## [1] "removing: surg_bun_t_uvula"
## [1] 33846
## [1] "removing: surg_bun_t_vag"
## [1] 33846
## [1] "removing: surg_bun_t_vagotomy"
## [1] 33846
## [1] "removing: surg_bun_t_valgus"
## [1] 33846
## [1] "removing: surg_bun_t_valve"
## [1] 33846
## [1] "removing: surg_bun_t_vascular"
## [1] 33846
## [1] "removing: surg_bun_t_vein"
## [1] 33846
## [1] "removing: surg_bun_t_venipuncture"
## [1] 33846
## [1] "removing: surg_bun_t_ventricular"
## [1] 33846
## [1] "removing: surg_bun_t_vertebra"
## [1] 33846
## [1] "removing: surg_bun_t_vessel"
## [1] 33846
## [1] "removing: surg_bun_t_vulva"
## [1] 33846
## [1] "removing: surg_bun_t_wall"
## [1] 33846
## [1] "removing: surg_bun_t_wax"
## [1] 33846
## [1] "removing: surg_bun_t_web"
## [1] 33846
## [1] "removing: surg_bun_t_wedge"
## [1] 33846
## [1] "removing: surg_bun_t_window"
## [1] 33846
## [1] "removing: surg_bun_t_windpipe"
## [1] 33846
## [1] "removing: surg_bun_t_wound"
## [1] 33846
## [1] "removing: surg_bun_t_wrinkle"
## [1] 33846
## [1] "removing: surg_bun_t_wrist"
## [1] 33846
## [1] "removing: surg_bun_t_x-ray"
## [1] 33846
## [1] "removing: surg_bun_t_xenograft"
## [1] 33846
## [1] "removing: surg_bun_t_xenogrft"
## [1] 33846
## [1] "removing: surg_bun_t_xgraft"
## [1] 33846
## [1] "removing: medi_bun_t_add-on"
## [1] 33710
## [1] "removing: medi_bun_t_analy"
## [1] 33699
## [1] "removing: medi_bun_t_applicat"
## [1] 33699
## [1] "removing: medi_bun_t_apply"
## [1] 33699
## [1] "removing: medi_bun_t_asses"
## [1] 33638
## [1] "removing: medi_bun_t_benzathine"
## [1] 33636
## [1] "removing: medi_bun_t_betamethasone"
## [1] 33623
## [1] "removing: medi_bun_t_boot"
## [1] 33597
## [1] "removing: medi_bun_t_candida"
## [1] 33597
## [1] "removing: medi_bun_t_cast"
## [1] 33597
## [1] "removing: medi_bun_t_cathet"
## [1] 33595
## [1] "removing: medi_bun_t_catheterisation"
## [1] 33595
## [1] "removing: medi_bun_t_choline"
## [1] 33592
## [1] "removing: medi_bun_t_clinic"
## [1] 33556
## [1] "removing: medi_bun_t_clostridium"
## [1] 33556
## [1] "removing: medi_bun_t_cochlear"
## [1] 33556
## [1] "removing: medi_bun_t_cold"
## [1] 33518
## [1] "removing: medi_bun_t_consult"
## [1] 33518
## [1] "removing: medi_bun_t_consultation"
## [1] 33518
## [1] "removing: medi_bun_t_crutch"
## [1] 33518
## [1] "removing: medi_bun_t_cyst"
## [1] 33518
## [1] "removing: medi_bun_t_depression"
## [1] 33384
## [1] "removing: medi_bun_t_derma"
## [1] 33384
## [1] "removing: medi_bun_t_dermal"
## [1] 33384
## [1] "removing: medi_bun_t_device"
## [1] 33362
## [1] "removing: medi_bun_t_dialys"
## [1] 33357
## [1] "removing: medi_bun_t_discharg"
## [1] 33357
## [1] "removing: medi_bun_t_drain"
## [1] 33357
## [1] "removing: medi_bun_t_draw"
## [1] 33343
## [1] "removing: medi_bun_t_dtap"
## [1] 33338
## [1] "removing: medi_bun_t_dtap-ipv"
## [1] 33338
## [1] "removing: medi_bun_t_emotional"
## [1] 33138
## [1] "removing: medi_bun_t_endocrine"
## [1] 33138
## [1] "removing: medi_bun_t_epstein-barr"
## [1] 33138
## [1] "removing: medi_bun_t_equipment"
## [1] 33138
## [1] "removing: medi_bun_t_establish"
## [1] 33108
## [1] "removing: medi_bun_t_estradiol"
## [1] 33107
## [1] "removing: medi_bun_t_evaluat"
## [1] 33070
## [1] "removing: medi_bun_t_flu"
## [1] 33022
## [1] "removing: medi_bun_t_follow-up"
## [1] 33022
## [1] "removing: medi_bun_t_gastroenterology"
## [1] 33022
## [1] "removing: medi_bun_t_gonorrhoeae"
## [1] 33022
## [1] "removing: medi_bun_t_haloperidol"
## [1] 33019
## [1] "removing: medi_bun_t_hearing"
## [1] 32997
## [1] "removing: medi_bun_t_hemorrhoid"
## [1] 32997
## [1] "removing: medi_bun_t_hepatitis"
## [1] 32995
## [1] "removing: medi_bun_t_hiv - 2"
## [1] 32995
## [1] "removing: medi_bun_t_hiv - 3"
## [1] 32995
## [1] "removing: medi_bun_t_home"
## [1] 32995
## [1] "removing: medi_bun_t_hormone"
## [1] 32995
## [1] "removing: medi_bun_t_hpv"
## [1] 32953
## [1] "removing: medi_bun_t_hydrate"
## [1] 32953
## [1] "removing: medi_bun_t_hydration"
## [1] 32922
## [1] "removing: medi_bun_t_iiv"
## [1] 32701
## [1] "removing: medi_bun_t_immun"
## [1] 32395
## [1] "removing: medi_bun_t_immune"
## [1] 32395
## [1] "removing: medi_bun_t_immunization"
## [1] 32395
## [1] "removing: medi_bun_t_immunoassay"
## [1] 32395
## [1] "removing: medi_bun_t_immunohisto"
## [1] 32395
## [1] "removing: medi_bun_t_immunotherapy"
## [1] 32395
## [1] "removing: medi_bun_t_influenza"
## [1] 32395
## [1] "removing: medi_bun_t_insulin"
## [1] 32393
## [1] "removing: medi_bun_t_intrauterine"
## [1] 32393
## [1] "removing: medi_bun_t_iodine"
## [1] 32393
## [1] "removing: medi_bun_t_laiv4"
## [1] 32393
## [1] "removing: medi_bun_t_manual"
## [1] 31641
## [1] "removing: medi_bun_t_medroxyprogesterone"
## [1] 31635
## [1] "removing: medi_bun_t_mental"
## [1] 31634
## [1] "removing: medi_bun_t_metabolic"
## [1] 31634
## [1] "removing: medi_bun_t_methylprednisolone"
## [1] 31603
## [1] "removing: medi_bun_t_mgmt"
## [1] 31598
## [1] "removing: medi_bun_t_mmr"
## [1] 31598
## [1] "removing: medi_bun_t_monitor"
## [1] 31592
## [1] "removing: medi_bun_t_muscle"
## [1] 31590
## [1] "removing: medi_bun_t_needle"
## [1] 31588
## [1] "removing: medi_bun_t_non-invasive"
## [1] 31588
## [1] "removing: medi_bun_t_non-stress"
## [1] 31588
## [1] "removing: medi_bun_t_normal"
## [1] 31570
## [1] "removing: medi_bun_t_nursing"
## [1] 31570
## [1] "removing: medi_bun_t_nutrition"
## [1] 31565
## [1] "removing: medi_bun_t_observ"
## [1] 31558
## [1] "removing: medi_bun_t_oral"
## [1] 31492
## [1] "removing: medi_bun_t_outpatient"
## [1] 31492
## [1] "removing: medi_bun_t_outpt"
## [1] 31492
## [1] "removing: medi_bun_t_pap"
## [1] 31485
## [1] "removing: medi_bun_t_physical"
## [1] 31475
## [1] "removing: medi_bun_t_portable"
## [1] 31475
## [1] "removing: medi_bun_t_prednisolone"
## [1] 31475
## [1] "removing: medi_bun_t_pressure"
## [1] 31475
## [1] "removing: medi_bun_t_progesterone"
## [1] 31475
## [1] "removing: medi_bun_t_psych"
## [1] 31464
## [1] "removing: medi_bun_t_psychiat"
## [1] 31464
## [1] "removing: medi_bun_t_psycho"
## [1] 31464
## [1] "removing: medi_bun_t_pulmonary"
## [1] 31464
## [1] "removing: medi_bun_t_rabies"
## [1] 31464
## [1] "removing: medi_bun_t_radiation"
## [1] 31464
## [1] "removing: medi_bun_t_radio"
## [1] 31463
## [1] "removing: medi_bun_t_reeducation"
## [1] 31355
## [1] "removing: medi_bun_t_reeval"
## [1] 31355
## [1] "removing: medi_bun_t_rehab"
## [1] 31355
## [1] "removing: medi_bun_t_remove"
## [1] 31355
## [1] "removing: medi_bun_t_rheumatoid"
## [1] 31355
## [1] "removing: medi_bun_t_riv4"
## [1] 31355
## [1] "removing: medi_bun_t_rubella"
## [1] 31355
## [1] "removing: medi_bun_t_saline"
## [1] 31355
## [1] "removing: medi_bun_t_sciatic"
## [1] 31355
## [1] "removing: medi_bun_t_serum"
## [1] 31355
## [1] "removing: medi_bun_t_services"
## [1] 31351
## [1] "removing: medi_bun_t_shave"
## [1] 31351
## [1] "removing: medi_bun_t_simple"
## [1] 31351
## [1] "removing: medi_bun_t_sleep"
## [1] 31347
## [1] "removing: medi_bun_t_smear"
## [1] 31347
## [1] "removing: medi_bun_t_sodium"
## [1] 31270
## [1] "removing: medi_bun_t_soft"
## [1] 31269
## [1] "removing: medi_bun_t_solution"
## [1] 31269
## [1] "removing: medi_bun_t_sound"
## [1] 31178
## [1] "removing: medi_bun_t_speech"
## [1] 31177
## [1] "removing: medi_bun_t_splint"
## [1] 31177
## [1] "removing: medi_bun_t_splt"
## [1] 31177
## [1] "removing: medi_bun_t_stent"
## [1] 31177
## [1] "removing: medi_bun_t_sterile"
## [1] 31177
## [1] "removing: medi_bun_t_stocking"
## [1] 31177
## [1] "removing: medi_bun_t_stone"
## [1] 31177
## [1] "removing: medi_bun_t_stool"
## [1] 31177
## [1] "removing: medi_bun_t_strep"
## [1] 31177
## [1] "removing: medi_bun_t_stress"
## [1] 31172
## [1] "removing: medi_bun_t_stud"
## [1] 31038
## [1] "removing: medi_bun_t_subseq"
## [1] 31037
## [1] "removing: medi_bun_t_substitute"
## [1] 31037
## [1] "removing: medi_bun_t_sulfate"
## [1] 30930
## [1] "removing: medi_bun_t_supplies"
## [1] 30930
## [1] "removing: medi_bun_t_suppls"
## [1] 30930
## [1] "removing: medi_bun_t_support"
## [1] 30929
## [1] "removing: medi_bun_t_susceptible"
## [1] 30929
## [1] "removing: medi_bun_t_swallow"
## [1] 30929
## [1] "removing: medi_bun_t_synth"
## [1] 30929
## [1] "removing: medi_bun_t_syphilis"
## [1] 30929
## [1] "removing: medi_bun_t_syringe"
## [1] 30929
## [1] "removing: medi_bun_t_system"
## [1] 30929
## [1] "removing: medi_bun_t_tdap"
## [1] 30926
## [1] "removing: medi_bun_t_test"
## [1] 30874
## [1] "removing: medi_bun_t_testing"
## [1] 30874
## [1] "removing: medi_bun_t_testosterone"
## [1] 30874
## [1] "removing: medi_bun_t_tests"
## [1] 30874
## [1] "removing: medi_bun_t_tetanus"
## [1] 30874
## [1] "removing: medi_bun_t_thyroxine"
## [1] 30874
## [1] "removing: medi_bun_t_tobacco"
## [1] 30874
## [1] "removing: medi_bun_t_toe"
## [1] 30874
## [1] "removing: medi_bun_t_tolerance"
## [1] 30874
## [1] "removing: medi_bun_t_topical"
## [1] 30874
## [1] "removing: medi_bun_t_toxin"
## [1] 30873
## [1] "removing: medi_bun_t_trach"
## [1] 30873
## [1] "removing: medi_bun_t_tracing"
## [1] 30668
## [1] "removing: medi_bun_t_traction"
## [1] 30084
## [1] "removing: medi_bun_t_training"
## [1] 30067
## [1] "removing: medi_bun_t_transfusion"
## [1] 30066
## [1] "removing: medi_bun_t_treat"
## [1] 30058
## [1] "removing: medi_bun_t_typhoid"
## [1] 30058
## [1] "removing: medi_bun_t_ultrasound"
## [1] 30058
## [1] "removing: medi_bun_t_urinalysis"
## [1] 30058
## [1] "removing: medi_bun_t_urinary"
## [1] 30058
## [1] "removing: medi_bun_t_urine"
## [1] 30058
## [1] "removing: medi_bun_t_uterine"
## [1] 30058
## [1] "removing: medi_bun_t_vacc"
## [1] 30046
## [1] "removing: medi_bun_t_valve"
## [1] 30046
## [1] "removing: medi_bun_t_vein"
## [1] 30045
## [1] "removing: medi_bun_t_virus"
## [1] 30045
## [1] "removing: medi_bun_t_visit"
## [1] 30043
## [1] "removing: medi_bun_t_visual"
## [1] 30039
## [1] "removing: medi_bun_t_vitamin"
## [1] 30023
## [1] "removing: medi_bun_t_wbc"
## [1] 30023
## [1] "removing: medi_bun_t_wheelchair"
## [1] 30023
## [1] "removing: medi_bun_t_wheezing"
## [1] 30023
## [1] "removing: medi_bun_t_whirlpool"
## [1] 30022
## [1] "removing: medi_bun_t_wire"
## [1] 30022
## [1] "removing: medi_bun_t_with"
## [1] 30003
## [1] "removing: medi_bun_t_without"
## [1] 30003
## [1] "removing: medi_bun_t_wkend"
## [1] 30003
## [1] "removing: medi_bun_t_wound"
## [1] 29923
## [1] "removing: medi_bun_t_menb"
## [1] 29923
## [1] "removing: medi_bun_t_hzv"
## [1] 29923
## [1] "removing: medi_bun_t_xgraft"
## [1] 29923
## [1] "removing: medi_bun_t_hepa"
## [1] 29923
## [1] "removing: medi_bun_t_hepb"
## [1] 29923
## [1] "removing: medi_bun_t_ppsv23"
## [1] 29923
## [1] "removing: medi_bun_t_psytx"
## [1] 29646
## [1] "removing: medi_bun_t_pcv13"
## [1] 29646
## [1] "removing: medi_bun_t_eye_exam"
## [1] 29573
## [1] "removing: medi_bun_t_var_vaccine_live_subq"
## [1] 29573
## [1] "removing: medi_bun_t_menacwy"
## [1] 29573
## [1] "removing: radi_bun_t_abdom"
## [1] 29431
## [1] "removing: radi_bun_t_adrenal"
## [1] 29431
## [1] "removing: radi_bun_t_angiography"
## [1] 29408
## [1] "removing: radi_bun_t_ankle"
## [1] 29252
## [1] "removing: radi_bun_t_aorta"
## [1] 29252
## [1] "removing: radi_bun_t_area"
## [1] 29252
## [1] "removing: radi_bun_t_arm"
## [1] 29210
## [1] "removing: radi_bun_t_artery"
## [1] 29210
## [1] "removing: radi_bun_t_bile"
## [1] 29210
## [1] "removing: radi_bun_t_bladder"
## [1] 29210
## [1] "removing: radi_bun_t_blade"
## [1] 29171
## [1] "removing: radi_bun_t_blockage"
## [1] 29171
## [1] "removing: radi_bun_t_body"
## [1] 29142
## [1] "removing: radi_bun_t_bone"
## [1] 28687
## [1] "removing: radi_bun_t_brain"
## [1] 28166
## [1] "removing: radi_bun_t_breast"
## [1] 28088
## [1] "removing: radi_bun_t_cardiac"
## [1] 28088
## [1] "removing: radi_bun_t_cast"
## [1] 28088
## [1] "removing: radi_bun_t_complex"
## [1] 28080
## [1] "removing: radi_bun_t_contrast"
## [1] 28028
## [1] "removing: radi_bun_t_cranial"
## [1] 28028
## [1] "removing: radi_bun_t_ct_scan"
## [1] 28027
## [1] "removing: radi_bun_t_delivery"
## [1] 28025
## [1] "removing: radi_bun_t_diagnostic"
## [1] 28025
## [1] "removing: radi_bun_t_dislocation"
## [1] 28025
## [1] "removing: radi_bun_t_duct"
## [1] 28025
## [1] "removing: radi_bun_t_echo"
## [1] 28025
## [1] "removing: radi_bun_t_elbow"
## [1] 27915
## [1] "removing: radi_bun_t_epidurograpy"
## [1] 27915
## [1] "removing: radi_bun_t_evaluation"
## [1] 27915
## [1] "removing: radi_bun_t_eye"
## [1] 27911
## [1] "removing: radi_bun_t_femur"
## [1] 27858
## [1] "removing: radi_bun_t_fetal"
## [1] 27858
## [1] "removing: radi_bun_t_finger"
## [1] 27834
## [1] "removing: radi_bun_t_fistula"
## [1] 27834
## [1] "removing: radi_bun_t_fluoroscopy"
## [1] 27834
## [1] "removing: radi_bun_t_fmri"
## [1] 27834
## [1] "removing: radi_bun_t_foot"
## [1] 27685
## [1] "removing: radi_bun_t_forearm"
## [1] 27685
## [1] "removing: radi_bun_t_fracture"
## [1] 27685
## [1] "removing: radi_bun_t_gadobutrol"
## [1] 27685
## [1] "removing: radi_bun_t_gastric"
## [1] 27685
## [1] "removing: radi_bun_t_general"
## [1] 27685
## [1] "removing: radi_bun_t_gential"
## [1] 27685
## [1] "removing: radi_bun_t_gestation"
## [1] 27685
## [1] "removing: radi_bun_t_gi_tract"
## [1] 27683
## [1] "removing: radi_bun_t_glands"
## [1] 27683
## [1] "removing: radi_bun_t_guide"
## [1] 27681
## [1] "removing: radi_bun_t_hand"
## [1] 27503
## [1] "removing: radi_bun_t_head"
## [1] 27479
## [1] "removing: radi_bun_t_heart"
## [1] 27479
## [1] "removing: radi_bun_t_heel"
## [1] 27476
## [1] "removing: radi_bun_t_humerus"
## [1] 27464
## [1] "removing: radi_bun_t_hyperthermia"
## [1] 27464
## [1] "removing: radi_bun_t_imag"
## [1] 27464
## [1] "removing: radi_bun_t_image"
## [1] 27464
## [1] "removing: radi_bun_t_imaging"
## [1] 27464
## [1] "removing: radi_bun_t_infant"
## [1] 27463
## [1] "removing: radi_bun_t_interpret"
## [1] 27463
## [1] "removing: radi_bun_t_jaw"
## [1] 27458
## [1] "removing: radi_bun_t_joint"
## [1] 27171
## [1] "removing: radi_bun_t_kidney"
## [1] 27171
## [1] "removing: radi_bun_t_leg"
## [1] 27148
## [1] "removing: radi_bun_t_lesion"
## [1] 27148
## [1] "removing: radi_bun_t_liver"
## [1] 27148
## [1] "removing: radi_bun_t_lower"
## [1] 27142
## [1] "removing: radi_bun_t_lung"
## [1] 27142
## [1] "removing: radi_bun_t_mammo"
## [1] 27110
## [1] "removing: radi_bun_t_mammogra"
## [1] 27110
## [1] "removing: radi_bun_t_management"
## [1] 27110
## [1] "removing: radi_bun_t_marrow"
## [1] 27110
## [1] "removing: radi_bun_t_mastoids"
## [1] 27110
## [1] "removing: radi_bun_t_mra"
## [1] 27110
## [1] "removing: radi_bun_t_myelography"
## [1] 27110
## [1] "removing: radi_bun_t_nuclear"
## [1] 27110
## [1] "removing: radi_bun_t_oncology"
## [1] 27110
## [1] "removing: radi_bun_t_pancreas"
## [1] 27110
## [1] "removing: radi_bun_t_penis"
## [1] 27110
## [1] "removing: radi_bun_t_pet"
## [1] 27110
## [1] "removing: radi_bun_t_pituitary"
## [1] 27110
## [1] "removing: radi_bun_t_plan"
## [1] 27110
## [1] "removing: radi_bun_t_plasma"
## [1] 27110
## [1] "removing: radi_bun_t_radiation"
## [1] 27110
## [1] "removing: radi_bun_t_radio"
## [1] 27109
## [1] "removing: radi_bun_t_repair"
## [1] 27109
## [1] "removing: radi_bun_t_rib"
## [1] 26936
## [1] "removing: radi_bun_t_scan"
## [1] 26936
## [1] "removing: radi_bun_t_sinuses"
## [1] 26933
## [1] "removing: radi_bun_t_skull"
## [1] 26927
## [1] "removing: radi_bun_t_socket"
## [1] 26927
## [1] "removing: radi_bun_t_speech"
## [1] 26927
## [1] "removing: radi_bun_t_spleen"
## [1] 26927
## [1] "removing: radi_bun_t_tear"
## [1] 26927
## [1] "removing: radi_bun_t_teeth"
## [1] 26927
## [1] "removing: radi_bun_t_therapy"
## [1] 26927
## [1] "removing: radi_bun_t_thigh"
## [1] 26927
## [1] "removing: radi_bun_t_thorax"
## [1] 26911
## [1] "removing: radi_bun_t_throat"
## [1] 26910
## [1] "removing: radi_bun_t_thrombosis"
## [1] 26910
## [1] "removing: radi_bun_t_thyroid"
## [1] 26910
## [1] "removing: radi_bun_t_toe"
## [1] 26899
## [1] "removing: radi_bun_t_tract"
## [1] 26899
## [1] "removing: radi_bun_t_treatment"
## [1] 26899
## [1] "removing: radi_bun_t_trunk"
## [1] 26897
## [1] "removing: radi_bun_t_tumor"
## [1] 26897
## [1] "removing: radi_bun_t_upper"
## [1] 26893
## [1] "removing: radi_bun_t_urethra"
## [1] 26893
## [1] "removing: radi_bun_t_urinary"
## [1] 26893
## [1] "removing: radi_bun_t_vein"
## [1] 26893
## [1] "removing: radi_bun_t_ventilation"
## [1] 26893
## [1] "removing: radi_bun_t_veous"
## [1] 26893
## [1] "removing: radi_bun_t_vessel"
## [1] 26893
## [1] "removing: radi_bun_t_wrist"
## [1] 26800
## [1] "removing: radi_bun_t_with_dye"
## [1] 26766
## [1] "removing: radi_bun_t_xray"
## [1] 26766
## [1] "removing: radi_bun_t_xtremity"
## [1] 26766
## [1] "removing: path_bun_t_acet"
## [1] 26765
## [1] "removing: path_bun_t_acid"
## [1] 26644
## [1] "removing: path_bun_t_acute"
## [1] 26635
## [1] "removing: path_bun_t_agent"
## [1] 26627
## [1] "removing: path_bun_t_alkaline"
## [1] 26626
## [1] "removing: path_bun_t_allerg"
## [1] 26621
## [1] "removing: path_bun_t_amikacin"
## [1] 26621
## [1] "removing: path_bun_t_amino"
## [1] 26597
## [1] "removing: path_bun_t_amitripytline"
## [1] 26597
## [1] "removing: path_bun_t_ammonia"
## [1] 26597
## [1] "removing: path_bun_t_amniotic"
## [1] 26597
## [1] "removing: path_bun_t_anabolic"
## [1] 26597
## [1] "removing: path_bun_t_analgesics"
## [1] 26596
## [1] "removing: path_bun_t_analysis"
## [1] 26596
## [1] "removing: path_bun_t_anti"
## [1] 26296
## [1] "removing: path_bun_t_antibod"
## [1] 26296
## [1] "removing: path_bun_t_antidepressant"
## [1] 26296
## [1] "removing: path_bun_t_antiepileptics"
## [1] 26296
## [1] "removing: path_bun_t_antipsychotic"
## [1] 26296
## [1] "removing: path_bun_t_aortic"
## [1] 26296
## [1] "removing: path_bun_t_autopsy"
## [1] 26296
## [1] "removing: path_bun_t_bact"
## [1] 26245
## [1] "removing: path_bun_t_bacteria"
## [1] 26245
## [1] "removing: path_bun_t_barbiturate"
## [1] 26245
## [1] "removing: path_bun_t_benzo"
## [1] 26245
## [1] "removing: path_bun_t_bile"
## [1] 26245
## [1] "removing: path_bun_t_bilirubin"
## [1] 26245
## [1] "removing: path_bun_t_biopsy"
## [1] 26245
## [1] "removing: path_bun_t_bleeding"
## [1] 26245
## [1] "removing: path_bun_t_blood"
## [1] 26189
## [1] "removing: path_bun_t_cadmium"
## [1] 26189
## [1] "removing: path_bun_t_calcium"
## [1] 26170
## [1] "removing: path_bun_t_cannabinoids"
## [1] 26170
## [1] "removing: path_bun_t_carbamazepine"
## [1] 26170
## [1] "removing: path_bun_t_carbon"
## [1] 26170
## [1] "removing: path_bun_t_catecholamine"
## [1] 26170
## [1] "removing: path_bun_t_chlamydia"
## [1] 26170
## [1] "removing: path_bun_t_chloride"
## [1] 26170
## [1] "removing: path_bun_t_chromosome"
## [1] 26170
## [1] "removing: path_bun_t_chromotograph"
## [1] 26170
## [1] "removing: path_bun_t_class"
## [1] 26170
## [1] "removing: path_bun_t_clot"
## [1] 26170
## [1] "removing: path_bun_t_clozapine"
## [1] 26170
## [1] "removing: path_bun_t_coagulation"
## [1] 26170
## [1] "removing: path_bun_t_cocaine"
## [1] 26170
## [1] "removing: path_bun_t_colon"
## [1] 26105
## [1] "removing: path_bun_t_count"
## [1] 26086
## [1] "removing: path_bun_t_cpk"
## [1] 26061
## [1] "removing: path_bun_t_creatinine"
## [1] 26045
## [1] "removing: path_bun_t_cryopreservd"
## [1] 26045
## [1] "removing: path_bun_t_culture"
## [1] 26021
## [1] "removing: path_bun_t_cyanide"
## [1] 26021
## [1] "removing: path_bun_t_cyclosporine"
## [1] 26021
## [1] "removing: path_bun_t_cystines"
## [1] 26021
## [1] "removing: path_bun_t_cytopath"
## [1] 25994
## [1] "removing: path_bun_t_cytopathol"
## [1] 25994
## [1] "removing: path_bun_t_dehydroepiandrosterone"
## [1] 25987
## [1] "removing: path_bun_t_desipramine"
## [1] 25987
## [1] "removing: path_bun_t_detect"
## [1] 25986
## [1] "removing: path_bun_t_digoxin"
## [1] 25986
## [1] "removing: path_bun_t_dihydro"
## [1] 25984
## [1] "removing: path_bun_t_dilation"
## [1] 25984
## [1] "removing: path_bun_t_dioxide"
## [1] 25984
## [1] "removing: path_bun_t_dipropylacetic"
## [1] 25984
## [1] "removing: path_bun_t_disease"
## [1] 25984
## [1] "removing: path_bun_t_dna"
## [1] 25941
## [1] "removing: path_bun_t_draw"
## [1] 25941
## [1] "removing: path_bun_t_drug"
## [1] 25908
## [1] "removing: path_bun_t_dysfunction"
## [1] 25908
## [1] "removing: path_bun_t_embryo"
## [1] 25908
## [1] "removing: path_bun_t_encephalitis"
## [1] 25908
## [1] "removing: path_bun_t_endocrine"
## [1] 25908
## [1] "removing: path_bun_t_enzyme"
## [1] 25900
## [1] "removing: path_bun_t_epstein-barr"
## [1] 25900
## [1] "removing: path_bun_t_estrogen"
## [1] 25897
## [1] "removing: path_bun_t_ethanol"
## [1] 25897
## [1] "removing: path_bun_t_evaluation"
## [1] 25897
## [1] "removing: path_bun_t_factor"
## [1] 25872
## [1] "removing: path_bun_t_fecal"
## [1] 25872
## [1] "removing: path_bun_t_fentanyl"
## [1] 25872
## [1] "removing: path_bun_t_ferritin"
## [1] 25819
## [1] "removing: path_bun_t_fetal"
## [1] 25819
## [1] "removing: path_bun_t_fibrinogen"
## [1] 25819
## [1] "removing: path_bun_t_fibronlytic"
## [1] 25819
## [1] "removing: path_bun_t_flouride"
## [1] 25819
## [1] "removing: path_bun_t_fluid"
## [1] 25818
## [1] "removing: path_bun_t_free"
## [1] 25700
## [1] "removing: path_bun_t_function"
## [1] 25671
## [1] "removing: path_bun_t_galacto"
## [1] 25671
## [1] "removing: path_bun_t_gases"
## [1] 25671
## [1] "removing: path_bun_t_gastrin"
## [1] 25671
## [1] "removing: path_bun_t_gene"
## [1] 25527
## [1] "removing: path_bun_t_general"
## [1] 25527
## [1] "removing: path_bun_t_genome"
## [1] 25527
## [1] "removing: path_bun_t_genotpying"
## [1] 25527
## [1] "removing: path_bun_t_giardia"
## [1] 25527
## [1] "removing: path_bun_t_glucagon"
## [1] 25527
## [1] "removing: path_bun_t_glucose"
## [1] 25525
## [1] "removing: path_bun_t_gonorr"
## [1] 25525
## [1] "removing: path_bun_t_growth"
## [1] 25525
## [1] "removing: path_bun_t_gynecologic"
## [1] 25525
## [1] "removing: path_bun_t_haloperidol"
## [1] 25525
## [1] "removing: path_bun_t_health"
## [1] 25525
## [1] "removing: path_bun_t_hearing"
## [1] 25525
## [1] "removing: path_bun_t_heavy"
## [1] 25521
## [1] "removing: path_bun_t_heinz"
## [1] 25521
## [1] "removing: path_bun_t_hematocrit"
## [1] 25519
## [1] "removing: path_bun_t_hemoglobin"
## [1] 25322
## [1] "removing: path_bun_t_hep"
## [1] 25316
## [1] "removing: path_bun_t_hepati"
## [1] 25316
## [1] "removing: path_bun_t_heridatary"
## [1] 25316
## [1] "removing: path_bun_t_heroin"
## [1] 25316
## [1] "removing: path_bun_t_herpes"
## [1] 25316
## [1] "removing: path_bun_t_hiv"
## [1] 25314
## [1] "removing: path_bun_t_hold"
## [1] 25314
## [1] "removing: path_bun_t_hormone"
## [1] 25253
## [1] "removing: path_bun_t_hpa"
## [1] 25253
## [1] "removing: path_bun_t_hsv"
## [1] 25253
## [1] "removing: path_bun_t_imipramine"
## [1] 25253
## [1] "removing: path_bun_t_immuno"
## [1] 25250
## [1] "removing: path_bun_t_immunoassay"
## [1] 25250
## [1] "removing: path_bun_t_immunohisto"
## [1] 25250
## [1] "removing: path_bun_t_influenza"
## [1] 25234
## [1] "removing: path_bun_t_insulin"
## [1] 25234
## [1] "removing: path_bun_t_iodine"
## [1] 25234
## [1] "removing: path_bun_t_iron"
## [1] 25232
## [1] "removing: path_bun_t_lab"
## [1] 25232
## [1] "removing: path_bun_t_lamotrigine"
## [1] 25232
## [1] "removing: path_bun_t_levetiracetam"
## [1] 25232
## [1] "removing: path_bun_t_lidocaine"
## [1] 25232
## [1] "removing: path_bun_t_lipid"
## [1] 25149
## [1] "removing: path_bun_t_lipoprotein"
## [1] 25148
## [1] "removing: path_bun_t_lithium"
## [1] 25148
## [1] "removing: path_bun_t_loss"
## [1] 25148
## [1] "removing: path_bun_t_low"
## [1] 25148
## [1] "removing: path_bun_t_lung"
## [1] 25148
## [1] "removing: path_bun_t_lyme"
## [1] 25148
## [1] "removing: path_bun_t_lympg"
## [1] 25148
## [1] "removing: path_bun_t_magnesium"
## [1] 25141
## [1] "removing: path_bun_t_malaria"
## [1] 25141
## [1] "removing: path_bun_t_mercury"
## [1] 25141
## [1] "removing: path_bun_t_metal"
## [1] 25141
## [1] "removing: path_bun_t_microbe"
## [1] 25141
## [1] "removing: path_bun_t_microscop"
## [1] 25141
## [1] "removing: path_bun_t_microscopic"
## [1] 25141
## [1] "removing: path_bun_t_microscopy"
## [1] 25141
## [1] "removing: path_bun_t_mitocondrial"
## [1] 25141
## [1] "removing: path_bun_t_mopath"
## [1] 25141
## [1] "removing: path_bun_t_mumps"
## [1] 25141
## [1] "removing: path_bun_t_mycophenolate"
## [1] 25141
## [1] "removing: path_bun_t_nicotine"
## [1] 25141
## [1] "removing: path_bun_t_nile"
## [1] 25141
## [1] "removing: path_bun_t_non-opioid"
## [1] 25141
## [1] "removing: path_bun_t_nonantibody"
## [1] 25141
## [1] "removing: path_bun_t_nucleic"
## [1] 25141
## [1] "removing: path_bun_t_obstretric"
## [1] 25141
## [1] "removing: path_bun_t_occult"
## [1] 25141
## [1] "removing: path_bun_t_onco"
## [1] 25141
## [1] "removing: path_bun_t_oocyte"
## [1] 25141
## [1] "removing: path_bun_t_opioid"
## [1] 25141
## [1] "removing: path_bun_t_ovulation"
## [1] 25141
## [1] "removing: path_bun_t_oxycodone"
## [1] 25141
## [1] "removing: path_bun_t_paper"
## [1] 25141
## [1] "removing: path_bun_t_patho"
## [1] 25130
## [1] "removing: path_bun_t_pathologist"
## [1] 25130
## [1] "removing: path_bun_t_period"
## [1] 25130
## [1] "removing: path_bun_t_phenobarbital"
## [1] 25130
## [1] "removing: path_bun_t_phenytoin"
## [1] 25130
## [1] "removing: path_bun_t_phosphatase"
## [1] 25130
## [1] "removing: path_bun_t_phosphorus"
## [1] 25129
## [1] "removing: path_bun_t_pituitary"
## [1] 25129
## [1] "removing: path_bun_t_plasminogen"
## [1] 25129
## [1] "removing: path_bun_t_poliovirus"
## [1] 25129
## [1] "removing: path_bun_t_potassium"
## [1] 25129
## [1] "removing: path_bun_t_pregabalin"
## [1] 25129
## [1] "removing: path_bun_t_pregnancy"
## [1] 25052
## [1] "removing: path_bun_t_primidone"
## [1] 25052
## [1] "removing: path_bun_t_procainamide"
## [1] 25052
## [1] "removing: path_bun_t_procedure"
## [1] 25052
## [1] "removing: path_bun_t_progesterone"
## [1] 25050
## [1] "removing: path_bun_t_prophobulinogen"
## [1] 25050
## [1] "removing: path_bun_t_prophyrins"
## [1] 25050
## [1] "removing: path_bun_t_propoxyphene"
## [1] 25050
## [1] "removing: path_bun_t_prostate"
## [1] 25050
## [1] "removing: path_bun_t_protein"
## [1] 24936
## [1] "removing: path_bun_t_psa"
## [1] 24921
## [1] "removing: path_bun_t_quantitative"
## [1] 24921
## [1] "removing: path_bun_t_quinide"
## [1] 24921
## [1] "removing: path_bun_t_quinine"
## [1] 24921
## [1] "removing: path_bun_t_reagent"
## [1] 24921
## [1] "removing: path_bun_t_renal"
## [1] 24921
## [1] "removing: path_bun_t_renin"
## [1] 24921
## [1] "removing: path_bun_t_retinal"
## [1] 24921
## [1] "removing: path_bun_t_rickettsia"
## [1] 24921
## [1] "removing: path_bun_t_rotovirus"
## [1] 24921
## [1] "removing: path_bun_t_rubella"
## [1] 24921
## [1] "removing: path_bun_t_salicylate"
## [1] 24921
## [1] "removing: path_bun_t_salmonella"
## [1] 24921
## [1] "removing: path_bun_t_scan"
## [1] 24921
## [1] "removing: path_bun_t_screen"
## [1] 24921
## [1] "removing: path_bun_t_sedative"
## [1] 24921
## [1] "removing: path_bun_t_semen"
## [1] 24920
## [1] "removing: path_bun_t_sequence"
## [1] 24920
## [1] "removing: path_bun_t_serotonin"
## [1] 24920
## [1] "removing: path_bun_t_sex"
## [1] 24920
## [1] "removing: path_bun_t_silica"
## [1] 24920
## [1] "removing: path_bun_t_sirolimus"
## [1] 24920
## [1] "removing: path_bun_t_skin"
## [1] 24920
## [1] "removing: path_bun_t_slenium"
## [1] 24920
## [1] "removing: path_bun_t_smear"
## [1] 24919
## [1] "removing: path_bun_t_specimen"
## [1] 24919
## [1] "removing: path_bun_t_sperm"
## [1] 24919
## [1] "removing: path_bun_t_steroid"
## [1] 24919
## [1] "removing: path_bun_t_stimulation"
## [1] 24919
## [1] "removing: path_bun_t_strep"
## [1] 24889
## [1] "removing: path_bun_t_syphilis"
## [1] 24889
## [1] "removing: path_bun_t_syringe"
## [1] 24889
## [1] "removing: path_bun_t_tacrolims"
## [1] 24889
## [1] "removing: path_bun_t_testosterone"
## [1] 24888
## [1] "removing: path_bun_t_thawing"
## [1] 24888
## [1] "removing: path_bun_t_theophylline"
## [1] 24888
## [1] "removing: path_bun_t_thyroid"
## [1] 24888
## [1] "removing: path_bun_t_time"
## [1] 24832
## [1] "removing: path_bun_t_tissue"
## [1] 24832
## [1] "removing: path_bun_t_tobramycin"
## [1] 24832
## [1] "removing: path_bun_t_tolerance"
## [1] 24832
## [1] "removing: path_bun_t_total"
## [1] 24763
## [1] "removing: path_bun_t_toxoplasma"
## [1] 24763
## [1] "removing: path_bun_t_tpiramate"
## [1] 24763
## [1] "removing: path_bun_t_transfusion"
## [1] 24763
## [1] "removing: path_bun_t_tumor"
## [1] 24763
## [1] "removing: path_bun_t_typing"
## [1] 24752
## [1] "removing: path_bun_t_urine"
## [1] 24752
## [1] "removing: path_bun_t_vancomycin"
## [1] 24752
## [1] "removing: path_bun_t_vein"
## [1] 24752
## [1] "removing: path_bun_t_venipuncture"
## [1] 24752
## [1] "removing: path_bun_t_virus"
## [1] 24751
## [1] "removing: path_bun_t_vitamin"
## [1] 24741
## [1] "removing: path_bun_t_west"
## [1] 24741
## [1] "removing: path_bun_t_with_scope"
## [1] 24741
## [1] "removing: path_bun_t_without_scope"
## [1] 24510
## [1] "removing: path_bun_t_zonisamide"
## [1] 24510
## [1] "removing: anes_bun_t_abdom"
## [1] 24510
## [1] "removing: anes_bun_t_add - on"
## [1] 24510
## [1] "removing: anes_bun_t_adrenal"
## [1] 24510
## [1] "removing: anes_bun_t_amputation"
## [1] 24510
## [1] "removing: anes_bun_t_analg"
## [1] 24510
## [1] "removing: anes_bun_t_anes"
## [1] 24501
## [1] "removing: anes_bun_t_anesth"
## [1] 24501
## [1] "removing: anes_bun_t_area"
## [1] 24501
## [1] "removing: anes_bun_t_arest"
## [1] 24501
## [1] "removing: anes_bun_t_arteriography"
## [1] 24501
## [1] "removing: anes_bun_t_artery"
## [1] 24501
## [1] "removing: anes_bun_t_biopsy"
## [1] 24501
## [1] "removing: anes_bun_t_bladder"
## [1] 24501
## [1] "removing: anes_bun_t_blood"
## [1] 24501
## [1] "removing: anes_bun_t_bone"
## [1] 24501
## [1] "removing: anes_bun_t_breast"
## [1] 24501
## [1] "removing: anes_bun_t_bronchi"
## [1] 24501
## [1] "removing: anes_bun_t_burn"
## [1] 24501
## [1] "removing: anes_bun_t_casting"
## [1] 24501
## [1] "removing: anes_bun_t_cervix"
## [1] 24501
## [1] "removing: anes_bun_t_chemonucleolysis"
## [1] 24501
## [1] "removing: anes_bun_t_chest"
## [1] 24501
## [1] "removing: anes_bun_t_cleft"
## [1] 24501
## [1] "removing: anes_bun_t_cord"
## [1] 24501
## [1] "removing: anes_bun_t_corneal"
## [1] 24501
## [1] "removing: anes_bun_t_cs"
## [1] 24501
## [1] "removing: anes_bun_t_debride"
## [1] 24501
## [1] "removing: anes_bun_t_deliver"
## [1] 24501
## [1] "removing: anes_bun_t_drainage"
## [1] 24501
## [1] "removing: anes_bun_t_duct"
## [1] 24501
## [1] "removing: anes_bun_t_ear"
## [1] 24501
## [1] "removing: anes_bun_t_elbow"
## [1] 24501
## [1] "removing: anes_bun_t_electroshock"
## [1] 24501
## [1] "removing: anes_bun_t_embolectomy"
## [1] 24501
## [1] "removing: anes_bun_t_endoscopy"
## [1] 24501
## [1] "removing: anes_bun_t_epidural"
## [1] 24501
## [1] "removing: anes_bun_t_esophag"
## [1] 24501
## [1] "removing: anes_bun_t_eyelid"
## [1] 24501
## [1] "removing: anes_bun_t_facial"
## [1] 24501
## [1] "removing: anes_bun_t_fat"
## [1] 24501
## [1] "removing: anes_bun_t_femoral"
## [1] 24501
## [1] "removing: anes_bun_t_fentanyl"
## [1] 24501
## [1] "removing: anes_bun_t_genitalia"
## [1] 24501
## [1] "removing: anes_bun_t_glands"
## [1] 24501
## [1] "removing: anes_bun_t_head"
## [1] 24501
## [1] "removing: anes_bun_t_hernia"
## [1] 24501
## [1] "removing: anes_bun_t_humer"
## [1] 24501
## [1] "removing: anes_bun_t_hydromorphone"
## [1] 24501
## [1] "removing: anes_bun_t_hysterectomy"
## [1] 24501
## [1] "removing: anes_bun_t_image"
## [1] 24501
## [1] "removing: anes_bun_t_insertion"
## [1] 24501
## [1] "removing: anes_bun_t_intracranial"
## [1] 24501
## [1] "removing: anes_bun_t_iredectomy"
## [1] 24501
## [1] "removing: anes_bun_t_joint"
## [1] 24501
## [1] "removing: anes_bun_t_ketorolac"
## [1] 24501
## [1] "removing: anes_bun_t_kidney"
## [1] 24501
## [1] "removing: anes_bun_t_knee"
## [1] 24501
## [1] "removing: anes_bun_t_larynx"
## [1] 24501
## [1] "removing: anes_bun_t_layer"
## [1] 24501
## [1] "removing: anes_bun_t_leg"
## [1] 24501
## [1] "removing: anes_bun_t_lip"
## [1] 24501
## [1] "removing: anes_bun_t_liver"
## [1] 24501
## [1] "removing: anes_bun_t_local"
## [1] 24501
## [1] "removing: anes_bun_t_lower"
## [1] 24501
## [1] "removing: anes_bun_t_lumbar"
## [1] 24501
## [1] "removing: anes_bun_t_lung"
## [1] 24501
## [1] "removing: anes_bun_t_major"
## [1] 24501
## [1] "removing: anes_bun_t_manipulation"
## [1] 24501
## [1] "removing: anes_bun_t_midazolam"
## [1] 24501
## [1] "removing: anes_bun_t_mod_sed"
## [1] 24499
## [1] "removing: anes_bun_t_neck"
## [1] 24499
## [1] "removing: anes_bun_t_nerve"
## [1] 24499
## [1] "removing: anes_bun_t_nose"
## [1] 24499
## [1] "removing: anes_bun_t_open"
## [1] 24499
## [1] "removing: anes_bun_t_organ"
## [1] 24499
## [1] "removing: anes_bun_t_pacemaker"
## [1] 24499
## [1] "removing: anes_bun_t_palate"
## [1] 24499
## [1] "removing: anes_bun_t_pelvic"
## [1] 24499
## [1] "removing: anes_bun_t_pelvis"
## [1] 24499
## [1] "removing: anes_bun_t_penis"
## [1] 24499
## [1] "removing: anes_bun_t_percutaneous"
## [1] 24499
## [1] "removing: anes_bun_t_plastic"
## [1] 24499
## [1] "removing: anes_bun_t_propofol"
## [1] 24499
## [1] "removing: anes_bun_t_prostate"
## [1] 24499
## [1] "removing: anes_bun_t_pump"
## [1] 24499
## [1] "removing: anes_bun_t_puncture"
## [1] 24499
## [1] "removing: anes_bun_t_radical"
## [1] 24499
## [1] "removing: anes_bun_t_release"
## [1] 24499
## [1] "removing: anes_bun_t_removal"
## [1] 24499
## [1] "removing: anes_bun_t_repair"
## [1] 24499
## [1] "removing: anes_bun_t_replace"
## [1] 24499
## [1] "removing: anes_bun_t_rib"
## [1] 24499
## [1] "removing: anes_bun_t_sedation"
## [1] 24499
## [1] "removing: anes_bun_t_shoulder"
## [1] 24499
## [1] "removing: anes_bun_t_sinus"
## [1] 24499
## [1] "removing: anes_bun_t_skull"
## [1] 24499
## [1] "removing: anes_bun_t_sperm"
## [1] 24499
## [1] "removing: anes_bun_t_spine"
## [1] 24499
## [1] "removing: anes_bun_t_spit"
## [1] 24499
## [1] "removing: anes_bun_t_sternal"
## [1] 24499
## [1] "removing: anes_bun_t_stone"
## [1] 24499
## [1] "removing: anes_bun_t_surg"
## [1] 24499
## [1] "removing: anes_bun_t_tendon"
## [1] 24499
## [1] "removing: anes_bun_t_testis"
## [1] 24499
## [1] "removing: anes_bun_t_thigh"
## [1] 24499
## [1] "removing: anes_bun_t_thyriod"
## [1] 24499
## [1] "removing: anes_bun_t_trachea"
## [1] 24499
## [1] "removing: anes_bun_t_transplant"
## [1] 24499
## [1] "removing: anes_bun_t_tromethamine"
## [1] 24499
## [1] "removing: anes_bun_t_tubal"
## [1] 24499
## [1] "removing: anes_bun_t_tympanotomy"
## [1] 24499
## [1] "removing: anes_bun_t_upper"
## [1] 24499
## [1] "removing: anes_bun_t_ureter"
## [1] 24499
## [1] "removing: anes_bun_t_urethra"
## [1] 24499
## [1] "removing: anes_bun_t_vag"
## [1] 24499
## [1] "removing: anes_bun_t_vein"
## [1] 24499
## [1] "removing: anes_bun_t_ventilation"
## [1] 24499
## [1] "removing: anes_bun_t_vesectomy"
## [1] 24499
## [1] "removing: anes_bun_t_vessel"
## [1] 24499
## [1] "removing: anes_bun_t_viteroretinal"
## [1] 24499
## [1] "removing: anes_bun_t_wall"
## [1] 24499
## [1] "removing: anes_bun_t_wrist"
## [1] 24499
## [1] "removing: faci_bun_t_advanced"
## [1] 24499
## [1] "removing: faci_bun_t_analysis"
## [1] 24499
## [1] "removing: faci_bun_t_annual"
## [1] 24499
## [1] "removing: faci_bun_t_assessmnt"
## [1] 24499
## [1] "removing: faci_bun_t_audit"
## [1] 24499
## [1] "removing: faci_bun_t_basic"
## [1] 24499
## [1] "removing: faci_bun_t_behav"
## [1] 24496
## [1] "removing: faci_bun_t_bus"
## [1] 24496
## [1] "removing: faci_bun_t_care"
## [1] 24471
## [1] "removing: faci_bun_t_chng"
## [1] 24471
## [1] "removing: faci_bun_t_computer"
## [1] 24471
## [1] "removing: faci_bun_t_conf"
## [1] 24471
## [1] "removing: faci_bun_t_consultation"
## [1] 24450
## [1] "removing: faci_bun_t_counseling"
## [1] 24450
## [1] "removing: faci_bun_t_critical"
## [1] 24450
## [1] "removing: faci_bun_t_data"
## [1] 24450
## [1] "removing: faci_bun_t_direct"
## [1] 24450
## [1] "removing: faci_bun_t_disability"
## [1] 24450
## [1] "removing: faci_bun_t_discharge"
## [1] 24449
## [1] "removing: faci_bun_t_dommicl"
## [1] 24449
## [1] "removing: faci_bun_t_exam"
## [1] 24449
## [1] "removing: faci_bun_t_fac"
## [1] 24447
## [1] "removing: faci_bun_t_facility"
## [1] 24447
## [1] "removing: faci_bun_t_follow-up"
## [1] 24410
## [1] "removing: faci_bun_t_handling"
## [1] 24409
## [1] "removing: faci_bun_t_home"
## [1] 24408
## [1] "removing: faci_bun_t_hospice"
## [1] 24408
## [1] "removing: faci_bun_t_hospital"
## [1] 24408
## [1] "removing: faci_bun_t_initial"
## [1] 24408
## [1] "removing: faci_bun_t_inpatient"
## [1] 24408
## [1] "removing: faci_bun_t_interfacility"
## [1] 24408
## [1] "removing: faci_bun_t_interpof"
## [1] 24408
## [1] "removing: faci_bun_t_lab"
## [1] 24408
## [1] "removing: faci_bun_t_life"
## [1] 24408
## [1] "removing: faci_bun_t_nonemergency"
## [1] 24408
## [1] "removing: faci_bun_t_noner"
## [1] 24408
## [1] "removing: faci_bun_t_nursing"
## [1] 24408
## [1] "removing: faci_bun_t_observation"
## [1] 24408
## [1] "removing: faci_bun_t_on_call"
## [1] 24408
## [1] "removing: faci_bun_t_online"
## [1] 24408
## [1] "removing: faci_bun_t_phlebotomy"
## [1] 24407
## [1] "removing: faci_bun_t_phone"
## [1] 24407
## [1] "removing: faci_bun_t_physician"
## [1] 24407
## [1] "removing: faci_bun_t_postop"
## [1] 24407
## [1] "removing: faci_bun_t_preventine"
## [1] 24407
## [1] "removing: faci_bun_t_prolong"
## [1] 24396
## [1] "removing: faci_bun_t_service"
## [1] 24394
## [1] "removing: faci_bun_t_smoking"
## [1] 24394
## [1] "removing: faci_bun_t_specimen"
## [1] 24394
## [1] "removing: faci_bun_t_standby"
## [1] 24394
## [1] "removing: faci_bun_t_subsequent"
## [1] 24394
## [1] "removing: faci_bun_t_support"
## [1] 24394
## [1] "removing: faci_bun_t_taxi"
## [1] 24394
## [1] "removing: faci_bun_t_team"
## [1] 24394
## [1] "removing: faci_bun_t_transort"
## [1] 24394
## [1] "removing: faci_bun_t_travel"
## [1] 24394
## [1] "removing: faci_bun_t_van"
## [1] 24394
xray_spine_c %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

xray_spine_c %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 32 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 radi_bun_t_without_dye      0.492 
##  2 faci_bun_t_dept             0.483 
##  3 faci_bun_t_emergency        0.483 
##  4 radi_bun_t_mri              0.420 
##  5 radi_bun_t_lumbar           0.394 
##  6 radi_bun_t_ct               0.264 
##  7 medi_bun_t_inject           0.236 
##  8 radi_bun_t_chest            0.206 
##  9 radi_bun_t_view             0.105 
## 10 radi_bun_t_pelvi            0.0901
## # ... with 22 more rows
xray_spine_c %>% get_tag_density_information("radi_bun_t_without_dye") %>% 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 ~ radi_bun_t_without_dye"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 102
## $dist_plots

## 
## $stat_tables

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

xray_spine_c %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 28 x 2
##    name                  correlation
##    <chr>                       <dbl>
##  1 faci_bun_t_dept            0.620 
##  2 faci_bun_t_emergency       0.620 
##  3 medi_bun_t_inject          0.324 
##  4 faci_bun_t_visit           0.210 
##  5 path_bun_t_urinalysis      0.117 
##  6 radi_bun_t_chest           0.109 
##  7 radi_bun_t_thorac          0.106 
##  8 radi_bun_t_view            0.0949
##  9 radi_bun_t_pelvi           0.0753
## 10 path_bun_t_assay           0.075 
## # ... with 18 more rows
xray_spine_c %>% 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 63.9
## $dist_plots

## 
## $stat_tables

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

xray_spine_c %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 26 x 2
##    name                correlation
##    <chr>                     <dbl>
##  1 medi_bun_t_inject        0.226 
##  2 faci_bun_t_visit         0.207 
##  3 faci_bun_t_office        0.201 
##  4 faci_bun_t_pat           0.201 
##  5 radi_bun_t_thorac        0.0829
##  6 radi_bun_t_view          0.0812
##  7 radi_bun_t_pelvi         0.0685
##  8 path_bun_t_assay         0.0641
##  9 faci_bun_t_est           0.0534
## 10 medi_bun_t_exercise      0.0528
## # ... with 16 more rows
xray_spine_c %>% 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 399
## $dist_plots

## 
## $stat_tables

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

xray_spine_c %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 25 x 2
##    name                correlation
##    <chr>                     <dbl>
##  1 faci_bun_t_visit         0.228 
##  2 faci_bun_t_office        0.222 
##  3 faci_bun_t_pat           0.222 
##  4 radi_bun_t_view          0.09  
##  5 radi_bun_t_thorac        0.0868
##  6 radi_bun_t_pelvi         0.0755
##  7 faci_bun_t_est           0.0603
##  8 medi_bun_t_exercise      0.0542
##  9 radi_bun_t_shoulder      0.053 
## 10 medi_bun_t_puncture      0.0507
## # ... with 15 more rows
xray_spine_c %>% 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 19.7
## $dist_plots

## 
## $stat_tables

xray_spine_w_visit <- xray_spine_c[faci_bun_t_visit==T]
xray_spine_n_visit <- xray_spine_c[faci_bun_t_visit==F]

xray_spine_w_visit

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 24 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 radi_bun_t_view             0.122 
##  2 faci_bun_t_est              0.0981
##  3 radi_bun_t_thorac           0.0711
##  4 radi_bun_t_pelvi            0.0669
##  5 radi_bun_t_chest            0.0663
##  6 medi_bun_t_exercise         0.0635
##  7 path_bun_t_assay            0.0569
##  8 medi_bun_t_chiropract       0.0548
##  9 medi_bun_t_therapeutic      0.0533
## 10 radi_bun_t_shoulder         0.053 
## # ... with 14 more rows
xray_spine_w_visit %>% get_tag_density_information("radi_bun_t_view") %>% 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 ~ radi_bun_t_view"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 31.4
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 22 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 faci_bun_t_est              0.107 
##  2 radi_bun_t_thorac           0.0783
##  3 radi_bun_t_pelvi            0.0716
##  4 medi_bun_t_exercise         0.0615
##  5 radi_bun_t_shoulder         0.0599
##  6 path_bun_t_assay            0.0578
##  7 medi_bun_t_chiropract       0.0529
##  8 medi_bun_t_therapeutic      0.052 
##  9 path_bun_t_comprehen        0.0347
## 10 path_bun_t_metabolic        0.0331
## # ... with 12 more rows
xray_spine_w_visit %>% get_tag_density_information("radi_bun_t_thorac") %>% 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 ~ radi_bun_t_thorac"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 20.4
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 21 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_est               0.124 
##  2 radi_bun_t_shoulder          0.0731
##  3 radi_bun_t_pelvi             0.0724
##  4 path_bun_t_assay             0.0709
##  5 medi_bun_t_exercise          0.0623
##  6 medi_bun_t_chiropract        0.0606
##  7 medi_bun_t_therapeutic       0.0502
##  8 medi_bun_t_puncture          0.0413
##  9 medi_bun_t_routine           0.0413
## 10 medi_bun_t_venipuncture      0.0413
## # ... with 11 more rows
xray_spine_w_visit %>% get_tag_density_information("radi_bun_t_shoulder") %>% 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 ~ radi_bun_t_shoulder"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 26.6
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 20 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_est              0.126  
##  2 path_bun_t_assay            0.0749 
##  3 radi_bun_t_pelvi            0.0724 
##  4 medi_bun_t_exercise         0.0654 
##  5 medi_bun_t_chiropract       0.0602 
##  6 medi_bun_t_therapeutic      0.053  
##  7 medi_bun_t_puncture         0.0438 
##  8 medi_bun_t_routine          0.0438 
##  9 medi_bun_t_venipuncture     0.0438 
## 10 path_bun_t_comprehen        0.0343 
## 11 path_bun_t_metabolic        0.0325 
## 12 path_bun_t_panel            0.0325 
## 13 radi_bun_t_neck             0.0311 
## 14 medi_bun_t_stimulation      0.0295 
## 15 path_bun_t_complete         0.0291 
## 16 faci_bun_t_office           0.0146 
## 17 faci_bun_t_pat              0.0146 
## 18 radi_bun_t_exam             0.00960
## 19 path_bun_t_urinalysis       0.0055 
## 20 path_bun_t_test             0.0038
xray_spine_w_visit %>% get_tag_density_information("path_bun_t_assay") %>% 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_assay"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 129
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 19 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_est               0.128 
##  2 radi_bun_t_pelvi             0.0703
##  3 medi_bun_t_exercise          0.0634
##  4 medi_bun_t_chiropract        0.0599
##  5 medi_bun_t_therapeutic       0.0512
##  6 medi_bun_t_puncture          0.034 
##  7 medi_bun_t_routine           0.034 
##  8 medi_bun_t_venipuncture      0.034 
##  9 radi_bun_t_neck              0.0313
## 10 medi_bun_t_stimulation       0.0294
## 11 path_bun_t_comprehen         0.0248
## 12 path_bun_t_metabolic         0.0232
## 13 path_bun_t_panel             0.0232
## 14 path_bun_t_complete          0.0198
## 15 faci_bun_t_office            0.0151
## 16 faci_bun_t_pat               0.0151
## 17 radi_bun_t_exam              0.0094
## 18 path_bun_t_test              0.0038
## 19 path_bun_t_urinalysis        0.0007
xray_spine_w_visit %>% get_tag_density_information("radi_bun_t_pelvi") %>% 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 ~ radi_bun_t_pelvi"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 24.2
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 18 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_est              0.132  
##  2 medi_bun_t_exercise         0.0647 
##  3 medi_bun_t_chiropract       0.0582 
##  4 medi_bun_t_therapeutic      0.0525 
##  5 medi_bun_t_puncture         0.0346 
##  6 medi_bun_t_routine          0.0346 
##  7 medi_bun_t_venipuncture     0.0346 
##  8 medi_bun_t_stimulation      0.0297 
##  9 radi_bun_t_neck             0.0263 
## 10 path_bun_t_comprehen        0.0226 
## 11 path_bun_t_metabolic        0.0209 
## 12 path_bun_t_panel            0.0209 
## 13 path_bun_t_complete         0.017  
## 14 faci_bun_t_office           0.0118 
## 15 faci_bun_t_pat              0.0118 
## 16 radi_bun_t_exam             0.0086 
## 17 path_bun_t_test             0.0037 
## 18 path_bun_t_urinalysis       0.00120
xray_spine_w_visit %>% get_tag_density_information("medi_bun_t_exercise") %>% 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_exercise"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 33.4
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 17 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_est               0.141 
##  2 medi_bun_t_chiropract        0.0723
##  3 medi_bun_t_stimulation       0.0415
##  4 medi_bun_t_puncture          0.0354
##  5 medi_bun_t_routine           0.0354
##  6 medi_bun_t_venipuncture      0.0354
##  7 radi_bun_t_neck              0.0238
##  8 path_bun_t_comprehen         0.0237
##  9 path_bun_t_metabolic         0.022 
## 10 path_bun_t_panel             0.022 
## 11 path_bun_t_complete          0.0181
## 12 faci_bun_t_office            0.013 
## 13 faci_bun_t_pat               0.013 
## 14 radi_bun_t_exam              0.0078
## 15 medi_bun_t_therapeutic       0.0073
## 16 path_bun_t_test              0.0036
## 17 path_bun_t_urinalysis        0.0008
xray_spine_w_visit %>% get_tag_density_information("medi_bun_t_chiropract") %>% 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_chiropract"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 16 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 faci_bun_t_est               0.167 
##  2 medi_bun_t_stimulation       0.0413
##  3 medi_bun_t_puncture          0.035 
##  4 medi_bun_t_routine           0.035 
##  5 medi_bun_t_venipuncture      0.035 
##  6 path_bun_t_comprehen         0.0233
##  7 path_bun_t_metabolic         0.0215
##  8 path_bun_t_panel             0.0215
##  9 radi_bun_t_neck              0.0209
## 10 path_bun_t_complete          0.0146
## 11 faci_bun_t_office            0.0117
## 12 faci_bun_t_pat               0.0117
## 13 radi_bun_t_exam              0.0094
## 14 medi_bun_t_therapeutic       0.0061
## 15 path_bun_t_test              0.004 
## 16 path_bun_t_urinalysis        0.0023
xray_spine_w_visit %>% get_tag_density_information("medi_bun_t_venipuncture") %>% 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_venipuncture"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 43.8
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 13 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 faci_bun_t_est              0.168 
##  2 medi_bun_t_stimulation      0.0413
##  3 radi_bun_t_neck             0.0213
##  4 path_bun_t_comprehen        0.0141
##  5 path_bun_t_metabolic        0.0125
##  6 path_bun_t_panel            0.0125
##  7 path_bun_t_complete         0.0111
##  8 radi_bun_t_exam             0.0092
##  9 faci_bun_t_office           0.0089
## 10 faci_bun_t_pat              0.0089
## 11 path_bun_t_urinalysis       0.0068
## 12 medi_bun_t_therapeutic      0.0063
## 13 path_bun_t_test             0.0039
xray_spine_w_visit %>% get_tag_density_information("medi_bun_t_stimulation") %>% 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_stimulation"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 16
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 12 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 faci_bun_t_est              0.174 
##  2 radi_bun_t_neck             0.0205
##  3 medi_bun_t_therapeutic      0.015 
##  4 path_bun_t_comprehen        0.0139
##  5 path_bun_t_metabolic        0.0124
##  6 path_bun_t_panel            0.0124
##  7 path_bun_t_complete         0.0109
##  8 radi_bun_t_exam             0.0095
##  9 faci_bun_t_office           0.0087
## 10 faci_bun_t_pat              0.0087
## 11 path_bun_t_urinalysis       0.007 
## 12 path_bun_t_test             0.004
xray_spine_w_visit %>% get_tag_density_information("path_bun_t_panel") %>% 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_panel"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 29.4
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 9 x 2
##   name                   correlation
##   <chr>                        <dbl>
## 1 faci_bun_t_est              0.175 
## 2 radi_bun_t_neck             0.0206
## 3 medi_bun_t_therapeutic      0.015 
## 4 radi_bun_t_exam             0.0095
## 5 path_bun_t_urinalysis       0.0094
## 6 faci_bun_t_office           0.0088
## 7 faci_bun_t_pat              0.0088
## 8 path_bun_t_complete         0.0078
## 9 path_bun_t_test             0.004
xray_spine_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 67
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 8 x 2
##   name                  correlation
##   <chr>                       <dbl>
## 1 faci_bun_t_est             0.176 
## 2 radi_bun_t_neck            0.0213
## 3 radi_bun_t_exam            0.0094
## 4 path_bun_t_urinalysis      0.0094
## 5 faci_bun_t_office          0.0088
## 6 faci_bun_t_pat             0.0088
## 7 path_bun_t_complete        0.0078
## 8 path_bun_t_test            0.004
xray_spine_w_visit %>% get_tag_density_information("path_bun_t_urinalysis") %>% 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_urinalysis"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 16.3
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 7 x 2
##   name                correlation
##   <chr>                     <dbl>
## 1 faci_bun_t_est           0.176 
## 2 radi_bun_t_neck          0.0216
## 3 radi_bun_t_exam          0.0095
## 4 faci_bun_t_office        0.0088
## 5 faci_bun_t_pat           0.0088
## 6 path_bun_t_complete      0.0078
## 7 path_bun_t_test          0.004
xray_spine_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 33.2
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 6 x 2
##   name              correlation
##   <chr>                   <dbl>
## 1 faci_bun_t_est         0.177 
## 2 radi_bun_t_neck        0.0217
## 3 radi_bun_t_exam        0.0095
## 4 faci_bun_t_office      0.0088
## 5 faci_bun_t_pat         0.0088
## 6 path_bun_t_test        0.004
xray_spine_w_visit %>% get_tag_density_information("path_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 ~ path_bun_t_test"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 30.9
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 5 x 2
##   name              correlation
##   <chr>                   <dbl>
## 1 faci_bun_t_est         0.177 
## 2 radi_bun_t_neck        0.0217
## 3 radi_bun_t_exam        0.0098
## 4 faci_bun_t_office      0.0092
## 5 faci_bun_t_pat         0.0092
xray_spine_w_visit %>% get_tag_density_information("radi_bun_t_exam") %>% 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 ~ radi_bun_t_exam"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.4
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 4 x 2
##   name              correlation
##   <chr>                   <dbl>
## 1 faci_bun_t_est         0.179 
## 2 radi_bun_t_neck        0.0227
## 3 faci_bun_t_office      0.0092
## 4 faci_bun_t_pat         0.0092
xray_spine_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 30.7
## $dist_plots

## 
## $stat_tables

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

xray_spine_w_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 2 x 2
##   name            correlation
##   <chr>                 <dbl>
## 1 faci_bun_t_est       0.180 
## 2 radi_bun_t_neck      0.0244
xray_spine_w_visit %>% get_tag_density_information("radi_bun_t_neck") %>% 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 ~ radi_bun_t_neck"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.2
## $dist_plots

## 
## $stat_tables

xray_spine_w_visit_w_neck <- xray_spine_w_visit[radi_bun_t_neck==T]
xray_spine_w_visit_n_neck <- xray_spine_w_visit[radi_bun_t_neck==F]

xray_spine_w_visit_w_neck

xray_spine_w_visit_w_neck <- xray_spine_w_visit_w_neck[,`:=`(
  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)
                           )]
xray_spine_w_visit_w_neck_btbv5 <- xray_spine_w_visit_w_neck %>% btbv5_radi_fac_only(
                                  osa_group_p=NA,
                                  osa_class_p = NA,
                                  osa_specialization_p = NA
                                  ) %>% 
  .[,`:=`(procedure_type=39, 
          procedure_modifier="X-ray Spine + Neck + Visit", 
          ingest_date =Sys.Date())]
xray_spine_w_visit_w_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    55.0   161.3   192.2   239.7   273.3  1695.4
spine_q1 <- xray_spine_w_visit_w_neck_btbv5$tp_med_med %>% quantile(.25)
spine_q3 <- xray_spine_w_visit_w_neck_btbv5$tp_med_med %>% quantile(.75)
spine_iqr <- xray_spine_w_visit_w_neck_btbv5$tp_med_med %>% IQR()
xray_spine_w_visit_w_neck_btbv5_final <- xray_spine_w_visit_w_neck_btbv5[tp_cnt_cnt > 1 ]
# mri_brain_clean_btbv5_final <- mri_brain_clean_btbv5[tp_med_med < 5000 ]
xray_spine_w_visit_w_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    55.0   161.3   192.2   239.7   273.3  1695.4
xray_spine_w_visit_w_neck_btbv5_final$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    55.0   161.3   192.2   239.7   273.3  1695.4
bq_table_upload(x=procedure_dev_table, values= xray_spine_w_visit_w_neck_btbv5_final, create_disposition='CREATE_IF_NEEDED', write_disposition='WRITE_APPEND')

xray_spine_w_visit_n_neck

xray_spine_w_visit_n_neck <- xray_spine_w_visit_n_neck[,`:=`(
  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)
                           )]
xray_spine_w_visit_n_neck_btbv5 <- xray_spine_w_visit_n_neck %>% btbv5_radi_fac_only(
                                  osa_group_p=NA,
                                  osa_class_p = NA,
                                  osa_specialization_p = NA
                                  ) %>% 
  .[,`:=`(procedure_type=39, 
          procedure_modifier="X-ray Spine + Visit", 
          ingest_date =Sys.Date())]
xray_spine_w_visit_n_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    55.0   170.9   206.8   241.8   291.6   774.2
spine_q1 <- xray_spine_w_visit_n_neck_btbv5$tp_med_med %>% quantile(.25)
spine_q3 <- xray_spine_w_visit_n_neck_btbv5$tp_med_med %>% quantile(.75)
spine_iqr <- xray_spine_w_visit_n_neck_btbv5$tp_med_med %>% IQR()
xray_spine_w_visit_n_neck_btbv5_final <- xray_spine_w_visit_n_neck_btbv5[tp_cnt_cnt > 1 ]
# mri_brain_clean_btbv5_final <- mri_brain_clean_btbv5[tp_med_med < 5000 ]
xray_spine_w_visit_n_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    55.0   170.9   206.8   241.8   291.6   774.2
xray_spine_w_visit_n_neck_btbv5_final$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    55.0   170.9   206.8   241.8   291.6   774.2
bq_table_upload(x=procedure_dev_table, values= xray_spine_w_visit_n_neck_btbv5_final, create_disposition='CREATE_IF_NEEDED', write_disposition='WRITE_APPEND')

xray_spine_n_visit

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 21 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 radi_bun_t_thorac            0.132 
##  2 radi_bun_t_pelvi             0.103 
##  3 medi_bun_t_puncture          0.101 
##  4 medi_bun_t_routine           0.101 
##  5 medi_bun_t_venipuncture      0.101 
##  6 radi_bun_t_view              0.0503
##  7 medi_bun_t_exercise          0.0447
##  8 medi_bun_t_therapeutic       0.0373
##  9 radi_bun_t_shoulder          0.0289
## 10 medi_bun_t_chiropract        0.0266
## # ... with 11 more rows
xray_spine_n_visit %>% get_tag_density_information("radi_bun_t_pelvi") %>% 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 ~ radi_bun_t_pelvi"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 64.1
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 20 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 radi_bun_t_thorac            0.130 
##  2 medi_bun_t_puncture          0.105 
##  3 medi_bun_t_routine           0.105 
##  4 medi_bun_t_venipuncture      0.105 
##  5 medi_bun_t_exercise          0.0491
##  6 radi_bun_t_view              0.045 
##  7 medi_bun_t_therapeutic       0.0419
##  8 medi_bun_t_chiropract        0.0404
##  9 radi_bun_t_shoulder          0.0305
## 10 path_bun_t_assay             0.0244
## 11 path_bun_t_comprehen         0.0237
## 12 path_bun_t_metabolic         0.0192
## 13 path_bun_t_panel             0.0192
## 14 path_bun_t_complete          0.0187
## 15 radi_bun_t_chest             0.0142
## 16 radi_bun_t_neck              0.0123
## 17 radi_bun_t_exam              0.0106
## 18 path_bun_t_urinalysis        0.006 
## 19 medi_bun_t_stimulation       0.0015
## 20 path_bun_t_test              0.0009
xray_spine_n_visit %>% get_tag_density_information("radi_bun_t_thorac") %>% 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 ~ radi_bun_t_thorac"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 31.5
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 19 x 2
##    name                    correlation
##    <chr>                         <dbl>
##  1 medi_bun_t_puncture         0.130  
##  2 medi_bun_t_routine          0.130  
##  3 medi_bun_t_venipuncture     0.130  
##  4 radi_bun_t_view             0.0441 
##  5 radi_bun_t_neck             0.0416 
##  6 medi_bun_t_chiropract       0.0403 
##  7 medi_bun_t_exercise         0.0364 
##  8 radi_bun_t_shoulder         0.0314 
##  9 medi_bun_t_therapeutic      0.0275 
## 10 path_bun_t_comprehen        0.0139 
## 11 path_bun_t_metabolic        0.0092 
## 12 path_bun_t_panel            0.0092 
## 13 path_bun_t_complete         0.0085 
## 14 radi_bun_t_exam             0.0066 
## 15 path_bun_t_urinalysis       0.0059 
## 16 path_bun_t_assay            0.00290
## 17 radi_bun_t_chest            0.00240
## 18 medi_bun_t_stimulation      0.0009 
## 19 path_bun_t_test             0.0001
xray_spine_n_visit %>% get_tag_density_information("medi_bun_t_venipuncture") %>% 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_venipuncture"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 266
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 16 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 radi_bun_t_view           0.0466  
##  2 radi_bun_t_neck           0.044   
##  3 medi_bun_t_chiropract     0.0402  
##  4 medi_bun_t_exercise       0.0382  
##  5 radi_bun_t_shoulder       0.0328  
##  6 medi_bun_t_therapeutic    0.0291  
##  7 path_bun_t_complete       0.0117  
##  8 radi_bun_t_exam           0.0062  
##  9 path_bun_t_urinalysis     0.006   
## 10 path_bun_t_comprehen      0.0047  
## 11 path_bun_t_assay          0.00290 
## 12 radi_bun_t_chest          0.0028  
## 13 path_bun_t_metabolic      0.0004  
## 14 path_bun_t_panel          0.0004  
## 15 medi_bun_t_stimulation    0.000300
## 16 path_bun_t_test           0.0002
xray_spine_n_visit %>% get_tag_density_information("radi_bun_t_view") %>% 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 ~ radi_bun_t_view"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 28.4
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 13 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 medi_bun_t_exercise         0.042 
##  2 medi_bun_t_chiropract       0.0389
##  3 radi_bun_t_neck             0.0382
##  4 radi_bun_t_shoulder         0.0342
##  5 medi_bun_t_therapeutic      0.0319
##  6 path_bun_t_complete         0.0182
##  7 path_bun_t_comprehen        0.0107
##  8 path_bun_t_urinalysis       0.0062
##  9 radi_bun_t_exam             0.0052
## 10 path_bun_t_metabolic        0.0043
## 11 path_bun_t_panel            0.0043
## 12 path_bun_t_test             0.0005
## 13 medi_bun_t_stimulation      0.0004
xray_spine_n_visit %>% get_tag_density_information("medi_bun_t_chiropract") %>% 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_chiropract"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.1
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 12 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 medi_bun_t_exercise       0.0691  
##  2 medi_bun_t_therapeutic    0.0601  
##  3 radi_bun_t_neck           0.0379  
##  4 radi_bun_t_shoulder       0.0344  
##  5 medi_bun_t_stimulation    0.0239  
##  6 path_bun_t_complete       0.0181  
##  7 path_bun_t_comprehen      0.0106  
##  8 path_bun_t_urinalysis     0.0063  
##  9 radi_bun_t_exam           0.0061  
## 10 path_bun_t_metabolic      0.0041  
## 11 path_bun_t_panel          0.0041  
## 12 path_bun_t_test           0.000300
xray_spine_n_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 42.8
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 11 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 medi_bun_t_exercise         0.0393
##  2 radi_bun_t_neck             0.0383
##  3 radi_bun_t_shoulder         0.0355
##  4 path_bun_t_complete         0.0184
##  5 path_bun_t_comprehen        0.0108
##  6 medi_bun_t_stimulation      0.0064
##  7 path_bun_t_urinalysis       0.0062
##  8 radi_bun_t_exam             0.0054
##  9 path_bun_t_metabolic        0.0043
## 10 path_bun_t_panel            0.0043
## 11 path_bun_t_test             0.0004
xray_spine_n_visit %>% get_tag_density_information("medi_bun_t_exercise") %>% 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_exercise"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 46.3
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 10 x 2
##    name                   correlation
##    <chr>                        <dbl>
##  1 radi_bun_t_neck             0.0372
##  2 radi_bun_t_shoulder         0.0357
##  3 path_bun_t_complete         0.0184
##  4 path_bun_t_comprehen        0.0108
##  5 medi_bun_t_stimulation      0.0065
##  6 path_bun_t_urinalysis       0.0062
##  7 radi_bun_t_exam             0.0053
##  8 path_bun_t_metabolic        0.0044
##  9 path_bun_t_panel            0.0044
## 10 path_bun_t_test             0.0005
# xray_spine_n_visit %>% get_tag_density_information("path_bun_t_urinalysis") %>% print()
xray_spine_n_visit <- xray_spine_n_visit[path_bun_t_urinalysis==F]
xray_spine_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 9 x 2
##   name                   correlation
##   <chr>                        <dbl>
## 1 radi_bun_t_neck             0.0371
## 2 radi_bun_t_shoulder         0.0357
## 3 path_bun_t_complete         0.0184
## 4 path_bun_t_comprehen        0.0108
## 5 medi_bun_t_stimulation      0.0065
## 6 radi_bun_t_exam             0.0053
## 7 path_bun_t_metabolic        0.0044
## 8 path_bun_t_panel            0.0044
## 9 path_bun_t_test             0.0005
xray_spine_n_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 43.7
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 7 x 2
##   name                   correlation
##   <chr>                        <dbl>
## 1 radi_bun_t_neck             0.0372
## 2 radi_bun_t_shoulder         0.0357
## 3 path_bun_t_metabolic        0.0077
## 4 path_bun_t_panel            0.0077
## 5 medi_bun_t_stimulation      0.0065
## 6 radi_bun_t_exam             0.0052
## 7 path_bun_t_test             0.0005
# xray_spine_n_visit %>% get_tag_density_information("path_bun_t_metabolic") %>% print()
xray_spine_n_visit <- xray_spine_n_visit[path_bun_t_metabolic==F]
xray_spine_n_visit %>% plot_med_density() %>% print()
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 5 x 2
##   name                   correlation
##   <chr>                        <dbl>
## 1 radi_bun_t_neck             0.0371
## 2 radi_bun_t_shoulder         0.0357
## 3 medi_bun_t_stimulation      0.0065
## 4 radi_bun_t_exam             0.0052
## 5 path_bun_t_test             0.0005
xray_spine_n_visit %>% get_tag_density_information("radi_bun_t_shoulder") %>% 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 ~ radi_bun_t_shoulder"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 34.7
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 4 x 2
##   name                   correlation
##   <chr>                        <dbl>
## 1 radi_bun_t_neck           0.0412  
## 2 medi_bun_t_stimulation    0.0068  
## 3 radi_bun_t_exam           0.00480 
## 4 path_bun_t_test           0.000600
xray_spine_n_visit %>% get_tag_density_information("medi_bun_t_stimulation") %>% 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_stimulation"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 19.3
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 3 x 2
##   name            correlation
##   <chr>                 <dbl>
## 1 radi_bun_t_neck    0.0421  
## 2 radi_bun_t_exam    0.00480 
## 3 path_bun_t_test    0.000600
xray_spine_n_visit %>% get_tag_density_information("radi_bun_t_exam") %>% 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 ~ radi_bun_t_exam"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 34.1
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 2 x 2
##   name            correlation
##   <chr>                 <dbl>
## 1 radi_bun_t_neck    0.0428  
## 2 path_bun_t_test    0.000600
xray_spine_n_visit %>% get_tag_density_information("path_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 ~ path_bun_t_test"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 43.8
## $dist_plots

## 
## $stat_tables

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

xray_spine_n_visit %>% get_tag_cor() %>% print()
## Warning in stats::cor(cor_data): the standard deviation is zero
## # A tibble: 1 x 2
##   name            correlation
##   <chr>                 <dbl>
## 1 radi_bun_t_neck      0.0428
xray_spine_n_visit %>% get_tag_density_information("radi_bun_t_neck") %>% 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 ~ radi_bun_t_neck"
## Warning in leveneTest.default(y = y, group = group, ...): group coerced to
## factor.
## Picking joint bandwidth of 26.7
## $dist_plots

## 
## $stat_tables

xray_spine_n_visit_w_neck <- xray_spine_n_visit[radi_bun_t_neck==T]
xray_spine_n_visit_n_neck <- xray_spine_n_visit[radi_bun_t_neck==F]

xray_spine_n_visit_w_neck

xray_spine_n_visit_w_neck <- xray_spine_n_visit_w_neck[,`:=`(
  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)
                           )]
xray_spine_n_visit_w_neck_btbv5 <- xray_spine_n_visit_w_neck %>% btbv5_radi_fac_only(
                                  osa_group_p=NA,
                                  osa_class_p = NA,
                                  osa_specialization_p = NA
                                  ) %>% 
  .[,`:=`(procedure_type=39, 
          procedure_modifier="X-ray Spine + Neck", 
          ingest_date =Sys.Date())]
xray_spine_n_visit_w_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   17.78   48.34   74.75  169.76  249.06 3262.49
spine_q1 <- xray_spine_n_visit_w_neck_btbv5$tp_med_med %>% quantile(.25)
spine_q3 <- xray_spine_n_visit_w_neck_btbv5$tp_med_med %>% quantile(.75)
spine_iqr <- xray_spine_n_visit_w_neck_btbv5$tp_med_med %>% IQR()
xray_spine_n_visit_w_neck_btbv5_final <- xray_spine_n_visit_w_neck_btbv5[tp_cnt_cnt > 1 ]
# mri_brain_clean_btbv5_final <- mri_brain_clean_btbv5[tp_med_med < 5000 ]
xray_spine_n_visit_w_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   17.78   48.34   74.75  169.76  249.06 3262.49
xray_spine_n_visit_w_neck_btbv5_final$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   17.78   55.00   93.39  184.37  256.25 3262.49
bq_table_upload(x=procedure_dev_table, values= xray_spine_n_visit_w_neck_btbv5_final, create_disposition='CREATE_IF_NEEDED', write_disposition='WRITE_APPEND')

xray_spine_n_visit_n_neck

xray_spine_n_visit_n_neck <- xray_spine_n_visit_n_neck[,`:=`(
  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)
                           )]
xray_spine_n_visit_n_neck_btbv5 <- xray_spine_n_visit_n_neck %>% btbv5_radi_fac_only(
                                  osa_group_p=NA,
                                  osa_class_p = NA,
                                  osa_specialization_p = NA
                                  ) %>% 
  .[,`:=`(procedure_type=39, 
          procedure_modifier="X-ray Spine", 
          ingest_date =Sys.Date())]
xray_spine_n_visit_n_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   11.29   50.25   76.10  145.73  231.12  594.46
spine_q1 <- xray_spine_n_visit_n_neck_btbv5$tp_med_med %>% quantile(.25)
spine_q3 <- xray_spine_n_visit_n_neck_btbv5$tp_med_med %>% quantile(.75)
spine_iqr <- xray_spine_n_visit_n_neck_btbv5$tp_med_med %>% IQR()
xray_spine_n_visit_n_neck_btbv5_final <- xray_spine_n_visit_n_neck_btbv5[tp_cnt_cnt > 1 ]
# mri_brain_clean_btbv5_final <- mri_brain_clean_btbv5[tp_med_med < 5000 ]
xray_spine_n_visit_n_neck_btbv5$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   11.29   50.25   76.10  145.73  231.12  594.46
xray_spine_n_visit_n_neck_btbv5_final$tp_med_med %>% summary()
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   11.29   52.37   93.43  157.37  243.45  594.46
bq_table_upload(x=procedure_dev_table, values= xray_spine_n_visit_n_neck_btbv5_final, create_disposition='CREATE_IF_NEEDED', write_disposition='WRITE_APPEND')