mobeets / nullSpaceControl

1 stars 0 forks source link

check out this shape #107

Closed mobeets closed 8 years ago

mobeets commented 8 years ago

i mean, check it out: it's got clear BOUNDARIES


D = io.quickLoadByDate('20120601');
Y = D.blocks(2).latents;
NB = D.blocks(2).fDecoder.NulM2;
[u,s,v] = svd(D.blocks(2).latents*NB);
Y2 = Y*NB*v;

Y2p = Y2(:,1:3);
plot3(Y2p(:,1), Y2p(:,2), Y2p(:,3), '.')

screen shot 2016-03-30 at 6 15 48 pm

mobeets commented 8 years ago

and here by thetaGrp

screen shot 2016-03-30 at 6 19 39 pm

mobeets commented 8 years ago

screen shot 2016-03-30 at 6 21 15 pm

mobeets commented 8 years ago

note that in all these cases, the colors are very clearly clustered. but the overall distribution still is in a sort of balloon.

and for 20120525, it has the same nub, but then from another view, two lobes: screen shot 2016-03-30 at 6 22 06 pm

and then the two 2013 days just have way less data:

screen shot 2016-03-30 at 6 26 35 pm

mobeets commented 8 years ago

"You could treat your data as a sample from a three-dimensional probability density, and estimate that density on a grid, e.g. via a 3d histogram, or better a 3d kernel density estimator. Then apply a threshold and extract the surface using isosurface."

mobeets commented 8 years ago

no isosurface, but i've evaluated by p-hat density estimation and found points where the p-hat is in the 75-76th percentile

screen shot 2016-04-01 at 10 19 54 am

mobeets commented 8 years ago

the black dots are data, and the red dots are the boundary defined by where p-hat is in the 90-91st percentile

screen shot 2016-04-01 at 11 38 28 am

mobeets commented 8 years ago

here's observed on the left, and cloud-hab on the right. note how cloud-hab doesn't fill the same space at all--it's less spread out. but at least it's not outside of it!

screen shot 2016-04-01 at 11 52 01 am

screen shot 2016-04-01 at 11 52 52 am

mobeets commented 8 years ago

and it turns out that left clump are the kinematics the cloud does poorly at:

screen shot 2016-04-01 at 11 56 47 am

screen shot 2016-04-01 at 11 57 12 am

(green are the 180+ directions)

mobeets commented 8 years ago

so this realization has led to two things: