SAS-R Integration Example: Transform in R, Analyze in SAS!

SAS-R integration example that I’m giving here is a small and simple one, but I think it represents a model as we proceed into the future of data science. It showcases the fast transformation capabilities of R, and the amazing analytic capabilities of SAS!

SAS-R Integration Example … Because I Was in a Hurry!

Normally, I would not do half a project in R and half a project in SAS! Here is a big picture overview of the problem I was faced with, and how I solved it with this SAS-R integration example:

I had been developing an analytic dataset from the Behavioral Risk Factor Surveillance Survey (BRFSS) for an analysis of just a subset of records (diabetes patients) – so I didn’t need the population weight variables from BRFSS. I was doing all this in R GUI.
While we were working on the diabetes paper, the person I was working with got an invitation in her e-mail to submit a paper to a journal on a different topic (healthcare quality) with a very close due date! I realized with just a tweak to my R GUI code, I could redo the current analytic dataset about diabetes that I had been working on, and set up a weighted analysis on the access to care topic!
I quickly created the access-to-care analytic dataset in R GUI! This accounts for the transformation steps (extract-transform-load, or ETL). But then when I tried to run the weighted analysis in R, it wouldn’t run. It just hung!
So, I exported a “skinny” dataset out of R GUI into *.csv format. Next, I uploaded it into the free online environment SAS OnDemand for Academics (ODA), which is new!
I was able to use SAS ODA to convert the *.csv to a SAS dataset, and run PROC SURVEYFREQ on it to get my weighted results!

I’m always trying to show learners a SAS-R integration example that showcases using each tool – SAS and R – optimally for what they are good at. That’s why took the effort of sharing this SAS-R integration example with you. I think it is a perfect illustration of the kind of thinking needed when developing modern SAS data warehouses, as I describe in my book, “Mastering SAS Programming for Data Warehousing”.

SAS-R Integration Example Starts with ETL in R

Transforming data is much simpler in R than SAS, because in SAS, you have to use data step sublanguage which is very complicated. As I recommend in my LinkedIn Learning course in data curation, I first developed a data dictionary for my project that specified the variables in my analytic dataset.

As you can see from the main data dictionary below, the only variables I planned to get from the native dataset were administrative variables having to do with the weighting I needed to do (SEQNO, X_PSU, X_STSTR, X_LLCPWT, and X_STATE), and the access to care variables for the analysis (HLTHPLN1 and MEDCOST).

This is a data dictionary main page for the SAS-R integration example given here. It shows both the native variables from the original dataset, as well as the planned transformed variables.

You can also see in the main dictionary that I plan to generate two flag variables in R for analysis: NOPLAN (flag for no health insurance plan), and MEDCOSTFLAG (flag for not having enough money to see a doctor in the last 12 months). Those tabs from the data dictionary are below.

This is the MEDCOST flag documentation. We see how MEDCOST is formatted, and we develop MEDCOST flag as a binary flag.

This is documentation of the native HLTHPLN1 variable from BRFSS, and how we plan to transform it into binary flag called NOPLAN.

SAS-R Integration Example Data and Code

The original dataset is from the BRFSS 2020 and is in SAS XPT format. You download it from this page, and need to put this huge file in a folder where your computer can access it.

The rest of the data and the code are located on my Github here. The resulting data files are big, so I only included the RDS so you could verify your work.

As you can see with the R GUI code below, the first step is to use the foreign package to convert that bad boy XPT file into an RDS so R can handle it with ease:

# call library(foreign) or load "foreign" package
#that way, we can load the SAS XPT file into R

library(foreign)

#Read XPT directly off of my computer
#and convert it into a dataframe named BRFSS_a.
#You will have to edit this code to add the name of the folder
#where the XPT file is on your computer.

BRFSS_a <- read.xport("C:/NAME OF FOLDER WHERE XPT FILE IS/LLCP2020.xpt")

#check number of columns and rows

colnames(BRFSS_a)
nrow(BRFSS_a)

#export as RDS

saveRDS(BRFSS_a, file = "BRFSS_a.rds")

You will see I import the XPT and export it using naming conventions from my SAS book, which is to advance an alpha suffix: BRFSS_a. Then, in the next code, I read in that BRFSS_a dataset, and start by making a vector called BRFSSVarList with the list of variables from my data dictionary that I know I want to keep:

#read in BRFSS_a

BRFSS_a <- readRDS(file = "BRFSS_a.rds")
colnames(BRFSS_a)

#define varlist
BRFSSVarList <- c("SEQNO", "HLTHPLN1", "MEDCOST", "X_STATE", 
"X_PSU", "X_STSTR", "X_LLCPWT")

