The invention relates to a method and a computer program for estimating a bilirubin level of a neonate, composed of the steps of:



Acquiring a series of bilirubin levels estimated at different time points from a sample obtained from a neonate,
Acquiring a plurality of covariates from the neonate, each composed of an information about a neonatal property,
Providing a pre-defined bilirubin model function, wherein the bilirubin model function is configured to describe a time course of a bilirubin level of a neonate,
Determining a plurality of model parameters of the bilirubin model function, wherein each model parameter is estimated from at least one covariate of the plurality of covariates and an associated population model parameter,
Determining from the series of acquired bilirubin levels and the bilirubin model function with the determined model parameters an expected bilirubin level of the neonate for a time particularly later than a lastly acquired bilirubin level of the series of bilirubin levels.