Abstract: The method for inferring epileptogenicity of a brain region not observed as recruited or not observed as not recruited, in a seizure activity of an epileptic patient brain, includes: providing a computerized model modelling various regions of a primate brain and connectivity between the regions; providing the computerized model with a model able to reproduce an epileptic seizure dynamic in the primate brain; providing structural data of the epileptic patient brain and personalizing the computerized model using the structural data to obtain a virtual epileptic patient (VEP) brain model; translating a state-space representation of the VEP brain model into a probabilistic programming language (PPL) using probabilistic state transitions to obtain a probabilistic virtual epileptic patient brain model (BVEP); and acquiring electro- or magneto-encephalographic data of the patient brain and fitting the probabilistic VEP brain model against the data to infer the epileptogenicity of the brain region that is not observe