/* This SAS program is to calculate the spectral density distribution
of the entropy, Factor Score Means/Variances profiles of bHLH domains 
through Finite Fourier transformation. The hamming window is applied
with the WEIGHT statement. The white noise tests are conducted with 
the whitetest option in PROC SPECTRA.*/
data bhlh;
input entropy f1v f2v f3v f4v f5v f1m f2m f3m f4m f5m @;
x=_N_;
cards;
2.5218	0.666645001	0.454495159	2.357253304	0.282970622	2.537528827	1.203664948	-0.064809278	0.692118557	0.207371134	1.713845361
1.8778	0.375205963	0.30686128	1.284978215	0.147260984	2.219134829	1.387570681	-0.124465969	1.01982199	0.304837696	2.072706806
3.775	1.228520073	0.67758172	3.958101509	0.565226346	2.613558913	0.192818182	-0.321767045	-0.100602273	0.404744318	0.122482955
3.5703	0.689390513	1.092351144	4.213843457	0.671576429	1.723467096	0.037505682	-0.145301136	-0.398102273	0.524454545	-0.040965909
2.8449	0.620110602	0.524144688	2.247920076	1.463985133	1.353562255	0.272655738	-0.413704918	-0.698961749	0.001765027	0.237775956
2.4388	0.536695991	0.516539091	3.212132422	0.299754885	1.743339375	0.557391753	0.679082474	0.467536082	0.127298969	0.31206701
3.4662	0.909382134	0.913589619	3.817955713	0.779669802	1.519328908	-0.188066327	-0.709040816	-0.144622449	0.591331633	-0.331994898
2.9641	1.41033437	0.591070628	3.994630024	0.345786311	3.642335336	0.33162234	-0.238446809	0.371388298	0.519888298	1.032654255
0.84337	0.303683387	0.346583475	0.776180256	0.126808091	0.285547599	1.204056122	-1.301306122	1.197102041	0.198204082	-0.751872449
1.1962	0.25541521	0.099287508	0.59557315	0.131875323	0.903346562	1.481020513	-0.170066667	1.207753846	0.292620513	2.472425641
3.2076	0.989712106	0.374986229	3.452149468	0.536739482	3.268650448	0.840459184	-0.043321429	0.301617347	0.100479592	1.078408163
0.84117	0.251381795	0.069610202	0.519655982	0.103768672	0.710000086	1.395069149	-0.057117021	1.423345745	0.398601064	2.652468085
2.2432	1.533413942	0.29519428	2.037688931	0.492428633	3.151544176	0.986837928	-0.380918367	0.778413265	0.464112245	1.349443878
3.4898	0.615566377	0.754235776	5.13883458	0.52765455	3.444173319	0.826112821	-0.218317949	-0.500902564	-0.022861538	-0.261138462
3.2771	0.704073382	0.659403948	5.149791379	0.359591855	4.468967028	0.949830769	0.044287179	-0.628671795	0.047338462	0.092789744
2.2486	0.465602095	0.162482979	2.697532517	0.498366764	0.940846746	-0.940545918	-0.77402551	0.580607143	0.612015306	-0.136780612
1.2538	0.180282389	0.207035797	1.364396849	0.100455447	0.761755118	0.958076923	0.658389744	1.067112821	-0.052933333	0.996810256
3.5331	0.84428974	1.049102657	5.920484947	0.476559389	3.072387429	0.51969898	-0.154734694	-0.574821429	0.29202551	-0.621193878
3.1301	0.782932615	1.307925623	4.149541919	1.098371891	2.131566039	-0.227204082	-0.264346939	-0.408933673	0.566270408	0.052586735
1.8109	0.143344076	0.177259191	2.026131748	0.473883818	0.446358968	-0.949963918	-0.576881443	1.368670103	0.029505155	0.307891753
3.5862	1.051954895	0.875399467	5.486517008	0.580537588	2.638898739	0.436545918	-0.392341837	-0.513867347	0.248760204	-0.647658163
3.2143	0.976100851	0.831449743	3.682739259	0.432595444	2.23355878	0.635102041	-0.606576531	0.281887755	0.400112245	-0.10222449
0.20632	0.030716803	0.00970303	0.125245266	0.010619532	0.017211365	-1.012867347	-0.976295918	-1.45652551	1.255295918	-0.908158163
2.2379	0.944091475	0.670352238	2.822631935	0.382465225	2.599709081	1.030780612	-0.122790816	0.39319898	0.40897449	1.55575
3.2845	0.683746353	0.692487973	6.150051822	0.366903386	4.963690488	0.890964286	-0.143647959	-0.761591837	0.140244898	-0.317765306
3.1021	0.733141789	0.451407413	3.246076276	1.080487009	1.434331465	-0.586994898	-0.628163265	-0.48044898	0.367908163	-0.28580102
2.2612	0.183652462	0.207386201	2.408970358	0.40505581	1.055903893	-1.077285714	-0.472862245	0.288030612	0.820071429	-0.245877551
1.6988	0.281733265	1.206122792	1.370749937	0.181396447	0.375737391	0.020517949	1.472466667	-1.598205128	0.457061538	-1.324097436
3.8851	0.702621036	0.946899758	4.824192253	0.617694525	2.188025882	0.080186047	-0.066156977	-0.362802326	0.422877907	-0.332889535
2.3729	0.315255011	0.322238007	1.199706239	0.513535999	0.486815661	-0.129564103	-0.210538462	-0.486064103	0.132230769	-0.138230769
3.6228	0.7159924	0.735122316	3.020740529	0.742498094	1.349121675	-0.37198895	-0.235955801	-0.014933702	0.483254144	-0.075718232
1.9197	0.384952588	0.462236088	1.354278143	0.21524877	1.679507422	1.491394872	-0.2392	0.412574359	-0.052225641	1.490666667
2.5369	0.419520919	0.749423889	1.739666989	0.722679869	0.689308803	-0.674302564	-0.63294359	-1.163128205	0.866210256	-0.752394872
2.5552	0.381342187	1.086418605	6.962512293	0.231179686	2.407182503	0.293132653	0.794377551	-2.252403061	0.408107143	-1.650811224
0.52621	0.138521139	0.029471797	0.315072276	0.061600099	0.27168206	1.746188776	-0.519209184	0.540969388	-0.220408163	1.611040816
2.9725	0.323156496	0.695201268	2.463141625	0.564367617	1.25968236	-0.855806122	-0.437362245	-0.15252551	0.919362245	-0.242755102
3.1283	0.967997342	1.297002714	6.178267966	0.386013566	2.385640405	0.341739796	-0.304061224	-0.553862245	0.476331633	-0.909984694
1.5697	0.310431615	0.188804974	1.13588793	0.132348969	0.585676904	-0.93452551	-0.328989796	1.733382653	0.60552551	0.664755102
0.42438	0.026792235	0.025887459	0.584209527	0.184399089	0.180241171	-1.005127551	-0.973709184	-1.349760204	1.176617347	-0.843346939
2.7717	0.451129661	0.400337593	3.918902196	0.326359318	4.06352021	1.202204082	-0.212045918	0.235433673	0.124494898	1.0275
3.3913	1.328166274	0.693244724	4.135306655	0.767978409	2.239584683	0.204415385	-0.527061538	-0.334123077	0.145871795	0.042861538
1.5868	0.142650317	0.838213501	2.885892766	0.134951809	0.992018616	-0.487663265	-0.714352041	-0.520183673	1.331459184	-0.163341837
2.7605	0.730403204	0.282424404	2.902459726	0.503857037	1.49748541	-0.774178571	-0.33652551	0.526545918	0.662433673	0.021214286
3.3021	0.696776322	1.102712949	6.72132881	0.389929076	3.819411093	0.67897449	-0.045214286	-0.908188776	0.258663265	-0.679673469
1.6571	0.275078404	0.32683462	3.982333296	0.377465451	1.093498853	0.101678571	0.525515306	1.994459184	-0.694520408	0.966908163
1.4756	0.040457561	0.147072069	2.698222373	0.256828044	0.77533231	-1.15244898	-0.628586735	1.140178571	0.580867347	0.320943878
3.5595	1.071703398	0.465709903	3.657944817	0.750636185	3.334088994	0.719676923	-0.349594872	-0.061353846	-0.001825641	0.397312821
3.9359	0.870935373	1.049168178	6.582358448	0.714130311	3.069616515	0.197561856	-0.031747423	-0.658494845	0.277701031	-0.33664433
1.6961	0.563498234	0.18905189	1.899810435	0.582669957	1.024742146	-0.743425641	-0.858676923	-0.899907692	0.887579487	-0.47745641
run;

/* Periodogram Analysis and White noise test*/
proc spectra data=bhlh whitetest out=specentropy p s adjmean ;
var entropy f1v f2v f3v f4v f5v f1m f2m f3m f4m f5m;
weights 0.08 0.54 1 0.54 0.08;  
run;

proc gplot data=specentropy;
plot s_01*period s_02*period s_03*period s_04*period s_05*period s_06*period s_07*period s_08*period s_09*period s_10*period s_11*period;
symbol1 i=join;
run;

PROC EXPORT DATA= WORK.SPECENTROPY
            OUTFILE= "D:\bhlh\Protein structure function bioinformatics\spectral_result.xls"
            DBMS=EXCEL2000 REPLACE;
RUN;

/* Decompose Entropy (dependent variable) into five 
   Factor Score Variance components (independent variables) */
proc reg data=bhlh;
model entropy=f1v f2v f3v f4v f5v;
run;