#subset by varlist
BRFSS_b <- BRFSS_a[BRFSSVarList]

#check columns
colnames(BRFSS_b)

#check rows
nrow(BRFSS_b)

#export as RDS
saveRDS(BRFSS_b, file = "BRFSS_b.rds")

Next, I trim off the variables I don’t need from BRFSS_a by using the BRFSSVarList vector, making the dataframe BRFSS_b – which is now very skinny! I export that as an RDS file.

Referring back to my data dictionary, I realize I need to transform two variables: MEDCOSTFLAG, and NOPLAN, which are both binary flags. I will use that value that equals “1” to designate the numerator of my weighted population percentage in my analysis.

You will see in my transformation code, I read in BRFSS_b, do the transformations (creating the flag variables), advance the transformed dataset to BRFSS_c. Then, to trim off unneeded variables and create analytic datasets in both RDS and CSV format, I first advance to BRFSS_d and then export.

#read in RDS of BRFSS_b

BRFSS_b <- readRDS(file = "BRFSS_b.rds")

#Create flag for no health insurance

BRFSS_b$NOPLAN <- 0
BRFSS_b$NOPLAN [BRFSS_b$HLTHPLN1 !=1 ] <- 1
table(BRFSS_b$NOPLAN, BRFSS_b$HLTHPLN1)

#Create flag for barrier due to cost 

