In [8]:
%matplotlib inline

In [9]:
import pandas as pd

In [10]:
background = pd.read_csv('/Users/user/Desktop/background.csv')
halos = pd.read_csv('/Users/user/Desktop/halos.csv')

In [11]:
background


Out[11]:
Unnamed: 0 ID RA Dec mag Mstar_obs z_obs eMod_int ePhi_int e1_int ... eMod ePhi kappa_halo gamma1_halo gamma2_halo g_halo e1_halo e2_halo eMod_halo ePhi_halo
0 0 0.0 0.024513 -0.024514 13.497907 3.834421e+11 1.3857 0.298210 35.522894 0.096862 ... 0.333865 38.316575 -0.000023 -0.067612 0.084970 (-0.0676106882294+0.0849684849947j) 0.019058 0.361238 0.361741 43.489989
1 1 1.0 0.024484 -0.024535 2.597448 7.917252e+11 1.3857 0.137932 -16.443335 0.115828 ... 0.098430 -0.174268 -0.000032 -0.047347 0.027848 (-0.0473451876871+0.0278475422844j) 0.068990 -0.047428 0.083720 -17.253507
2 2 2.0 0.024504 -0.024496 0.872094 5.288954e+11 1.3857 0.065941 61.708902 -0.036317 ... 0.104823 50.778468 -0.000024 -0.012544 0.106721 (-0.0125439131722+0.10671855771j) -0.048044 0.160893 0.167913 53.313015

3 rows × 28 columns


In [12]:
background.columns


Out[12]:
Index([u'Unnamed: 0', u'ID', u'RA', u'Dec', u'mag', u'Mstar_obs', u'z_obs',
       u'eMod_int', u'ePhi_int', u'e1_int', u'e2_int', u'strong_flag',
       u'kappa', u'gamma1', u'gamma2', u'g', u'e1', u'e2', u'eMod', u'ePhi',
       u'kappa_halo', u'gamma1_halo', u'gamma2_halo', u'g_halo', u'e1_halo',
       u'e2_halo', u'eMod_halo', u'ePhi_halo'],
      dtype='object')

In [13]:
halos


Out[13]:
Unnamed: 0 112009306000027 37000046000024 37014522009504 456000089002589 65000614000013 log-likelihood
0 0 1.368440e+12 1.876400e+13 1.782420e+12 5.606320e+12 8.288130e+11 2.340629

In [14]:
background['RA']


Out[14]:
0    0.024513
1    0.024484
2    0.024504
Name: RA, dtype: float64

In [15]:
background['e1_int']


Out[15]:
0    0.096862
1    0.115828
2   -0.036317
Name: e1_int, dtype: float64

In [16]:
background


Out[16]:
Unnamed: 0 ID RA Dec mag Mstar_obs z_obs eMod_int ePhi_int e1_int ... eMod ePhi kappa_halo gamma1_halo gamma2_halo g_halo e1_halo e2_halo eMod_halo ePhi_halo
0 0 0.0 0.024513 -0.024514 13.497907 3.834421e+11 1.3857 0.298210 35.522894 0.096862 ... 0.333865 38.316575 -0.000023 -0.067612 0.084970 (-0.0676106882294+0.0849684849947j) 0.019058 0.361238 0.361741 43.489989
1 1 1.0 0.024484 -0.024535 2.597448 7.917252e+11 1.3857 0.137932 -16.443335 0.115828 ... 0.098430 -0.174268 -0.000032 -0.047347 0.027848 (-0.0473451876871+0.0278475422844j) 0.068990 -0.047428 0.083720 -17.253507
2 2 2.0 0.024504 -0.024496 0.872094 5.288954e+11 1.3857 0.065941 61.708902 -0.036317 ... 0.104823 50.778468 -0.000024 -0.012544 0.106721 (-0.0125439131722+0.10671855771j) -0.048044 0.160893 0.167913 53.313015

3 rows × 28 columns


In [44]:
l0 = pd.read_csv('/Users/user/Desktop/lightcone_0.csv')

In [45]:
len(l0)


Out[45]:
472

In [40]:
l0


Out[40]:
Unnamed: 0 RA Dec z_obs Mhalo_obs Type GalID mag_SDSS_i Mh x y r phi z Mh_obs
0 15442 0.024418 -0.023977 2.096700 1.290990e+11 0 233002149000031 25.8436 1.290990e+11 0.327373 1.845184 1.874001 1.395203 2.096700 11.110923
1 15501 0.024783 -0.024626 2.048570 7.048780e+11 0 22000121000031 25.8783 7.048780e+11 -0.926342 -0.384194 1.002853 -2.748441 2.048570 11.848114
2 15502 0.024523 -0.024725 2.079910 5.267220e+11 0 233000373000031 24.0311 5.267220e+11 -0.032455 -0.724875 0.725601 -1.615539 2.079910 11.721581
3 15503 0.024137 -0.024825 2.134700 2.082790e+11 0 279000084035274 25.1953 2.082790e+11 1.292045 -1.069681 1.677378 -0.691522 2.134700 11.318645
4 15504 0.025030 -0.024506 2.029570 1.833200e+11 0 23003340009344 25.0323 1.833200e+11 -1.775900 0.027304 1.776110 3.126219 2.029570 11.263210
5 15505 0.025034 -0.024673 2.032430 2.547540e+11 0 23003340005709 25.3176 2.547540e+11 -1.790671 -0.545425 1.871895 -2.845928 2.032430 11.406121
6 15507 0.024600 -0.024896 2.103080 1.308200e+11 0 489021309000031 25.9926 1.308200e+11 -0.298797 -1.314105 1.347647 -1.794372 2.103080 11.116674
7 15508 0.024890 -0.024472 2.101870 2.495910e+11 0 489000062000242 24.5315 2.495910e+11 -1.294074 0.145562 1.302235 3.029580 2.101870 11.397229
8 15510 0.024678 -0.024242 2.030510 1.222130e+11 0 23005429000031 25.6129 1.222130e+11 -0.567555 0.934525 1.093369 2.116580 2.030510 11.087117
9 15511 0.024776 -0.024615 2.095920 4.200010e+11 0 22004065000031 24.3258 4.200010e+11 -0.902629 -0.348786 0.967672 -2.772856 2.095920 11.623250
10 15512 0.024787 -0.024349 2.102610 2.022540e+11 0 278017644000031 23.3518 2.022540e+11 -0.940782 0.569093 1.099517 2.597567 2.102610 11.305897
11 15514 0.024287 -0.024971 2.110740 1.729920e+11 0 278017142000031 25.3238 1.729920e+11 0.778257 -1.570561 1.752811 -1.110733 2.110740 11.238026
12 15516 0.025074 -0.024664 2.092310 3.743860e+11 0 22004094000031 25.7742 3.743860e+11 -1.926764 -0.516548 1.994804 -2.879661 2.092310 11.573320
13 15518 0.024848 -0.024501 2.095910 1.919270e+11 0 22004065000093 25.5163 1.919270e+11 -1.149731 0.044836 1.150605 3.102615 2.095910 11.283136
14 15519 0.024850 -0.024216 2.144380 1.179100e+11 0 279000203000157 25.7667 1.179100e+11 -1.157643 1.024594 1.545940 2.417088 2.144380 11.071551
15 15833 0.025082 -0.024452 2.101920 2.065580e+11 0 489000062000174 24.9722 2.065580e+11 -1.953237 0.213974 1.964922 3.032480 2.101920 11.315042
16 15834 0.024798 -0.024517 2.080750 2.291070e+12 0 233000140000568 24.6521 2.291070e+12 -0.977895 -0.010167 0.977948 -3.131196 2.080750 12.360038
17 15835 0.024796 -0.024495 2.014490 2.323770e+11 0 23003446000776 24.7610 2.323770e+11 -0.969991 0.065119 0.972174 3.074559 2.014490 11.366193
18 15836 0.024181 -0.024708 2.037520 1.704100e+11 0 233022813004200 25.9590 1.704100e+11 1.141521 -0.667808 1.322512 -0.529329 2.037520 11.231495
19 15837 0.024980 -0.024685 2.145380 1.075820e+11 0 279000153000589 25.7114 1.075820e+11 -1.603026 -0.587021 1.707128 -2.790563 2.145380 11.031740
20 15838 0.024312 -0.024850 2.078860 2.805740e+11 0 233000373000206 25.6916 2.805740e+11 0.691655 -1.155281 1.346499 -1.031341 2.078860 11.448047
21 15839 0.024313 -0.024784 2.135660 9.983620e+10 0 279000084031265 25.8863 9.983620e+10 0.688907 -0.928390 1.156071 -0.932406 2.135660 10.999288
22 15840 0.024260 -0.024850 2.147000 4.440990e+11 0 279000457000031 24.8868 4.440990e+11 0.871050 -1.155281 1.446859 -0.924758 2.147000 11.647480
23 15841 0.024426 -0.024725 2.135740 4.621730e+11 0 279000084031172 25.3267 4.621730e+11 0.299186 -0.724875 0.784192 -1.179354 2.135740 11.664805
24 15843 0.024377 -0.024754 2.146010 2.366810e+11 0 279000153000663 24.9864 2.366810e+11 0.469646 -0.825257 0.949535 -1.053414 2.146010 11.374163
25 15844 0.024916 -0.024194 2.145480 3.124180e+11 0 279000203000031 25.4718 3.124180e+11 -1.385500 1.100225 1.769210 2.470460 2.145480 11.494736
26 15845 0.024440 -0.024393 2.004500 1.686890e+11 0 23005379000031 25.5336 1.686890e+11 0.251075 0.416113 0.485992 1.027894 2.004500 11.227087
27 15846 0.025007 -0.024262 2.134340 1.299590e+11 0 279000084041477 25.7212 1.299590e+11 -1.695491 0.866458 1.904059 2.669155 2.134340 11.113806
28 15847 0.024354 -0.024721 2.135120 4.819680e+11 0 279000084035430 25.1659 4.819680e+11 0.548347 -0.710093 0.897171 -0.913226 2.135120 11.683018
29 15849 0.024673 -0.024465 2.113540 4.217220e+11 0 279000304000031 25.7822 4.217220e+11 -0.549337 0.167220 0.574225 2.846101 2.113540 11.625026
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
874 276123 0.024434 -0.024314 3.109770 1.101640e+11 0 131003593003333 25.8618 1.101640e+11 0.271352 0.686320 0.738016 1.194286 3.109770 11.042040
875 276124 0.024054 -0.024597 3.146770 4.147510e+12 0 130003564001180 24.9906 4.147510e+12 1.578330 -0.285875 1.604010 -0.179182 3.146770 12.617787
876 276125 0.024807 -0.024048 3.166930 1.213530e+11 0 130004391000036 25.8914 1.213530e+11 -1.010555 1.601104 1.893345 2.133814 3.166930 11.084051
877 276126 0.024787 -0.024609 2.915100 1.471720e+11 0 383004539000036 25.5683 1.471720e+11 -0.939401 -0.327128 0.994730 -2.806496 2.915100 11.167825
878 276127 0.024598 -0.024426 3.063190 9.983620e+10 0 128008291001021 25.9441 9.983620e+10 -0.290553 0.303011 0.419806 2.335208 3.063190 10.999288
879 276134 0.024553 -0.023950 3.117030 1.359840e+11 0 131003681000721 25.9289 1.359840e+11 -0.137277 1.938004 1.942859 1.641512 3.117030 11.133488
880 276150 0.024571 -0.024772 3.130880 1.179100e+11 1 131026328000119 25.9025 1.179100e+11 -0.198447 -0.886105 0.908055 -1.791115 3.130880 11.071551
881 276151 0.024009 -0.024413 3.139040 2.177460e+11 0 130003679000604 25.7894 2.177460e+11 1.732646 0.348046 1.767257 0.198237 3.139040 11.337950
882 276152 0.024636 -0.025040 3.159070 2.272130e+11 0 130004049000036 25.5852 2.272130e+11 -0.420798 -1.807422 1.855760 -1.799538 3.159070 11.356433
883 276153 0.024890 -0.024500 3.062630 2.194680e+11 1 128008291000109 25.4513 2.194680e+11 -1.292698 0.048274 1.293599 3.104266 3.062630 11.341371
884 276154 0.024148 -0.024394 3.136120 2.022540e+11 0 130003966000036 25.9644 2.022540e+11 1.253911 0.410956 1.319537 0.316708 3.136120 11.305897
885 276155 0.024376 -0.024363 3.150660 3.451230e+11 0 130003564010421 24.9168 3.451230e+11 0.470338 0.519933 0.701105 0.835439 3.150660 11.537974
886 276158 0.024667 -0.024662 3.121680 1.256560e+11 0 131026312003423 25.7723 1.256560e+11 -0.528714 -0.509328 0.734135 -2.374868 3.121680 11.099183
887 276159 0.024077 -0.024568 3.117620 1.170490e+11 0 131003661001261 25.9378 1.170490e+11 1.500661 -0.185149 1.512040 -0.122758 3.117620 11.068368
888 276160 0.024118 -0.024918 3.048830 1.377050e+11 0 128009641000036 25.4321 1.377050e+11 1.359744 -1.388016 1.943063 -0.795687 3.048830 11.138950
889 276161 0.024703 -0.024520 2.932030 1.437300e+11 0 383004144000036 25.7251 1.437300e+11 -0.653125 -0.019449 0.653415 -3.111822 2.932030 11.157547
890 276162 0.024392 -0.024777 3.117990 2.693860e+11 0 131003909000074 25.5801 2.693860e+11 0.418439 -0.904325 0.996442 -1.137424 3.117990 11.430375
891 276163 0.024383 -0.024662 3.148930 8.606570e+11 0 130003564008606 24.8742 8.606570e+11 0.449027 -0.510016 0.679516 -0.848906 3.148930 11.934830
892 276164 0.024556 -0.024736 3.129790 1.590490e+12 0 131026328000036 25.3305 1.590490e+12 -0.145866 -0.764409 0.778202 -1.759351 3.129790 12.201531
893 276165 0.024383 -0.024969 3.114140 8.993870e+11 0 131003679000036 24.5913 8.993870e+11 0.448680 -1.565748 1.628766 -1.291715 3.114140 11.953947
894 276166 0.024384 -0.024935 3.118380 3.537300e+11 0 131003909000036 24.2076 3.537300e+11 0.445931 -1.448177 1.515279 -1.272084 3.118380 11.548672
895 276167 0.024697 -0.024125 3.108300 1.988120e+11 0 131003593002003 25.5640 1.988120e+11 -0.631136 1.338461 1.479801 2.011417 3.108300 11.298443
896 276168 0.024757 -0.024606 3.122630 9.897560e+10 0 131026312003455 25.9330 9.897560e+10 -0.838019 -0.317158 0.896027 -2.779790 3.122630 10.995528
897 276170 0.024682 -0.024660 3.129570 1.112830e+12 0 131026328000079 24.7004 1.112830e+12 -0.579921 -0.503484 0.767988 -2.426630 3.129570 12.046429
898 276171 0.024770 -0.024287 3.149050 1.185120e+12 0 130003564008704 25.6769 1.185120e+12 -0.882703 0.782233 1.179429 2.416466 3.149050 12.073762
899 276172 0.024734 -0.024204 3.108560 1.669670e+11 0 131003593001481 25.6267 1.669670e+11 -0.757951 1.064472 1.306748 2.189560 3.108560 11.222631
900 279113 0.024932 -0.024443 0.097020 2.065580e+10 0 501030332000004 24.8241 2.065580e+10 -1.437730 0.243538 1.458210 2.973795 0.097020 10.315042
901 281978 0.023986 -0.024453 0.272747 2.151640e+10 0 54017844000010 25.8890 2.151640e+10 1.813063 0.208129 1.824970 0.114294 0.272747 10.332770
902 282450 0.023991 -0.024400 0.279547 2.840170e+10 1 457000143001683 24.7406 2.840170e+10 1.795882 0.393768 1.838544 0.215846 0.279547 10.453344
903 285742 0.024393 -0.024914 0.160179 7.057390e+10 0 449008977000006 22.1766 7.057390e+10 0.412939 -1.376328 1.436940 -1.279312 0.160179 10.848644

