二级分组百分比构成图

2011/08/25
By
二级分组百分比构成图1

引子: 有朋友在QQ群上问道的一个问题,这种图有个名称叫份额图(share charts),比较有代表性。注意这个是等宽的,还有种不等宽的叫mekko图。 已有临床试验数据集bylw.dose6631,一共19893条观察值, 变量有 pid、   time、  group1 、dosegroup, 分别表示病人的id号、时间点(4,8,12)、组类型(对照和干预组)和剂量类型(1,2,3)。 这时很想知道知道使用样本的构成情况,很显然用下面代码很好搞定: libname bylw "D:\SAStem"; proc gchart data=bylw.dose6631; vbar group1/ discrete type=pct group=time subgroup=dosegroup inside=subpct g100 width=12; run; quit; 得到的图如下:

Read more »

线性规划之合金原料配比问题

2011/08/16
By
合金 打铁

前言:虽然这是个比较简单的线性规划问题,但是怎么把问题转化为SAS代码,中间需要一定的小技巧,特叙此例。 问题: Metalco公司希望得到一种新的合金,其中锡40%、锌35%、铅25%。原料合金的成分如下表: 合金 原料代号 成    分 价格 锡含量 / % 锌含量 / % 铅含量 / % 美元/磅 1 60 10 30 77 2 25 15 60 70 3 45 45 10 88 4 20 50 30 84 5 50 40 10 94 工厂怎样从这几种合金原料选择来制得需要的合金,要求生产成本最小。 分析: 显然,这是个线性规划问题。 首先分析下目标,目标是生产成本最小,也就是生产单位重量(磅)的新合金用到的原料合金总价格最低。

Read more »

学习SAS的小故事2——手摇式计算机

2011/05/25
By
手摇式计算器

大学时代,最喜欢两个老师,他们都是已经退休后被返聘的老教授。其中一个就是教多元统计的余老师,他教我们多元统计,还教SAS编程。在他的一次统计课上,提到过当年农科院需要统计农作物数据的时候,用到了手摇式计算机,据说算一个路径分析,算几天才算出一个结果。手摇式计算机,也许80后可能没看过,90后可能没听说过,00后。。。

Read more »

学习SAS的小故事1——调时间

2011/05/23
By
穿越

从今天起,慢慢更新这个个人博客,讲讲学习SAS过程的一些小事,与大家分享。 记得第一次去上多元统计的SAS上机课, 那时用的SAS版本是6.12。刚坐下开好电脑, 那位和蔼的老教授扶着黑框大眼镜说:”请同学们把电脑的时间调到1996年……。”

Read more »

读书人偷书不算偷 from SFG2011

2011/05/21
By
孔乙己

这是件发生在今年SAS全球用户大会上的真事,当然我是道听途说的。 在这次的SAS大会上,SAS公司开发员Sanjay Matange and Dan Heath把他们写的一本关于SGplot画图程序的草稿,放在大会的展示室中。后来有人想去看看,结果被告知被偷了。据说这是SAS大会第一次有人偷书的事。可见这本书对这位SASor有多么大的吸引力啊,都迫不及待的到直接偷了。

Read more »

用SAS实现不等比例扇形的南丁格尔玫瑰图

2011/04/28
By
rose324

有朋友看到上一篇文章: 用SAS实现堆积式南丁格尔玫瑰图Nightingale Rose Diagram (上),然后就自己改造了一下,因为他的数据中并不是等比例的,因此对应的饼图并不是等比例的,不过由于对那段SAS代码理解不深,因此出现了错误。 后来我帮忙修正了下代码,照样也很好用。

Read more »

用SAS实现堆积式南丁格尔玫瑰图Nightingale Rose Diagram (上)

2011/02/26
By
rose283

ExcelPro的图表博客的作者刘老师最近用EXCEL做出了一个非常炫的图叫堆积式南丁格尔玫瑰图,还叫半径不等的扇形图、极区图 南丁格尔玫瑰图、风玫瑰图等等,据说最早是一个叫南丁格尔的护士最早使用的,她自己叫这个图为“鸡冠花图(coxcomb)”,一种圆形的直方图,除了讲述用EXCEL的具体做法,还对此图从专业方面给出了评价,指出了该图的优缺点和适用场合。

Read more »

SAS JMP9和IML/Studio3.2开始兼容R

2010/11/05
By

SAS公司产品JMP9和IML/Studio3.2开始兼容R,  注意不是SAS系统。 2010.10.8日Rick Wicklin在他博客中提到在SAS公司的独立模块IML/Studio3.2(可以调用SAS IML模块和使用IMPLUS语言),用IMLPLUS语言(非SAS原系统语言,在IML模块语言上进行了扩展)也可以调用R语句,当然前提是你要两者都安装了。继IML/stuidio模块开始兼容R以后,2010.10.18日,SAS公司新推出的另一个工业界颇有建树的JMP9也可以调用R软件。

Read more »

Some discriminat methodes and SAS codes

2010/10/29
By

NOTICE: The following text is one part of my published paper. Do not distribute it! All right reserved. 介绍线性判别方法、针对强共线性数据的两种降维判别方法及SAS实现代码。 1, linear discriminant analysis (LDA) Because of its simplicity and robustness, LDA has been one of the most frequently used classification techniques since 1936. LDA: proc discrim data=ex1 testdata=ex2; class...

Read more »