The paper presents a model for the evolution of an infectious disease in a population with individual-specific immunity. The immune state of an individual varies with time according to its own dynamics, depending on whether the individual is infected or not. The model involves a system of size-structured (first-order) PDEs that capture both the dynamics of the immune states and the transition between compartments consisting of infected, susceptible, etc. Due to the unavailability of precise data about the immune states of the individuals, the main focus in the paper is on developing a technique for set-membership estimations of aggregated quantities of interest. The technique involves solving specific optimization problems for the underlying PDE system and is developed up to a numerical method. Results of numerical simulations are presented for a benchmark model of SIS-type, potentially applicable to diseases like influenza and to various sexually transmitted diseases.