904 rows × 15 columns


In [26]:
l1 = pd.read_csv('/Users/user/Desktop/lightcone_1.csv')

In [27]:
l1


Out[27]:
Unnamed: 0 RA Dec z_obs Mhalo_obs Type GalID mag_SDSS_i Mh x y r phi z Mh_obs
0 15442 0.024418 -0.023977 2.096700 1.290990e+11 0 233002149000031 25.8436 1.290990e+11 0.227577 1.917921 1.931376 1.452691 2.096700 11.110923
1 15501 0.024783 -0.024626 2.048570 7.048780e+11 0 22000121000031 25.8783 7.048780e+11 -1.026137 -0.311457 1.072363 -2.846906 2.048570 11.848114
2 15502 0.024523 -0.024725 2.079910 5.267220e+11 0 233000373000031 24.0311 5.267220e+11 -0.132250 -0.652138 0.665413 -1.770877 2.079910 11.721581
3 15503 0.024137 -0.024825 2.134700 2.082790e+11 0 279000084035274 25.1953 2.082790e+11 1.192250 -0.996944 1.554142 -0.696420 2.134700 11.318645
4 15504 0.025030 -0.024506 2.029570 1.833200e+11 0 23003340009344 25.0323 1.833200e+11 -1.875696 0.100041 1.878362 3.088308 2.029570 11.263210
5 15505 0.025034 -0.024673 2.032430 2.547540e+11 0 23003340005709 25.3176 2.547540e+11 -1.890466 -0.472688 1.948665 -2.896579 2.032430 11.406121
6 15507 0.024600 -0.024896 2.103080 1.308200e+11 0 489021309000031 25.9926 1.308200e+11 -0.398592 -1.241368 1.303790 -1.881489 2.103080 11.116674
7 15508 0.024890 -0.024472 2.101870 2.495910e+11 0 489000062000242 24.5315 2.495910e+11 -1.393869 0.218299 1.410860 2.986241 2.101870 11.397229
8 15510 0.024678 -0.024242 2.030510 1.222130e+11 0 23005429000031 25.6129 1.222130e+11 -0.667351 1.007262 1.208278 2.155936 2.030510 11.087117
9 15511 0.024776 -0.024615 2.095920 4.200010e+11 0 22004065000031 24.3258 4.200010e+11 -1.002424 -0.276049 1.039739 -2.872872 2.095920 11.623250
10 15512 0.024787 -0.024349 2.102610 2.022540e+11 0 278017644000031 23.3518 2.022540e+11 -1.040578 0.641830 1.222599 2.588911 2.102610 11.305897
11 15514 0.024287 -0.024971 2.110740 1.729920e+11 0 278017142000031 25.3238 1.729920e+11 0.678463 -1.497824 1.644320 -1.145479 2.110740 11.238026
12 15518 0.024848 -0.024501 2.095910 1.919270e+11 0 22004065000093 25.5163 1.919270e+11 -1.249526 0.117573 1.255046 3.047774 2.095910 11.283136
13 15519 0.024850 -0.024216 2.144380 1.179100e+11 0 279000203000157 25.7667 1.179100e+11 -1.257440 1.097331 1.668919 2.424083 2.144380 11.071551
14 15834 0.024798 -0.024517 2.080750 2.291070e+12 0 233000140000568 24.6521 2.291070e+12 -1.077690 0.062570 1.079505 3.083599 2.080750 12.360038
15 15835 0.024796 -0.024495 2.014490 2.323770e+11 0 23003446000776 24.7610 2.323770e+11 -1.069786 0.137856 1.078632 3.013436 2.014490 11.366193
16 15836 0.024181 -0.024708 2.037520 1.704100e+11 0 233022813004200 25.9590 1.704100e+11 1.041726 -0.595071 1.199710 -0.519001 2.037520 11.231495
17 15837 0.024980 -0.024685 2.145380 1.075820e+11 0 279000153000589 25.7114 1.075820e+11 -1.702821 -0.514284 1.778789 -2.848285 2.145380 11.031740
18 15838 0.024312 -0.024850 2.078860 2.805740e+11 0 233000373000206 25.6916 2.805740e+11 0.591860 -1.082544 1.233775 -1.070466 2.078860 11.448047
19 15839 0.024313 -0.024784 2.135660 9.983620e+10 0 279000084031265 25.8863 9.983620e+10 0.589112 -0.855653 1.038843 -0.967834 2.135660 10.999288
20 15840 0.024260 -0.024850 2.147000 4.440990e+11 0 279000457000031 24.8868 4.440990e+11 0.771255 -1.082544 1.329186 -0.951765 2.147000 11.647480
21 15841 0.024426 -0.024725 2.135740 4.621730e+11 0 279000084031172 25.3267 4.621730e+11 0.199391 -0.652138 0.681939 -1.274072 2.135740 11.664805
22 15843 0.024377 -0.024754 2.146010 2.366810e+11 0 279000153000663 24.9864 2.366810e+11 0.369851 -0.752520 0.838497 -1.113985 2.146010 11.374163
23 15844 0.024916 -0.024194 2.145480 3.124180e+11 0 279000203000031 25.4718 3.124180e+11 -1.485296 1.172962 1.892602 2.473154 2.145480 11.494736
24 15845 0.024440 -0.024393 2.004500 1.686890e+11 0 23005379000031 25.5336 1.686890e+11 0.151279 0.488850 0.511722 1.270685 2.004500 11.227087
25 15847 0.024354 -0.024721 2.135120 4.819680e+11 0 279000084035430 25.1659 4.819680e+11 0.448552 -0.637356 0.779372 -0.957544 2.135120 11.683018
26 15849 0.024673 -0.024465 2.113540 4.217220e+11 0 279000304000031 25.7822 4.217220e+11 -0.649133 0.239957 0.692064 2.787513 2.113540 11.625026
27 15850 0.024037 -0.024256 2.096260 1.979510e+11 0 233000736000031 25.3426 1.979510e+11 1.538002 0.959821 1.812928 0.557930 2.096260 11.296558
28 21016 0.023918 -0.024634 2.198190 2.969270e+11 0 276013253001046 24.5082 2.969270e+11 1.946268 -0.338616 1.975505 -0.172258 2.198190 11.472650
29 25176 0.023980 -0.024335 0.678461 9.295100e+10 0 82003711000018 23.9289 9.295100e+10 1.733205 0.687552 1.864599 0.377653 0.678461 10.968254
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
880 276122 0.024235 -0.024129 2.913160 1.454510e+11 0 383004761000036 25.5790 1.454510e+11 0.858219 1.397790 1.640231 1.020159 2.913160 11.162717
881 276123 0.024434 -0.024314 3.109770 1.101640e+11 0 131003593003333 25.8618 1.101640e+11 0.171556 0.759057 0.778202 1.348519 3.109770 11.042040
882 276124 0.024054 -0.024597 3.146770 4.147510e+12 0 130003564001180 24.9906 4.147510e+12 1.478534 -0.213138 1.493818 -0.143168 3.146770 12.617787
883 276126 0.024787 -0.024609 2.915100 1.471720e+11 0 383004539000036 25.5683 1.471720e+11 -1.039197 -0.254391 1.069881 -2.901518 2.915100 11.167825
884 276127 0.024598 -0.024426 3.063190 9.983620e+10 0 128008291001021 25.9441 9.983620e+10 -0.390348 0.375748 0.541811 2.375250 3.063190 10.999288
885 276150 0.024571 -0.024772 3.130880 1.179100e+11 1 131026328000119 25.9025 1.179100e+11 -0.298242 -0.813368 0.866323 -1.922248 3.130880 11.071551
886 276151 0.024009 -0.024413 3.139040 2.177460e+11 0 130003679000604 25.7894 2.177460e+11 1.632850 0.420783 1.686196 0.252211 3.139040 11.337950
887 276152 0.024636 -0.025040 3.159070 2.272130e+11 0 130004049000036 25.5852 2.272130e+11 -0.520592 -1.734684 1.811117 -1.862352 3.159070 11.356433
888 276153 0.024890 -0.024500 3.062630 2.194680e+11 1 128008291000109 25.4513 2.194680e+11 -1.392494 0.121011 1.397742 3.054908 3.062630 11.341371
889 276154 0.024148 -0.024394 3.136120 2.022540e+11 0 130003966000036 25.9644 2.022540e+11 1.154115 0.483694 1.251376 0.396866 3.136120 11.305897
890 276155 0.024376 -0.024363 3.150660 3.451230e+11 0 130003564010421 24.9168 3.451230e+11 0.370542 0.592670 0.698970 1.012047 3.150660 11.537974
891 276158 0.024667 -0.024662 3.121680 1.256560e+11 0 131026312003423 25.7723 1.256560e+11 -0.628510 -0.436591 0.765269 -2.534469 3.121680 11.099183
892 276159 0.024077 -0.024568 3.117620 1.170490e+11 0 131003661001261 25.9378 1.170490e+11 1.400866 -0.112412 1.405369 -0.080073 3.117620 11.068368
893 276160 0.024118 -0.024918 3.048830 1.377050e+11 0 128009641000036 25.4321 1.377050e+11 1.259950 -1.315279 1.821382 -0.806880 3.048830 11.138950
894 276161 0.024703 -0.024520 2.932030 1.437300e+11 0 383004144000036 25.7251 1.437300e+11 -0.752921 0.053288 0.754804 3.070936 2.932030 11.157547
895 276162 0.024392 -0.024777 3.117990 2.693860e+11 0 131003909000074 25.5801 2.693860e+11 0.318644 -0.831588 0.890547 -1.204877 3.117990 11.430375
896 276163 0.024383 -0.024662 3.148930 8.606570e+11 0 130003564008606 24.8742 8.606570e+11 0.349232 -0.437279 0.559621 -0.896880 3.148930 11.934830
897 276164 0.024556 -0.024736 3.129790 1.590490e+12 0 131026328000036 25.3305 1.590490e+12 -0.245661 -0.691672 0.734002 -1.912069 3.129790 12.201531
898 276165 0.024383 -0.024969 3.114140 8.993870e+11 0 131003679000036 24.5913 8.993870e+11 0.348886 -1.493011 1.533233 -1.341236 3.114140 11.953947
899 276166 0.024384 -0.024935 3.118380 3.537300e+11 0 131003909000036 24.2076 3.537300e+11 0.346137 -1.375440 1.418325 -1.324260 3.118380 11.548672
900 276167 0.024697 -0.024125 3.108300 1.988120e+11 0 131003593002003 25.5640 1.988120e+11 -0.730933 1.411198 1.589258 2.048702 3.108300 11.298443
901 276168 0.024757 -0.024606 3.122630 9.897560e+10 0 131026312003455 25.9330 9.897560e+10 -0.937814 -0.244421 0.969142 -2.886636 3.122630 10.995528
902 276169 0.024270 -0.025066 3.113350 2.203280e+11 0 131003679000163 25.5162 2.203280e+11 0.737915 -1.826129 1.969585 -1.186771 3.113350 11.343070
903 276170 0.024682 -0.024660 3.129570 1.112830e+12 0 131026328000079 24.7004 1.112830e+12 -0.679716 -0.430747 0.804710 -2.576750 3.129570 12.046429
904 276171 0.024770 -0.024287 3.149050 1.185120e+12 0 130003564008704 25.6769 1.185120e+12 -0.982499 0.854970 1.302412 2.425488 3.149050 12.073762
905 276172 0.024734 -0.024204 3.108560 1.669670e+11 0 131003593001481 25.6267 1.669670e+11 -0.857747 1.137209 1.424421 2.217016 3.108560 11.222631
906 279113 0.024932 -0.024443 0.097020 2.065580e+10 0 501030332000004 24.8241 2.065580e+10 -1.537525 0.316275 1.569718 2.938718 0.097020 10.315042
907 281978 0.023986 -0.024453 0.272747 2.151640e+10 0 54017844000010 25.8890 2.151640e+10 1.713268 0.280866 1.736137 0.162491 0.272747 10.332770
908 282450 0.023991 -0.024400 0.279547 2.840170e+10 1 457000143001683 24.7406 2.840170e+10 1.696086 0.466505 1.759072 0.268411 0.279547 10.453344
909 285742 0.024393 -0.024914 0.160179 7.057390e+10 0 449008977000006 22.1766 7.057390e+10 0.313145 -1.303591 1.340675 -1.335046 0.160179 10.848644

