diff --git a/homework/AthishThiruvengadam/06/HW06+7.ipynb b/homework/AthishThiruvengadam/06/HW06+7.ipynb new file mode 100644 index 0000000..0360ead --- /dev/null +++ b/homework/AthishThiruvengadam/06/HW06+7.ipynb @@ -0,0 +1,1512 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# ASTR 596: FDS Homework 6+7: Gaussian Processes (200 pts)\n", + "\n", + "### This is a double HW set so you get extra time - until reading day (May 4th, 2023) at noon to do it. \n", + "### After that, it's finals time. \n", + "\n", + "\n", + "# P1. Gaussian Processes\n", + "\n", + "### Last HW, you worked on finding periodic planet signals in the light curve of Kepler-90, a star that is photometrically stable. The periodogram worked nicely because \n", + "\n", + "### a) we cleaned the light curve to squelch red noise\n", + "### b) the signals really were periodic and we could implictly make a strong assumption about the covariance between points.\n", + "\n", + "### Life gets harder when the star itself has quasi-periodic variations because it has a magnetic field and is rotating (ruh oh...) " + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "import matplotlib.pyplot as plt\n", + "import numpy as np\n", + "\n", + "from astropy.table import Table\n", + "import scipy.stats as st\n", + "import sklearn\n", + "import sklearn.ensemble" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
| time | flux | error | quarter |
|---|---|---|---|
| float64 | float64 | float64 | int64 |
| 539.4710179205795 | 5221.16455078125 | 4.891251087188721 | 6 |
| 539.4914521464307 | 5216.07958984375 | 4.892819404602051 | 6 |
| 539.5118864718097 | 5225.43359375 | 4.891888618469238 | 6 |
| 539.5323205971945 | 5233.111328125 | 4.892074108123779 | 6 |
| 539.5527548221144 | 5221.17333984375 | 4.891964435577393 | 6 |
| 539.573189147035 | 5220.09521484375 | 4.891523361206055 | 6 |
| ... | ... | ... | ... |
| 905.8044143491425 | 5862.43798828125 | 4.119598388671875 | 9 |
| 905.824848745775 | 5863.95458984375 | 4.118742942810059 | 9 |
| 905.845282941933 | 5870.349609375 | 4.119640827178955 | 9 |
| 905.865717337867 | 5867.27099609375 | 4.119102478027344 | 9 |
| 905.8861516335674 | 5868.67822265625 | 4.119019031524658 | 9 |
| 905.9065858290269 | 5864.173828125 | 4.119871616363525 | 9 |
| 905.9270202242624 | 5861.61572265625 | 4.119407653808594 | 9 |
| \n", + " | imageID | \n", + "specObjID | \n", + "bestObjID | \n", + "class | \n", + "subClass | \n", + "z | \n", + "zErr | \n", + "objID | \n", + "type | \n", + "flags | \n", + "... | \n", + "p_dk | \n", + "P_disk | \n", + "P_edge_on | \n", + "P_bar_GZ2 | \n", + "P_bulge | \n", + "P_cigar | \n", + "P_merg | \n", + "P_bar_Nair10 | \n", + "TT_pred | \n", + "P_S0 | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", + "604197 | \n", + "3109758158738843648 | \n", + "1237668333638516913 | \n", + "GALAXY | \n", + "NaN | \n", + "0.115187 | \n", + "0.000022 | \n", + "1237668333638516913 | \n", + "3 | \n", + "35253360136720 | \n", + "... | \n", + "0.094 | \n", + "2.002438e-06 | \n", + "0.000882 | \n", + "0.006050 | \n", + "0.881993 | \n", + "0.002938 | \n", + "0.114075 | \n", + "0.348058 | \n", + "-1.870139 | \n", + "0.917128 | \n", + "
| 1 | \n", + "410863 | \n", + "1559448938690406400 | \n", + "1237662664284438746 | \n", + "GALAXY | \n", + "NaN | \n", + "0.130501 | \n", + "0.000023 | \n", + "1237662664284438746 | \n", + "3 | \n", + "103347785728 | \n", + "... | \n", + "NaN | \n", + "1.579445e-07 | \n", + "0.434940 | \n", + "0.001006 | \n", + "0.848765 | \n", + "0.352186 | \n", + "0.049307 | \n", + "0.214841 | \n", + "-2.130526 | \n", + "0.864439 | \n", + "
| 2 | \n", + "352194 | \n", + "1888165532671174656 | \n", + "1237662499464675448 | \n", + "GALAXY | \n", + "STARBURST | \n", + "0.135090 | \n", + "0.000006 | \n", + "1237662499464675448 | \n", + "3 | \n", + "68987912704 | \n", + "... | \n", + "0.080 | \n", + "3.629266e-05 | \n", + "0.393671 | \n", + "0.015625 | \n", + "0.003221 | \n", + "0.961169 | \n", + "0.009676 | \n", + "0.255024 | \n", + "3.664025 | \n", + "0.996161 | \n", + "
| 3 | \n", + "391105 | \n", + "1109072471920764928 | \n", + "1237663543139959665 | \n", + "GALAXY | \n", + "NaN | \n", + "0.224838 | \n", + "0.000050 | \n", + "1237663543139959665 | \n", + "3 | \n", + "35253360132368 | \n", + "... | \n", + "0.158 | \n", + "1.203650e-03 | \n", + "0.171578 | \n", + "0.001327 | \n", + "0.235339 | \n", + "0.477329 | \n", + "0.727791 | \n", + "0.370840 | \n", + "-0.272969 | \n", + "0.815209 | \n", + "
| 4 | \n", + "33174 | \n", + "507800681206278144 | \n", + "1237651249878597774 | \n", + "GALAXY | \n", + "STARFORMING | \n", + "0.180001 | \n", + "0.000012 | \n", + "1237651249878597774 | \n", + "3 | \n", + "2450547277824 | \n", + "... | \n", + "0.109 | \n", + "4.230215e-04 | \n", + "0.000493 | \n", + "0.008908 | \n", + "0.051623 | \n", + "0.074105 | \n", + "0.031973 | \n", + "0.249713 | \n", + "0.258085 | \n", + "0.956782 | \n", + "
| ... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "
| 413215 | \n", + "163794 | \n", + "933472216362805248 | \n", + "1237657606427115859 | \n", + "GALAXY | \n", + "NaN | \n", + "0.068819 | \n", + "0.000029 | \n", + "1237657606427115859 | \n", + "3 | \n", + "68988043776 | \n", + "... | \n", + "0.143 | \n", + "8.809507e-04 | \n", + "0.004825 | \n", + "0.001832 | \n", + "0.341333 | \n", + "0.000163 | \n", + "0.062108 | \n", + "0.240102 | \n", + "-1.975903 | \n", + "0.929764 | \n", + "
| 413216 | \n", + "2942 | \n", + "345812102876784640 | \n", + "1237648704593330407 | \n", + "GALAXY | \n", + "NaN | \n", + "0.102469 | \n", + "0.000021 | \n", + "1237648704593330407 | \n", + "3 | \n", + "35253360136208 | \n", + "... | \n", + "0.069 | \n", + "2.408648e-07 | \n", + "0.102996 | \n", + "0.001373 | \n", + "0.967661 | \n", + "0.000316 | \n", + "0.573898 | \n", + "0.102444 | \n", + "-2.340441 | \n", + "0.818810 | \n", + "
| 413217 | \n", + "10230 | \n", + "301908603076372480 | \n", + "1237648722828198080 | \n", + "GALAXY | \n", + "NaN | \n", + "0.097740 | \n", + "0.000022 | \n", + "1237648722828198080 | \n", + "3 | \n", + "35255507619856 | \n", + "... | \n", + "0.091 | \n", + "5.508738e-02 | \n", + "0.011437 | \n", + "0.074087 | \n", + "0.732646 | \n", + "0.000677 | \n", + "0.983481 | \n", + "0.325126 | \n", + "-1.425996 | \n", + "0.677585 | \n", + "
| 413218 | \n", + "301804 | \n", + "1937748846054500352 | \n", + "1237661949200302413 | \n", + "GALAXY | \n", + "NaN | \n", + "0.033959 | \n", + "0.000014 | \n", + "1237661949200302413 | \n", + "3 | \n", + "35253360136720 | \n", + "... | \n", + "0.039 | \n", + "6.966625e-05 | \n", + "0.079195 | \n", + "0.126760 | \n", + "0.561732 | \n", + "0.523639 | \n", + "0.251242 | \n", + "0.263267 | \n", + "-1.844626 | \n", + "0.847763 | \n", + "
| 413219 | \n", + "160515 | \n", + "1006650218101368832 | \n", + "1237657595688059046 | \n", + "GALAXY | \n", + "NaN | \n", + "0.040637 | \n", + "0.000010 | \n", + "1237657595688059046 | \n", + "3 | \n", + "35253360136208 | \n", + "... | \n", + "0.000 | \n", + "2.975413e-02 | \n", + "0.021821 | \n", + "0.235166 | \n", + "0.962262 | \n", + "0.552131 | \n", + "0.389143 | \n", + "0.454696 | \n", + "-1.686107 | \n", + "0.839427 | \n", + "
413220 rows × 121 columns
\n", + "| ug | gr | gi | gz | petroRad_i | petroR50_i | deVAB_i | fracDev_i |
|---|---|---|---|---|---|---|---|
| float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 |
| 1.999718 | 0.9592915 | 1.370582 | 1.692535 | 5.454633 | 2.344662 | 0.8708095 | 1.0 |
| 1.98884 | 0.9332962 | 1.361919 | 1.668918 | 5.177937 | 2.424093 | 0.9946066 | 0.8865207 |
| 1.020407 | 1.601521 | 2.396175 | 2.827263 | 5.846272 | 1.158836 | 0.6898293 | 1.0 |
| 1.83742 | 1.211882 | 1.688545 | 1.989878 | 7.145427 | 2.795447 | 0.7493948 | 1.0 |
| 4.55904 | 1.536175 | 2.29286 | 2.712597 | 4.741546 | 2.09032 | 0.5173162 | 0.8244159 |
| 4.532835 | 1.433081 | 2.251518 | 2.669952 | 8.930706 | 1.656738 | 0.6116931 | 1.0 |
| ... | ... | ... | ... | ... | ... | ... | ... |
| 4.392231 | 1.8193 | 2.470142 | 2.87352 | 3.710293 | 1.530089 | 0.6766507 | 1.0 |
| 2.502085 | 1.588123 | 2.214443 | 2.619152 | 6.212802 | 2.041965 | 0.7738834 | 1.0 |
| 1.215888 | 0.4952621 | 0.69736 | 0.8714724 | 15.58729 | 7.597391 | 0.4803397 | 0.2104624 |
| 1.750438 | 0.8658667 | 1.292953 | 1.691803 | 13.22304 | 5.658463 | 0.2999924 | 0.55682 |
| 0.8426914 | 0.2660866 | 0.3770752 | 0.4642525 | 6.997887 | 3.126914 | 0.8816415 | 0.07251444 |
| 0.8361015 | 0.4111938 | 0.701189 | 0.788826 | 3.364222 | 1.517709 | 0.9998366 | 0.6170198 |
| 0.9144573 | 1.080809 | 1.867247 | 2.458876 | 2.523635 | 1.270804 | 0.5280833 | 0.4554564 |
| ug | gr | gi | gz | petroRad_i | petroR50_i | deVAB_i | fracDev_i |
|---|---|---|---|---|---|---|---|
| float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 |
| 1.999718 | 0.9592915 | 1.370582 | 1.692535 | 5.454633 | 2.344662 | 0.8708095 | 1.0 |
| 1.98884 | 0.9332962 | 1.361919 | 1.668918 | 5.177937 | 2.424093 | 0.9946066 | 0.8865207 |
| 1.020407 | 1.601521 | 2.396175 | 2.827263 | 5.846272 | 1.158836 | 0.6898293 | 1.0 |
| 1.83742 | 1.211882 | 1.688545 | 1.989878 | 7.145427 | 2.795447 | 0.7493948 | 1.0 |
| 4.55904 | 1.536175 | 2.29286 | 2.712597 | 4.741546 | 2.09032 | 0.5173162 | 0.8244159 |
| 4.532835 | 1.433081 | 2.251518 | 2.669952 | 8.930706 | 1.656738 | 0.6116931 | 1.0 |
| ... | ... | ... | ... | ... | ... | ... | ... |
| 0.3703461 | 0.7442894 | 1.246332 | 1.428238 | 2.177634 | 1.006152 | 0.7540503 | 1.0 |
| 4.392231 | 1.8193 | 2.470142 | 2.87352 | 3.710293 | 1.530089 | 0.6766507 | 1.0 |
| 2.502085 | 1.588123 | 2.214443 | 2.619152 | 6.212802 | 2.041965 | 0.7738834 | 1.0 |
| 1.215888 | 0.4952621 | 0.69736 | 0.8714724 | 15.58729 | 7.597391 | 0.4803397 | 0.2104624 |
| 1.750438 | 0.8658667 | 1.292953 | 1.691803 | 13.22304 | 5.658463 | 0.2999924 | 0.55682 |
| 0.8361015 | 0.4111938 | 0.701189 | 0.788826 | 3.364222 | 1.517709 | 0.9998366 | 0.6170198 |
| 0.9144573 | 1.080809 | 1.867247 | 2.458876 | 2.523635 | 1.270804 | 0.5280833 | 0.4554564 |
| attributes | mean | label 0 mean | label 1 mean | label 2 mean | std | label 0 std | label 1 std | label 2 std |
|---|---|---|---|---|---|---|---|---|
| str10 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 |
| ug | 1.9736512245600564 | 1.6904043577463796 | 2.4520162178109417 | 1.5776578690333711 | 0.8961525007326261 | 0.45550165450762264 | 1.1695704572872982 | 0.4518459673774183 |
| gr | 1.1834039754149408 | 0.9048157888218405 | 1.661776047727539 | 0.7730855675839932 | 0.4481004969196309 | 0.2255762352281746 | 0.19378145612386463 | 0.24118402346210072 |
| gi | 1.7106376014767795 | 1.3288281979595375 | 2.3584627985444335 | 1.1689089788766394 | 0.60344629124241 | 0.28702561715780883 | 0.27149001077776136 | 0.32273227402652105 |
| gz | 2.0498659006526823 | 1.6371854860110822 | 2.758415652947522 | 1.4422511732347703 | 0.6801238040911088 | 0.3607720904778782 | 0.30906975710387347 | 0.42775409949934656 |
| petroRad_i | 4.701980754413328 | 4.0778936018991265 | 4.112893787031611 | 8.177807855146403 | 2.2714444460805217 | 1.2478232890698362 | 1.622579164584602 | 2.97375887958622 |
| petroR50_i | 2.0079861804822747 | 1.786073761482041 | 1.7069565242486928 | 3.4862432415753695 | 0.8816581020395255 | 0.4853879460342915 | 0.4895153759103631 | 1.1395754092047046 |
| deVAB_i | 0.6710848352143614 | 0.6514132790595476 | 0.7173066029052122 | 0.609171389219638 | 0.20204103779651814 | 0.2108634452326184 | 0.1743299984311314 | 0.2217801464439574 |
| fracDev_i | 0.8456755631637496 | 0.8602413748941022 | 0.9175590730199601 | 0.6107622598601209 | 0.23205721294445394 | 0.20887139756394518 | 0.15680123405999022 | 0.30924868938226646 |
| attributes | mean | label 0 mean | label 1 mean | label 2 mean | std | label 0 std | label 1 std | label 2 std |
|---|---|---|---|---|---|---|---|---|
| str10 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 |
| ug | 1.9736516957311474 | 1.6904051193682015 | 2.4520173817157955 | 1.5776611672921523 | 0.8933883784815817 | 0.451842000758611 | 1.1548427995102553 | 0.4432670222084905 |
| gr | 1.1834043192038515 | 0.9048161860536302 | 1.6617772108925426 | 0.7730863869980574 | 0.44311439091916527 | 0.22356338060653957 | 0.19132014266353436 | 0.23733112079231797 |
| gi | 1.7106379324178813 | 1.3288284834707818 | 2.35846428141515 | 1.1689100435314685 | 0.596731290031454 | 0.2846615010483446 | 0.2681093703204747 | 0.3180687439521963 |
| gz | 2.0498660911028224 | 1.6371856494257526 | 2.7584171610393424 | 1.4422521375252524 | 0.672556364172342 | 0.3580222185065494 | 0.30482849496180736 | 0.4221408217876327 |
| petroRad_i | 4.7019837488433955 | 4.077895955097395 | 4.112905038608338 | 8.177835447163947 | 2.257857797273874 | 1.2290869507502458 | 1.6067297827940823 | 2.9290545213767945 |
| petroR50_i | 2.007986912659033 | 1.786074649127245 | 1.7069578058719908 | 3.486249784032634 | 0.876650115916793 | 0.47804444461645734 | 0.4837541173445275 | 1.1248448901039718 |
| deVAB_i | 0.6711212001839 | 0.6514468011902352 | 0.717359646936289 | 0.6091823886713287 | 0.19977073347981858 | 0.2076591524652485 | 0.17174315647199279 | 0.21658687648835012 |
| fracDev_i | 0.8464330839232015 | 0.8609780381608905 | 0.9184447151644157 | 0.6111822524475524 | 0.2294430300042191 | 0.20572872049824822 | 0.15396450652079668 | 0.3020113921230028 |
| HJD | MAG_0 | MAG_1 | MAG_2 | MAG_3 | MAG_4 | MER_0 | MER_1 | MER_2 | MER_3 | MER_4 | GRADE | FRAME |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | str1 | int64 |
| 1877.72747 | 12.174 | 12.255 | 12.266 | 12.272 | 12.28 | 0.198 | 0.153 | 0.136 | 0.134 | 0.131 | B | 1525 |
| 1879.71766 | 11.964 | 11.948 | 11.963 | 11.956 | 11.958 | 0.174 | 0.14 | 0.122 | 0.124 | 0.129 | B | 1835 |
| 1880.71166 | 12.272 | 12.277 | 12.265 | 12.259 | 12.257 | 0.194 | 0.153 | 0.135 | 0.128 | 0.139 | B | 1988 |
| 1881.70621 | 11.433 | 11.471 | 11.464 | 11.445 | 11.452 | 0.209 | 0.17 | 0.142 | 0.141 | 0.143 | B | 2144 |
| 1885.70256 | 12.032 | 12.005 | 12.025 | 12.031 | 12.038 | 0.109 | 0.085 | 0.086 | 0.085 | 0.085 | B | 2451 |
| 1889.69087 | 12.274 | 12.279 | 12.292 | 12.293 | 12.349 | 0.195 | 0.154 | 0.125 | 0.118 | 0.118 | B | 3087 |
| ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
| 5086.87177 | 11.944 | 11.912 | 11.894 | 11.878 | 11.871 | 0.05 | 0.047 | 0.042 | 0.046 | 0.051 | A | 326589 |
| 5098.88785 | 12.248 | 12.212 | 12.179 | 12.133 | 12.142 | 0.048 | 0.044 | 0.041 | 0.046 | 0.053 | A | 328055 |
| 5123.83111 | 12.032 | 12.054 | 12.072 | 12.06 | 12.071 | 0.049 | 0.047 | 0.047 | 0.051 | 0.055 | B | 330804 |
| 5131.8245 | 11.83 | 11.858 | 11.863 | 11.879 | 11.947 | 0.056 | 0.048 | 0.043 | 0.046 | 0.05 | B | 331799 |
| 5154.83688 | 12.223 | 12.232 | 12.243 | 12.251 | 12.302 | 0.048 | 0.047 | 0.044 | 0.047 | 0.05 | B | 333419 |
| 5162.78817 | 11.888 | 11.871 | 11.86 | 11.871 | 11.879 | 0.047 | 0.046 | 0.045 | 0.048 | 0.052 | B | 334245 |
| 5165.75062 | 11.92 | 11.923 | 11.906 | 11.904 | 11.864 | 0.048 | 0.045 | 0.042 | 0.048 | 0.055 | A | 334609 |
| \n", + " | ID | \n", + "ID_Star | \n", + "Star | \n", + "AltName | \n", + "Karmn | \n", + "RA_J2000 | \n", + "DE_J2000 | \n", + "SpT | \n", + "SpTnum | \n", + "d_pc | \n", + "... | \n", + "eurho_gcm-3 | \n", + "edrho_gcm-3 | \n", + "Teq_K | \n", + "euTeq_K | \n", + "edTeq_K | \n", + "TSM | \n", + "ESM | \n", + "Method | \n", + "DiscoveryRef | \n", + "ParameterRef | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", + "1 | \n", + "1 | \n", + "GJ 3053 | \n", + "LHS 1140 | \n", + "J00449-152 | \n", + "00:44:59.33 | \n", + "-15:16:17.5 | \n", + "M4.5 V | \n", + "4.5 | \n", + "14.964 | \n", + "... | \n", + "0.84 | \n", + "0.80 | \n", + "378.9 | \n", + "4.3 | \n", + "4.2 | \n", + "8.88 | \n", + "0.19 | \n", + "RV | \n", + "Dit17 | \n", + "Dit17 Men19 Lil20 | \n", + "
| 1 | \n", + "2 | \n", + "1 | \n", + "GJ 3053 | \n", + "LHS 1140 | \n", + "J00449-152 | \n", + "00:44:59.33 | \n", + "-15:16:17.5 | \n", + "M4.5 V | \n", + "4.5 | \n", + "14.964 | \n", + "... | \n", + "0.81 | \n", + "0.74 | \n", + "708.9 | \n", + "8.0 | \n", + "7.8 | \n", + "25.90 | \n", + "3.45 | \n", + "RV | \n", + "Dit17 | \n", + "Dit17 Men19 Lil20 | \n", + "
| 2 | \n", + "3 | \n", + "2 | \n", + "TOI-1468 | \n", + "LSPM J0106+1913 | \n", + "J01066+192 | \n", + "01:06:36.98 | \n", + "+19:13:33.2 | \n", + "M3 V | \n", + "3.0 | \n", + "24.721 | \n", + "... | \n", + "1.05 | \n", + "0.92 | \n", + "682.2 | \n", + "7.4 | \n", + "6.9 | \n", + "9.00 | \n", + "6.30 | \n", + "RV | \n", + "Cha22 | \n", + "Cha22 | \n", + "
| 3 | \n", + "4 | \n", + "2 | \n", + "TOI-1468 | \n", + "LSPM J0106+1913 | \n", + "J01066+192 | \n", + "01:06:36.98 | \n", + "+19:13:33.2 | \n", + "M3 V | \n", + "3.0 | \n", + "24.721 | \n", + "... | \n", + "0.21 | \n", + "0.19 | \n", + "337.5 | \n", + "3.7 | \n", + "3.4 | \n", + "59.00 | \n", + "NaN | \n", + "RV | \n", + "Cha22 | \n", + "Cha22 | \n", + "
| 4 | \n", + "5 | \n", + "3 | \n", + "GJ 3090 | \n", + "TOI-177 | \n", + "NaN | \n", + "01:21:45.39 | \n", + "-46:42:51.8 | \n", + "M2 V | \n", + "2.0 | \n", + "22.454 | \n", + "... | \n", + "0.52 | \n", + "0.45 | \n", + "693.0 | \n", + "18.0 | \n", + "18.0 | \n", + "221.00 | \n", + "NaN | \n", + "RV | \n", + "Alm22a | \n", + "Alm22a | \n", + "
5 rows × 44 columns
\n", + "| \n", + " | mean | \n", + "sd | \n", + "hdi_3% | \n", + "hdi_97% | \n", + "mcse_mean | \n", + "mcse_sd | \n", + "ess_bulk | \n", + "ess_tail | \n", + "r_hat | \n", + "
|---|---|---|---|---|---|---|---|---|---|
| a | \n", + "1.023 | \n", + "0.006 | \n", + "1.012 | \n", + "1.035 | \n", + "0.0 | \n", + "0.0 | \n", + "2490.0 | \n", + "2676.0 | \n", + "1.0 | \n", + "
| b | \n", + "0.317 | \n", + "0.012 | \n", + "0.293 | \n", + "0.339 | \n", + "0.0 | \n", + "0.0 | \n", + "2481.0 | \n", + "2697.0 | \n", + "1.0 | \n", + "
| \n", + " | mean | \n", + "sd | \n", + "hdi_3% | \n", + "hdi_97% | \n", + "mcse_mean | \n", + "mcse_sd | \n", + "ess_bulk | \n", + "ess_tail | \n", + "r_hat | \n", + "
|---|---|---|---|---|---|---|---|---|---|
| z[33] | \n", + "0.502 | \n", + "0.500 | \n", + "0.000 | \n", + "1.000 | \n", + "0.245 | \n", + "0.187 | \n", + "4.0 | \n", + "4.0 | \n", + "5.58 | \n", + "
| z[34] | \n", + "0.484 | \n", + "0.500 | \n", + "0.000 | \n", + "1.000 | \n", + "0.239 | \n", + "0.182 | \n", + "4.0 | \n", + "4.0 | \n", + "3.74 | \n", + "
| z[35] | \n", + "0.314 | \n", + "0.464 | \n", + "0.000 | \n", + "1.000 | \n", + "0.205 | \n", + "0.154 | \n", + "5.0 | \n", + "5.0 | \n", + "2.23 | \n", + "
| z[36] | \n", + "0.676 | \n", + "0.468 | \n", + "0.000 | \n", + "1.000 | \n", + "0.150 | \n", + "0.109 | \n", + "10.0 | \n", + "10.0 | \n", + "1.33 | \n", + "
| a[0] | \n", + "1.084 | \n", + "0.045 | \n", + "1.015 | \n", + "1.142 | \n", + "0.021 | \n", + "0.016 | \n", + "5.0 | \n", + "33.0 | \n", + "2.06 | \n", + "
| a[1] | \n", + "1.037 | \n", + "0.055 | \n", + "0.960 | \n", + "1.134 | \n", + "0.027 | \n", + "0.021 | \n", + "5.0 | \n", + "28.0 | \n", + "2.41 | \n", + "
| b[0] | \n", + "0.346 | \n", + "0.024 | \n", + "0.295 | \n", + "0.378 | \n", + "0.011 | \n", + "0.008 | \n", + "6.0 | \n", + "32.0 | \n", + "1.88 | \n", + "
| b[1] | \n", + "0.365 | \n", + "0.064 | \n", + "0.287 | \n", + "0.477 | \n", + "0.031 | \n", + "0.024 | \n", + "5.0 | \n", + "32.0 | \n", + "2.08 | \n", + "
| p[0] | \n", + "0.486 | \n", + "0.101 | \n", + "0.303 | \n", + "0.679 | \n", + "0.016 | \n", + "0.012 | \n", + "39.0 | \n", + "1616.0 | \n", + "1.07 | \n", + "
| p[1] | \n", + "0.514 | \n", + "0.101 | \n", + "0.321 | \n", + "0.697 | \n", + "0.016 | \n", + "0.012 | \n", + "39.0 | \n", + "1616.0 | \n", + "1.07 | \n", + "
<xarray.DataArray 'b' (chain: 4, draw: 1000, b_dim_0: 2)>\n", + "array([[[0.2966302 , 0.3682176 ],\n", + " [0.29150693, 0.35987771],\n", + " [0.30444869, 0.35571687],\n", + " ...,\n", + " [0.32787692, 0.36133043],\n", + " [0.33148199, 0.36345292],\n", + " [0.34565882, 0.3666832 ]],\n", + "\n", + " [[0.37185024, 0.32911081],\n", + " [0.37669095, 0.34146443],\n", + " [0.36533733, 0.34656975],\n", + " ...,\n", + " [0.36158203, 0.33904753],\n", + " [0.36008476, 0.32199052],\n", + " [0.35990644, 0.31166547]],\n", + "\n", + " [[0.35191582, 0.45181092],\n", + " [0.33252437, 0.47010269],\n", + " [0.32975338, 0.46711415],\n", + " ...,\n", + " [0.33840952, 0.47062463],\n", + " [0.34040357, 0.47464959],\n", + " [0.33163708, 0.46341786]],\n", + "\n", + " [[0.36864732, 0.33666572],\n", + " [0.3651342 , 0.32946099],\n", + " [0.36438941, 0.32528019],\n", + " ...,\n", + " [0.36866852, 0.31601492],\n", + " [0.36751868, 0.31890441],\n", + " [0.36840261, 0.29711909]]])\n", + "Coordinates:\n", + " * chain (chain) int64 0 1 2 3\n", + " * draw (draw) int64 0 1 2 3 4 5 6 7 8 ... 992 993 994 995 996 997 998 999\n", + " * b_dim_0 (b_dim_0) int64 0 1