In [42]:

    

#import necessary python libraries
import json
import matplotlib.pyplot as plt
import numpy as np
from scipy.optimize import curve_fit, minimize
import pandas as pd
import math
from JSON_to_DF import JSON_to_DataFrame
from Lightcurve_class import *
import celerite
import pickle
import seaborn as sns
%matplotlib notebook


    

In [56]:

    

#Create Kernels for Gaussian Process

#Real term parameter initialization
a = 1e-4
c = 1
#Matern term parameter initialization
sig = 1e-2
rho = 100

#Bounds on parameters 
bounds = dict(log_a = (-15,15), log_c = (-15,15))
bounds = dict(log_sigma = (-15, 15), log_rho = (-15, 15))

#Create Kernels
Real_Kernel = celerite.terms.RealTerm(log_a = np.log(a), log_c = np.log(c), bounds=bounds)
Matern_Kernel = celerite.terms.Matern32Term(log_sigma = np.log(sig), log_rho = np.log(rho))


    

In [57]:

    

poly_4 = []
poly_6 = []
poly_8 = []
GP_real = []
GP_matern = []
Bazin  = []
Kapernka = []
directory = "../../../OSC_data/pickled_data/"
for root, dirs, files in os.walk(directory):
    for file in files:
        file_name, ext = os.path.splitext(os.path.basename(file))
     
        if(ext == ".pickle"):
            SN = deserialize(directory + file)
           
            for key in SN.Lightcurves.keys():
                if(SN.Lightcurves[key].n_good_obs <= 6):
                    continue
                for m_fit in SN.Lightcurves[key].medians.keys():
                    if(m_fit[0:4] == 'poly'):
                        degree = m_fit[-1]
                        R_fit = 'polynomial_' + str(degree)
                    else:
                        R_fit = m_fit
                    tup = (SN.Lightcurves[key].Rchi2[R_fit], SN.Lightcurves[key].medians[m_fit])
                    if(np.isnan(tup[0]) or np.isnan(tup[1])):
                        print(file_name)
                        print(key)
                        print(m_fit)
                    if(m_fit == "poly_4"):
                        poly_4.append(tup)
                    elif(m_fit == "poly_6"):
                        poly_6.append(tup)
                    elif(m_fit == "poly_8"):
                        poly_8.append(tup)
                    elif(m_fit == "GP_Real"):
                        GP_real.append(tup)
                    elif(m_fit == "GP_Matern"):
                        GP_matern.append(tup)
                    elif(m_fit == "Bazin"):
                        Bazin.append(tup)
                    elif(m_fit == "Kapernka"):
                        Kapernka.append(tup)


    

    




