Methods, systems, and computer readable media for predicting early onset glaucoma via structural parameters obtained via time and spectral-domain optical coherence tomography (SD-OCT) are disclosed. SD-OCT technology measures the thicknesses of the layers of the retina that contain the cells that die in glaucoma before there is any functional loss. A multivariable predictive model may be used for the diagnosis of early glaucoma using a combination of optic nerve head (ONH), peripapillary retinal nerve fiber layer (RNFL), and macular ganglion cell-inner plexiform layer (GCIPL) parameters measured via SD-OCT. The model generates a glaucoma probability score for individual patient, which is valuable information in assisting clinicians in glaucoma diagnosis and treatment.