解决实际问题之画热力图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. /*自定义标签*/
   
2. proc format;
   
3. value pdb
   
4. 1="1W2X"
   
5. 2="3H0J"
   
6. 3="3H0Q"
   
7. 4="3TV5"
   
8. 5="3TVU"
   
9. 6="3TVW"
   
10. 7="3TZ3"
    
11. 8="4WYO"
    
12. 9="4WZ8"
    
13. 10="5CTB"
    
14. 11="5CTC"
    
15. 12="5CTE";
    
16. 
    
17. data rmsd1;
    
18. format RowLab pdb.;
    
19. format ColLab pdb.;
    
20. do RowLab=1 to 12;
    
21. do ColLab=1 to 12;
    
22. input Value @@;
    
23. output;
    
24. end;
    
25. end;
    
26. cards;
    
27. 0        0.224        0.25        0.226        0.236        0.248        0.267        0.407        0.418        0.239        0.27        0.363
    
28. 0.224        0        0.246        0.23        0.242        0.225        0.303        0.442        0.459        0.245        0.257        0.41
    
29. 0.25        0.246        0        0.305        0.275        0.242        0.324        0.394        0.386        0.262        0.289        0.378
    
30. 0.226        0.23        0.305        0        0.145        0.129        0.216        0.398        0.41        0.149        0.162        0.331
    
31. 0.236        0.242        0.275        0.145        0        0.137        0.194        0.324        0.306        0.127        0.137        0.301
    
32. 0.248        0.225        0.242        0.129        0.137        0        0.224        0.34        0.326        0.122        0.137        0.301
    
33. 0.267        0.303        0.324        0.216        0.194        0.224        0        0.39        0.392        0.214        0.257        0.345
    
34. 0.407        0.442        0.394        0.398        0.324        0.34        0.39        0        0.132        0.351        0.309        0.205
    
35. 0.418        0.459        0.386        0.41        0.306        0.326        0.392        0.132        0        0.328        0.294        0.181
    
36. 0.239        0.245        0.262        0.149        0.127        0.122        0.214        0.351        0.328        0        0.14        0.296
    
37. 0.27        0.257        0.289        0.162        0.137        0.137        0.257        0.309        0.294        0.14        0        0.273
    
38. 0.363        0.41        0.378        0.331        0.273        0.301        0.345        0.205        0.181        0.296        0.273        0
    
39. ;
    
40. run;
    
41. 
    
42. data rmsd2;
    
43.    if _n_ = 1 then /* Scale matrix size by the number or rows/columns */
    
44.    call symputx('size',ceil(18 * sqrt(n)));   
    
45.    set rmsd1 nobs=n;
    
46.    Outline = ifn(Value > 1, Value, .); /* Outline the nonmissing values (Value > 1) */
    
47.    output;
    
48. run;
    
49. 
    
50. proc print data=rmsd2(obs=13);
    
51. id RowLab ColLab;
    
52. run;
    
53. 
    
54. %put &size;
    
55. 
    
56. proc template;
    
57.    define statgraph matrix;
    
58.       begingraph / designheight=%eval(250+&size) /* Size: a bit higher than wider */
    
59.                    designwidth =%eval(200+&size);
    
60.          entrytitle "RMSD Matrix";
    
61.          layout overlay / yaxisopts=(discreteopts=(tickvaluefitpolicy=none)
    
62.                                      display=(tickvalues) reverse=true)
    
63.                           xaxisopts=(discreteopts=(tickvaluefitpolicy=rotate)
    
64.                                      display=(tickvalues));
    
65.             * Heat map provides the background color for each cell;
    
66.             heatmapparm y=RowLab x=ColLab colorresponse=Value /
    
67.                ColorModel=(GraphWalls:color ThreeColorRamp:StartColor);
    
68.             * Heat map provides the outlines;
    
69.             heatmapparm y=RowLab x=ColLab colorresponse=Outline /
    
70.                ColorModel=(GraphWalls:Color ThreeColorRamp:StartColor)
    
71.                display=all includemissingcolor=false fillattrs=(transparency=1)
    
72.                outlineattrs=graphdata2(thickness=1);
    
73.             * Textplot provides the values;
    
74.             textplot y=RowLab x=ColLab text=eval(put(Value,5.3)) /
    
75.                      textattrs=(size=12px) position=center;
    
76.          endlayout;
    
77.       endgraph;
    
78.    end;
    
79. run;
    
80. 
    
81. proc sgrender data=rmsd2 template=matrix;
    
82. run;

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：标签