悬赏：用SAS画出调和曲线图来！

最近看了一本 统计建模与Ｒ软件这本书（第186页3.5.3），其中一个调和曲线图不知是否可以用SAS来实现。

各位高手有兴趣的可以考虑下！


最好能将您的程序的运行结果贴上！
画出该书中的图形3.17的程序方可获得奖励！

奖励在2010.12.12日左右发放到获奖者手中，之后的回复就不算了哈~~·


统计建模与Ｒ软件 下载地址：

http://www.pinggu.org/bbs/a-344921.html

数据：




====================
ATT:以后各位牛人有好的题目可以跟我联系，各位牛人出题，本版主出论坛币，共同悬赏好的思路和程序！

知识创造价值！论坛币神马都是浮云！

程序：
data aa;
input x1 x2 x3 x4 x5;
n = _n_;
cards;
99 94 93 100 100
99 88 96 99 97   
100 98 81 96 100
93 88 88 99 96
100 91 72 96 78
90 78 82 75 97
75 73 88 97 89
93 84 83 68 88
87 73 60 76 84
95 82 90 62 39
76 72 43 67 78
85 75 50 34 37
;
run;
data aa;
set aa;
do t=-3.14 to 3.14 by 0.0001;
f=x1/sqrt(2)+x2*sin(t)+x3*cos(t)+x4*sin(2*t)+x5*cos(2*t);
output;
end;
run;
proc sgplot data =aa noautolegend;
series x = t y = f/group =n;
xaxis values =(-3.14 to 3.14 by 0.0001);
run;

大家跑跑吧！！

1# peijiamei

用SAS当然可以啊 我还做过呢 呵呵

囧 程序记得以前保留了 怎么找不到了

楼主是要这样的图吗？

等着 我来做 先去吃饭 下午来做  呵呵

呵呵 饭后我师弟会帮版主写一下这个程序

1. data had;
   
2. input politics chinese english math physics ;  
   
3. cards;
   
4. 99 94 93 100 100
   
5. 99 88 96 99 97
   
6. 100 98 81 96 100
   
7. 93 88 88 99 96
   
8. 100 91 72 96 78
   
9. 90 78 82 75 97
   
10. 75 73 88 97 89
    
11. 93 84 83 68 88
    
12. 87 73 60 76 84
    
13. 95 82 90 62 39
    
14. 76 72 43 67 78
    
15. 85 75 50 34 37
    
16. ;
    
17. data have;
    
18. set had;
    
19. array v{*} politics--physics;  
    
20. n = _n_; pi = constant('pi'); y0 = v[1] / sqrt(2);
    
21. do i = -1 to 1 by 1e-3;
    
22.   y = y0; t = pi*i;  
    
23.   do j = 2 to dim(v);
    
24.    k = int(j/2);
    
25.      y = sum(y, (mod(j,2)*cos(k*t)+(1-mod(j,2))*sin(k*t) )*v[j]);     
    
26.   end;
    
27.   output;
    
28.   end;
    
29. run;
    
30. title "Andrews' Curve";
    
31. proc sgplot data = have;
    
32. series x = t y = y/group =n;
    
33. xaxis values =(-3.2 to 3.2 by 0.8);
    
34. yaxis label = 'f(t)';
    
35. keylegend/title = 'observation number';
    
36. run;

复制代码

互相学习为上。
关于Andrews'curve的想法：
从公式来看，这条曲线是根据不同权重的转换曲线。因为第一个向量的weight很大，所以很多人建议把最重要的向量放在第一位。它的功用似乎可以鉴别观测之间的不同。因此有点factor分析的味道。不过显而易见，当观测值很多导致曲线拥挤而此图变得没有信息载量。说实话，我想不出来怎么把它用到自己的工作中去。不过楼上楼下的，应该有更好的诠释。大家都愿意聆听学习。京剧

难道jinju大大就是楼楼上的亲师弟?
阿门......

我瞎画的。不知道何谓调和曲线？也没有效率可言。或许公式也理解错了。虽然赶不上hopewell，sxlion等人，sas画图是我常用的东西。