Update Issues with covaraince matrix authored by Soumen Roy's avatar Soumen Roy
......@@ -38,13 +38,13 @@ The above plot shows the array of the eigenvalues. This plot indicates that the
## GW190814
First, we demonstrate an example from the GW190814 event to find an appropriate cut-off on the low-rank approximation of the covariance matrix, and its effect on $`\beta`$. We choose top 5 [PE samples](https://ldas-jobs.ligo.caltech.edu/~charlie.hoy/PE/O3/S190814bv/C01/SEOBNRv4HM/samples/SEOBNRv4HM_pesummary.dat) ( largest likelihood values ) and plot $`\beta`$ as a function of the rank of the covariance matrix.
| m1 | m2 | s1z | s2z | logLikelihood |
| ------- | ------- | -------- | ------- | -------------- |
| 24.2080 | 2.7402 | -0.0031 | -0.1011 | -45762.6371 |
| 24.2444 | 2.7392 | -0.0162 | 0.1414 | -45762.8520 |
| 23.4810 | 2.7965 | -0.0026 | -0.5558 | -45763.1355 |
| 24.0278 | 2.7539 | -0.0236 | 0.0809 | -45763.1803 |
| 23.5563 | 2.7887 | 0.0045 | -0.6138 | -45763.2090 |
| m1 | m2 | s1z | s2z | logLikelihood | beta |
| ------- | ------- | -------- | ------- | -------------- | ------ |
| 24.2080 | 2.7402 | -0.0031 | -0.1011 | -45762.6371 | 8.7 |
| 24.2444 | 2.7392 | -0.0162 | 0.1414 | -45762.8520 | 8.5 |
| 23.4810 | 2.7965 | -0.0026 | -0.5558 | -45763.1355 | 8.1 |
| 24.0278 | 2.7539 | -0.0236 | 0.0809 | -45763.1803 | 8.2 |
| 23.5563 | 2.7887 | 0.0045 | -0.6138 | -45763.2090 | 8.1 |
The spectrum of the eigenvalues is also plotted.
......@@ -57,12 +57,12 @@ We can see that the median $`\beta`$ (black trace) saturates after a certain num
A similar plot for the GW190412 event (top 5 [PE samples](https://git.ligo.org/pe/O3/S190412m/-/blob/master/Preferred/dat/M4-PhenomHM_pesummary.dat)) is given below, and we see that in this case, the plateau occurs once we use the top-350 eigenvectors to approximate the covariance matrix. In this case, the Frobenius norm of the residuals is 3.5e-10.
| m1 | m2 | s1z | s2z | logLikelihood |
| m1 | m2 | s1z | s2z | logLikelihood |
| ------- | ------- | -------- | ------- | -------------- |
| 30.7206 | 10.4021 | 0.0465 | 0.7059 | -48629.6293 |
| 33.7896 | 9.7169 | 0.3550 | -0.2997 | -48629.9648 |
| 33.0223 | 9.8901 | 0.2510 | 0.1500 | -48630.0465 |
| 31.1000 | 10.2668 | 0.0572 | 0.6907 | -48630.0520 |
| 33.7223 | 9.7035 | 0.3080 | -0.0899 | -48630.0727 |
| 30.7206 | 10.4021 | 0.0465 | 0.7059 | -48629.6293 |
| 33.7896 | 9.7169 | 0.3550 | -0.2997 | -48629.9648 |
| 33.0223 | 9.8901 | 0.2510 | 0.1500 | -48630.0465 |
| 31.1000 | 10.2668 | 0.0572 | 0.6907 | -48630.0520 |
| 33.7223 | 9.7035 | 0.3080 | -0.0899 | -48630.0727 |
<img src="uploads/b950108077a19ca098f0b79b9949e030/beta_lambda_all.png" width="440" >
\ No newline at end of file