Site icon R-bloggers

Create SAS Code from R ‘tree’ Objects

[This article was first published on Mario Segal - Professional Site » R, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

I recently was faced with the desire to port some tree models developed in R to SAS so I could score a large database. To me this makes sense as SAS is better with large files (or at least so as not to offend anyone, I am better with large files in SAS).

I started with a web search and identified the following post:

Model decision tree in R, score in Base SAS

which basically talked about the same situation, and created code to write the SAS code. Unfortunately for me they used package party::ctree for the trees, while I use package tree.

Inspired by the great idea, I developed a code that generates the if statements to insert into  the Data Step for ‘tree’ objects developed by package tree. It is likely not the most efficient way to move through the tree, but it works (I have not dealt with binary trees in code since programming 2 in college – in Pascal of all languages – over 20 years ago). I welcome any suggestions.

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    <http://www.gnu.org/licenses/>.

tree2sas <- function (treeobj,predictor_name,filename) {

#Developed by Mario Segal. 

#initialize;
frame <- treeobj$frame
frame$printed <- 0;
frame$close <- 0;

code <- data.frame(line=””,stringsAsFactors =F)
j <- 1
count=0;
end <- row.names(frame[dim(frame)[1],])
mymax <- 0;
while (row.names(frame[i,]) != end) {
i <- mymax+1
#cycle the process until you get to the last leaf on the left from current node;
while (frame[i,]$var != ‘<leaf>’) {
code[j,] = paste(‘if’,frame[i,]$var,frame[i,]$splits[1,’cutleft’],’then do;’)
i<- i+1;
j<-j+1;
}

mymax <- max(i,mymax)
# we are at the farthest left possible, so write the assignment;
code[j,] = paste(predictor_name,’ = ‘,”‘”,frame[i,]$yval,”‘”,”;”,sep=””)
#since I printed the prediction mark as printed;
frame[i,7] <- 1;
j<- j+1

#move up until you find the first not printed. To account for the way back from right;
#check also if you are back at node #1 if not stop
while (frame[i,]$printed==1) {
if (i==1) {break}
i<-i-1
if ( frame[i,]$close==0) {
#if not not closed ,close it, I am closing leafs as I prinmt them so otherwise I duplicate;
code[j,] = ‘end;’
j <- j+1;
frame[i,8] <- 1;
}
}

if (i == 1 & mymax==dim(frame)[1]) {break}
code[j,] = paste(‘if’,frame[i,]$var,frame[i,]$splits[1,’cutright’],’then do;’)
j<- j+1;
frame[i,7] <- 1;
#the next right one here will be max+1;

}

for (k in 1:dim(code)[1]) {
cat(code[k,],sep=”\n”,append=T,file=filename)
}

}


To leave a comment for the author, please follow the link and comment on their blog: Mario Segal - Professional Site » R.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.