setwd("~/Desktop/ar_2013")
library(haven)
library(readr)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
CCR <- read_csv("cc2013NIS.csv")
## Parsed with column specification:
## cols(
## `'HOSP_NIS'` = col_character(),
## `'YEAR'` = col_double(),
## `'CCR_NIS'` = col_double(),
## `'WAGEINDEX'` = col_double()
## )
ccr<-data.frame(as.numeric(substr(CCR[[1]],2,6)),CCR[[3]])
colnames(ccr)<-c("HOSPID","ratio")
head(ccr)
## HOSPID ratio
## 1 80235 0.241
## 2 80256 0.288
## 3 80283 0.177
## 4 80291 0.939
## 5 80069 0.409
## 6 80320 0.531
head(CCR)
## # A tibble: 6 x 4
## `'HOSP_NIS'` `'YEAR'` `'CCR_NIS'` `'WAGEINDEX'`
## <chr> <dbl> <dbl> <dbl>
## 1 '80235' 2013 0.241 0.903
## 2 '80256' 2013 0.288 0.899
## 3 '80283' 2013 0.177 0.935
## 4 '80291' 2013 0.939 0.892
## 5 '80069' 2013 0.409 0.933
## 6 '80320' 2013 0.531 0.898
ar_ahal_13 <- read_sas("ar_sidc_ahal_2013.sas7bdat")
head(ar_ahal_13)
## # A tibble: 6 x 7
## AHAID DSHOSPID HOSPID COMMUNITY_NONREHAB_NONLTAC HFIPSSTCO HOSPST YEAR
## <chr> <chr> <dbl> <dbl> <dbl> <chr> <dbl>
## 1 6710023 ADFS 5018 0 5131 AR 2013
## 2 6710530 ADOC 5117 NA NA AR 2013
## 3 6710034 ADSR 5118 0 5145 AR 2013
## 4 6710011 ADVA 5020 0 5051 AR 2013
## 5 6710033 ALLE 5119 0 5119 AR 2013
## 6 6710010 AMER 5021 0 5119 AR 2013
Core
Ar_core_13<- read_sas("ar_sidc_core_2013.sas7bdat")
Ar_core_13$DSHOSPID[1:20]
## [1] "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS"
## [11] "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS" "ADFS"
Charges
ar_chgs_13<- read_sas("ar_sidc_chgs_2013.sas7bdat")
head(ar_chgs_13)
## # A tibble: 6 x 5
## CHARGE KEY RATE REVCODE UNITS
## <dbl> <dbl> <dbl> <chr> <dbl>
## 1 31198 5.20e13 821 0111 38
## 2 12146. 5.20e13 25.5 0250 477
## 3 320. 5.20e13 3.95 0250 81
## 4 5571. 5.20e13 143. 0250 39
## 5 1291. 5.20e13 430. 0250 3
## 6 161. 5.20e13 5.55 0250 29