910 rows × 15 columns


In [28]:
l2 = pd.read_csv('/Users/user/Desktop/lightcone_2.csv')
l2


Out[28]:
Unnamed: 0 RA Dec z_obs Mhalo_obs Type GalID mag_SDSS_i Mh x y r phi z Mh_obs
0 15442 0.024418 -0.023977 2.096700 1.290990e+11 0 233002149000031 25.8436 1.290990e+11 0.295013 1.784578 1.808798 1.406966 2.096700 11.110923
1 15501 0.024783 -0.024626 2.048570 7.048780e+11 0 22000121000031 25.8783 7.048780e+11 -0.958702 -0.444801 1.056862 -2.707189 2.048570 11.848114
2 15502 0.024523 -0.024725 2.079910 5.267220e+11 0 233000373000031 24.0311 5.267220e+11 -0.064815 -0.785482 0.788151 -1.653126 2.079910 11.721581
3 15503 0.024137 -0.024825 2.134700 2.082790e+11 0 279000084035274 25.1953 2.082790e+11 1.259685 -1.130288 1.692441 -0.731309 2.134700 11.318645
4 15504 0.025030 -0.024506 2.029570 1.833200e+11 0 23003340009344 25.0323 1.833200e+11 -1.808260 -0.033303 1.808567 -3.123178 2.029570 11.263210
5 15505 0.025034 -0.024673 2.032430 2.547540e+11 0 23003340005709 25.3176 2.547540e+11 -1.823031 -0.606031 1.921123 -2.820655 2.032430 11.406121
6 15507 0.024600 -0.024896 2.103080 1.308200e+11 0 489021309000031 25.9926 1.308200e+11 -0.331157 -1.374711 1.414036 -1.807185 2.103080 11.116674
7 15508 0.024890 -0.024472 2.101870 2.495910e+11 0 489000062000242 24.5315 2.495910e+11 -1.326434 0.084956 1.329152 3.077632 2.101870 11.397229
8 15510 0.024678 -0.024242 2.030510 1.222130e+11 0 23005429000031 25.6129 1.222130e+11 -0.599915 0.873919 1.060015 2.172381 2.030510 11.087117
9 15511 0.024776 -0.024615 2.095920 4.200010e+11 0 22004065000031 24.3258 4.200010e+11 -0.934989 -0.409392 1.020689 -2.728882 2.095920 11.623250
10 15512 0.024787 -0.024349 2.102610 2.022540e+11 0 278017644000031 23.3518 2.022540e+11 -0.973142 0.508486 1.097982 2.660092 2.102610 11.305897
11 15514 0.024287 -0.024971 2.110740 1.729920e+11 0 278017142000031 25.3238 1.729920e+11 0.745897 -1.631167 1.793619 -1.141906 2.110740 11.238026
12 15518 0.024848 -0.024501 2.095910 1.919270e+11 0 22004065000093 25.5163 1.919270e+11 -1.182091 -0.015770 1.182196 -3.128253 2.095910 11.283136
13 15519 0.024850 -0.024216 2.144380 1.179100e+11 0 279000203000157 25.7667 1.179100e+11 -1.190004 0.963988 1.531464 2.460741 2.144380 11.071551
14 15833 0.025082 -0.024452 2.101920 2.065580e+11 0 489000062000174 24.9722 2.065580e+11 -1.985597 0.153367 1.991512 3.064506 2.101920 11.315042
15 15834 0.024798 -0.024517 2.080750 2.291070e+12 0 233000140000568 24.6521 2.291070e+12 -1.010255 -0.070774 1.012731 -3.071651 2.080750 12.360038
16 15835 0.024796 -0.024495 2.014490 2.323770e+11 0 23003446000776 24.7610 2.323770e+11 -1.002351 0.004513 1.002361 3.137091 2.014490 11.366193
17 15836 0.024181 -0.024708 2.037520 1.704100e+11 0 233022813004200 25.9590 1.704100e+11 1.109161 -0.728415 1.326961 -0.581089 2.037520 11.231495
18 15837 0.024980 -0.024685 2.145380 1.075820e+11 0 279000153000589 25.7114 1.075820e+11 -1.635386 -0.647628 1.758952 -2.764531 2.145380 11.031740
19 15838 0.024312 -0.024850 2.078860 2.805740e+11 0 233000373000206 25.6916 2.805740e+11 0.659295 -1.215887 1.383131 -1.073935 2.078860 11.448047
20 15839 0.024313 -0.024784 2.135660 9.983620e+10 0 279000084031265 25.8863 9.983620e+10 0.656547 -0.988996 1.187083 -0.984745 2.135660 10.999288
21 15840 0.024260 -0.024850 2.147000 4.440990e+11 0 279000457000031 24.8868 4.440990e+11 0.838690 -1.215887 1.477086 -0.966965 2.147000 11.647480
22 15841 0.024426 -0.024725 2.135740 4.621730e+11 0 279000084031172 25.3267 4.621730e+11 0.266826 -0.785482 0.829565 -1.243329 2.135740 11.664805
23 15843 0.024377 -0.024754 2.146010 2.366810e+11 0 279000153000663 24.9864 2.366810e+11 0.437286 -0.885864 0.987914 -1.112260 2.146010 11.374163
24 15844 0.024916 -0.024194 2.145480 3.124180e+11 0 279000203000031 25.4718 3.124180e+11 -1.417860 1.039618 1.758162 2.508911 2.145480 11.494736
25 15845 0.024440 -0.024393 2.004500 1.686890e+11 0 23005379000031 25.5336 1.686890e+11 0.218714 0.355507 0.417398 1.019261 2.004500 11.227087
26 15846 0.025007 -0.024262 2.134340 1.299590e+11 0 279000084041477 25.7212 1.299590e+11 -1.727852 0.805851 1.906533 2.705193 2.134340 11.113806
27 15847 0.024354 -0.024721 2.135120 4.819680e+11 0 279000084035430 25.1659 4.819680e+11 0.515987 -0.770699 0.927480 -0.980831 2.135120 11.683018
28 15849 0.024673 -0.024465 2.113540 4.217220e+11 0 279000304000031 25.7822 4.217220e+11 -0.581698 0.106614 0.591387 2.960324 2.113540 11.625026
29 15850 0.024037 -0.024256 2.096260 1.979510e+11 0 233000736000031 25.3426 1.979510e+11 1.605438 0.826478 1.805685 0.475417 2.096260 11.296558
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
865 276123 0.024434 -0.024314 3.109770 1.101640e+11 0 131003593003333 25.8618 1.101640e+11 0.238992 0.625713 0.669802 1.205946 3.109770 11.042040
866 276124 0.024054 -0.024597 3.146770 4.147510e+12 0 130003564001180 24.9906 4.147510e+12 1.545970 -0.346481 1.584320 -0.220476 3.146770 12.617787
867 276125 0.024807 -0.024048 3.166930 1.213530e+11 0 130004391000036 25.8914 1.213530e+11 -1.042916 1.540498 1.860324 2.165918 3.166930 11.084051
868 276126 0.024787 -0.024609 2.915100 1.471720e+11 0 383004539000036 25.5683 1.471720e+11 -0.971762 -0.387734 1.046259 -2.761947 2.915100 11.167825
869 276127 0.024598 -0.024426 3.063190 9.983620e+10 0 128008291001021 25.9441 9.983620e+10 -0.322913 0.242405 0.403773 2.497656 3.063190 10.999288
870 276134 0.024553 -0.023950 3.117030 1.359840e+11 0 131003681000721 25.9289 1.359840e+11 -0.169637 1.877397 1.885045 1.660909 3.117030 11.133488
871 276150 0.024571 -0.024772 3.130880 1.179100e+11 1 131026328000119 25.9025 1.179100e+11 -0.230807 -0.946712 0.974441 -1.809930 3.130880 11.071551
872 276151 0.024009 -0.024413 3.139040 2.177460e+11 0 130003679000604 25.7894 2.177460e+11 1.700285 0.287439 1.724411 0.167470 3.139040 11.337950
873 276152 0.024636 -0.025040 3.159070 2.272130e+11 0 130004049000036 25.5852 2.272130e+11 -0.453158 -1.868028 1.922207 -1.808785 3.159070 11.356433
874 276153 0.024890 -0.024500 3.062630 2.194680e+11 1 128008291000109 25.4513 2.194680e+11 -1.325058 -0.012332 1.325116 -3.132286 3.062630 11.341371
875 276154 0.024148 -0.024394 3.136120 2.022540e+11 0 130003966000036 25.9644 2.022540e+11 1.221551 0.350350 1.270799 0.279310 3.136120 11.305897
876 276155 0.024376 -0.024363 3.150660 3.451230e+11 0 130003564010421 24.9168 3.451230e+11 0.437978 0.459327 0.634670 0.809186 3.150660 11.537974
877 276158 0.024667 -0.024662 3.121680 1.256560e+11 0 131026312003423 25.7723 1.256560e+11 -0.561074 -0.569935 0.799769 -2.348361 3.121680 11.099183
878 276159 0.024077 -0.024568 3.117620 1.170490e+11 0 131003661001261 25.9378 1.170490e+11 1.468301 -0.245755 1.488726 -0.165837 3.117620 11.068368
879 276160 0.024118 -0.024918 3.048830 1.377050e+11 0 128009641000036 25.4321 1.377050e+11 1.327384 -1.448623 1.964805 -0.829044 3.048830 11.138950
880 276161 0.024703 -0.024520 2.932030 1.437300e+11 0 383004144000036 25.7251 1.437300e+11 -0.685485 -0.080056 0.690144 -3.025332 2.932030 11.157547
881 276162 0.024392 -0.024777 3.117990 2.693860e+11 0 131003909000074 25.5801 2.693860e+11 0.386079 -0.964932 1.039303 -1.190195 3.117990 11.430375
882 276163 0.024383 -0.024662 3.148930 8.606570e+11 0 130003564008606 24.8742 8.606570e+11 0.416667 -0.570622 0.706556 -0.940090 3.148930 11.934830
883 276164 0.024556 -0.024736 3.129790 1.590490e+12 0 131026328000036 25.3305 1.590490e+12 -0.178226 -0.825016 0.844047 -1.783554 3.129790 12.201531
884 276165 0.024383 -0.024969 3.114140 8.993870e+11 0 131003679000036 24.5913 8.993870e+11 0.416320 -1.626354 1.678795 -1.320194 3.114140 11.953947
885 276166 0.024384 -0.024935 3.118380 3.537300e+11 0 131003909000036 24.2076 3.537300e+11 0.413571 -1.508784 1.564439 -1.303259 3.118380 11.548672
886 276167 0.024697 -0.024125 3.108300 1.988120e+11 0 131003593002003 25.5640 1.988120e+11 -0.663497 1.277854 1.439840 2.049707 3.108300 11.298443
887 276168 0.024757 -0.024606 3.122630 9.897560e+10 0 131026312003455 25.9330 9.897560e+10 -0.870379 -0.377765 0.948823 -2.732104 3.122630 10.995528
888 276170 0.024682 -0.024660 3.129570 1.112830e+12 0 131026328000079 24.7004 1.112830e+12 -0.612281 -0.564091 0.832518 -2.397137 3.129570 12.046429
889 276171 0.024770 -0.024287 3.149050 1.185120e+12 0 130003564008704 25.6769 1.185120e+12 -0.915063 0.721627 1.165369 2.473836 3.149050 12.073762
890 276172 0.024734 -0.024204 3.108560 1.669670e+11 0 131003593001481 25.6267 1.669670e+11 -0.790311 1.003866 1.277630 2.237726 3.108560 11.222631
891 279113 0.024932 -0.024443 0.097020 2.065580e+10 0 501030332000004 24.8241 2.065580e+10 -1.470090 0.182932 1.481428 3.017793 0.097020 10.315042
892 281978 0.023986 -0.024453 0.272747 2.151640e+10 0 54017844000010 25.8890 2.151640e+10 1.780703 0.147523 1.786803 0.082657 0.272747 10.332770
893 282450 0.023991 -0.024400 0.279547 2.840170e+10 1 457000143001683 24.7406 2.840170e+10 1.763522 0.333161 1.794716 0.186718 0.279547 10.453344
894 285742 0.024393 -0.024914 0.160179 7.057390e+10 0 449008977000006 22.1766 7.057390e+10 0.380579 -1.436935 1.486480 -1.311886 0.160179 10.848644

