开始学SAS，第一个小程序
SAS for Windows (v8) 统计分析系统教程新编 洪楠, 侯军编著 北京：清华大学出版社：北京交通大学出版社，2004 558页 CNY55.00 (含光盘) C819/28
SAS应用统计实验 邵建利主编 上海：上海财经大学出版社，2006 180页 CNY19.00 (含光盘) C819/38
SAS编程技术与金融数据处理 朱世武著 北京：清华大学出版社，2003 21, 491页 CNY52.00 (含光盘) F830.49/30
下面这个程序是第二章P12页的代码，其中
PLOT EDUC*OCU1=SEX;
这一句没看懂，书上的解释是：以每人的性别为图点，画出其教育年限（纵轴）与职业（横轴）交点图。这什么意思？是不是分性别的教育年限与职业曲线？没弄懂。下面运行出来的结果中也没有这一项。
另外发现，SAS中屏幕中拷贝过来的内容贴到word或者gmail编辑窗口都变形了，各列不能对齐，stata的运行结果也有这个问题，吴晓刚老师教了一招，把字体改成current new就ok了，另外字体可以弄小点，这点经验在sas中也适用。
另外，sas输出的交互表，统计图全是这种文本格式的，非常难看，虽然可以直接用一个记事本容纳所有内容，还是不够美观，不爽，不知道sas更高版本是否有更漂亮的图形。stata可以做出很漂亮的统计图形来，尤其事件史分析中的图形，看起来很专业。
第一个SAS小程序，呵呵：
data older;
input ID1 1-2 CASEID 3-5 NAME 6 SEX AGE 8-9 EDC 10 OCU1 11 OCU2 12 SAL1 13-15 SAL2 16-18 V1 19 V2 20 V3 21 V4 22 V5 23;
LIST;
CARD;
11001116031015012021214
11002116542218018011210
11030227023013012021310
;
PROC PRINT;
TITLE '老龄人社会保障抽样调查';
PROC PLOT;
PLOT EDUC*OCU1=SEX;
PROC CHART;
VBAR OCU1;
HBAR SEX;
PROC FREQ;
TABLES EDC;
TABLES EDC*SAL1;
PROC MEANS;
PROC SORT;
BY SEX EDC;
PROC PRINT;
PROC MEANS;
BY SEX EDC;
Run;
运行结果：
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 7
Obs ID1 CASEID NAME SEX AGE EDC OCU1 OCU2 SAL1 SAL2 V1 V2 V3 V4 V5
1 11 1 1 1 60 3 1 0 150 120 2 1 2 1 4
2 11 2 1 1 65 4 2 2 180 180 1 1 2 1 0
3 11 30 2 2 70 2 3 0 130 120 2 1 3 1 0
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 8
Frequency
1 ? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
? ***** ***** *****
妰儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍
1 2 3
OCU1 Midpoint
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 9
SEX Cum. Cum.
Midpoint Freq Freq Percent Percent
?
1.0 ?*************************************** 2 2 66.67 66.67
?
1.5 ? 0 2 0.00 66.67
?
2.0 ?******************* 1 3 33.33 100.00
?
妰儍儍儍儍儍儍儍儍儍垉儍儍儍儍儍儍儍儍儍?
1 2
Frequency
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 10
The FREQ Procedure
Cumulative Cumulative
EDC Frequency Percent Frequency Percent
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍
2 1 33.33 1 33.33
3 1 33.33 2 66.67
4 1 33.33 3 100.00
Table of EDC by SAL1
EDC SAL1
Frequency?
Percent ?
Row Pct ?
Col Pct ? 130? 150? 180? Total
儍儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍?
2 ? 1 ? 0 ? 0 ? 1
? 33.33 ? 0.00 ? 0.00 ? 33.33
?100.00 ? 0.00 ? 0.00 ?
?100.00 ? 0.00 ? 0.00 ?
儍儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍?
3 ? 0 ? 1 ? 0 ? 1
? 0.00 ? 33.33 ? 0.00 ? 33.33
? 0.00 ?100.00 ? 0.00 ?
? 0.00 ?100.00 ? 0.00 ?
儍儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍?
4 ? 0 ? 0 ? 1 ? 1
? 0.00 ? 0.00 ? 33.33 ? 33.33
? 0.00 ? 0.00 ?100.00 ?
? 0.00 ? 0.00 ?100.00 ?
儍儍儍儍儓儍儍儍儍垉儍儍儍儓儍儍儍儍?
Total 1 1 1 3
33.33 33.33 33.33 100.00
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 11
The MEANS Procedure
Variable N Mean Std Dev Minimum Maximum
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
ID1 3 11.0000000 0 11.0000000 11.0000000
CASEID 3 11.0000000 16.4620776 1.0000000 30.0000000
NAME 3 1.3333333 0.5773503 1.0000000 2.0000000
SEX 3 1.3333333 0.5773503 1.0000000 2.0000000
AGE 3 65.0000000 5.0000000 60.0000000 70.0000000
EDC 3 3.0000000 1.0000000 2.0000000 4.0000000
OCU1 3 2.0000000 1.0000000 1.0000000 3.0000000
OCU2 3 0.6666667 1.1547005 0 2.0000000
SAL1 3 153.3333333 25.1661148 130.0000000 180.0000000
SAL2 3 140.0000000 34.6410162 120.0000000 180.0000000
V1 3 1.6666667 0.5773503 1.0000000 2.0000000
V2 3 1.0000000 0 1.0000000 1.0000000
V3 3 2.3333333 0.5773503 2.0000000 3.0000000
V4 3 1.0000000 0 1.0000000 1.0000000
V5 3 1.3333333 2.3094011 0 4.0000000
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 12
Obs ID1 CASEID NAME SEX AGE EDC OCU1 OCU2 SAL1 SAL2 V1 V2 V3 V4 V5
1 11 1 1 1 60 3 1 0 150 120 2 1 2 1 4
2 11 2 1 1 65 4 2 2 180 180 1 1 2 1 0
3 11 30 2 2 70 2 3 0 130 120 2 1 3 1 0
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 13
-------------------------------------------- SEX=1 EDC=3 ---------------------------------------------
The MEANS Procedure
Variable N Mean Std Dev Minimum Maximum
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
ID1 1 11.0000000 . 11.0000000 11.0000000
CASEID 1 1.0000000 . 1.0000000 1.0000000
NAME 1 1.0000000 . 1.0000000 1.0000000
AGE 1 60.0000000 . 60.0000000 60.0000000
OCU1 1 1.0000000 . 1.0000000 1.0000000
OCU2 1 0 . 0 0
SAL1 1 150.0000000 . 150.0000000 150.0000000
SAL2 1 120.0000000 . 120.0000000 120.0000000
V1 1 2.0000000 . 2.0000000 2.0000000
V2 1 1.0000000 . 1.0000000 1.0000000
V3 1 2.0000000 . 2.0000000 2.0000000
V4 1 1.0000000 . 1.0000000 1.0000000
V5 1 4.0000000 . 4.0000000 4.0000000
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
-------------------------------------------- SEX=1 EDC=4 ---------------------------------------------
Variable N Mean Std Dev Minimum Maximum
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
ID1 1 11.0000000 . 11.0000000 11.0000000
CASEID 1 2.0000000 . 2.0000000 2.0000000
NAME 1 1.0000000 . 1.0000000 1.0000000
AGE 1 65.0000000 . 65.0000000 65.0000000
OCU1 1 2.0000000 . 2.0000000 2.0000000
OCU2 1 2.0000000 . 2.0000000 2.0000000
SAL1 1 180.0000000 . 180.0000000 180.0000000
SAL2 1 180.0000000 . 180.0000000 180.0000000
V1 1 1.0000000 . 1.0000000 1.0000000
V2 1 1.0000000 . 1.0000000 1.0000000
V3 1 2.0000000 . 2.0000000 2.0000000
V4 1 1.0000000 . 1.0000000 1.0000000
V5 1 0 . 0 0
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
-------------------------------------------- SEX=2 EDC=2 ---------------------------------------------
Variable N Mean Std Dev Minimum Maximum
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
ID1 1 11.0000000 . 11.0000000 11.0000000
CASEID 1 30.0000000 . 30.0000000 30.0000000
NAME 1 2.0000000 . 2.0000000 2.0000000
AGE 1 70.0000000 . 70.0000000 70.0000000
OCU1 1 3.0000000 . 3.0000000 3.0000000
OCU2 1 0 . 0 0
SAL1 1 130.0000000 . 130.0000000 130.0000000
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
老龄人社会保障抽样调查 10:58 Thursday, September 1, 2007 14
-------------------------------------------- SEX=2 EDC=2 ---------------------------------------------
The MEANS Procedure
Variable N Mean Std Dev Minimum Maximum
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
SAL2 1 120.0000000 . 120.0000000 120.0000000
V1 1 2.0000000 . 2.0000000 2.0000000
V2 1 1.0000000 . 1.0000000 1.0000000
V3 1 3.0000000 . 3.0000000 3.0000000
V4 1 1.0000000 . 1.0000000 1.0000000
V5 1 0 . 0 0
儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍儍?
--
最新奇有趣的新技术、科幻、社会学信息：
奇想录 http://www.qixianglu.cn
订 阅 http://feed.feedsky.com/woodphone
forcode科幻小说《抽水马桶的秘密》：
http://www.cmfu.com/showbook.asp?Bl_id=60008
手机阅读
0 条评论:
发表评论
订阅 博文评论 [Atom]
<< 主页
============================================================
格兰仕微波炉报价单 英语六级历年真题听力下载 大陆身份证生成器
東方神起的所有反轉劇!!.[含东方剧场](會繼續更新以後的)(已可覲看)
电视剧《靠近你温暖我》全集下载（BT/迅雷/电驴/剧照）
精彩的洞庭湖人鼠大战（4视频+forcode点评）
一百多个电影字幕下载网站，精心收集整理！
(视频)（CCTV10“走进科学”－科幻之旅专题－克隆人 8.14）
国外BT站点和BT种子搜索站（国外完整bt搜索列表）
============================================================
forcode科幻小说《抽水马桶的秘密》正在起点中文网连载
《抽水马桶的秘密》读者评论：
（1）你的书很好看，比大刘，王晋康的创意好太多了，努力吧将来出实体书我一定会买的。（2）很有想象力的作者啊！！估计是看了不少科幻小说的人，希望不要浪费你非凡的想象力。 （3） 读者在自己博客或论坛对《抽水马桶的秘密》的评论。（4）点击此处查看全部的读者评论（18页，1000条以上）。
《抽水马桶的秘密》相关帖子：
《远程面包机》提纲|| 《进化论危机》提纲大家一起来设计
抽水马桶是外星人的试管|| 读者推荐超一万票
《抽水马桶的秘密》内容简介：
地球哺乳界正在发生的一次大规模跨物种升级，DNA机制并非人们所想像的那样是决定生命的最终遗传载体，而是类似浏览器这样的转译机制，真正的遗传物质存储在弥漫整个宇宙的光子数据库中，DNA机制实际上是一种设定了进化路径的文明压缩包的解压机制，数十亿年前灭绝的三栖人发明了光子数据库和DNA机制，目的是为了让这个机制最终复活三栖人文明，而人类（裸猿）这一物种在三栖人社会里其实是一种宠物，但是DNA机制似乎出现了点问题，或者说不知道什么原因裸猿突然变得太过聪明了，在播撒了始祖菌（DNA种子）的所有星球，进化路径发展到裸猿阶段，并没有继续演化出最后一步：三栖人，而强大的光子数据库一旦意识到DNA进化机制的这个漏洞，立刻关闭了这些星球对光子数据库的访问权限，这样，这些星球的生态系统都面临着灭绝的危险，因为他们脱离了光子数据库的遗传支持再也无法自然繁殖，只能靠遗传工程来复制现有的基因，或者做些小打小闹的修改，整个宇宙各星球上的基于DNA机制的生态系统都面临崩溃的危险。最终在13世纪，裸猿一族在银河边缘一个不起眼的小星系发现了地球这个由于某种原因至今还刚进化到裸猿初级阶段的星球，为了催熟地球的进化速度，外星裸猿文明开始介入地球的发展，为了防止光子数据库察觉到非地球文明的介入并关闭地球的权限，这种介入始终是暗中进行，因为介入方式的分歧，银河系裸猿文明分裂为两大集团，这两大集团的争斗伴随着人类近现代的发展，于是，文艺复兴开始了、三次科技革命出现了、两次世界大战也来了，直到今天，地球人类为自己的技术进步而沾沾自喜，丝毫不知道技术迅速发展的真正原因以及潜藏的危机。
============================================================
forcode2003年以前的习作：未来的婚姻、远程面包机
forcode朗诵《蜀道难》||forcode的一百多个科幻构思
奇想录：最新奇有趣的新技术和新闻点评|| 订阅“奇想录”