In [3]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

In [6]:
path = "/Users/adrian/Downloads/R_12675046899699465570"
d = np.loadtxt(path, usecols=[0,1,2])

# trim out points with big uncertainties
d = d[d[:,2] < 0.15]

In [12]:
hjd0 = 55763.907589
period = 0.6043698

phase = ((d[:,0] - hjd0) % period) / period
plt.errorbar(phase, d[:,1], d[:,2], linestyle='none', 
             ecolor='#666666', marker='o', alpha=0.75)
plt.ylim(17.65, 16.75)


Out[12]:
(17.65, 16.75)

In [ ]: