强大的SAS/GRAPH 到来了，以前的质疑和摇头都可以暂时停一停。
先示一例子，见代码和得到的图形。
data ex;
input day t h;
cards;
0 2 455
0 4 488
0 4 499
5 8 388
5 9 345
5 9 322
10 16 255
10 15 256
10 15 264
15 20 100
15 22 112
15 21 123
;
run;
proc sgplot;
vbar day / response=t stat=mean limits=both LIMITSTAT=STDERR BARWIDTH=0.4 ;
vline day / response=h stat=mean limits=both LIMITSTAT=STDERR y2axis;
yaxis min=0 max=30 ;
y2axis min=0 max=600 ;
run;
ANOTHER PIC:
Functions and CALL Routines
Returns a SAS date value of a Chinese specified traditional holiday for a specified year.
Category:
Date and Time
CNHOLIDAY(‘Chinese_holiday’, year)
is a character constant, variable, or expression that specifies one of the values listed in the following table.
Values for Chinese holiday can be in uppercase or lowercase
Chinese Traditional Holidays Recognized By SAS NOW
| CNHoliday Value | Description | Date Celebrated |
| CHUNJIE | 春节 | 正月初一 |
| YUANXIAO | 元宵节 | 正月十五 |
| ERYUEER | 二月二 龙抬头 | 农历二月二 |
| SANYUESAN | 三月三 | 农历三月三 |
| DUANWU | 端午节 | 农历五月五 |
| QIXI | 七夕 中国情人节 | 农历七月七 |
| ZHONGQIU | 中秋节 | 农历八月十五 |
| KONGZI | 孔子诞辰日 | 农历八月二十七 |
| CHONGYANG | 重阳节 | 农历九月九 |
| LABA | 腊八节 | 农历十二月八 |
| ZAOWANG | 灶王 | 农历十二月二十三 |
| XIAONIAN | 小年 | 农历十二月二十四 |
| CHUXI | 除夕 | 农历一年的最后一天 |
is a numeric constant, variable, or expression that specifies a four-digit year.
Range: 1968-2031
The CNHOLIDAY function computes the date on which a specific holiday occurs in a specified year. Only certain Chinese holidays are defined for use with this function (see Chinese Traditional Holidays Recognized By SAS NOW for a list of valid holidays) (字母拼音相对应的中文节日名也行).
The CNHOLIDAY function returns a SAS date value. To convert the SAS date value to a calendar date, use any valid SAS date format, such as the DATE9. format.
和HOLIDAY函数相比，CNHOLIDAY函数可以返回中国的传统节日。
和函数CNTOG函数相比，CNHOLIDAY函数可以返回相同的结果。
如，CNHOLIDAY(‘CHUNJIE’,2009)或者CNHOLIDAY(‘春节’,2009);返回DNTOG(2009,1,1,0)同样的结果。
options cmplib=(sasuser.datetr sasuser.cnh);
DATA _NULL_;
Fday= CNHOLIDAY("chuxi",2008);
put Fday= date9.;
run;
LOG OUTPUT:
Fday=25JAN2009
options cmplib=(sasuser.datetr sasuser.cnh);
DATA _NULL_;
Fday= CNHOLIDAY("除夕",2008);
put Fday= date9.;
run;
LOG OUTPUT:
Fday=25JAN2009
Functions:
DNTOG Function
HOLIDAY Function
Functions and CALL Routines
Category: 日期函数 Date
DNTOG(NYEAR,NMONTH,NDAY,YUNYUE);
NYEAR 中国农历年份，赋值范围：1968-2031
NONTH 中国农历月份，赋值范围：1-12
NDAY 中国农历日期，赋值范围：1-30
YUNYUE 是否农历闰月，赋值范围：0 or1 ，0表示否，1表示是。
…
代码示例1，把中国农历日期转化为公历日期
options cmplib=(sasuser.datetr);
proc fcmp;
Gday= dNtoG(1968,7,11,1);
Gday2=dNtoG(1968,7,11,0);
put Gday= date9.;
put Gday2= date9.;
run;
quit;
结果：
The FCMP Procedure
Transfered date : 03SEP1968
Transfered date : 04AUG1968
Gday=03SEP1968
Gday2=04AUG1968
代码例2，想知道2009年中国农历春节是哪一天吗？
options cmplib=(sasuser.datetr);
data _null_;
Gday=dNtoG(2009,1,1,0);
put "2009 Chinese Spring Festival day is " Gday date9.;
put "GONG XI FA CAI !";
run;
LOG中显示：
Transfered date : 26JAN2009
2009 Chinese Spring Festival day is 26JAN2009
GONG XI FA CAI !
注：如对此函数感兴趣，请邮件索取自定义函数DNTOG（）。
求解方程，在matlab中简直是小菜一碟，并且还提供数值解和解析解，功能很强大。当然，R也不错的功能强大的求解备选软件之一。
以前幻想着SAS里面也可以做得到，曾经在SASor发帖咨询，未果！【当然IML和用类似于C语言一样的算法求解除外，但不具备第四代语言的风格。】
作为号称第四代语言的SAS，现在9.2版本可以做到了：
代码示例：解方程：y=x+x**2,y=25，求x值。
proc fcmp;
/* define the function */
function inversesqrt(x);
return(x+x**2);
endsub;
y=25;
x=solve("inversesqrt", {.}, y, .);
put x;
run;
最后得到输出值：
4.5249378106
注：1.不可解方程组；2.不可解虚根，一律missing。