ASASSN-14az
C_5
GP_Real
ASASSN-14az
C_5
GP_Matern
ASASSN-14az
g_7,10
GP_Real
ASASSN-14az
g_7,10
GP_Matern
ASASSN-14bb
i_7,8
GP_Matern
ASASSN-14bb
g_7,8
GP_Matern
ASASSN-14bb
I_7,8
GP_Matern
ASASSN-14bd
W1_5
GP_Real
ASASSN-14bd
W1_5
GP_Matern
ASASSN-14bd
U_5
GP_Real
ASASSN-14bd
U_5
GP_Matern
ASASSN-14bd
B_5
GP_Real
ASASSN-14bd
B_5
GP_Matern
ASASSN-14bd
W2_5
GP_Matern
ASASSN-14bd
M2_5
GP_Matern
ASASSN-14bd
V_5
GP_Real
ASASSN-14bd
V_5
GP_Matern
ASASSN-14dc
W1_5
GP_Real
ASASSN-14dc
W1_5
GP_Matern
ASASSN-14dc
U_5
GP_Real
ASASSN-14dc
U_5
GP_Matern
ASASSN-14dc
B_5
GP_Real
ASASSN-14dc
B_5
GP_Matern
ASASSN-14dc
W2_5
GP_Real
ASASSN-14dc
W2_5
GP_Matern
ASASSN-14dc
V_5
GP_Real
ASASSN-14dc
V_5
GP_Matern
ASASSN-14dc
M2_5
GP_Real
ASASSN-14dc
M2_5
GP_Matern
ASASSN-14fa
W1_4
GP_Real
ASASSN-14fa
U_4
GP_Real
ASASSN-14fa
W2_4
GP_Real
ASASSN-14fa
V_4
GP_Real
ASASSN-14my
w_2,3,5
GP_Matern
CSS111118:074033+240900
C_3
GP_Matern
CSS130315:131135-042408
C_5
GP_Real
CSS130912:025702-001844
R_1
GP_Real
CSS130912:025702-001844
r_1
GP_Real
CSS130912:025702-001844
g_1
GP_Real
CSS130912:025702-001844
i_1
GP_Real
CSS130925:010249-004431
C_7
GP_Real
CSS130925:010249-004431
C_7
GP_Matern
CSS130925:010249-004431
R_1
GP_Real
CSS130925:010249-004431
r_1
GP_Real
CSS131126:032030+211829
C_3
GP_Real
CSS131126:032030+211829
C_3
GP_Matern
CSS131213:112237+122931
C_3
GP_Real
CSS131213:112237+122931
C_3
GP_Matern
CSS140218:095739+123318
C_4
GP_Real
CSS140218:095739+123318
C_4
GP_Matern
CSS140227:114929+083432
C_2
GP_Real
CSS140227:114929+083432
C_2
GP_Matern
CSS140306:094630+061119
C_2
GP_Real
CSS140306:094630+061119
C_2
GP_Matern
CSS140326:122257+282955
R_1
GP_Real
CSS140326:122257+282955
R_1
GP_Matern
CSS140326:122257+282955
r_1
GP_Real
CSS140326:122257+282955
r_1
GP_Matern
CSS140914:010107-101840
W1_3
GP_Real
CSS140914:010107-101840
W1_3
GP_Matern
CSS140914:010107-101840
U_3
GP_Real
CSS140914:010107-101840
U_3
GP_Matern
CSS140914:010107-101840
B_3
GP_Real
CSS140914:010107-101840
B_3
GP_Matern
CSS140914:010107-101840
V_3
GP_Real
CSS140914:010107-101840
V_3
GP_Matern
CSS140925:011731+054106
C_4
GP_Real
CSS140925:011731+054106
C_4
GP_Matern
CSS141025:075616+301945
w_3
GP_Real
CSS141025:075616+301945
w_3
GP_Matern
CSS141117:030348+002415
w_5
GP_Real
Gaia14aag
g_4,5
GP_Real
Gaia14aag
g_4,5
GP_Matern
Gaia14aag
r_4,5
GP_Real
Gaia14aag
r_4,5
GP_Matern
Gaia14aan
B_4,5
GP_Real
Gaia14aan
B_4,5
GP_Matern
Gaia14aan
V_4,5
GP_Real
Gaia14aan
V_4,5
GP_Matern
Gaia14aan
r_4,5
GP_Real
Gaia14aan
r_4,5
GP_Matern
Gaia14aan
i_4,5
GP_Real
Gaia14aan
i_4,5
GP_Matern
Gaia14abh
r_4,5
GP_Real
Gaia14abh
r_4,5
GP_Matern
Gaia14abh
i_4,5
GP_Real
Gaia14abh
i_4,5
GP_Matern
Gaia14abh
V_4,5
GP_Real
Gaia14abh
V_4,5
GP_Matern
Gaia14acf
B_3,4
GP_Real
Gaia14acf
B_3,4
GP_Matern
Gaia14acf
r_3,4
GP_Real
Gaia14acf
r_3,4
GP_Matern
Gaia14acf
g_3,4
GP_Real
Gaia14acf
g_3,4
GP_Matern
Gaia14acq
V_4,5
GP_Real
Gaia14acq
V_4,5
GP_Matern
Gaia14acq
r_4,5
GP_Real
Gaia14acq
r_4,5
GP_Matern
Gaia14acq
i_4,5
GP_Real
Gaia14acq
i_4,5
GP_Matern
Gaia14ada
g_4,5
GP_Real
Gaia14ada
g_4,5
GP_Matern
Gaia14ada
r_4,5
GP_Real
Gaia14ada
r_4,5
GP_Matern
Gaia14adl
r_3,4
GP_Real
Gaia14adl
r_3,4
GP_Matern
Gaia14adm
r_4,5
GP_Real
Gaia14adm
r_4,5
GP_Matern
Gaia14adm
g_4,5
GP_Real
Gaia14adm
g_4,5
GP_Matern
Gaia14adm
i_4,5
GP_Real
Gaia14adm
i_4,5
GP_Matern
Gaia14adm
u_4,5
GP_Real
Gaia14adm
u_4,5
GP_Matern
Gaia14ado
V_4,5
GP_Real
Gaia14ado
V_4,5
GP_Matern
Gaia14ado
r_4,5
GP_Real
Gaia14ado
r_4,5
GP_Matern
Gaia14ado
g_4,5
GP_Real
Gaia14ado
g_4,5
GP_Matern
Gaia14ado
i_4,5
GP_Real
Gaia14ado
i_4,5
GP_Matern
Gaia14adq
g_5,6
GP_Real
Gaia14adq
g_5,6
GP_Matern
Gaia14adq
r_5,6
GP_Real
Gaia14adq
r_5,6
GP_Matern
Gaia14adq
i_5,6
GP_Real
Gaia14adq
i_5,6
GP_Matern
Gaia14adr
V_4,5
GP_Real
Gaia14adr
V_4,5
GP_Matern
Gaia14adr
B_4,5
GP_Real
Gaia14adr
B_4,5
GP_Matern
Gaia14adr
r_4,5
GP_Real
Gaia14adr
r_4,5
GP_Matern
Gaia14adr
i_4,5
GP_Real
Gaia14adr
i_4,5
GP_Matern
Gaia14adz
g_4,5
GP_Real
Gaia14adz
g_4,5
GP_Matern
Gaia14adz
r_4,5
GP_Real
Gaia14adz
r_4,5
GP_Matern
Gaia14adz
u_4,5
GP_Real
Gaia14adz
u_4,5
GP_Matern
Gaia15aab
r_4,5
GP_Real
Gaia15aab
r_4,5
GP_Matern
Gaia15aab
B_4,5
GP_Real
Gaia15aab
B_4,5
GP_Matern
Gaia15aab
V_4,5
GP_Real
Gaia15aab
V_4,5
GP_Matern
Gaia15aab
i_4,5
GP_Real
Gaia15aab
i_4,5
GP_Matern
Gaia15abh
r_3,4
GP_Matern
Gaia15abk
g_3,4
GP_Real
Gaia15abk
g_3,4
GP_Matern
iPTF13ajg
R_2
GP_Real
iPTF13ajg
R_2
GP_Matern
iPTF13ajg
i_2
GP_Real
iPTF13ajg
i_2
GP_Matern
iPTF13ajg
r_2
GP_Real
iPTF13ajg
r_2
GP_Matern
iPTF13ajg
B_2
GP_Real
iPTF13ajg
B_2
GP_Matern
iPTF13ajg
g_2
GP_Real
iPTF13ajg
g_2
GP_Matern
iPTF13bjx
R_1
GP_Real
iPTF13bjx
R_1
GP_Matern
iPTF13bjx
r_1
GP_Real
iPTF13bjx
r_1
GP_Matern
iPTF13bld
R_1
GP_Real
iPTF13bld
R_1
GP_Matern
iPTF13bsg
R_1
GP_Real
iPTF13bsg
R_1
GP_Matern
iPTF13bsg
r_1
GP_Real
iPTF13bsg
r_1
GP_Matern
LSQ12btn
g_2
GP_Real
LSQ12btn
r_2
GP_Real
LSQ12btn
i_2
GP_Real
LSQ12btw
G_2
GP_Real
LSQ12btw
V_2
GP_Real
LSQ12ca
g_2
GP_Real
LSQ12ca
r_2
GP_Real
LSQ12ca
i_2
GP_Real
LSQ12ca
B_2
GP_Real
LSQ12ca
V_2
GP_Real
LSQ12cdl
g_2
GP_Real
LSQ12cdl
r_2
GP_Real
LSQ12cdl
i_2
GP_Real
LSQ12cdl
V_2
GP_Real
LSQ12cdl
B_2
GP_Real
LSQ12dlf
V_3
GP_Real
LSQ12dlf
B_3
GP_Real
LSQ12dlf
I_3
GP_Real
LSQ12dlf
R_3
GP_Real
LSQ12fuk
B_2
GP_Real
LSQ12fuk
V_2
GP_Real
LSQ12fuk
r_2
GP_Real
LSQ12fuk
i_2
GP_Real
LSQ12fxd
V_2
GP_Real
LSQ12fxd
B_2
GP_Real
LSQ12fxd
i_2
GP_Real
LSQ12fxd
u_2
GP_Real
LSQ12fxd
r_2
GP_Real
LSQ12fxd
G_3
GP_Real
LSQ12fxd
I_3
GP_Real
LSQ12fxd
R_3
GP_Real
LSQ12gdj
i_3
GP_Real
LSQ12gdj
r_3
GP_Real
LSQ12gdj
V_3
GP_Real
LSQ12gdj
B_3
GP_Real
LSQ12gdj
B_4
GP_Real
LSQ12gdj
I_4
GP_Real
LSQ12gdj
R_4
GP_Real
LSQ12gdj
V_4
GP_Real
LSQ12gdj
W1_6
GP_Real
LSQ12gdj
M2_4
GP_Real
LSQ12gdj
U_4
GP_Real
LSQ12gdj
W1_4
GP_Real
LSQ12gdj
W2_4
GP_Real
LSQ12gdj
U_6
GP_Real
LSQ12gdj
B_6
GP_Real
LSQ12gdj
W2_6
GP_Real
LSQ12gdj
V_6
GP_Real
LSQ12gdj
M2_6
GP_Real
LSQ12gdj
G_4
GP_Real
LSQ12gdj
Z_4
GP_Real
LSQ12gef
V_2
GP_Real
LSQ12gef
B_2
GP_Real
LSQ12gef
r_2
GP_Real
LSQ12gef
i_2
GP_Real
LSQ12gln
i_2
GP_Real
LSQ12gln
r_2
GP_Real
LSQ12gln
V_2
GP_Real
LSQ12gln
B_2
GP_Real
LSQ12gpw
B_2
GP_Real
LSQ12gpw
V_2
GP_Real
LSQ12gpw
r_2
GP_Real
LSQ12gxj
r_2
GP_Real
LSQ12gxj
i_2
GP_Real
LSQ12gxj
V_2
GP_Real
LSQ12gxj
B_2
GP_Real
LSQ12gyc
r_2
GP_Real
LSQ12gyc
i_2
GP_Real
LSQ12gyc
V_2
GP_Real
LSQ12gyc
B_2
GP_Real
LSQ12gzm
r_2
GP_Real
LSQ12gzm
V_2
GP_Real
LSQ12gzm
B_2
GP_Real
LSQ12hjm
r_2
GP_Real
LSQ12hjm
i_2
GP_Real
LSQ12hjm
V_2
GP_Real
LSQ12hjm
B_2
GP_Real
LSQ12hno
V_2
GP_Real
LSQ12hno
B_2
GP_Real
LSQ12hno
r_2
GP_Real
LSQ12hno
i_2
GP_Real
LSQ12hvj
B_2
GP_Real
LSQ12hvj
V_2
GP_Real
LSQ12hvj
r_2
GP_Real
LSQ12hvj
i_2
GP_Real
LSQ12hxx
i_2
GP_Real
LSQ12hxx
r_2
GP_Real
LSQ12hxx
V_2
GP_Real
LSQ12hxx
B_2
GP_Real
LSQ12hzj
r_2
GP_Real
LSQ12hzj
i_2
GP_Real
LSQ12hzj
V_2
GP_Real
LSQ12hzj
B_2
GP_Real
LSQ14ahb
C_4
GP_Real
LSQ14ahb
C_4
GP_Matern
LSQ14an
r_1
GP_Matern
LSQ14an
V_1
GP_Matern
LSQ14an
g_1
GP_Matern
LSQ14an
i_1
GP_Matern
LSQ14bcj
C_5
GP_Real
LSQ14bcj
C_5
GP_Matern
LSQ14wp
C_4
GP_Real
LSQ14wp
C_4
GP_Matern
MASTER OT J075929.72+194610.7
C_3
GP_Real
MASTER OT J075929.72+194610.7
C_3
GP_Matern
MASTER OT J140536.83-012210.6
C_3
GP_Real
MASTER OT J140536.83-012210.6
C_3
GP_Matern
MLS100215:150114-123607
C_3
GP_Real
MLS100215:150114-123607
C_3
GP_Matern
MLS100219:120858+034250
C_2
GP_Real
MLS100219:120858+034250
C_2
GP_Matern
MLS100313:131349-062327
C_3
GP_Real
MLS100313:131349-062327
C_3
GP_Matern
MLS100315:073816+195257
C_3
GP_Matern
MLS100315:114953-002513
C_3
GP_Real
MLS100315:114953-002513
C_3
GP_Matern
MLS100407:093835+142119
C_3
GP_Real
MLS100409:104506+081910
C_3
GP_Real
MLS100505:133811-093052
C_3
GP_Real
MLS100610:221229-092622
C_3
GP_Real
MLS100610:221229-092622
C_3
GP_Matern
MLS100613:134627-095129
C_3
GP_Real
MLS100613:134627-095129
C_3
GP_Matern
MLS100617:214000-143660
C_3
GP_Matern
MLS100619:221052-103113
C_3
GP_Real
MLS100619:221052-103113
C_3
GP_Matern
MLS100901:224523-070002
C_3
GP_Matern
MLS100902:004233+042524
C_3
GP_Real
MLS100902:004233+042524
C_3
GP_Matern
MLS100905:035122+181501
C_3
GP_Real
MLS100905:035122+181501
C_3
GP_Matern
MLS100911:033634+193727
C_3
GP_Matern
MLS100918:001802+144556
C_3
GP_Real
MLS100918:022749+162612
C_3
GP_Real
MLS100918:022749+162612
C_3
GP_Matern
MLS100929:004423-013814
C_3
GP_Matern
MLS100929:232031-043249
C_3
GP_Real
MLS100929:232031-043249
C_3
GP_Matern
MLS101009:004045+040316
C_2
GP_Real
MLS101009:004045+040316
C_2
GP_Matern
MLS101009:004857+053147
C_3
GP_Real
MLS101009:004857+053147
C_3
GP_Matern
MLS101011:021500+192736
C_3
GP_Real
MLS101114:033154+290505
C_2
GP_Real
MLS101114:033154+290505
C_2
GP_Matern
MLS101127:013210+081548
C_2
GP_Real
MLS101127:013210+081548
C_2
GP_Matern
MLS101201:041720+252926
C_3
GP_Real
MLS101201:041720+252926
C_3
GP_Matern
MLS101205:031059+235346
C_3
GP_Real
MLS101205:031059+235346
C_3
GP_Matern
MLS101205:084407+180717
C_3
GP_Real
MLS101205:084407+180717
C_3
GP_Matern
MLS101208:085152+155343
C_3
GP_Real
MLS101208:085152+155343
C_3
GP_Matern
MLS101212:005238+053040
C_3
GP_Real
MLS101212:005238+053040
C_3
GP_Matern
MLS110102:071845+174038
C_3
GP_Real
MLS110102:071845+174038
C_3
GP_Matern
MLS110104:075613+271148
C_3
GP_Real
MLS110104:075613+271148
C_3
GP_Matern
MLS110108:105331+010911
C_3
GP_Real
MLS110108:105331+010911
C_3
GP_Matern
MLS110110:105658+111548
C_3
GP_Real
MLS110110:105658+111548
C_3
GP_Matern
MLS110113:233254+005651
C_3
GP_Real
MLS110113:233254+005651
C_3
GP_Matern
MLS110114:074706+233205
C_2
GP_Real
MLS110114:074706+233205
C_2
GP_Matern
MLS110114:080505+292140
C_3
GP_Real
MLS110114:080505+292140
C_3
GP_Matern
MLS110128:091122+131254
C_2
GP_Real
MLS110128:091122+131254
C_2
GP_Matern
MLS110129:035607+192755
C_3
GP_Real
MLS110129:035607+192755
C_3
GP_Matern
MLS110208:115629-003653
C_3
GP_Real
MLS110208:115629-003653
C_3
GP_Matern
MLS110208:120417-012415
C_2
GP_Real
MLS110208:120417-012415
C_2
GP_Matern
MLS110210:101045+074134
C_3
GP_Matern
MLS110211:091006+131404
C_3
GP_Real
MLS110211:091006+131404
C_3
GP_Matern
MLS110212:084550+170838
C_3
GP_Matern
MLS110213:014203+112145
C_3
GP_Real
MLS110213:014203+112145
C_3
GP_Matern
MLS110225:093527+131432
C_3
GP_Matern
MLS110225:094836+114017
C_3
GP_Matern
MLS110225:104818+085623
C_3
GP_Matern
MLS110226:110325+070631
C_3
GP_Real
MLS110226:110325+070631
C_3
GP_Matern
MLS110301:092035+111205
C_3
GP_Real
MLS110301:092035+111205
C_3
GP_Matern
MLS110301:144813-120124
C_3
GP_Matern
MLS110302:091540+095610
C_3
GP_Real
MLS110302:091540+095610
C_3
GP_Matern
MLS110302:094259+074050
C_3
GP_Matern
MLS110304:105205+030524
C_3
GP_Real
MLS110304:105205+030524
C_3
GP_Matern
MLS110305:101539+104656
C_3
GP_Real
MLS110305:101539+104656
C_3
GP_Matern
MLS110306:100601+072628
C_3
GP_Real
MLS110306:100601+072628
C_3
GP_Matern
MLS110309:112638+134509
C_3
GP_Real
MLS110309:112638+134509
C_3
GP_Matern
MLS110310:112922+295114
C_3
GP_Real
MLS110310:112922+295114
C_3
GP_Matern
MLS110312:110709+154832
C_3
GP_Matern
MLS110314:112615+105157
C_3
GP_Matern
MLS110328:115422+001547
C_3
GP_Real
MLS110328:115422+001547
C_3
GP_Matern
MLS110329:090700+152330
C_3
GP_Real
MLS110329:090700+152330
C_3
GP_Matern
MLS110329:115737+021302
C_3
GP_Matern
MLS110330:124640-085116
C_3
GP_Real
MLS110330:124640-085116
C_3
GP_Matern
MLS110330:125334-053143
C_3
GP_Matern
MLS110404:110052+101203
C_3
GP_Real
MLS110404:110052+101203
C_3
GP_Matern
MLS110406:122742-034057
C_3
GP_Real
MLS110406:122742-034057
C_3
GP_Matern
MLS110429:114433-001908
C_3
GP_Real
MLS110429:114433-001908
C_3
GP_Matern
MLS110430:132739-075944
C_3
GP_Real
MLS110430:132739-075944
C_3
GP_Matern
MLS110430:132800-102852
C_2
GP_Matern
MLS110503:134315-070246
C_2
GP_Matern
MLS110512:151647-072202
C_2
GP_Real
MLS110512:151647-072202
C_2
GP_Matern
MLS110525:134153+011557
C_3
GP_Real
MLS110525:134153+011557
C_3
GP_Matern
MLS110525:134316+004749
C_3
GP_Real
MLS110525:134316+004749
C_3
GP_Matern
MLS110526:093516+131102
C_3
GP_Matern
MLS110526:094440+135936
C_3
GP_Matern
MLS110526:104421+082012
C_3
GP_Matern
MLS110603:135215-001421
C_3
GP_Real
MLS110603:135215-001421
C_3
GP_Matern
MLS110604:095908+115727
C_3
GP_Real
MLS110604:095908+115727
C_3
GP_Matern
MLS110605:105619+104444
C_3
GP_Real
MLS110605:105619+104444
C_3
GP_Matern
MLS110607:110336+093819
C_3
GP_Real
MLS110607:110336+093819
C_3
GP_Matern
MLS110607:140351-113108
C_3
GP_Real
MLS110607:140351-113108
C_3
GP_Matern
MLS110623:220331-095729
C_3
GP_Real
MLS110623:220331-095729
C_3
GP_Matern
MLS110624:125414-014544
C_3
GP_Matern
MLS110626:164208-002929
C_3
GP_Real
MLS110626:164208-002929
C_3
GP_Matern
MLS110918:233252-015838
C_2
GP_Matern
MLS110919:235442+025617
C_3
GP_Real
MLS110919:235442+025617
C_3
GP_Matern
MLS110921:231509-073149
C_3
GP_Real
MLS110921:231509-073149
C_3
GP_Matern
MLS110922:205621-144117
C_3
GP_Real
MLS110922:205621-144117
C_3
GP_Matern
MLS110923:040547+191222
C_2
GP_Matern
MLS110923:230624-074618
C_3
GP_Real
MLS110923:230624-074618
C_3
GP_Matern
MLS110928:003522+010256
C_3
GP_Real
MLS110928:003522+010256
C_3
GP_Matern
MLS110928:003544+013947
C_3
GP_Real
MLS110928:003544+013947
C_3
GP_Matern
MLS110928:013335+091325
C_3
GP_Real
MLS110928:013335+091325
C_3
GP_Matern
MLS110928:014044+101718
C_3
GP_Real
MLS110928:014044+101718
C_3
GP_Matern
MLS110928:025340+200040
C_3
GP_Real
MLS110928:025340+200040
C_3
GP_Matern
MLS110929:032550+210325
C_3
GP_Matern
MLS111001:001444+065819
C_3
GP_Real
MLS111001:001444+065819
C_3
GP_Matern
MLS111001:002408+073155
C_2
GP_Real
MLS111001:002408+073155
C_2
GP_Matern
MLS111001:235057-001403
C_3
GP_Matern
MLS111003:001735+015233
C_2
GP_Real
MLS111003:001735+015233
C_2
GP_Matern
MLS111003:002437-012122
C_3
GP_Real
MLS111003:002437-012122
C_3
GP_Matern
MLS111018:035322+190349
C_3
GP_Real
MLS111018:035322+190349
C_3
GP_Matern
MLS111018:230236-053124
C_3
GP_Real
MLS111018:230236-053124
C_3
GP_Matern
MLS111018:230738-051313
C_3
GP_Real
MLS111018:230738-051313
C_3
GP_Matern
MLS111019:014622+145108
C_3
GP_Real
MLS111019:014622+145108
C_3
GP_Matern
MLS111019:230046-081411
C_3
GP_Real
MLS111019:230046-081411
C_3
GP_Matern
MLS111019:234959+002444
C_3
GP_Real
MLS111019:234959+002444
C_3
GP_Matern
MLS111020:020314+121721
C_3
GP_Real
MLS111020:020314+121721
C_3
GP_Matern
MLS111020:222136-082829
C_3
GP_Real
MLS111020:222136-082829
C_3
GP_Matern
MLS111022:081449+192315
C_3
GP_Real
MLS111022:081449+192315
C_3
GP_Matern
MLS111024:015106+091908
C_3
GP_Real
MLS111024:015106+091908
C_3
GP_Matern
MLS111024:081839+174654
C_3
GP_Real
MLS111024:081839+174654
C_3
GP_Matern
MLS111028:092534+134312
C_2
GP_Real
MLS111028:092534+134312
C_2
GP_Matern
MLS111029:230618-013254
C_3
GP_Real
MLS111029:230618-013254
C_3
GP_Matern
MLS111030:041520+222014
C_3
GP_Real
MLS111030:041520+222014
C_3
GP_Matern
MLS111031:015238+013231
C_3
GP_Real
MLS111031:015238+013231
C_3
GP_Matern
MLS111031:015916+084815
C_3
GP_Real
MLS111031:015916+084815
C_3
GP_Matern
MLS111031:020050+010720
C_3
GP_Real
MLS111031:020050+010720
C_3
GP_Matern
MLS111101:012707+003017
C_3
GP_Real
MLS111101:012707+003017
C_3
GP_Matern
MLS111102:105504+081547
C_3
GP_Real
MLS111102:105504+081547
C_3
GP_Matern
MLS111103:105634+041319
C_3
GP_Real
MLS111103:105634+041319
C_3
GP_Matern
MLS111116:004252+031300
C_3
GP_Real
MLS111116:004252+031300
C_3
GP_Matern
MLS111123:000235+001552
C_3
GP_Real
MLS111123:000235+001552
C_3
GP_Matern
MLS111126:071406+271059
C_3
GP_Real
MLS111126:071406+271059
C_3
GP_Matern
MLS111126:092534+134312
C_3
GP_Real
MLS111126:092534+134312
C_3
GP_Matern
MLS111127:022116+125551
C_3
GP_Real
MLS111127:022116+125551
C_3
GP_Matern
MLS111127:033342+131632
C_3
GP_Real
MLS111127:033342+131632
C_3
GP_Matern
MLS111128:100109+124221
C_3
GP_Real
MLS111128:100109+124221
C_3
GP_Matern
MLS111129:012905+114913
C_3
GP_Real
MLS111129:012905+114913
C_3
GP_Matern
MLS111225:113846+035506
C_3
GP_Real
MLS111225:113846+035506
C_3
GP_Matern
MLS111225:125106-000447
C_3
GP_Real
MLS111225:125106-000447
C_3
GP_Matern
MLS111227:071633+293258
C_3
GP_Real
MLS111227:071633+293258
C_3
GP_Matern
MLS111227:080448+173647
C_3
GP_Real
MLS111227:080448+173647
C_3
GP_Matern
MLS111227:124651-055601
C_3
GP_Real
MLS111227:124651-055601
C_3
GP_Matern
MLS111228:035911+152547
C_3
GP_Real
MLS111228:035911+152547
C_3
GP_Matern