895 rows × 15 columns


In [29]:
count = 287822
ra_max = -0.0174533
ra_min = -0.034898
dec_max = -0.0174533
dec_min = -0.0348981

xi = -0.0251909188704
xf = -0.0250163859452
yi = -0.0251909188704
yf = -0.0250163859452

In [30]:
area1 = (ra_max - ra_min)*(dec_max-dec_min)
area2 = (xf - xi)* (yf - yi)

In [31]:
area1


Out[31]:
0.00030431930255999995

In [32]:
area2


Out[32]:
3.0461741978869316e-08

In [33]:
frac = area2 / area1
ex_count = count * frac

In [34]:
ex_count


Out[34]:
28.81039561436791

In [37]:
out = pd.read_csv('/Users/user/Desktop/halos.csv')

In [38]:
out


Out[38]:
Unnamed: 0 111006726000565 111006726004073 111007167000028 111007406000027 111007703000027 111008078000027 112009253001465 112009253001583 112009306000027 ... 97001163000025 97020576000025 98000008000404 98003134000025 98016480000024 98016480000159 98017091000024 98018070000024 98018361000024 log-likelihood
0 0 2.174880e+12 1.729920e+11 4.595910e+11 3.347960e+11 2.168860e+11 1.798770e+11 1.093030e+11 1.351230e+11 1.368440e+12 ... 3.236070e+11 1.256560e+11 8.004110e+10 7.831980e+10 5.800830e+11 1.093030e+11 2.168860e+11 1.179100e+11 9.295100e+10 2.389705

