n = 20000
ACSSub = ACSdata.sample(n)
ACSmag = ACSSub["Vvega"]
ACScolor = ACSSub["Vvega"] - ACSSub["Ivega"]

HUGSSub = HUGSdata[1].sample(n)
HUGSmag = HUGSSub["F606W"]
HUGScolor = HUGSSub["F606W"] - HUGSSub["F814W"]