In [58]:

    

#Create new supernovae object
SN = Supernovae('../../../OSC_data/JSON_data/sne-2010-2014-master/Gaia14adz.json')
#Load in lightcurve data
SN.load_LightCurves()


    

    




/Users/ethanmarx/anaconda/envs/NUREU17/lib/python3.6/site-packages/pandas/core/generic.py:1299: PerformanceWarning: 
your performance may suffer as PyTables will pickle object types that it cannot
map directly to c-types [inferred_type->mixed,key->block1_values] [items->['band', 'source', 'telescope']]

  return pytables.to_hdf(path_or_buf, key, self, **kwargs)

In [59]:

    

print(SN.Lightcurves['g_4,5'].Gaussian_process(Real_Kernel))
print(SN.Lightcurves['g_4,5'].Gaussian_process(Matern_Kernel))


    

    




Real
[ 0.96725273  0.96725273  0.96725273  0.96725273  0.96725388  0.96725296
  0.96725273  0.96725273  0.96725267  0.96725372  0.96725273  0.96725273
  0.96725273  0.96725273  0.96725273  0.96725273  0.96725273  0.96725273
  0.96725273  0.96725273  0.96725273  0.96725273  0.96726084  0.96718958
  0.96742634  0.96714646  0.96725568  0.96720949  0.96724471  0.9672425
  0.96723466  0.96726056  0.96724774  0.96725142  0.96725679  0.96719228
  0.96725977  0.96724126  0.96724098  0.96736279  0.96725435         nan
         nan  0.96725966  0.96726509  0.96726033  0.96724565  0.96724018
  0.96725225  0.96725273  0.96725273]