1 rows × 546 columns


In [46]:
f = pd.read_csv('/Users/user/Desktop/foo.csv')
t = pd.read_csv('/Users/user/Desktop/foo2.csv')

In [48]:
t['r'].mean()


Out[48]:
1.3756341326660508

In [49]:
f['r'].mean()


Out[49]:
1.3756341326638073

In [50]:
print f['x'].mean()
print t['x'].mean()


0.0732333796174
0.0732333796496

In [51]:
print len(f)
print len(t)


472
472

In [54]:
diff = []
for i in f['x']:
    if i not in t['x']:
        diff.append(i)

In [55]:
diff2 = []
for i in t['x']:
    if i not in f['x']:
        diff2.append(i)

In [56]:
len(diff2)


Out[56]:
91

In [57]:
len(diff)


Out[57]:
91

In [74]:
diff[5]


Out[74]:
-1.2796314305700001

In [73]:
diff2[5]


Out[73]:
-1.2796314305399998

In [66]:
import numpy as np

In [67]:
np.cos(1000)


Out[67]:
0.56237907629070294

In [68]:
np.random.seed(111)
np.cos(1000)


Out[68]:
0.56237907629070294

In [85]:
np.cos(1002) * ((180.0 * 60.0) / np.pi)


Out[85]:
-3389.314804341921

In [86]:
np.cos(1002) * 180.0 * 60.0 / np.pi


Out[86]:
-3389.3148043419214

In [88]:
arcmin2rad = (1.0/60.0)*np.pi/180.0
arcmin2deg = (1.0/60.0)
rad2arcmin = 1.0/arcmin2rad

In [89]:
rad2arcmin


Out[89]:
3437.7467707849396

In [90]:
Angle.radian_to_arcmin(1)


---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-90-415a5b133522> in <module>()
----> 1 Angle.radian_to_arcmin(1)

NameError: name 'Angle' is not defined

In [91]:
diff[0:10]


Out[91]:
[-1.4813612131600002,
 -1.1878732729600001,
 -1.1483619360500001,
 -1.1452715354799998,
 -1.22773356119,
 -1.2796314305700001,
 -1.2469857679799998,
 -1.2878785685900001,
 -1.33324855801,
 -1.5831085207100002]

In [93]:
np.pi


Out[93]:
3.141592653589793

In [97]:
np.pi*10**15 + 0.2


Out[97]:
3141592653589793.0

In [100]:
pred = pd.read_csv('/Users/user/Desktop/pred.csv')

In [104]:
pred.columns


Out[104]:
Index([u'Unnamed: 0', u'0.0', u'0.0.1', u'1.0', u'1.0.1', u'2.0', u'2.0.1'], dtype='object')

In [106]:
np.zeros((5,5))


Out[106]:
array([[ 0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.],
       [ 0.,  0.,  0.,  0.,  0.]])

In [109]:
d = pd.read_table('/Users/user/Code/Pangloss/data/GGL_los_8_0_0_0_0_N_4096_ang_4_Guo_galaxies_on_plane_27_to_63.images.txt')

In [110]:
d


Out[110]:
GalID HaloID SubhaloID Type PlaneNumber z_spec pos_0[rad] pos_1[rad] Dc_los[Mpc/h] M_Halo[M_sol/h] M_Subhalo[M_sol/h] M_Stellar[M_sol/h] mag_SDSS_u mag_SDSS_g mag_SDSS_r mag_SDSS_i mag_SDSS_z mag_J mag_H mag_K
0 22005657000031 22005657000031 32002200007540 0 32 2.077090 -0.031597 -0.029582 3878.500 1.153280e+11 1.377050e+11 2.538560e+09 26.4515 25.9714 25.9779 25.9212 25.6947 24.8963 24.5919 24.3665
1 23003396000121 23003396000094 32002300002152 1 32 2.004520 -0.030872 -0.030504 3797.800 5.766400e+11 1.368440e+11 3.423400e+09 26.5567 26.0925 26.0241 25.9720 25.7234 24.7098 24.3246 24.0779
2 23003339006170 23003339002325 32002300002416 0 32 1.992640 -0.030766 -0.029958 3783.120 7.083200e+11 5.404930e+11 1.212330e+10 26.8170 25.8597 25.8092 25.5507 24.9120 23.9694 23.5292 23.0121
3 23003385000769 23003385000232 32002300005650 0 32 1.995730 -0.030627 -0.029594 3790.110 1.609430e+11 1.962300e+11 1.414890e+09 25.7759 25.6043 25.6101 25.6653 25.5959 24.7451 24.5413 24.4404
4 22000106000031 22000106000031 32002200006033 0 32 2.041820 -0.030769 -0.029931 3842.240 1.833200e+11 1.841810e+11 4.642430e+09 26.8309 25.9563 25.9070 25.6981 25.2301 24.3321 23.9834 23.6662
5 489020438001227 489020438001916 32049000000387 0 32 2.124030 -0.031282 -0.030422 3933.170 4.541690e+12 4.195700e+12 4.624150e+10 26.8694 26.1084 25.7770 25.5346 24.8123 23.2087 22.4581 22.0212
6 23005411000031 23005411000031 32002300010212 0 32 2.010650 -0.031274 -0.029770 3805.330 7.918040e+10 8.864770e+10 2.814590e+08 25.6097 25.5477 25.6489 25.7497 25.8087 25.3611 25.4079 25.4507
7 23003339000771 23003339001608 32002300000469 2 32 1.994780 -0.030579 -0.030002 3783.620 1.155350e+13 0.000000e+00 8.061260e+09 26.2996 25.5214 25.5081 25.3058 24.7935 23.8448 23.4391 23.0691
8 489020448004718 489020448001872 32049000005538 0 32 2.149090 -0.031530 -0.029740 3958.590 1.609430e+11 1.729920e+11 1.798840e+09 25.9459 25.5596 25.5285 25.5160 25.3935 24.6104 24.4181 24.2959
9 23003339001406 23003339001608 32002300000469 2 32 1.991430 -0.030633 -0.029939 3783.720 1.155350e+13 0.000000e+00 6.010610e+09 26.1708 25.5023 25.4905 25.3534 24.9863 24.0450 23.6678 23.3826
10 23003339001615 23003339001608 32002300000469 2 32 1.992220 -0.030625 -0.029962 3783.820 1.155350e+13 0.000000e+00 2.007910e+10 26.4384 25.4217 25.4436 25.1767 24.5282 23.4385 22.9097 22.4723
11 22000106000306 22000106000211 32002200007581 0 32 2.042320 -0.031023 -0.030046 3842.080 1.316800e+11 1.334020e+11 7.680040e+08 25.3776 25.1787 25.2130 25.2700 25.2612 24.6590 24.6098 24.5887
12 23003339001583 23003339001608 32002300000469 2 32 1.989740 -0.030629 -0.030004 3783.910 1.155350e+13 0.000000e+00 2.421100e+10 25.9685 24.9399 24.8657 24.6463 24.0885 22.7949 22.3892 22.0492
13 489020438013554 489020438005952 32048900007480 0 32 2.120850 -0.031017 -0.030232 3927.670 1.247950e+11 1.385660e+11 1.394700e+09 26.2199 25.7756 25.7657 25.7372 25.5883 24.8544 24.6581 24.5277
14 23003331032688 23003331013722 32002300001792 1 32 1.971940 -0.030582 -0.029902 3765.250 5.112300e+11 5.077880e+10 3.707350e+09 26.3215 25.6971 25.6555 25.5467 25.2308 24.2751 23.9523 23.7188
15 23003339009419 23003339005019 32002300001978 0 32 1.994530 -0.030581 -0.029466 3789.140 7.444680e+11 1.029350e+12 1.778600e+10 25.4687 24.9303 24.8125 24.7482 24.4083 23.2439 22.7254 22.4115
16 23003339012562 23003339005892 32002300004917 0 32 1.993080 -0.031255 -0.029557 3787.370 2.272130e+11 2.556150e+11 8.383700e+08 25.0957 24.9691 25.0335 25.1134 25.1321 24.5687 24.5333 24.5218
17 23003334013157 23003334005866 32002300005900 0 32 2.032420 -0.031379 -0.029523 3828.360 1.084430e+11 1.652460e+11 2.200240e+09 26.2523 25.8627 25.8435 25.8194 25.6261 24.7220 24.4238 24.2381
18 23003339014535 23003339008344 32002300006402 0 32 1.989440 -0.030716 -0.029958 3780.720 1.273770e+11 1.540580e+11 2.152270e+09 26.5568 25.8670 25.8247 25.7014 25.3951 24.5318 24.2884 24.0947
19 23003331037591 23003331016100 32002300002534 0 32 1.973850 -0.030712 -0.030123 3766.950 6.670090e+11 7.238130e+11 1.368020e+10 26.6578 25.6650 25.5853 25.3258 24.6909 23.5801 23.1579 22.7169
20 489020517000263 489020517000087 32048900005332 0 32 2.108280 -0.030969 -0.029490 3915.340 1.815990e+11 2.177460e+11 1.149430e+09 25.5293 25.3716 25.4112 25.4786 25.4504 24.7320 24.5819 24.5024
21 22004397000031 22004397000031 32002200005357 0 32 2.085230 -0.030965 -0.029882 3890.970 1.902050e+11 2.263530e+11 3.317590e+09 25.7862 25.3178 25.3180 25.2685 25.0621 24.2932 24.0094 23.8099
22 23003675000174 23003675000129 32002300007028 0 32 2.011700 -0.031436 -0.029787 3805.750 1.308200e+11 1.497540e+11 1.425220e+09 25.9897 25.4913 25.4997 25.4430 25.2314 24.5637 24.3676 24.2174
23 23003334016733 23003334007881 32002300003460 2 32 2.026620 -0.031426 -0.029794 3822.900 4.053690e+11 0.000000e+00 2.523990e+09 26.0858 25.5850 25.5670 25.5147 25.2943 24.4585 24.1877 24.0055
24 23005110000031 23005110000031 32002300006710 0 32 2.017230 -0.031429 -0.029611 3811.990 1.454510e+11 1.643850e+11 2.166410e+09 25.3981 25.0073 25.0154 25.0035 24.8634 24.1969 24.0014 23.8727
25 489020448001768 489020448001313 32049000007687 0 32 2.149080 -0.031427 -0.029542 3958.740 9.122960e+10 1.093030e+11 5.895810e+08 25.9614 25.7999 25.8532 25.9218 25.9145 25.3170 25.2468 25.2075
26 23003331013748 23003331006124 32002300010232 0 32 1.974910 -0.031631 -0.030205 3765.880 8.692630e+10 8.778700e+10 6.598210e+08 25.6942 25.5582 25.6234 25.7003 25.7103 25.1631 25.0820 25.0304
27 23004504000031 23004504000031 32002300003286 0 32 1.987380 -0.031396 -0.029872 3778.890 4.595910e+11 3.941810e+11 1.289410e+10 25.9597 25.2405 25.2330 25.0445 24.5721 23.5079 22.9810 22.6258
28 23003334023784 23003334009885 32002300007680 0 32 2.029520 -0.031415 -0.030000 3827.530 1.239350e+11 1.316810e+11 1.416000e+09 26.0890 25.5263 25.5436 25.4559 25.2036 24.5548 24.3543 24.1896
29 22004259000031 22004259000031 32002200004127 0 32 2.091940 -0.031419 -0.029809 3896.620 3.012300e+11 3.261890e+11 6.895720e+09 25.7880 25.1205 25.0227 24.9207 24.6135 23.6237 23.2907 23.0669
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
287792 54000085004173 54000085000010 53005300001690 2 53 0.265379 -0.023552 -0.020424 754.168 4.081150e+13 0.000000e+00 3.731290e+08 26.0746 24.9619 23.9155 23.5462 23.3109 22.9955 22.9507 22.8749
287793 54000085009403 54000085004389 53005300001697 1 53 0.265488 -0.023087 -0.020302 753.888 4.081150e+13 8.434440e+10 4.602520e+09 23.7237 22.4532 21.2601 20.8569 20.5986 20.1582 20.0627 19.9173
287794 54000085003876 54000085000010 53005300001690 2 53 0.264504 -0.023093 -0.020336 753.636 4.081150e+13 0.000000e+00 1.182640e+09 26.8623 24.9796 23.5520 23.0434 22.7194 22.1560 22.0467 21.8781
287795 54000085003919 54000085000010 53005300001690 2 53 0.265127 -0.022943 -0.019665 753.750 4.081150e+13 0.000000e+00 1.104930e+08 27.7590 26.4028 25.3185 24.9215 24.6756 24.3886 24.3598 24.3060
287796 54000085004438 54000085000010 53005300001690 2 53 0.263735 -0.022935 -0.019660 753.875 4.081150e+13 0.000000e+00 6.067410e+07 27.1930 26.3225 25.4770 25.1771 24.9857 24.7882 24.8050 24.7905
287797 54000085009247 54000085004172 53005300001702 1 53 0.265274 -0.023159 -0.020644 753.667 4.081150e+13 6.454930e+10 2.269620e+09 25.8597 23.9850 22.5747 22.0840 21.7793 21.2488 21.1408 20.9645
287798 54000085006332 54000085000010 53005300001690 2 53 0.264283 -0.023170 -0.020092 752.960 4.081150e+13 0.000000e+00 1.104430e+08 27.8233 26.3958 25.3016 24.8991 24.6511 24.3648 24.3376 24.2853
287799 457000346000365 457000346000207 53045700011743 0 53 0.286548 -0.023335 -0.020531 812.607 7.573780e+10 8.090180e+10 8.023770e+08 24.4522 23.7478 22.8918 22.6151 22.4281 22.1248 22.0862 21.9950
287800 457000346000440 457000346000248 53045700011403 2 53 0.287163 -0.023350 -0.020257 812.410 6.713120e+10 0.000000e+00 5.568940e+07 27.7170 26.7733 25.8946 25.5826 25.3826 25.1708 25.1847 25.1609
287801 457004670000010 457004670000010 53045700011589 0 53 0.276481 -0.023075 -0.018786 790.167 7.745910e+10 1.006970e+11 1.244510e+09 25.0292 24.0311 22.9056 22.5074 22.2434 21.7842 21.6885 21.5417
287802 54000085007551 54000085003741 53005300001693 1 53 0.264856 -0.022983 -0.019044 753.997 4.081150e+13 1.213530e+11 3.205690e+09 24.2812 22.7480 21.4984 21.0838 20.8244 20.3947 20.2950 20.1564
287803 54000085004646 54000085000010 53005300001690 2 53 0.264549 -0.023552 -0.019489 753.792 4.081150e+13 0.000000e+00 3.568230e+09 25.7899 23.8095 22.3399 21.8196 21.5043 20.9397 20.8198 20.6335
287804 54000085005972 54000085000010 53005300001690 2 53 0.264315 -0.023070 -0.019126 753.709 4.081150e+13 0.000000e+00 5.675970e+08 26.1596 24.8363 23.6788 23.2693 23.0098 22.6399 22.5778 22.4749
287805 54000085007008 54000085002627 53005300001712 1 53 0.263421 -0.023079 -0.019136 753.869 4.081150e+13 2.495910e+10 8.519830e+09 25.5278 23.5494 21.9491 21.3733 20.9948 20.1723 19.9616 19.7322
287806 54000085008022 54000085003989 53005300001700 1 53 0.266278 -0.022937 -0.019344 753.638 4.081150e+13 6.971320e+10 1.143980e+10 24.6611 22.7485 21.2091 20.6616 20.3249 19.6327 19.4495 19.2376
287807 54000085006972 54000085002598 53005300022819 0 53 0.264629 -0.023847 -0.020025 753.959 3.184430e+10 3.442630e+10 1.372670e+08 26.1674 25.3926 24.5591 24.2892 24.1078 23.8735 23.8671 23.8258
287808 54000085006694 54000085000010 53005300001690 2 53 0.264358 -0.022953 -0.019958 754.100 4.081150e+13 0.000000e+00 3.424900e+08 27.1035 25.4949 24.2675 23.8367 23.5692 23.1995 23.1352 23.0315
287809 54000085005932 54000085000010 53005300001690 2 53 0.265545 -0.023123 -0.020075 753.456 4.081150e+13 0.000000e+00 3.730640e+08 27.2614 25.6441 24.3750 23.9243 23.6397 23.2231 23.1522 23.0356
287810 54000085006399 54000085000010 53005300001690 2 53 0.265313 -0.023650 -0.019863 753.761 4.081150e+13 0.000000e+00 4.678540e+07 28.6472 27.2249 26.1719 25.7774 25.5368 25.2861 25.2774 25.2438
287811 54000085006523 54000085000010 53005300001690 2 53 0.264455 -0.023143 -0.019876 753.572 4.081150e+13 0.000000e+00 3.283740e+07 27.7890 26.9215 26.0849 25.7885 25.5999 25.4093 25.4294 25.4177
287812 54000085005895 54000085000010 53005300001690 2 53 0.263553 -0.023167 -0.020071 753.613 4.081150e+13 0.000000e+00 1.809330e+08 27.7774 26.2038 25.0119 24.5788 24.3088 23.9591 23.9107 23.8281
287813 54000085006430 54000085000010 53005300001690 2 53 0.265252 -0.023170 -0.020068 753.612 4.081150e+13 0.000000e+00 3.765210e+09 26.2415 24.2583 22.7058 22.1522 21.7854 21.0274 20.8486 20.6302
287814 54000085003305 54000085000010 53005300001690 2 53 0.263759 -0.023483 -0.019933 754.285 4.081150e+13 0.000000e+00 1.663730e+09 25.7149 23.9564 22.6682 22.2126 21.9327 21.5217 21.4477 21.3212
287815 54000085008629 54000085002856 53005300001692 1 53 0.265385 -0.022909 -0.019994 754.075 4.081150e+13 2.590580e+11 2.478550e+10 23.3025 21.9629 20.5587 20.0711 19.7165 18.8286 18.5747 18.3333
287816 457009347000010 457009347000010 53045700019544 0 53 0.296568 -0.023214 -0.019661 835.886 2.409840e+10 4.303290e+10 3.848050e+07 28.2130 27.0106 26.0863 25.7454 25.5397 25.3548 25.3689 25.3581
287817 54000085010042 54000085005212 53005300001713 1 53 0.265297 -0.023519 -0.019925 754.508 4.081150e+13 2.237710e+10 8.211960e+06 24.6312 24.6693 24.6080 24.7590 24.8474 25.0462 25.3615 25.5697
287818 54000085005428 54000085000010 53005300001690 2 53 0.265640 -0.023061 -0.020333 753.557 4.081150e+13 0.000000e+00 1.696080e+08 28.2060 26.4894 25.2284 24.7757 24.4936 24.1068 24.0501 23.9438
287819 54000085006736 54000085000010 53005300001690 2 53 0.263647 -0.022992 -0.020081 753.910 4.081150e+13 0.000000e+00 9.158140e+09 24.4597 22.5756 21.1567 20.6581 20.3506 19.8168 19.7106 19.5351
287820 54000085009704 54000085004656 53005300001709 1 53 0.262917 -0.023282 -0.020388 753.552 4.081150e+13 3.098370e+10 2.649590e+08 26.6341 25.2964 24.2185 23.8350 23.5980 23.3134 23.2757 23.2177
287821 54000085002923 54000085000010 53005300001690 2 53 0.264561 -0.023019 -0.020047 754.000 4.081150e+13 0.000000e+00 6.758000e+09 23.8866 22.5520 21.2798 20.8278 20.5386 20.0222 19.9098 19.7401

