goptions reset=all; Title " "; Footnote " "; axis1 label=(font=swiss h=1.5 angle=90); *******************************************; data fabric; input Y @@; location+.25; label location = "Distance from fabric end in inches"; pi = constant("pi"); s1 = sin(2*pi/32); c1=cos(2*pi/32); s2 = sin(2*pi/10); c1=cos(2*pi/10); cards; 9.6522 9.2729 10.2347 10.9434 11.1482 11.0472 10.2253 10.9372 10.7710 10.8765 10.4798 10.9468 9.4090 9.0172 9.1186 9.9570 9.2555 9.0580 9.6437 8.1358 9.7593 7.8261 9.1227 10.0401 11.1420 9.4384 10.0978 10.8139 8.9794 12.9167 13.2704 10.0762 7.7174 10.6704 10.7988 9.1012 10.5630 9.4283 11.6136 10.6398 9.0910 9.8009 9.4473 9.3971 9.4226 9.1219 9.2276 11.7425 9.6363 8.3284 10.9591 10.7108 10.6995 10.0938 10.2316 9.3275 10.2522 7.7255 9.3808 7.6761 10.3801 11.7105 7.8500 8.4237 9.5786 9.3868 9.5714 11.6126 9.7259 10.6922 11.1485 10.7860 10.2519 11.6800 9.4545 8.4145 9.9007 8.2529 9.8603 8.6640 10.6281 8.7543 9.6294 10.0341 10.3447 11.1238 10.0826 10.7192 11.8742 8.3195 10.0983 10.1450 9.8657 11.7137 10.3360 10.5882 10.7266 8.3389 9.6894 9.3531 8.6141 8.8263 10.1913 9.8914 10.0110 10.2276 10.5815 10.4273 13.2196 10.3786 10.9648 9.5913 10.6570 7.5938 6.9137 10.3719 8.4523 6.4303 9.1602 10.3286 8.6877 10.3651 8.7841 9.8956 8.2851 9.5538 8.9536 9.9543 11.6016 11.6509 10.2554 11.0539 10.7543 9.2623 10.0903 11.3440 9.5806 11.8309 10.7510 8.8843 8.9535 9.9128 9.0117 11.0192 10.3225 9.9886 11.2521 10.2255 11.2754 10.8223 10.6520 10.1185 10.1605 10.6420 8.2162 8.9028 11.2344 10.9330 8.1114 10.1532 8.7317 9.8679 8.9660 10.8063 10.0091 10.5798 9.8049 12.5082 11.5374 8.9669 10.5180 12.5141 10.6083 10.8392 10.3839 10.0463 9.6511 11.5636 8.5517 9.0417 10.6874 10.1072 9.1569 9.5618 8.5188 10.2297 8.7976 9.7678 9.9599 11.0236 9.9843 13.0321 9.4554 8.8452 9.9885 10.0929 9.6221 10.6091 10.0960 11.1333 9.5918 9.2733 8.9020 9.3018 11.5969 11.4310 11.6446 10.6139 10.2448 10.1650 10.6837 9.1776 10.0802 9.3467 10.7734 8.8215 10.5727 9.1634 8.2465 8.7384 7.2661 9.0559 10.9496 9.0294 11.7795 10.4594 9.6836 10.6138 8.2314 12.1232 11.3560 10.3465 11.2749 11.0889 10.0561 9.7273 9.0754 8.5864 9.4325 7.6567 9.7680 9.1855 9.4966 9.2977 10.3489 8.5471 9.8378 9.6257 11.9055 9.0645 11.0092 9.9594 10.1208 10.1482 9.4279 9.4084 9.4182 8.3468 11.2017 9.8725 9.1876 11.6496 11.4896 9.7804 11.4930 10.7894 10.9897 10.8671 11.2747 11.0493 8.0933 10.8397 10.6629 9.6363 9.6731 8.2106 8.8985 9.5183 9.6916 8.0730 8.3481 9.0504 8.6459 8.9209 9.9356 9.8471 10.0948 9.9655 10.1336 10.5903 10.7610 11.2472 9.0259 9.4563 10.7706 9.1227 10.8688 11.2718 11.6350 9.0087 11.4075 11.1843 9.8823 10.0783 10.5793 9.4220 10.3141 10.0179 9.6195 10.2122 9.9894 8.0515 9.6481 9.8099 9.9093 10.5282 11.4239 8.3169 11.2628 8.2139 9.9884 10.5121 9.7995 10.1926 11.4225 9.4547 11.4619 8.9644 12.6083 12.0276 11.8491 10.3545 11.6757 9.9730 11.6692 8.4360 10.7959 9.5530 9.1348 8.2133 9.9552 9.3619 10.5835 7.6814 9.8116 10.4108 11.0747 8.9451 9.0697 7.6755 11.8451 11.4056 10.4704 10.8396 11.4497 10.6049 9.2203 9.1703 9.5912 11.0472 9.6238 11.0481 8.5500 11.1568 8.1351 10.2309 9.9181 9.0875 12.0132 11.0703 11.3421 9.0859 6.7166 10.3290 8.6440 8.4814 8.1888 6.7693 8.7888 10.4132 10.9575 10.0273 11.1589 10.2656 11.5069 10.8283 10.6386 11.4215 11.6061 13.0207 11.2677 11.1443 11.4897 11.3800 10.6470 11.0611 9.6274 10.6027 10.0972 9.8935 8.8031 10.6574 10.3456 9.3204 10.7040 11.0988 10.4414 11.4296 11.7952 9.4463 10.7339 10.3838 6.7520 10.3326 10.1236 8.4882 8.9324 9.9690 11.1028 9.6834 10.6232 9.5139 9.8925 7.6972 10.4602 11.9742 9.0514 10.3034 9.1859 11.6082 12.0778 9.9100 10.3802 10.2470 10.2981 8.2597 9.3090 9.1055 8.6080 10.2145 8.4912 10.4154 10.1289 10.4772 9.5580 9.5547 11.8105 9.9851 11.0071 12.5906 10.9930 10.0571 11.2286 11.0106 10.8763 9.1213 10.4495 9.7276 9.5743 9.0943 9.4269 9.8780 10.3364 11.1778 11.4280 12.4034 10.4509 8.9858 10.2729 9.9455 10.3513 11.0260 10.0941 10.8615 9.9538 8.6196 10.0294 9.0198 10.5916 10.4572 8.1046 11.6617 9.3167 9.8054 11.1742 10.8962 9.6408 10.9001 10.5413 10.6256 9.2619 12.3197 8.3492 10.4393 9.6345 10.0560 8.3832 8.7877 10.6046 9.6061 9.3789 8.5419 8.5403 7.6492 9.3288 7.5269 11.2788 10.8388 10.4793 8.8586 10.5429 8.9308 ; Title "Fabric Thickness Data"; proc gplot data=fabric; plot Y*location; symbol1 v=none i=join c=green; proc univariate normal plot; var Y; proc arima data=fabric; i var=y; run; proc spectra p s adjmean out=outspec; var Y; proc plot data=outspec; plot P_01*freq;run; proc print data=outspec; where P_01>30; run; proc gplot; plot P_01*freq S_01*freq/vaxis=axis1; symbol1 v=circle i=needle c=red; DATA Waves; set fabric; s5 = sin(2*pi*location/5); c5 = cos(2*pi*location/5); s6 = sin(2*pi*location/6); c6 = cos(2*pi*location/6); s8 = sin(2*pi*location/8); c8 = cos(2*pi*location/8); s11 = sin(2*pi*location/11); c11 = cos(2*pi*location/11); s16 = sin(2*pi*location/16); c16 = cos(2*pi*location/16); s20 = sin(2*pi*location/20); c20 = cos(2*pi*location/20); proc reg; model Y = s5--c20/ss1 ss2; T_5: test s5=0, c5=0; T_6: test s6=0, c6=0; T_8: test s8=0, c8=0; T_11: test s11=0, c11=0; T_16: test s16=0, c16=0; T_20: test s20=0, c20=0; proc reg data=waves; model Y = s5 c5 s8 c8; output out=out1 residual=r; proc spectra p s adjmean whitetest data=out1 out=outspec; var r; weight 1 2 3 4 5 6 7 6 5 4 3 2 1; proc gplot data=outspec; plot (P_01 P_01)*freq/overlay vaxis=axis1; plot (S_01 S_01)*freq/overlay vaxis=axis1; symbol1 v=circle i=needle c=red; symbol2 v=none i=smooth40 c=black w=3; run;