坊间关于SAS和R以及其他各种工具的比较和争论已经很多了,本文在此不再赘述。我们今天聊一聊一个更实际的问题:如何在SAS里调用R?
如何在SAS里调用R?简单来说,目前我知道的有三种方式,具体如下:
SAS/IML interface
%Proc_R
Base SAS Java Object
这三种方式,在使用方式、易用性以及稳定性上都有差别。第一种是借助SAS 的IML过程调用R;第二种是SAS爱好者个人Xin Wei写的一个SAS 宏程序;第三种是SAS官方发布的在Base SAS中直接调用R、Python等开源软件的开源工具。
本文示例如何通过SAS/IML interface 调用R,步骤如下:
1.检查SAS许可是否限制了系统选项。
proc options restrict;
run;
一般结果都会提示SAS options没有收到限制。
2.检查是否禁止了SAS调用R选项,一般默认是禁止的。
proc options option=Rlang;
run;
3.在SAS配置文件里添加「-Rlang」。注意默认的配置文件是英文版本,若想各个版本都能调用R,请在nls文件下的所有CFG文件里都添加「-Rlang」。
4.开启SAS调用R选项。
proc options option=RLANG;
run;
5. 可以用SAS愉快的调戏R了。
比如,玩一下R的抽样函数sample。注意其中蓝色的字是固定的语法,就像一个夹板,夹住R语句。
proc iml;
submit / R;
sample(c("Prof. Zhang","Dr. Wang" , "Brother Gu"), 1)
endsubmit;
quit;
SAS的output直接反馈的结果如下:
再比如,玩一下R的结构函数 str。
proc iml;
run ExportDataSetToR("Sashelp.Class", "df" );
submit / R;
str(df)
endsubmit;
quit;
SAS output的结果如下:
注意,其中我们用了 run ExportDataSetToR进行SAS和R之间的数据交换。类型的还有:
不过,需要留意SAS和R有版本的对应的关系。
对于在SAS调用R,我能想到的使用场景就是:
调用R里一些独特的包和函数,方便快捷的实现某些功能。
比如,你可以在SAS里试试ggplot2的画图。
proc iml;
submit/R;
library(ggplot2)
qplot(carat, price, data = diamonds, color= color)
endsubmit;
*quit;
如此,你就可以在SAS里看到酷炫的ggplot2画图效果。
2.做一些SAS和R的教学,演示和比较。
接着上面的,赶紧补一个SAS版的图,代码和效果可以有一个比较。
proc sgplot data=diamonds noborder;
scatter x=carat y=price/group=color markerattrs=(size=3);
keylegend / position=right across=1 noopaque noborder;
run;
SAS默认效果确实是丑,SAS的伙伴不要恨我。如果想画SAS 美图,可以翻翻我的「你没见过的SAS绘图系列教程」
SAS绘图是不真的又难又丑啊?
如何用SAS画正态、t、卡方以及F分布曲线?
多图预警:如何又快又美地用SAS画各种统计图形?
一图看尽丑小鸭变白天鹅全过程,手把手教你SAS画图!
怎么画一个丰满又美丽的生存分析K-M曲线图?
另外两种:请google
%PROC_R: A SAS Macro that Enables Native R Programming in the Base SAS Environment
Open Source Integration Using the Base SAS? Java Object
也许有童鞋会问,那么,有没有可以在R调用SAS的方法?好像是有的,R有一个包叫「R2SAS」,好像可以实现。但是,作者在2011年PO了一篇文章后,一直没看到包的下载地址和后续消息。
最后,按常规,关注「StatsThinking」,回复「SASCR」获取彩蛋。
如果你有问题,想获得及时回复,可以扫描下面二维码,来「分答」问我。
版权说明:
感谢支持鼓励,形式包括不限于转发,转载,推荐,赞赏。
微信公众平台可直接转载,欲获取修改白名单权限,请后台回复「51BMD」。
其他平台转载,不得省略作者信息,包括公众号二维码。
↓长按二维码,关注StatsThinking↓
联系客服