Setup

Load Packages

library( irs990efile )  # functions for parsing efile XMLs
library( dplyr )        # data wrangling
library( purrr )        # data wrangling 
library( pander )       # formatting
library( knitr )        # formatting
library( DT )           # table printing
library( httr )         # web requests 
library( xmltools )     # xml utilities
library( xml2 )         # xml utilities
library( XML )          # xml utilities 
library( data.tree )    # network visualization
library( networkD3 )    # network visualization
library( igraph )       # network visualization

source( "R/rdb-functions-v2.R")
source( "R/utils.R")

Data

head( index )

The Concordance File

The concordance is an RDB to XML Crosswalk that maps all xpaths onto a specified set of variables. It is located in the irs990efile package.


The full data dictionary is available on GitHub.

DATA DICTIONARY


concordance %>% 
  filter( rdb_table == table.name )  








xpath id
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/AddressLine1 1
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/AddressLine2 2
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/City 3
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/Country 4
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/PostalCode 5
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/ProvinceOrState 6
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/AddressLine1 7
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/AddressLine2 8
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/City 9
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/State 10
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/ZIPCode 11
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/AddressLine1 12
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/AddressLine2 13
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/City 14
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/Country 15
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/PostalCode 16
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/ProvinceOrState 17
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/AddressLine1 18
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/AddressLine2 19
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/City 20
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/State 21
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/ZIPCode 22
xpaths <- as.character( t.xpaths$xpath )
el <- create_edgelist_v1( xpaths )
nd <- FromDataFrameNetwork( network=el )
print( nd )
##                                    levelName
## 1  Return                                   
## 2   °--ReturnData                           
## 3       °--IRS990EZ                         
## 4           ¦--CompensationOfHighestPaidEmpl
## 5           ¦   ¦--AddressForeign           
## 6           ¦   ¦   ¦--AddressLine1         
## 7           ¦   ¦   ¦--AddressLine2         
## 8           ¦   ¦   ¦--City                 
## 9           ¦   ¦   ¦--Country              
## 10          ¦   ¦   ¦--PostalCode           
## 11          ¦   ¦   °--ProvinceOrState      
## 12          ¦   °--AddressUS                
## 13          ¦       ¦--AddressLine1         
## 14          ¦       ¦--AddressLine2         
## 15          ¦       ¦--City                 
## 16          ¦       ¦--State                
## 17          ¦       °--ZIPCode              
## 18          °--OfficerDirectorTrusteeKeyEmpl
## 19              ¦--AddressForeign           
## 20              ¦   ¦--AddressLine1         
## 21              ¦   ¦--AddressLine2         
## 22              ¦   ¦--City                 
## 23              ¦   ¦--Country              
## 24              ¦   ¦--PostalCode           
## 25              ¦   °--ProvinceOrState      
## 26              °--AddressUS                
## 27                  ¦--AddressLine1         
## 28                  ¦--AddressLine2         
## 29                  ¦--City                 
## 30                  ¦--State                
## 31                  °--ZIPCode
SetGraphStyle( nd, rankdir = "LR")
SetEdgeStyle( nd, arrowhead = "vee", color = "grey20", penwidth = 2 )
SetNodeStyle( nd, 
              style = "filled,rounded", 
              shape = "box", 
              fillcolor = "LightBlue", 
              fontname = "helvetica", 
              fontcolor="black",
              tooltip = GetDefaultTooltip )

# SetNodeStyle(acme$IT, fillcolor = "LightBlue", penwidth = "5px")
plot( nd )
# old version
find_group_names( table.name )
## [1] "//CompensationOfHighestPaidEmpl" "//OfficerDirectorTrusteeKeyEmpl"
# new version
find_table_headers( table.name )
##  [1] "//CompensationOfHighestPaidEmpl/AddressForeign"                
##  [2] "//CompensationOfHighestPaidEmpl/AddressUS"                     
##  [3] "//CompensationOfHighestPaidEmpl/AverageHoursPerWeek"           
##  [4] "//CompensationOfHighestPaidEmpl/Compensation"                  
##  [5] "//CompensationOfHighestPaidEmpl/EmployeeBenefits"              
##  [6] "//CompensationOfHighestPaidEmpl/Name"                          
##  [7] "//CompensationOfHighestPaidEmpl/Title"                         
##  [8] "//OfficerDirectorTrusteeKeyEmpl/AddressForeign"                
##  [9] "//OfficerDirectorTrusteeKeyEmpl/AddressUS"                     
## [10] "//OfficerDirectorTrusteeKeyEmpl/AvgHoursPerWkDevotedToPosition"
## [11] "//OfficerDirectorTrusteeKeyEmpl/BusinessName"                  
## [12] "//OfficerDirectorTrusteeKeyEmpl/Compensation"                  
## [13] "//OfficerDirectorTrusteeKeyEmpl/ContriToEmplBenefitPlansEtc"   
## [14] "//OfficerDirectorTrusteeKeyEmpl/Name"                          
## [15] "//OfficerDirectorTrusteeKeyEmpl/PersonName"                    
## [16] "//OfficerDirectorTrusteeKeyEmpl/Title"
# manual
if( ! is.null(table.headers) )
{table.headers}