None
Matern






    




/Users/ethanmarx/anaconda/envs/NUREU17/lib/python3.6/site-packages/numpy/lib/function_base.py:3858: RuntimeWarning: Invalid value encountered in median
  r = func(a, **kwargs)






    




[ 0.96677612  0.96560035  0.96502048  0.96515965  0.96507628  0.96482924
  0.96474374         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan  0.96725273]
None

In [60]:

    

#Create new supernovae object
SN = Supernovae('../../../OSC_data/JSON_data/sne-2010-2014-master/LSQ12ca.json')
#Load in lightcurve data
SN.load_LightCurves()


    

    




/Users/ethanmarx/anaconda/envs/NUREU17/lib/python3.6/site-packages/pandas/core/generic.py:1299: PerformanceWarning: 
your performance may suffer as PyTables will pickle object types that it cannot
map directly to c-types [inferred_type->mixed,key->block1_values] [items->['band', 'source', 'telescope']]

  return pytables.to_hdf(path_or_buf, key, self, **kwargs)

In [61]:

    

print(SN.Lightcurves['B_2'].Gaussian_process(Real_Kernel))
print(SN.Lightcurves['B_2'].Gaussian_process(Matern_Kernel))


    

    




/Users/ethanmarx/anaconda/envs/NUREU17/lib/python3.6/site-packages/numpy/lib/function_base.py:3858: RuntimeWarning: Invalid value encountered in median
  r = func(a, **kwargs)






    




Real
[ 0.62395343         nan         nan         nan         nan         nan
         nan         nan         nan         nan  0.62395343]
None
Matern
[ 0.62395343         nan         nan         nan         nan         nan
         nan         nan         nan         nan  0.62395343]
None

In [66]:

    

#Create new supernovae object
SN = Supernovae('../../../OSC_data/JSON_data/sne-2010-2014-master/iPTF13bjx.json')
#Load in lightcurve data
SN.load_LightCurves()


    

In [68]:

    

SN.Lightcurves['R_1'].Gaussian_process(Real_Kernel)
SN.Lightcurves['R_1'].Gaussian_process(Matern_Kernel)


    

    




Real
[ 0.84455021         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan  0.84455021]
Matern
[ 0.84455021         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan         nan
         nan         nan         nan         nan         nan  0.84455021]






    




/Users/ethanmarx/anaconda/envs/NUREU17/lib/python3.6/site-packages/numpy/lib/function_base.py:3858: RuntimeWarning: Invalid value encountered in median
  r = func(a, **kwargs)

In [ ]: