请教高手如何用SAS对曲线作光滑处理？

问题：附件中X和Y的关系如图所示，曲线既有波峰，也有波谷。现要求曲线呈单调曲线或单峰曲线，不希望有双峰出现，同时光滑之后，要保证曲线的整体走向（比如上升或下降）。请教高手，此问题如何用SAS编程处理？多谢！

DATA ONE;
  INPUT X Y;
CARDS;
10 20
20 30
25 30
50 40
55 70
70 60
70 50
80 40
90 40
92 39
;
PROC IML;
  USE ONE;
  READ ALL INTO MAT;
  SMVAL=45;   
  Z=MAT(|,1|);
  S=(Z-Z(|+|)/NROW(Z))(|##|);
  SM=S##1.5 / (10##(9-(SMVAL*.1)));
   CALL SPLINEC(FITTED,C,ENDVAL,MAT,SM,,200);
  CREATE TWO VAR {NEWX NEWY };
  APPEND FROM FITTED;

DATA THREE;
  MERGE ONE TWO;
TITLE "SM45 AND SPLINE FROM IML OVERLAID";
PROC GPLOT;
   PLOT Y*X NEWY*NEWX/OVERLAY;
      SYMBOL1 I=SM45 C=BLUE V=NONE;
      SYMBOL2 I=JOIN C=RED V=NONE;
RUN;

NEWX and NEWY are  what you want.
Try in 9.1!
Bingo in sas 9.2.

1. proc sgplot;
   
2. loess x = x y = y/degree = 1 interpolation = cubic smooth = 0.75 lineattrs = (color = blue);
   
3. run;

复制代码

Jingju11玩9.2版本的很熟啊^_^.

2# jingju11
过程 proc sgplot在SAS9.1.3中没有啊，请教高手给个能在9.1.3版本下处理的程序？多谢！

symbol    i=SM25p  color=blue    ;
proc gplot;
  plot  y*x;
  run;
  quit;

注：SM表示平滑，SM加数字0到99，需调整。效果和jingju给的SAS9.2的sgplot一样。

5# sxlion
高手给出的这种方法，比较好。但是我想从平滑后的曲线中提取出平滑后的数据，又该如何做呢？谢谢！

tj0412ymy 发表于 2010-8-4 20:15
2# jingju11
过程 proc sgplot在SAS9.1.3中没有啊，请教高手给个能在9.1.3版本下处理的程序？多谢！

1. 
   
2. proc loess data=had;
   
3.       ods output OutputStatistics = HadFit(keep = x depVar Pred rename = (depVar = y Pred = smoothedY));
   
4.       model y = x / degree=2 direct smooth = 1;
   
5.    run;
   
6.    symbol1 color=black value=dot ;
   
7.    symbol2 color=blue interpol=spline w = 2 value=none;  
   
8.    proc gplot data=HadFit;
   
9.       plot (y SmoothedY) *x/overlay;
   
10.    run; quit;

复制代码

IML总令人望而生畏。一点说明：产生的光滑曲线取决于你所使用的方法

graph里面的sm（可控）和spline选项用来光滑曲线都是基于样条插值的。

其实BASE里面也有个函数可以做到，不过是需要9.2的，楼主没有，只好作罢。

sxlion 发表于 2010-8-7 00:03
DATA ONE;
  INPUT X Y;
CARDS;
10 20
20 30
25 30
50 40
55 70
70 60
70 50
80 40
90 40
92 39
;
PROC IML;
  USE ONE;
  READ ALL INTO MAT;
  SMVAL=45;   
  Z=MAT(|,1|);
  S=(Z-Z(|+|)/NROW(Z))(|##|);
  SM=S##1.5 / (10##(9-(SMVAL*.1)));
   CALL SPLINEC(FITTED,C,ENDVAL,MAT,SM,,200);
  CREATE TWO VAR {NEWX NEWY };
  APPEND FROM FITTED;

DATA THREE;
  MERGE ONE TWO;
TITLE "SM45 AND SPLINE FROM IML OVERLAID";
PROC GPLOT;
   PLOT Y*X NEWY*NEWX/OVERLAY;
      SYMBOL1 I=SM45 C=BLUE V=NONE;
      SYMBOL2 I=JOIN C=RED V=NONE;
RUN;

NEWX and NEWY are  what you want.
Try in 9.1!
Bingo in sas 9.2.

请教高手，要保证平滑后的数据为非负数，不能出现负数！谢谢！