Build the Table

# create new folder
dir.create( folder.name )
setwd( folder.name )

# create sample
set.seed( 1234 )
sample.index <- sample_n( index, sample.size )
write.csv( index, "sample-index.csv", row.names=F )
sample.urls <- sample.index$URL

# erase existing log files
file.create("XPATH-LOG.txt") 
file.create("FAIL-LOG.txt")


start.build.time <- Sys.time()    # --------------------

results.list <- list()
xpath.list <- list()


for( i in 1:length( sample.urls ) )
{
  
  url <- sample.urls[i]
  
  results.list[[i]] <- 
    build_rdb_table_v2( url, 
                        table.name, 
                        table.headers=table.headers )
  
  xpath.list[[i]]   <- 
    get_table_xpaths( url, 
                      table.name, 
                      table.headers=table.headers )
  
  if( i %% 100 == 0 ){ print(i) }
  
}


end.build.time <- Sys.time()      # --------------------


df <- dplyr::bind_rows( results.list )


fail.log <- readLines( "FAIL-LOG.txt" )
all.xpaths <- xpath.list %>% unlist() %>% unique() %>% sort()
cat( all.xpaths, sep="\n", file="XPATH-LOG.txt" )

TABLE: F9-P99-T01-COMPENSATION

## [1] "NUMBER OF RETURNS PROCESSED:   1000"
## [1] "TOTAL NUMBER OF TABLE ROWS:   606"

TOTAL RUN TIME:

## Time difference of 23.98762 mins
## Time difference of -0.3997936 hours

Write to File

write.csv( df, paste0( folder.name, "/", table.name, ".csv" ), row.names=F )
head(df) 

Are all names standardized?

names(df)
##  [1] "OBJECT_ID"                      "EIN"                           
##  [3] "NAME"                           "TAXYR"                         
##  [5] "FORMTYPE"                       "URL"                           
##  [7] "PersonName"                     "Title"                         
##  [9] "AvgHoursPerWkDevotedToPosition" "Compensation"                  
## [11] "F9_99_COMP_DTK_ADDR_L1_HCE"     "F9_99_COMP_DTK_ADDR_CITY_HCE"  
## [13] "F9_99_COMP_DTK_ADDR_STATE_HCE"  "F9_99_COMP_DTK_ADDR_ZIP_HCE"   
## [15] "ContriToEmplBenefitPlansEtc"    "ExpenseAccountOtherAllowances" 
## [17] "F9_99_COMP_DTK_ADDR_L2_HCE"








Diagnostics

Missing or Extra XPATHS

Check for xpaths that were used in the tables but do not exist in the Concordance file.

c.table <- dplyr::filter( concordance, rdb_table == table.name ) 
xpath.concordance <- c.table$xpath
setdiff( all.xpaths, xpath.concordance )
## [1] ""                                                                                        
## [2] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS"                     
## [3] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AvgHoursPerWkDevotedToPosition"
## [4] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/BusinessName"                  
## [5] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/BusinessName/BusinessNameLine1"
## [6] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/Compensation"                  
## [7] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/ContriToEmplBenefitPlansEtc"   
## [8] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/PersonName"                    
## [9] "/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/Title"

Capture to update the Concordance file.

these <- setdiff( all.xpaths, xpath.concordance )
cat( these, sep="\n", file=paste0( folder.name, "/MISSING-XPATHS.txt" ) )

Failes URLs

Check the status of failed URLs.

if( length(fail.log) > 0 )
{ 
  sapply( fail.log, get_message ) %>% 
  table() %>% 
  pander() 
} else {
  print("ALL URLs ARE VALID")
}
Client error: (404) Not Found
1








XML Table Structure

Create Network Edgelists

The create_edgelist() functions v1 and v2 are in R/utils.R.

V1 retains the root nodes Return and ReturnData.

V2 drops the root nodes and uses form versions “990” and “990EZ” as the root to visualize the sets of tables separately.

# KEEP ROOT NOTES Return & ReturnData
el <- create_edgelist_v1( all.xpaths )



Edgelist format:

V1 V2
Return ReturnData
ReturnData IRS990
IRS990 Form990PartVIISectionAGrp
Form990PartVIISectionAGrp PersonNm
Form990PartVIISectionAGrp TitleTxt
Form990PartVIISectionAGrp AverageHoursPerWeekRt
Form990PartVIISectionAGrp IndividualTrusteeOrDirectorInd
Form990PartVIISectionAGrp OfficerInd
Form990PartVIISectionAGrp ReportableCompFromOrgAmt
Form990PartVIISectionAGrp ReportableCompFromRltdOrgAmt








Visualize Tree Structure

