Abstract: A method of reconstructing a causation process from time varying data describing an event, the data consisting in a certain number of entities each having a position in a space, and each of the entities being characterized by at least a quantity or value relatively to at least one feature and in the quantity or value relatively to at least one of the features of the entities at least at two different times or at each time instant of a sequence of time instants; the method describing the higher likelihood transition of all entities i, j from the time n to the time n+1 as a function of the position coordinate of the entity I and of the entity j and the quantity of the at least one feature of the entity I and of the entity j at the time n and at the time n+1: Mi,j[n,n+1]=?(xi,yi,qi[n],qi[n+1],xj,yj,qj[n],qj[n+1]).