287822 rows × 20 columns


In [111]:
d = pd.read_csv('/Users/user/Desktop/output.csv')

In [112]:
d


Out[112]:
Unnamed: 0 37000046000024 37014522009504 456000089002589 65000614000013 82020319006828 log-likelihood
0 0 1.734220e+12 1.662620e+13 5.545210e+12 1.214470e+13 3.414230e+12 -0.346273
1 1 1.999310e+12 2.239430e+12 3.510880e+13 2.235130e+12 2.379720e+12 -0.322402
2 2 2.366810e+12 2.495040e+13 3.309230e+12 4.417750e+12 8.873370e+11 -0.342119
3 3 1.562090e+12 1.082790e+13 1.000080e+12 9.647970e+12 1.151820e+13 -0.347341
4 4 4.222380e+12 1.774670e+12 8.224440e+12 2.154220e+12 1.779580e+13 -0.310862
5 5 2.863410e+12 8.348370e+12 1.746270e+12 2.657710e+12 2.889230e+12 -0.338001
6 6 2.143900e+12 6.251810e+12 9.777060e+11 1.187710e+12 1.183490e+13 -0.331005
7 7 7.900830e+12 7.739890e+12 1.855580e+12 3.579470e+12 5.261200e+12 -0.337661
8 8 1.230740e+12 8.480050e+12 5.666570e+12 2.214300e+13 2.009630e+12 -0.359540
9 9 7.891360e+12 1.000080e+12 9.540380e+12 2.777340e+13 1.090450e+12 -0.363210
10 10 1.775110e+13 2.026850e+12 9.678090e+12 1.215160e+13 4.502100e+12 -0.338870
11 11 3.309230e+12 4.260250e+12 7.409830e+13 1.302170e+12 1.047590e+13 -0.349466
12 12 8.985260e+11 8.107390e+12 9.531780e+12 2.071600e+12 7.315580e+11 -0.322878
13 13 6.523780e+11 8.159030e+11 1.222820e+13 1.601680e+12 6.092590e+12 -0.313117
14 14 4.920380e+12 1.319390e+13 9.153090e+12 2.068330e+13 1.310950e+13 -0.349005
15 15 5.690660e+12 8.592800e+12 1.949390e+12 2.418790e+13 1.093030e+11 -0.372727
16 16 1.213530e+12 1.017300e+13 1.522500e+12 4.489190e+12 1.062050e+12 -0.344490
17 17 1.117650e+13 1.530250e+12 3.750050e+13 1.833200e+11 5.573610e+12 -0.317041
18 18 6.713120e+10 7.670180e+12 8.964600e+12 4.755990e+12 2.925370e+12 -0.326864
19 19 6.337880e+12 2.235130e+12 1.274630e+12 1.517340e+12 1.445470e+13 -0.327558

In [ ]: