
cOpts.Type='gauss';cOpts.NormalizationType='beltrami';cOpts.Delta=0.1;cOpts.Precision=1e-2;cOpts.MaxEigenVals=30;cOpts.WhichEigenVals='small';
cXNorm=zeros(size(cX));for lk = 1:size(cX,1);cXNorm(lk,:)=cX(lk,:)/norm(cX(lk,:));end;
lW=pdist(cXNorm);
help pdist
lW=pdist(cXNorm);
lW=1-cXNorm*cXNorm';
diary analysis1
cXNorm=zeros(size(cX));for lk = 1:size(cX,1);cXNorm(lk,:)=cX(lk,:)/norm(cX(lk,:));end;
clear cXNorm
lDelta=1.0;lWW=exp(-(lW/lDelta).^2);
figure;imagesc(lWW);
lDD=sum(lWW,2);
cOpts.NormalizationType='graph';cOpts.WhichEigenVals='small';
cOpts
[vEigenVectors,vEigenValues,vLaplacian,vDiffusion] = FastLaplacianEigs( cOpts, lDD, lWW, 30, 1e-8 );
vEigenVals
vEigenValues
lDelta=0.2;lWW=exp(-(lW/lDelta).^2);lDD=sum(lWW,2);[vEigenVectors,vEigenValues,vLaplacian,vDiffusion] = FastLaplacianEigs( cOpts, lDD, lWW, 30, 1e-8 );
vEIgenValues
vEigenValues
lDelta=0.35;lWW=exp(-(lW/lDelta).^2);lDD=sum(lWW,2);[vEigenVectors,vEigenValues,vLaplacian,vDiffusion] = FastLaplacianEigs( cOpts, lDD, lWW, 30, 1e-8 );
vEigenValues
figure;scatter3(vEigenVectors(2,:),vEigenVectors(3,:),vEigenVectors(4,:),10,cLabels(:,2),'filled');
lDelta=0.5;lWW=exp(-(lW/lDelta).^2);lDD=sum(lWW,2);[vEigenVectors,vEigenValues,vLaplacian,vDiffusion] = FastLaplacianEigs( cOpts, lDD, lWW, 30, 1e-8 );
figure;scatter3(vEigenVectors(2,:),vEigenVectors(3,:),vEigenVectors(4,:),10,cLabels(:,2),'filled');
figure;scatter3(vEigenVectors(5,:),vEigenVectors(6,:),vEigenVectors(7,:),10,cLabels(:,2),'filled');
figure;scatter3(vEigenVectors(8,:),vEigenVectors(9,:),vEigenVectors(10,:),10,cLabels(:,2),'filled');
length(find(cLabels(:,2)))
length(find(cLabels(:,2)==8))
lChiPhysics=zeros(size(lW,1));lChiPhysics(find(cLabels(:,2)==8))=1;
lChiAstronomy=zeros(size(lW,1));lChiAstronomy(find(cLabels(:,2)==2))=1;
figure;plot(sum(vDiffusion,2))
figure;plot(sum(vDiffusion,1))
lChiAstronomy=zeros(1,size(lW,1));lChiAstronomy(find(cLabels(:,2)==2))=1;
lChiPhysics=zeros(1,size(lW,1));lChiPhysics(find(cLabels(:,2)==8))=1;
lTime=10;lDiffChiPhysics=lChiPhysics*(vDiffusion^lTime);lDiffChiAstronomy=lChiAstronomy*(vDiffusion^lTime);
lPhyAstroIdxs=find((lDiffChiPhysics>0.5)&(lDiffChiAstronomy>0.5));
lPhyAstroIdxs=find((lDiffChiPhysics>0.2)&(lDiffChiAstronomy>0.2));
lPhyAstroIdxs=find((lDiffChiPhysics>0.1)&(lDiffChiAstronomy>0.1));
lPhyAstroIdxs=find((lDiffChiPhysics>0.2)&(lDiffChiAstronomy>0.2));
figure;imagesc(vDiffusion);
vDiffusion=sparse(vDiffusion);
lTime=10;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;
lTime=30;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;
lPhyAstroIdxs=find((lDiffChiPhysics>0.2)&(lDiffChiAstronomy>0.2));
lPhyAstroIdxs=find((lDiffChiPhysics>0.1)&(lDiffChiAstronomy>0.1));
L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);
figure;plot(sort(L1dist));
lTime=10;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;
L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);
figure;plot(sort(L1dist));
figure;imagesc(lPowerD);
svd(lPowerD)
svds(lPowerD)
figure;imagesc(lPowerD);
svds(vDiffusion)
lTime=3;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;
L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);
figure;plot(sort(L1dist));
L1dist=abs(lChiPhysics-lChiAstronomy);
figure;plot(sort(L1dist));
L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);
figure;plot(sort(L1dist));
[lSorted lSortedIdxs]=sort(L1dist);
help sort
lSortedIdxs(1:5)
intersect(lSortedIdxs(1:5),find(lChiPhysics>0))
intersect(lSortedIdxs(1:10),find(lChiPhysics>0))
setdiff(lSortedIdxs(1:10),find(lChiPhysics>0))
lChiAstronomy(78)
lChiPhysics(78)
cLabels(78)
cLabels(78,:)
lTime=1;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;
lTime=1;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);
lTime=1;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);figure;plot(sort(L1dist));
lTime=1;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);figure;plot(sort(L1dist));[lSorted lSortedIdxs]=sort(L1dist);
lSortedIdxs(1:10)
lSortedIdxs(1:10)'
lChiPhysics(lSortedIdxs(1:10))
lChiAstronomy(lSortedIdxs(1:10))
lDiffChiAstronomy(lSortedIdxs(1:10))
lDiffChiPhysics(lSortedIdxs(1:10))
cLabels(lSortedIdxs,2)
cLabels(lSortedIdxs(1:10),2)'
lTime=1/2;lPowerD=vDiffusion^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);figure;plot(sort(L1dist));[lSorted lSortedIdxs]=sort(L1dist);
lTime=1/2;lPowerD=full(vDiffusion)^lTime;lDiffChiPhysics=lChiPhysics*lPowerD;lDiffChiAstronomy=lChiAstronomy*lPowerD;L1dist=abs(lDiffChiPhysics-lDiffChiAstronomy);figure;plot(sort(L1dist));[lSorted lSortedIdxs]=sort(L1dist);
close all;
help pdist
cXNorm=zeros(size(cX));for lk = 1:size(cX,1);cXNorm(lk,:)=signum(cX(lk,:));end;
cXNorm=zeros(size(cX));for lk = 1:size(cX,1);cXNorm(lk,:)=sign(cX(lk,:));end;
figure;plot(cXNorm(1,:))
figure;plot(cXNorm(2,:))
figure;plot(cXNorm(3,:))
lW=pdist(cXNorm,'seuclidean');