%---- compute and sort eigenvalues/eigenvectors ----- hold on A = scAGCT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scATGC; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scAGCT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scCGCG; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scAGCT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scCGGC; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scCGTA; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scTAAT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scTATA; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) A = scTCGA; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(vpA) hold off line(158*ones(1,160),[1:160]) line(159*ones(1,5),[1:5]) line((158+9)*ones(1,160),[1:160]) %---- same, log scale ----- hold on A = scAGCT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scATGC; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scAGCT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scCGCG; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scAGCT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scCGGC; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scCGTA; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scTAAT; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scTATA; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) A = scTCGA; [vA, vpA] = eig(A); [vpA, s] = sort(diag(vpA)); %sort eigenvalues plot(log(vpA)) hold off line(158*ones(1,5),[1:5]) line(159*ones(1,5),[1:5]) line((158+9)*ones(1,5),[1:5])