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