BRFSS_b$MEDCOSTFLAG <- 0
BRFSS_b$MEDCOSTFLAG[BRFSS_b$MEDCOST == 1] <- 1
table(BRFSS_b$MEDCOSTFLAG, BRFSS_b$MEDCOST, useNA = c("always")

#advance suffix of data frame
BRFSS_c <- BRFSS_b

#trim off unneeded variables

keep_in_analytic <- c("SEQNO", "X_PSU", "X_STSTR", "X_LLCPWT", "X_STATE",
"NOPLAN", "MEDCOSTFLAG")

BRFSS_d <- BRFSS_c[keep_in_analytic]

#export as an analytic dataset
#in CSV and in RDS

colnames(BRFSS_d)
write.csv(BRFSS_d, file = "access_analytic.csv", row.names = FALSE)
saveRDS(BRFSS_d, file = "access_analytic.rds")

You can post your code on GitHub so everyone can share.

Watch our SAS integration recorded livestreams for use-cases you can apply in your own context!

Doing this strategic renaming of dataframes as we transform creates rollback points in ETL code. The final analytic dataset in CSV format is the one we will put in SAS ODA.

SAS-R Integration Example – Moving on to Analytics in SAS

Now we turn to SAS! First, please log into SAS ODA, and upload the files to the right places, as shown in the diagram below. As you can see, you will need to make a folder named X, and we will be mapping a LIBNAME to it. You will need to use the import utility to get the CSV into SAS ODA and converted into a SAS file.

This graphic shows you where to put the files you upload into SAS ODA. You will need to create folders as directories and develop LIBNAMES to map the directory. You can use the upload utility to upload CSV data into the SAS ODA environment.

If you are having trouble with this part, please take my free course in using SAS ODA which will help you get set up in SAS ODA. To learn how to get a CSV into SAS, take the second course in the series, “Getting Data into SAS ODA”.

Let’s start by looking at our import code. Because we are transforming data in R and making the analytic dataset into a CSV, we are able to make it light and skinny. However, SAS can’t use the CSV until it converts it with PROC IMPORT. This is a little confusing, in that my course “Getting Data into SAS ODA” teaches you how to get the CSV into the environment, and this code shows you how to convert it to a SAS dataset so SAS can use it.

/*Here, I set LIBNAME to point to folder in SAS ODA*/
/*I made earlier, where I will put data and code for this*/

LIBNAME X "/home/dethwench0/X";
run;

/*Earlier, I used the upload utility to get the CSV into SAS.*/
/*Here, I identify uploaded CSV file as REFFILE to SAS.*/

FILENAME REFFILE '/home/dethwench0/X/access_analytic.csv';

/*PROC IMPORT code: I used the SAS ODA import utility to help me write this code.*/
/*This code converts the REFFILE to a SAS dataset called "analytic"*/
/*and puts it in the X folder mapped to the LIBNAME.*/

PROC IMPORT DATAFILE=REFFILE
DBMS=CSV
OUT=x.analytic;
GETNAMES=YES;
RUN;

/*Here we can review contents of imported SAS dataset*/

PROC CONTENTS DATA=x.analytic;
RUN;

This code outputs the converted SAS dataset to the folder designated by the X LIBNAME, and names it simply analytic. If you are following along, please make sure you edit the example code to point to the folder in your instance of SAS ODA. We also run a PROC CONTENTS to verify that our skinny dataset was uploaded and converted to SAS format successfully. Here is some of our PROC CONTENTS output:

This PROC CONTENTS output shows the transformed and native variables in the dataset. This is the analytic dataset we made in R that is now ready for analysis in SAS.

Also, after you edit and run that code successfully, you will notice that the converted SAS file analytic.sas7bdat will show up in your X folder.

Leveraging SAS’s Analytics Capabilities in SAS-R Integration Example

Now we get to the part where we showcase SAS’s analytics abilities. As I said before, the CDC offers guidance on how to do a weighted analysis of BRFSS data in both R and SAS. When I originally tried this in R GUI, I used the code provided by the CDC, and it just was not efficient enough on my system to run. It hung!

That’s the problem with R – it can’t handle big data analytics like SAS can. In fact, I was lucky to be able to do this SAS-R integration example, because the source BRFSS dataset was big, but not too big. If it had been too big to get into R GUI, I would definitely have not been able to put it in SAS ODA. I would have probably needed a SQL server instance at that point!

But now that we could provide SAS a manicured analytic dataset, it handled our weighted frequency analysis with ease in PROC SURVEYFREQ:

/*just to make sure LIBNAME is mapped*/
LIBNAME X "/home/dethwench0/X";
run;

/*Running PROC SURVEYFREQ on access to care variables*/
/*in "skinny" analytic dataset prepared in R and uploaded as CSV*/
/*Code provided by BRFSS/CDC for modification. That way, I knew*
/*I put the stratum variables in the right options.*/

proc surveyfreq
	data = x.analytic;
/* Use common weight because both questions are from the core section */
	weight x_llcpwt;
	strata x_ststr;
	cluster x_psu;
table noplan;
table medcostflag;
run;

I ran this code and got the output I needed from PROC SURVEYFREQ – which was the weighted estimates. I show you below:

Here are the results from PROC SURVEYFREQ in SAS when doing a weighted analysis. You will see the unweighted and weighted frequencies, and the weighted percentages, along with other calculations and estimates.

Remember, if you are getting started with either R or SAS, or just data science in general, take my courses on LinkedIn Learning. They are great for beginners new to data science who have experience in other fields.

Updated December 26, 2021. Added videos, updated links, and added banners on March 4, 2024.

Read all of our data science blog posts!

You can make mean, median and mode with a continuous variable.

Make Categorical Variable Out of Continuous Variable

Make categorical variables by cutting up continuous ones. But where to put the boundaries? Get [...]

04
Nov

Remove Rows in R with the Subset Command

Remove rows by criteria is a common ETL operation – and my blog post shows [...]

28
Oct

The open government data movements in many countries have resulted in government data being available online.

CDC Wonder for Studying Vaccine Adverse Events: The Shameful State of US Open Government Data

CDC Wonder is an online query portal that serves as a gateway to many government [...]

24
Jun

If you are a data scientist, you may want to do statistics, but you may also be interested in machine learning and artificial intelligence.

AI Careers: Riding the Bubble

AI careers are not easy to navigate. Read my blog post for foolproof advice for [...]

16
Jun

The Black Friday Death Count database has a list of news reports of deaths or injuries on Black Friday.

Descriptive Analysis of Black Friday Death Count Database: Creative Classification

Descriptive analysis of Black Friday Death Count Database provides an example of how creative classification [...]

27
Nov

What if you have too many categories in a categorical variable? Your cardinality is too high for a chi-square analysis.

Classification Crosswalks: Strategies in Data Transformation

Classification crosswalks are easy to make, and can help you reduce cardinality in categorical variables, [...]

24
Nov

Want to learn more about pharmacy data? You can use adverse event data in a data science portfolio project.

FAERS Data: Getting Creative with an Adverse Event Surveillance Dashboard

FAERS data are like any post-market surveillance pharmacy data – notoriously messy. But if you [...]

22
Nov

If you work on a big data project with multiple source datasets, you run the risk of forgetting exactly how you blended them together.

Dataset Source Documentation: Necessary for Data Science Projects with Multiple Data Sources

Dataset source documentation is good to keep when you are doing an analysis with data [...]

09
Nov

In base R, you can execute SQL-like joins, as long as you use the correct code syntax.

Joins in Base R: Alternative to SQL-like dplyr

Joins in base R must be executed properly or you will lose data. Read my [...]

07
Nov

If you are interested in population-level surveillance data, you might have thought about using NHANES data in portfolio projects.

NHANES Data: Pitfalls, Pranks, Possibilities, and Practical Advice

NHANES data piqued your interest? It’s not all sunshine and roses. Read my blog post [...]

02
Nov

When using big data, you will want to make visualizations. How do you use color to the greatest communicative advantage?

Color in Visualizations: Using it to its Full Communicative Advantage

Color in visualizations of data curation and other data science documentation can be used to [...]

09
Oct

The defaults in PowerPoint are really set up for making presentations, not data visualizations.

Defaults in PowerPoint: Setting Them Up for Data Visualizations

Defaults in PowerPoint are set up for slides – not data visualizations. Read my blog [...]

06
Oct

Adding text and arrows to diagrams can help your audience navigate the image, and understand what you are trying to communicate.

Text and Arrows in Dataviz Can Greatly Improve Understanding

Text and arrows in dataviz, if used wisely, can help your audience understand something very [...]

04
Oct

If you use good judgment in choosing chapes and images to add to your data visualizations, your audience will be enlightened.

Shapes and Images in Dataviz: Making Choices for Optimal Communication

Shapes and images in dataviz, if chosen wisely, can greatly enhance the communicative value of [...]

01
Oct

When you use a data analysis program like R or SAS, you often have to do some data editing. It can be difficult because the software was intended for calculations, not transformation.

Table Editing in R is Easy! Here Are a Few Tricks…

Table editing in R is easier than in SAS, because you can refer to columns, [...]

24
Aug

Logistic regression calculate the log odds of the probability of the outcome. Many people are used to using SAS for logistic regression, but you can also use R.

R for Logistic Regression: Example from Epidemiology and Biostatistics

R for logistic regression in health data analytics is a reasonable choice, if you know [...]

1 Comments

17
Aug

Did you know it is possible to integrate SAS with other data environments, like Microsoft SQL or Excel?

Connecting SAS to Other Applications: Different Strategies

Connecting SAS to other applications is often necessary, and there are many ways to do [...]

30
Jul

Are you a data scientist who is interested in doing independent portfolio projects to sharpen your skills? Then I strongly suggest you get a coach or a mentor.

Portfolio Project Examples for Independent Data Science Projects

Portfolio project examples are sometimes needed for newbies in data science who are looking to [...]

13
Jul

If you are a health data analyst or a biostatistician, we might find computer programmers and application developers use different terminology for the same ideas and concepts.

Project Management Terminology for Public Health Data Scientists

Project management terminology is often used around epidemiologists, biostatisticians, and health data scientists, and it’s [...]

25
Jun

If you work on front-ends or back-ends of health applications, you are probably already familiar with the concepts of Agile and rapid application development.

Rapid Application Development Public Health Style

“Rapid application development” (RAD) refers to an approach to designing and developing computer applications. In [...]

22
Jun

Data systems started being in use in the 1960s and 1970s, but these were flat systems, usually using IBM mainframes.

Understanding Legacy Data in a Relational World

Understanding legacy data is necessary if you want to analyze datasets that are extracted from [...]

19
Jun

How the front-end and back-end are connected can impact how data are stored in the application. So if you extract the data, you can have data quality problems caused by the front-end.

Front-end Decisions Impact Back-end Data (and Your Data Science Experience!)

Front-end decisions are made when applications are designed. They are even made when you design [...]

16
Jun

Slow queries can happen in SAS, R, Python, SQL or any database language. These slow queries have a cost.

Reducing Query Cost (and Making Better Use of Your Time)

Reducing query cost is especially important in SAS – but do you know how to [...]

13
Jun

If you need data to do a project, read this blog post for information.

Curated Datasets: Great for Data Science Portfolio Projects!

Curated datasets are useful to know about if you want to do a data science [...]

31
May

Public health, artificial intelligence, and data science trivia! Fun! Educational! Test your knowledge!

Statistics Trivia for Data Scientists

Statistics trivia for data scientists will refresh your memory from the courses you’ve taken – [...]

12
Apr

When working in data science, there are some tips and tricks to managing your communication and relationship with superiors that can help you advance in your career.

Management Tips for Data Scientists

Management tips for data scientists can be used by anyone – at work and in [...]

23
Mar

REDCap mess on your hands? The REDCap designers made the application so loosey goosey, you can really program yourself into a messy corner if you don't plan well.

REDCap Mess: How it Got There, and How to Clean it Up

REDCap mess happens often in research shops, and it’s an analysis showstopper! Read my blog [...]

07
Mar

If you are an aspiring data scientist, you will need to know how GitHub works. You will probably want to use it for your projects.

GitHub Beginners in Data Science: Here’s an Easy Way to Start!

GitHub beginners – even in data science – often feel intimidated when starting their GitHub [...]

28
Feb

This blog post shows you how to properly document your extract, transform, and load code.

ETL Pipeline Documentation: Here are my Tips and Tricks!

ETL pipeline documentation is great for team communication as well as data stewardship! Read my [...]

22
Feb

How do you measure how long it takes for code to run in different programs? And why would you want to measure something like that? Mainly, the reason to benchmark runtime is so that you can figure out how to optimize your code.

Benchmarking Runtime is Different in SAS Compared to Other Programs

Benchmarking runtime is different in SAS compared to other programs, where you have to request [...]

13
Dec

What is an end-to-end AI pipeline? And why are academics so bad at making one? These are different ideas we will examine in this blog post.

End-to-End AI Pipelines: Can Academics Be Taught How to Do Them?

End-to-end AI pipelines are being created routinely in industry, and one complaint is that academics [...]

25
Nov

There are different ways to refer to variables in R dataframes. You can use a field names, and you can also use field numbers.

Referring to Columns in R by Name Rather than Number has Pros and Cons

Referring to columns in R can be done using both number and field name syntax. [...]

17
Oct

The paste command is used to concatenate strings in R. You can use it different ways, which is what I demonstrate in my blog and videos.

The Paste Command in R is Great for Labels on Plots and Reports

The paste command in R is used to concatenate strings. You can leverage the paste [...]

12
Oct

You do not need to use the default R colors on your plot. You don't even need to limit yourself to named colors on cheat sheets.

Coloring Plots in R using Hexadecimal Codes Makes Them Fabulous!

Recoloring plots in R? Want to learn how to use an image to inspire R [...]

11
Oct

Error bars on plots can provide the audience an estimate of the amount of certainty you have with your estimates.

Adding Error Bars to ggplot2 Plots Can be Made Easy Through Dataframe Structure

Adding error bars to ggplot2 in R plots is easiest if you include the width [...]

11
Oct

AI on the edge refers to doing the AI processing and equations at the site of the object collecting the data.

AI on the Edge: What it is, and Data Storage Challenges it Poses

“AI on the edge” was a new term for me that I learned from Marc [...]

06
Jun

How do you make a pie chart in ggplot2 package in R? It's not that obvious

Pie Chart ggplot Style is Surprisingly Hard! Here’s How I Did it

Pie chart ggplot style is surprisingly hard to make, mainly because ggplot2 did not give [...]

20
Apr

Time series plots can be customized if you use package ggplot2 in R. You can place labels and configure axes.

Time Series Plots in R Using ggplot2 Are Ultimately Customizable

Time series plots in R are totally customizable using the ggplot2 package, and can come [...]

11
Apr

It is possible to use data curation to solve the problem of a confusion vector containing options.

Data Curation Solution to Confusing Options in R Package UpSetR

Data curation solution that I posted recently with my blog post showing how to do [...]

08
Apr

If you are having trouble setting options using R making plots, then you should read this blog post.

Making Upset Plots with R Package UpSetR Helps Visualize Patterns of Attributes

Making upset plots with R package UpSetR is an easy way to visualize patterns of [...]

4 Comments

03
Apr

There are two main ways to make box plots in R, and this blog post shows you how, and explains the differences.

20
Sep

This is the forrest plot resulting from analysis with open source statistical software R using package rmeta.

14
Oct

The Likert package in R can visualize categorical data.

Doing Surveys? Try my R Likert Plot Data Hack!

I love the Likert package in R, and use it often to visualize data. The [...]

2 Comments

11
Oct

Epidemiologic Curve of 2015 Middle East Respiratory Virus Outbreak Using R EpiCurve Package

I Used the R Package EpiCurve to Make an Epidemiologic Curve. Here’s How It Turned Out.

With all this talk about “flattening the curve” of the coronavirus, I thought I would [...]

14
Mar

Little table of X Y Pairs with a Regresion Diagram with Least Squares Line

Which Independent Variables Belong in a Regression Equation? We Don’t All Agree, But Here’s What I Do.

During my failed attempt to get a PhD from the University of South Florida, my [...]

24
Aug

SAS-R Integration Example … Because I Was in a Hurry!

got questions

about data science?

SAS-R Integration Example Starts with ETL in R

make gorgeous

spreadsheets!

SAS-R Integration Example Data and Code

learn the basics

of analyzing data from applications

Watch our SAS integration recorded livestreams for use-cases you can apply in your own context!

SAS-R Integration Example – Moving on to Analytics in SAS

make a

data science portfolio project

Leveraging SAS’s Analytics Capabilities in SAS-R Integration Example

Monika Wahi

Leave a Reply Cancel reply