# library( data.tree )
nd <- FromDataFrameNetwork( network=el )
print( nd )
##                                         levelName
## 1  Return                                        
## 2   °--ReturnData                                
## 3       °--IRS990EZ                              
## 4           °--OfficerDirectorTrusteeKeyEmpl     
## 5               ¦--AddressUS                     
## 6               ¦   ¦--AddressLine1              
## 7               ¦   ¦--AddressLine2              
## 8               ¦   ¦--City                      
## 9               ¦   ¦--State                     
## 10              ¦   °--ZIPCode                   
## 11              ¦--AvgHoursPerWkDevotedToPosition
## 12              ¦--BusinessName                  
## 13              ¦   °--BusinessNameLine1         
## 14              ¦--Compensation                  
## 15              ¦--ContriToEmplBenefitPlansEtc   
## 16              ¦--PersonName                    
## 17              °--Title
SetGraphStyle( nd, rankdir = "LR")
SetEdgeStyle( nd, arrowhead = "vee", color = "grey20", penwidth = 2 )
SetNodeStyle( nd, 
              style = "filled,rounded", 
              shape = "box", 
              fillcolor = "LightBlue", 
              fontname = "helvetica", 
              fontcolor="black",
              tooltip = GetDefaultTooltip )

# SetNodeStyle(acme$IT, fillcolor = "LightBlue", penwidth = "5px")
plot( nd )
# library( networkD3 )
# net <- ToDataFrameNetwork( nd, "name" )
# simpleNetwork( net, fontSize = 12 )
simpleNetwork( el, fontSize = 12 )
#plot with networkD3
radnet <- ToListExplicit( nd, unname = TRUE )
radialNetwork( radnet, fontSize = 12 )








XPATH Summary

Unique Xpaths for this RDB Table

t.xpaths <- 
  concordance %>%
  filter( rdb_table == table.name ) %>% 
  select( xpath ) %>% 
  unique() %>% 
  arrange( xpath ) %>% 
  mutate( id=row_number() ) 

t.xpaths %>% 
  kable( align="l" )
xpath id
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/AddressLine1 1
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/AddressLine2 2
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/City 3
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/Country 4
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/PostalCode 5
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressForeign/ProvinceOrState 6
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/AddressLine1 7
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/AddressLine2 8
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/City 9
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/State 10
/Return/ReturnData/IRS990EZ/CompensationOfHighestPaidEmpl/AddressUS/ZIPCode 11
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/AddressLine1 12
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/AddressLine2 13
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/City 14
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/Country 15
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/PostalCode 16
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressForeign/ProvinceOrState 17
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/AddressLine1 18
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/AddressLine2 19
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/City 20
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/State 21
/Return/ReturnData/IRS990EZ/OfficerDirectorTrusteeKeyEmpl/AddressUS/ZIPCode 22

Grouping Variables

Identify all unique groups used in xpaths.

if( is.null(table.headers) )
{ find_table_headers( table.name ) } else
table.headers
##  [1] "//CompensationOfHighestPaidEmpl/AddressForeign"                
##  [2] "//CompensationOfHighestPaidEmpl/AddressUS"                     
##  [3] "//CompensationOfHighestPaidEmpl/AverageHoursPerWeek"           
##  [4] "//CompensationOfHighestPaidEmpl/Compensation"                  
##  [5] "//CompensationOfHighestPaidEmpl/EmployeeBenefits"              
##  [6] "//CompensationOfHighestPaidEmpl/Name"                          
##  [7] "//CompensationOfHighestPaidEmpl/Title"                         
##  [8] "//OfficerDirectorTrusteeKeyEmpl/AddressForeign"                
##  [9] "//OfficerDirectorTrusteeKeyEmpl/AddressUS"                     
## [10] "//OfficerDirectorTrusteeKeyEmpl/AvgHoursPerWkDevotedToPosition"
## [11] "//OfficerDirectorTrusteeKeyEmpl/BusinessName"                  
## [12] "//OfficerDirectorTrusteeKeyEmpl/Compensation"                  
## [13] "//OfficerDirectorTrusteeKeyEmpl/ContriToEmplBenefitPlansEtc"   
## [14] "//OfficerDirectorTrusteeKeyEmpl/Name"                          
## [15] "//OfficerDirectorTrusteeKeyEmpl/PersonName"                    
## [16] "//OfficerDirectorTrusteeKeyEmpl/Title"

Groups can then be used to extract the 1:M RDB tables:

if( is.null(table.headers) )
{ table.headers <- find_table_headers( table.name ) }
demo.url <- "https://nccs-efile.s3.us-east-1.amazonaws.com/xml/201913079349300206_public.xml"
doc <- xml2::read_xml( demo.url )
xml2::xml_ns_strip( doc )
df.demo <- get_table_v2( doc, table.name, table.headers )
df.demo 
## NULL

Standardize Variable Names

Mapping of xml variable names to concordance variable names:

get_var_map( table.name=table.name )

A table with proper names:

v.map <- get_var_map( table.name=table.name )

if( ! is.null(df.demo) )
{re_name( df.demo, v.map ) %>% as.data.frame()}