Pierre-Alexandre Murena
email: murena@telecom-paristech.fr

Antoine Cornuéjols
email: antoine.cornuejols@agroparistech.fr

Jean-Louis Dessalles
email: dessalles@telecom-paristech.fr

Incremental Learning with the Minimum Description Length Principle

Whereas a large number of machine learning methods focus on offline learning over a single batch of data called training data set, the increasing number of automatically generated data leads to the emergence of new issues that offline learning cannot cope with. Incremental learning designates online learning of a model from streaming data. In non-stationary environments, the process generating these data may change over time, hence the learned concept becomes invalid. Adaptation to this non-stationary nature, called concept drift, is an intensively studied topic and can be reached algorithmically by two opposite approaches: active or passive approaches. We propose a formal framework to deal with concept drift, both in active and passive ways. Our framework is derived from the Minimum Description Length principle and exploits the algorithmic theory of information to quantify the model adaptation. We show that this approach is consistent with state of the art techniques and has a valid probabilistic counterpart. We propose two simple algorithms to use our framework in practice and tested both of them on real and simulated data.

I. INTRODUCTION

The current development of sensors and of the Internet of Things leads to a dramatically increasing number of automatically generated data: this volume was estimated to 2.8 ZB in 2012 by the Digital Universe Study [START_REF] Gantz | The digital universe in 2020: Big data, bigger digital shadows, and biggest growth in the far east[END_REF]. Many of these data appear in a stream fashion and thus cannot be stored in the machine's memory.

This new way of collecting data, as well as the gathering of data in ever-growing social media, induces two major problems that traditional machine learning algorithms cannot cope with. First, older data in streams have to be removed from the system memory and thus cannot be used anymore by the learner. Secondly, the data generating process may vary over time, in particular because of ageing effect [START_REF] Akhtar | Adaptive Facial Recognition Under Ageing Effect[END_REF] or of changes in the environment (e.g. in fraud detection [START_REF] Pozzolo | Credit card fraud detection and concept-drift adaptation with delayed supervised information[END_REF], [START_REF] Malekian | An adaptive profile based fraud detection framework for handling concept drift[END_REF]). This nonstationary nature of environments is called concept drift.

There exists several ways a concept drift can occur in a stream. If the data distribution changes permanently at a given time step, the drift is called abrupt. On the contrary, an incremental concept drift will be characterized by a continuous change in the distribution which can possibly end up in a new stationary state. In recurrent concept drift, the system alternates between several stationary distributions.

To deal with concept drift, two strategies can be adopted: active algorithms and passive algorithms. Active algorithms aim at identifying the distribution changes (in particular by using a triggering mechanism) and adapt the model only once a change is detected. On the contrary, passive algorithms adapt the model at each time step without detecting a change.

Many methods have been proposed to deal with concept drifts and state of the art algorithms are known to obtain good performances in streaming environments [START_REF] Gama | A survey on concept drift adaptation[END_REF], [START_REF] Ditzler | Learning in nonstationary environments: A survey[END_REF]. However, a formal approach to the problem is less often addressed.

In this paper, we propose a formal approach based on the Minimum Description Length (MDL) Principle. The MDL principle is an inductive principle assuming that the best theory to describe a system is the theory the programming of which has the shortest expression on a universal Turing machine. In particular, we will see that several existing methods are particular cases of our proposed framework. The MDL principle has already been applied to a wide variety of problems, among which model selection [START_REF] Mark | Model selection and the principle of minimum description length[END_REF], analogy reasoning [START_REF] Cornuéjols | Analogy and induction : which (missing) link ?[END_REF] and transfer learning [START_REF] Murena | Minimum description length principle applied to structure adaptation for classification under concept drift[END_REF].

In section II, we will expose the general framework and the Minimum Description Length Principle. We will show that our approach is robust to all concept drift types we previously exposed. After developing a general theory, we will provide simplifications to deal with problems of different natures. In section III, we will then propose two general algorithms which can be used to deal directly with the MDL principle. In section IV, we provide three different classes of models that can be used in our framework: probabilistic models, prototype-based models (inspired by Learning Vector Quantization and Self Organizing Maps [START_REF]Self-Organizing Maps[END_REF]) and models on finite alphabets. Using the prototype-based class of models, we give some experimental results we obtained using the proposed algorithms both for active and passive incremental learning.

II. A FRAMEWORK FOR INCREMENTAL LEARNING A. Notations for Incremental Learning

Let P be a problem space and S a solution space. At a time step t, the system receives a problem X t ∈ P and aims at predicting the solution Y t ∈ S. After predicting, the system may receive the actual solution Y t to the given problem.

This formalism is consistent with several usual situations. In particular it covers both cases of instance-incremental and batch-incremental learning [START_REF] Read | Batch-incremental versus instance-incremental learning in dynamic and evolving data[END_REF].

In instance-incremental learning, the system receives data one by one. At a step t, the learner receives a point x ∈ X where X is an input space (typically, X = R d ) and has to predict an output y ∈ Y (where the output space Y can be either continuous in regression, or finite in classification).

In batch-incremental learning, the learner receives a batch of data x 1 , . . . , x p and has to attribute a label y 1 , . . . , y p to each of the input points.

Using our notation is direct in both cases. In instanceincremental learning, the problem space and the input space are the same. In batch-incremental learning, a problem consists of a batch of instances of the input space X .

In both cases, the problem at time t is denoted by X t . In instance-incremental learning, X t ∈ X is directly an element of the input space. If X is a vector space, we denote by X i t the i-th coordinate of the vector X t . More generally, we will use the upper-script index to designate the coordinate of a vector. In batch-incremental learning, X t is given in form of a list and we will designate by X t,n the n-th element of X t . In particular, X t,n is an element of the input space X .

B. Minimum Description Length Principle

The idea of our approach is to use the Minimum Description Length Principle as an inductive principle in incremental learning. Firstly introduced by Solomonoff's theory of induction [START_REF] Solomonoff | A Formal Theory of Inductive Inference: Parts 1 & 2[END_REF], this inductive principle is exposed formally by Rissanen [START_REF] Rissanen | Modeling by shortest data description[END_REF]. The MDL principle is defined as follows:

The best theory chosen to describe observed data is the one which minimizes the sum of the description length (in bits) of:

• the theory description • the data encoded using the theory The description lengths are calculated using the notion of Kolmogorov complexity [START_REF] Li | An Introduction to Kolmogorov Complexity and Its Applications[END_REF]. Given a string x, the Kolmogorov complexity of x, denoted by K(x) is defined as the length (in bits) of the shortest program on a prefix Universal Turing Machine M (i.e. a Turing machine which produces decodable codes) generating the string x:

K M (x) = min p∈P M {l(p); p() = x} (1)
In equation 1, the term P M designates the set of all programs on the Universal Turing Machine M; l(p) is the length of the program p and p() is its output. A similar definition is proposed for the conditional Kolmogorov complexity of x given a string y:

K M (x|y) = min p∈P M {l(p); p(y) = x} (2)
where p(y) is the output of the program p taking y as input parameter.

It is important to notice in equations 1 and 2 that the quantity is defined for a given Universal Turing Machine M. The invariance theorem guarantees that the complexities for two different machines are equal up to a constant. Based on this theorem, a machine independent definition of Kolmogorov complexity can be defined but will not be used in this article. This choice is motivated by two reasons: first, the exact Kolmogorov complexity is not calculable, as it requires an evaluation over all programs on all Universal Turing Machines. On the other hand, the choice of M can be used and exploited as an inductive bias: by selecting a specific Turing machine, we will control how the compression is operated. This choice may seem arbitrary at first sight; but it is of a same nature as the choice of a restricted class of hypotheses which is used in Machine Learning for both theoretical and practical purpose. In the following, we will abusively omit to mention the machine in the expression of the complexity.

The ideal MDL can formulated using the notion of Kolmogorov complexity. Given data D, the model M chosen to describe an observation is the model minimizing the objective K(M ) + K(D|M ).

C. Simplifying Assumptions for Incremental Learning

The following assumptions correspond to a choice for Turing machine M.

To deal with incremental learning, we will consider that each time step t will be associated to a model M t . This model is used to describe both the problem X t and the solution Y t . The problem X t is described directly with the model M t . The solution Y t is described with the help of a decision function β t (the classifier in a classification problem or regressor in a regression problem). The decision function β t is described directly with the help of the model M t . In order to apply this description (summarized graphically in figure 1) to the mathematical expression of MDL principle, we use the chain rule: if A and B are two strings, the following equality is verified:

K(A, B) = K(A) + K(B|A) + O(1).
Applying this equality to our description, we obtain:

K(X t , Y t |M t ) = K(X t |M t ) + K(β t |M t , X t ) + K(Y t |M t , X t , β t ) (3)
This value of complexity expressed in equation 3 applies only to one time step. Further hypotheses over model transfer have to be expressed to deal with the whole streaming process.

In the following, we will use the notation X :T to designate the vector (X 1 , . . . , X T ) of the T first values of X. We will use the equivalent notation for the solutions (Y :T ).

We decompose the global model as a description of all models up to time T . This decomposition is denoted M :T . Besides, we suppose that the data (X t , Y t ) at any time step t are described only by the corresponding model M t . In particular, they are not described by any previous data. As a consequence, at a time step T , the global description of the data is given by:

K(D|M ) = K(X :T , Y :T |M :T ) = T t=1 K(X t , Y t |M t ) (4)

D. Model Complexity

The term K(M :T ) = K(M 1 , . . . , M T ) is of major importance in the sense that it encodes the whole transfer process between time steps and will be responsible for the adaptation to the concept drift.

We will adopt a markovian point of view on the model description. At time step T , we consider that all the models M t with t ≤ T are described using previously defined models. For any t ≤ T , we define an association function

∆ t : {1, . . . , t -1} → {0, 1} (5)
such that ∆ i (j) = 1 if model i can be described with the help of model j, and ∆ i (j) = 0 otherwise. Using these association functions, we can express the complexity of the models up to time T as:

K(M 1 , . . . , M T ) = T t=1 K(M t |M ∆ -1 t ({1}) ) (6)
In this equation, the notation ∆ -1 t ({1}) designates the set of indices i such that ∆ t (i) = 1.

In practice, the association functions ∆ t can be either fixed by the system or learned online.

The choice of the functions ∆ t is of major importance in theory and in practice, because it offers to the system the possibility to store previously acquired knowledge and thus to memorize states of interest. A constant effort for obtaining this property has been deployed in recent techniques [START_REF] Hosseini | Using a classifier pool in accuracy based tracking of recurring concepts in data stream classification[END_REF]. Several choice scenarios can be considered in our case, which all correspond to state of the art methods (figure 2):

• When ∆ i (j) = 0 for all i, j: all models are a priori independent in terms of description. The model is learned completely at each time step. • When ∆ i (j) = 1 for all i, j: the whole past models are taken into account to describe the present model. • When ∆ i (j) = 1 only for j ≥ i -h with a fixed h: the present model can be described with the last h models, which correspond to a sliding window of fixed size. The fixed sliding window is used in several algorithms like FLORA [START_REF] Widmer | Learning in the presence of concept drift and hidden contexts[END_REF]. • When ∆ i (j) = 1 only for j ≥ i -h with a size h estimated by the system: the size of the sliding window is not fixed anymore but heuristically adapted to the current problem. Such techniques are used in ADWIN [START_REF] Bifet | Mining adaptively frequent closed unlabeled rooted trees in data streams[END_REF]. Besides the choice of the association functions ∆ t have consequences over the minimization objective at each time step. Indeed, 1 values of these functions impose to add terms relative to previous observations in the complete Kolmogorov complexity. In practice, this is not necessarily possible, in particular because old data are not stored in memory anymore. Hence, considering old states in the complexity is possible only with simplifying assumptions or with the use of a stochastic gradient descent for the minimization.

E. Optimization problems

On the contrary to offline learning which proceeds to model learning in only one step, incremental learning has to operate in multiple phases: At each time step, the system has to answer to the new occurring problem X t and to adapt to the actual solution Y t .

In a strict use of the MDL principle, it would be possible to optimize the model M when observing the new problem X t : Indeed the data D can be considered as the delimited concatenation of all previously observed problem and solution couples and of the new observation:

D = (X 1 , Y 1 ), . . . , (X t-1 , Y t-1 ), X t (7)
Using the previously defined notations and complexity terms, the global complexity of the system at time step T is then defined as:

K T = K(M :T ) + T -1 t=1 K(X t , Y t |M t ) (8)
The learning problem that the system is facing, from the Minimum Description Length point of view, is the minimization of this global complexity K T term over the whole history of models M :T up to time T . This minimization is complex for various reasons. One of them is the interdependency of the M t which appear in several terms of the sum (because of the ∆ t functions). In particular, all models have to be learned again at each time step. In the perspective of data stream mining, this global update is obviously impossible and algorithmic hypotheses have to be chosen in order to overcome this difficulty.

At time step T , we choose to optimize the complexity over M T -1 and M T only. The past models M :T -2 are supposed to be fixed and can be used for the description of the last two models. This assumption is coherent with the incremental paradigm in which data are not stored and computations must be fast. Re-optimizing the whole modeling process would require to store all data X :T and would be highly time-consuming.

Hence, the optimization problem to solve at each time step is the following: minimize

M T -1 ,M T K T -K T -1 + C(X T -1 |M T -1 ) + C(M T -1 |M ∆ -1 T -1 ({1})) (9)

F. Probabilistic interpretation

A use of Kolmogorov complexity to measure the quality of a learning is not the most frequent paradigm in Machine Learning: A large range of methods rely on statistics and probability theory.

Our method is not in opposition to those techniques; on the contrary, they can be shown to be limit cases of or complexity-based framework. Basically there exists a strong correlation between Kolmogorov complexity and probabilities. This correlation is justified intuitively by the Shannon-Fano code. Over a finite alphabet, the Shannon-Fano code consists in associating each character to a complexity based on its occurrence frequency with the rule P r(X) = 2 -K(X) (or equivalently K(X) = -log P (X), where log designates the logarithm with base 2 as always in this paper).

The precise theoretical conditions to substitute the terms K(M ) and K(X|M ) by their probabilistic counterparts are discussed in [START_REF] Vitanyi | Minimum description length induction, bayesianism, and kolmogorov complexity[END_REF].

Focusing on this interpretation, the minimization objective in equation 9 corresponds to Bayes rule and the minimization over the model (hence the probability distribution or the classifier) corresponds to the Bayesian criterion used in Maximium A Posteriori and Maximum Likelihood Principles.

The very nature of the objective remains unclear. Two approaches may be adopted, which correspond to the notions of discriminative and generative learning [START_REF] Ng | On discriminative vs. generative classifiers: A comparison of logistic regression and naive bayes[END_REF]. The purpose of discriminative learning is to establish a separation between classes without estimating the density of the input points. On the contrary, generative learning focuses on the estimation of the generating process of the points.

In the terms of our framework, generative learning will consist in choosing a class of models describing the data: in a probabilistic setting, this class may be the class of joint distributions over the input and output space. Discriminative learning corresponds to a situation where the model M cannot be used to describe data, which means in terms of complexity that K(X|M ) = K(X). In generative learning, the model will thus correspond directly to the decision function β.

G. A general principle

As stated above, the proposed framework offers a common description to a large variety of existing methods. We propose to examine how state of the art methods can be expressed in the terms of our framework.

1) Active methods: Active approaches of incremental learning aim to detect change time steps t explicitly. Previous knowledge is used while no change has been detected.

Adaptive sliding window approach (ADWIN) [START_REF] Bifet | Mining adaptively frequent closed unlabeled rooted trees in data streams[END_REF] stores all previous data points into a sliding window. At each step, the sliding window is split in two sub-windows: the method detects a drift if the difference of the errors for each subwindow is higher than a parameter. ADWIN framework can be interpreted in terms of Kolmogorov complexity as follows. The model M t at time t corresponds to the concatenation of data memorized in the sliding window and of the corresponding classifier. Considering constant sizes for the description of data and description functions, and assuming that the decision function β is a constant during the whole process, the objective to minimize over β is

T t=1 K(Y t |X t , β).
R n+1:T (β) -R n+1:T (β 1 ) > K(β 1 ) T -N (10)
Similar arguments may be given to other methods based on sliding windows. In Probabilistic Adaptive Windowing (PAW) [START_REF] Bifet | Efficient data stream classification via probabilistic adaptive windows[END_REF], the probabilistic selection of samples to remove from the window corresponds to a non-deterministic way to minimize the global complexity.

All methods based on a statistical test are closely related to our model: they all check a regularity in data, hence a way to compress them. For instance, the validity of the Sequential Probability Ratio Test (SPRT) [START_REF] Wald | Sequential Analysis[END_REF] in our context can be established in the same way as done for ADWIN.

2) Passive methods: Passive approaches of incremental learning do not consider abrupt changes but continuously adapt the learned decision functions to new incoming data.

Most passive methods rely on ensemble learning, and in particular on bagging. The key point of bagging is that at any time t the system relies on a pool of base learners {h i t } 1≤i≤N . We consider that this pool corresponds to the model M t . In order to deal efficiently with concept drift, ensemble methods have to remove outdated experts from the pool. Several strategies are used to select learners to eliminate, including systematic removal of experts with performance lower than a threshold, or elimination of the worst base learner. The elimination of a base learner corresponds to a change in the model, which has a cost in terms of complexity. The higher value of K(M t |M ∆ -1 t ({1}) ) has to be compensated by a lower value of K(Y t |M t , X t , β t ), hence a better performance. To express the value of Y t using M t and β t , a method based on a majority vote of all expert learners is used, thus the correction term K(Y t |M t , X t , β t ) can depend on all experts.

III. ALGORITHMS

In this section, we will develop methods to solve optimization problem 9 in a context of incremental learning (i.e. with low memory and high speed).

A. Dealing with previous models

We propose to classify algorithms depending on the way they deal with previously acquired models. As mentioned earlier, the dependency on the past is given by the complexity term K(M t |M :t-1 ) which encodes the description length of model M t at step t and the previous models M :t-1 up to step t. Using the previously defined association functions ∆, the expression has already been simplified into K(M t |M ∆ -1 t ({1}) ). In order to describe model M t with the help of the set M ∆ -1 t ({1}) , two strategies may be chosen: either use many models in the set or select one single model.

The first strategy is employed in all ensemble learning methods [START_REF] Street | A streaming ensemble algorithm (sea) for large-scale classification[END_REF], [START_REF] Elwell | Incremental Learning of Variable Rate Concept Drift[END_REF], [START_REF] Brzezinski | Combining block-based and online methods in learning ensembles from concept drifting data streams[END_REF]. Such as in classical machine learning, ensemble learning methods construct the solution to a new problem by considering a weighted sum of the predictions of previous models.

In the second strategy, the key idea is to select the optimal model inside the set of predecessors M ∆ -1 t ({1}) . The selected predecessor is the best model in the sense of MDL principle.

In the perspective of selecting one single predecessor for each model, the total objective of equation 9 can be divided in two parts. The first part corresponds to the description of completed data at time step t once the solution has been given to the system:

φ 1 (M t-1 , M ) = K(M t-1 |M ) + K(X t-1 |M t-1 ) + K(β t-1 |M t-1 , X t-1 ) + K(Y t-1 |M t-1 , X t-1 , β t-1 ) (11)
The second part corresponds to the description of incomplete data at time step t:

φ 2 (M t , M ) = K(M t |M ) + K(X t |M t ) (12)
Equation 9 can be reformulated as:

minimize Mt,Mt-1, M1, M1 φ 1 (M t-1 , M1 ) + φ 2 (M t , M2 ) (13)
In the following, we propose two algorithms to solve problem 13 with a general class of models. The first algorithm calculates model transformations at each step; The second algorithm exploits a dynamic memory of previously calculated models.

B. First Algorithm: Continuous Adaptation

In Continuous Adaptation Incremental Learning, the system infers a new model at each time step t for both model M t-1 and M t . The system has access to all previously learned models M :t and chooses a predecessor among the models M ∆ -1 t-1 ({1}) and M ∆ -1 t ({1}) . In practice, we separate the choice of the predecessor for M t-1 and for M t . Such a separation is motivated by the fact that the description of data at time t depend on model M t-1 only by the transfer term in the case where the predecessor of

M t is M t-1 .
Consequently, and using the notations introduced previously, we can describe the learning algorithm in three steps:

1) Minimize the objective φ 2 (M t , M ) over predecessor

M ∈ M ∆ -1 t ({1})\{M t-1 } and M t 2) Minimize the objective φ 1 (M t-1 , M ) over predecessor M ∈ M ∆ -1 t-1 ({1}) and M t-1 3) Minimize the objective φ 1 (M t-1 , M ) + φ 2 (M t , M t-1 )
over predecessor M ∈ M ∆ -1 t-1 ({1}) , models M t-1 and M t In practice, this algorithm can have interesting properties in terms of comprehension of the underlying process. We propose to represent the dependency between two models by a vertex in a graph of models. Such a graphical representation makes the dependencies obvious and would enable a user interpret the decision, for example by detect easily periodic behaviours.

C. Second Algorithm: Memory Based Learning

The second algorithm we propose uses an ordered memory of models. These models are used as a simplification in the minimization.

At time step t, the system has access to a memory of m t models denoted M m 0 , . . . , M m mt . The stored models are used as an approximation for the description models: Each model M t is associated to a model in memory. On contrary to the previous algorithm, a perfect adaptation to the data is not required at any step, which implies that the algorithm is based on approximations and is less precise than the first algorithm. However, it is more convenient for high speed streaming data. The algorithm relies on several approximations. First, a model is chosen only among the previously learned models. Once a model is chosen in the list, it can be updated using the new data. If the complexity difference between the base model and the updated version is high, the updated model may either replace the base model or be inserted as a new memorized model (figure 3). Another approximation consists in using the difference of position in ordered memory as complexity for model transfer:

K(M m k |M m l ) = log(1 + |l -k|) ( 14
)
This algorithm is more efficient than the first proposed algorithm in terms of computation time, in the sense that it avoids optimizing the model at each step. The model choice is made only by comparing the values of the complexity with different accessible models; the optimization of the model will be done only in a case where no model fits the observed data.

Hence, this algorithm corresponds to a process of active learning: the change of model has to be detected before the system actually calculates the new model.

IV. APPLICATIONS A. Symbolic incremental learning

A first application of our model concerns non-continuous data. The use of MDL principle has been introduced naturally to analogy reasoning [START_REF] Cornuéjols | Analogy and induction : which (missing) link ?[END_REF]. Analogy reasoning can be defined as a form of reasoning in which one entity is inferred to be similar to another entity in a certain respect, on the basis of the known similarity between the entities in other respects [START_REF] Indurkhya | Modes of analogy[END_REF]. In particular, proportional analogy concerns situations of the form "A is to B as C is to D" (which will be denoted A : B ::

C : D).
An interesting example of proportional analogies have been suggested by D. Hofstadter [START_REF] Hofstadter | Fluid Concepts and Creative Analogies: Computer Models of the Fundamental Mechanisms of Thought[END_REF]. In these analogies, the terms are strings made up of alphabet letters. In this micro-world, it is possible to define a model as a generic construction rule for both words and the associated transformation β.

The framework we proposed for incremental learning can be applied directly to sequences of proportional analogies. We define a sequence of proportional analogy as a sequence of the form X 1 : Y 1 :: X 2 : Y 2 :: . . . :: X T : Y T . Some empirical methods have been proposed to solve one single problem of proportional analogy (for example the Copycat project [START_REF] Mitchell | Copycat: A Computer Model of High-level Perception and Conceptual Slippage in Analogy-making[END_REF]) but the issue of sequences of analogies has never been addressed, whereas it is of major importance in human cognition. Our framework offers a direct way to deal with this issue.

B. Probabilistic Models

Probabilistic models form another class of models which can be used in the Minimum Description Length Principle and especially in our approach of incremental learning. We define a probabilistic model as a joint distribution P X ×Y over the input space X and the output space Y.

Using the probabilistic model M = P X ×Y to describe data, we obtain:

K(X|M ) = n i=1 -log P X [X i ] (15)
The classifier β is obtained from M as follows:

β(x) = arg max y∈Y P Y|X [Y = y|X = x] (16)
and thus K(β|X, M ) = 0.

If the models M t and M u correspond to the respective distributions P t and P u , then the model transfer term K(M t |M u ) can be interpreted as the cross-entropy H betwen P u and P t :

K(M t |M u ) = H(P u , P t ) (17)
The properties of such models in relation with information theory opens new theoretical perspectives regarding incremental learning.

C. Prototype-based models

In order to test our framework, we have selected a simple visual class of models. Our model is inspired by the Nearest Prototype Classifiers such as Learning Vector Quantification (LVQ) [START_REF]Self-Organizing Maps[END_REF]. Such algorithms rely on a set of virtual points called prototypes used to approximate the classification locally: the classifier maps a point in the input space to the class of its closest prototype.

Consider a prototype set P of size K. A prototype P k ∈ P is a tuple P k = (p k , c k ) where p k ∈ X and c k ∈ Y. Using the prototype set P , a vector x ∈ X is linked to the closest prototype in P and its complete description includes the specification of the index of the prototype and the encoding of the relative position between the two vectors. With this modeling, the complexity of the data point x is:

K(x|P ) = min k=1,...,K K(x -p k ) (18)
We choose to work on a basic Universal Turing Machine which does not compress over the structure of a set of data: A set of n vectors is encoded by the individual description of each of the vectors. Besides, a vector is given by the description of all his coordinates one by one. Hence, if X is a design matrix, then:

K(X) = n i=1 K(X i ) = n i=1 d j=1 K(X d i ) (19)
Given a prototype model P , a classifier β is built by selecting the class of the closest prototype. The method used to build β simply enumerates the prototypes, hence K(β|M, X) = 0.

Given a classifier β and input data X ∈ X n , the complexity K(Y |P, X, β) corresponds to a correction term between the predictions Y = β(X) and the actual outputs Y . The correction method enumerates a list of correction terms, displaying both the identifier of the point (of description length log n) and the actual class (of description length log |Y|). Hence: [START_REF] Bifet | Efficient data stream classification via probabilistic adaptive windows[END_REF] where I designates the identity function the value of which is 1 if its argument is true and 0 otherwise.

K(Y |P, X, β) = (log n + log |Y|) n i=1 I(Y i = β(X i ))
In order to define the complexity K(M t |M ∆ -1 t ({1}) ), we propose the following process: Each model M t is described as a transformation Φ(M ) of a single model M ∈ ∆ -1 t ({1}). Thus, the transfer algorithm is divided in two parts:

1) For all models M ∈ ∆ -1 t ({1}), calculate the optimal transformation Φ 2) Select the past model M and the transformation Φ with lower density. The first step is obviously the longest and its time complexity depends on the number of previous models considered by the system (i.e. the cardinality of ∆ -1 t ({1})). In case of a firstorder markovian process, only the model M t-1 is considered by the algorithm.

In our prototype-based framework, we choose basic transformations for the models. Three transformations are considered:

1) Prototype addition: a prototype is added to the model (for example when a new class appears) 2) Prototype deletion: a prototype is removed because it has become useless. 3) Prototype shifting: the prototype position is changed. The methods suggested for both prototype addition and removal are highly similar to the methods proposed by [START_REF] Grbovic | Learning vector quantization with adaptive prototype addition and removal[END_REF] for Online Learning Vector Quantization. Our method just like theirs is based on a score for both operations, however the score we propose is more direct as it originates directly from the theory itself.

We consider that a prototype model P (of cardinal K) can be extended into a prototype model P + of cardinal K + 1 if adding one extra prototype to the model makes the global complexity decrease. Besides, adding a prototype is compulsory when a new class appears. When adding a new prototype, an additional term K add has to be added to the complexity K(M t |M :t )): this term corresponds to the complexity of adding a prototype p K+1 to P :

K add = K(p K+1 ) + K log(K + 1) (21)
In a similar way, a prototype model P of cardinal K can be reduced into a prototype model P -of cardinality K -1 if removing one particular prototype from the model makes the global complexity decrease. This is the case in particular when no input point is linked to a prototype in a new configuration for a batch incremental learning problem. Removing prototype p i from P requires to add an extra complexity term:

K rem = K(i) + K log(K -1) (22)
The term K(i) can be chosen either independent of i or not. If K(i) is constant, the cost to remove any prototype is the same. On the contrary if K(i) = log(1 + i), the cost to remove a recently added prototype is higher than the cost to remove an older prototype.

As explained before, we choose to apply a global transformation function Φ to the whole set of prototypes in order to describe the prototype shifting. In our proof of concept, we focus on linear functions of the form Φ(p) = p+δ+δ(p) where δ is a global translation vector and δ(p) a local translation vector. The translation vector δ applies to all prototypes in the model. The vector δ(p) depends on the prototype p and is a correction term. The complexity of the shifting process is:

K shif t = K(∆) + K k=1 C(δ(p)) (23)
The algorithm used to infer the prototype-based models at time step T is a variant of the Expectation-Maximization (EM) algorithm [START_REF] Dempster | Maximum likelihood from incomplete data via the em algorithm[END_REF]. EM algorithm alternates two steps: the Expectation step (E step) and the Maximization step (M step).

In the E step, data points are associated to their closest prototype, both in data for time T -1 and for time T . Using this fixed point-prototype association, the position of the prototypes is modified by minimizing the global complexity. This new position is then supposed to be fixed.

In the M step, the classes of the prototypes are changed if necessary using their new positions. The decision to insert or remove prototypes from the considered model is also made in the M step.

D. Experimental results

In order to illustrate the pertinence of our framework in practice, we have tested its performances on classical data sets with the naive prototype-based model. We considered three datasets: SEA [START_REF] Street | A streaming ensemble algorithm (sea) for large-scale classification[END_REF] (50000 instances, 3 attributes, artificial), Weather [START_REF] Elwell | Incremental learning of concept drift in nonstationary environments[END_REF] (18159 instances, 8 attributes, real) and Electricity Market [START_REF] Harries | Splice-2 comparative evaluation: Electricity pricing[END_REF] (45312 instances, 3 attributes, real).

For all datasets, we tested the instance-incremental version of our algorithm (by streaming directly over the data) and the batch-incremental version (by grouping successive data into a same batch). The experiments were all done with a fix sliding window size: |∆ -1 t ({1})| = 3. We tested the two proposed algorithms.

Table I presents the performances for different size S of batches. When S = 1, the situation corresponds to a problem of instance-incremental learning. Otherwise, the situation corresponds to batch-incremental learning.

As the purpose of this paper is not to establish a competitive performance for the suggested algorithms with the prototypebased model, we do not propose any comparison to state of the art algorithms. The key idea is that the obtained results are not necessarily better but similar to existing methods. Developing more accurate algorithms will be an improvement perspective to this paper.

In practice, the calculation time with the passive approach makes impossible to use this algorithm directly for a realtime process: the optimization algorithms take too much time even in the case of instance-incremental and makes the system unable to deal with a real data stream. We notice that in general, the active algorithm has lower accuracy than passive algorithm. This is coherent with the fact that this algorithm relies on strong simplifying hypotheses.

The obtained results are not competitive with state of the art algorithms, which was expected: Our method is not specifically designed to perform well and the class of models is very basic. We tested a very direct application of the equations presented above regardless of time complexity nor performance of the method. An effort has to be made in this direction in future works. The results are good enough to validate our framework yet.

V. CONCLUSION

In this paper, we have presented a new modeling the problem of incremental learning. This modeling offers a common description for several existing algorithms and proposes a unified approach of passive and active incremental learning by the mean of Kolmogorov complexity. In particular, it has been shown that our modeling is consistent with the sliding window methods, ensemble learning and Markovian models. We have proposed two elementary algorithms based on the chosen inductive principle and we have tested our framework with a simple class of models. By the mean of this tests, we have shown its efficiency to deal both with batch incremental and instance incremental learning in a context of concept drift.

Our work opens new perspectives in terms of theoretical approach of incremental learning. Besides, the formalism we developed can address a similar yet different problem: education. Instead of adapting models from step to step, education aims at finding the right steps to transfer a source concept into a target concept. Finally, the link with information geometry cannot be neglected: using probabilistic models, we have enlightened the use of well-known quantities to infer a discrete trajectory inside a manifold of probability distributions.

Fig. 1 .

 1 Fig. 1. Simplified framework for the description of problem Xt and solution Yt using model Mt. On this graphic representation, an arrow from A to B means that string B is described using string A.

Fig. 2 .

 2 Fig. 2. Possible choices for ∆ function. The color of the square at line i and column j indicates the value of ∆ i (j). If the square is dark, ∆ i (j) = 0; If the square is white, ∆ i (j) = 1

  If a drift is detected at time step n, the last T -n terms in the sum are changed in take into account the introduction of a new decision function β 1 and of the new data description using β 1 . As the complexity term K(Y t |X t , β) corresponds to a correction term, a link can be established with empirical risk. Denoting R n+1:T (β) the empirical risk of binary classifier β on the last T -n data, the detection of a break point obeys the criterion:

Fig. 3 .

 3 Fig. 3. General schema of Memory based algorithm

ACKNOWLEDGMENT

This research is supported by the programme Futur et Ruptures (Institut Mines-Telecom).


Le Stum

How high is herbage intake of organic rabbits grazing fescue or sainfoin ?

.02

Organic rabbit breeding Conventional Organic Photo Lebas F. Photo Gidenne T.

Potential production unknown

Most relevant management practices yet to be defined (herbage allowance, stocking rate, regrowth) Grazing management Roinsard et al (2016) Lack of knowledge on herbage intake / nutrition

.03

The moving cages Herbage intake measurements .08

Results Over the whole grazing period .014

Thank you for your attention

•

  Photo Legendre H.

  Photo Legendre H.Photo Duprat A.

Acknowledgements

Financial support from INRA division PHASE (MarkPast project)+ INRA AgriBio4 (CuniPat project) + GISA metaprogram of INRA (PROF project)

The authors thank the AVEM association for providing organic sainfoin's seeds, and all of the colleagues involved in data collection and analysis


Anna Bastet

Baptiste Lederer

Nathalie Giovinazzo

Joan Estevan

Christophe Robaglia

Jean-Luc Gallois

A synthetic eIF4E1 resistance allele allows resistance pyramiding and broad resistance spectrum to potyviruses in Arabidopsis thaliana

des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Developing genetic resistance is an important alternative to the use of pesticides in order to avoid viral infections in crops. However, this approach can be limited by the lack of natural resistant alleles available for plant breeding. Biotechnologies can supply for this lack by altering susceptibility host factors used by viruses. However, strategies relying on knocking out plant host factors can be undermined by associated developmental default, or gene redundancy (Gauffier et al., 2016). Synthetic biology could offer the opportunity to develop functional synthetic alleles associated with plant resistance at no developmental cost. The translation initiation factors (eIF4E), that are associated with resistance to most single strand positive RNA (including the major group of potyviruses), are the ideal targets to validate such approaches. We focused on the Arabidopsis thaliana eIF4E1 gene that encodes a susceptibility factor to a potyvirus, the Clover yellow vein virus. Artificial resistance alleles were created by directed mutagenesis and used to complement loss-of-function eif4e1 plants for assessing in planta both their functionality and resistance efficiency. By combining these synthetic resistance alleles with others eIF4E factors-mediated resistance, a broad resistance spectrum to potyviruses is expected as well as increased resistance durability. This study will make proof of the feasibility of this approach to obtaining large resistance spectrum in plants with the perspective of extending it to crops.

Reference : Gauffier, C., Lebaron, C., Moretti, A., Constant, C., Moquet, F., Bonnet, G., Caranta, C. and Gallois, J.L. (2016) A TILLING approach to generate broad-spectrum resistance to potyviruses in tomato is hampered by eIF4E gene redundancy. Plant J. 85: 717-729. .

synthetic eIF4E1 resistance allele allows resistance pyramiding and broad resistance spectrum to potyviruses in Arabidopsis thaliana

  16eme rencontres de Virologie Végétale, du 15 au 19 janvier 2017, AUSSOIS, France. Communication Orale, Anna Bastet A Bastet Anna 1-2 , Lederer Baptiste 1 , Giovinazzo Nathalie 1 , Estevan Joan 1 , Robaglia Christophe 2 and

	Gallois Jean-Luc 1
	1 Unité Génétique et Amélioration des Fruits et Légumes, INRA Avignon, BP 94, F-84143
	Montfavet, France
	2 Laboratoire de Génétique et Biophysique des Plantes, CEA-CNRS-Université Aix-Marseille
	II, Faculté des Sciences de Luminy, F-13009 Marseille, France
	Mots-clés : Potyvirus, resistance, Arabidopsis thaliana, eIF4E


B Wenden

T Barreneche

M Meland

M M Blanke

Harmonisation of phenology stages and selected cherry cultivars as bioindicators for climate change

Keywords: cherry (Prunus avium L., Prunus cerasus), climate change, flowering, global warming, low chill, phenology, sustainable horticulture

Perennial fruit crops phenology such as cherry is an ideal bio-indicator of climate change due to their long-lasting features, in particular, dates of flower opening and full bloom. This implies i) the use of several generations of cherry trees/orchards and ii) the use of the same original cherry cultivars, which existed as bearing trees and were replanted after the orchard had been grubbed. A comparison of available definitions of phenological stages in cherry previously used independently throughout Europe showed overlaps and shortcomings; hence, harmonisation was reached in this respect in the COST Cherry FA 1104 working group 2 (cherry phenology and climate change) based largely on the acceptance of the BBCH scale. This contribution presents the agreed phenology stages in both visual and wording evidence. Similarly, this contribution presents the agreed cultivars to be monitored in future for phenology and climate change effects for harmonisation. For sweet cherry, this EU-wide harmonisation includes 'Burlat', 'Cristobalina' and 'Rita' as early, 'Stella' and 'Van' as medium flowering and 'Sweetheart', 'Regina' and 'Bigarreau Noire de Meched/Germersdorfer' for late flowering cultivars for climate change effects. For sour cherry, this harmonisation resulted in 'Meteor korai' and 'Anglaise Hative' for early flowering, 'Chrisana Pandy' and 'Erdibotermo' for medium flowering and 'Schattemorelle', 'Iiva, Ujfehrtoifurtos (Balaton)' for late flowering.

INTRODUCTION

Cherry is an ideal candidate for climate change studies (Luedeling et al., 2013a,b), since distinct phenological stages were defined decades ago (Baggiolini, 1952;Lichou, 1990;[START_REF] Meier | Phänologische Entwicklungsstadien des Kernobstes (Malus domestica Borkh. und Pyrus communis L.), des Steinobstes (Prunus-Arten), der Johannisbeere (Ribes-Arten) und der Erdbeere (Fragaria x ananassa Duch.)[END_REF] and have been recorded by simple means like filing cards for a long time as part of phenology research in the 1970s before computers and climate change emerged. In the 1990s, these records were often classified as old-fashioned research and no longer continued or research sites closed down or re-located. At the few locations where the records were maintained, it was important that the records continued on the same original cherry cultivars to show any long-term changes in both climate and phenology. In addition, these cultivars are often used as reference in breeding and evaluation programs, thus maintaining observation campaigns. Four phenological definitions or scales have been used to assess cherry flowering in Europe since the 1970s (Fleckinger, Baggiolini, 1952;Lichou, 1990;[START_REF] Meier | Phänologische Entwicklungsstadien des Kernobstes (Malus domestica Borkh. und Pyrus communis L.), des Steinobstes (Prunus-Arten), der Johannisbeere (Ribes-Arten) und der Erdbeere (Fragaria x ananassa Duch.)[END_REF]. The combination of such historic meteorological and phenological data from the same site are now considered a valuable tool and treasure in assessing the impacts of climate change on a local scale.

The objective of this work was to harmonise terminology for phenological stages and provide a guideline for future European phenology records in cherry and beyond.

MATERIAL AND METHODS

Phenological data on filing cards included calendar dates of first flower opening (BBCH 61), full bloom (BBCH 65 or F1040 DD), harvest (BBCH 81) and leaf drop (BBCH 95) as well as late frost and consequent yield loss for a range of cherry cultivars, using only fruit bearing trees at the time. BBCH stages harmonisation for observation, together with reference cultivars for sweet and sour cherries were agreed within the meeting held by the COST Cherry FA 1104 working group 2 (cherry phenology and climate change) in Ullensvang (April 2014).

RESULTS AND DISCUSSION

Harmonisation of phenology stages in cherry

Figure 1a-f shows the descriptions for the agreed phenology stages for cherry, largely based on the BBCH scale [START_REF] Meier | Phänologische Entwicklungsstadien des Kernobstes (Malus domestica Borkh. und Pyrus communis L.), des Steinobstes (Prunus-Arten), der Johannisbeere (Ribes-Arten) und der Erdbeere (Fragaria x ananassa Duch.)[END_REF]. Table 1. Phenology harmonisation for future flowering records of sweet/sour cherry [START_REF] Meier | Phänologische Entwicklungsstadien des Kernobstes (Malus domestica Borkh. und Pyrus communis L.), des Steinobstes (Prunus-Arten), der Johannisbeere (Ribes-Arten) und der Erdbeere (Fragaria x ananassa Duch.)[END_REF].

Category

Harmonisation of cherry cultivars to be monitored in the EU for climate change effects

A major result and outcome of the EU COST Cherry FA 1104 WG 2 (Cherry phenology and climate change) meeting was the acceptance of one universal scale for all European countries (Table 1) and agreed standard cultivars (Table 2) for phenology monitoring. Cultivars were selected on the basis of early, medium and late flowering and most widely grown throughout Europe.

The BBCH-scale is a scale used to identify the phenological development stages of a plant. A series of BBCH-scales have been developed for a range of crop species. Phenological development stages of plants are used in a number of scientific disciplines (crop physiology, phytopathology, entomology and plant breeding) and in the agriculture industry (timing of pesticide application, fertilisation, agricultural insurance). The BBCH-scale uses a decimal code system, which is divided into principal and secondary growth stages, and is based on the cereal code system (Zadoks scale) . BBCH officially stands for "Biologische Bundesanstalt, Bundessortenamt und CHemische industrie.

CONCLUSIONS

The experts of the EU COST Cherry agreed to use the detailed and more precise German BBCH throughout Europe in the future for cherry flower monitoring and provide detailed descriptions in all European languages plus pictorial digital evidence in this chapter of this Acta Horticulturae. Continued data acquisition throughout Europe will support the development of predictive models for assessing the effects of climate change on cherry phenology but also as useful tools for breeding strategies.

Figure 1 .

 1 Figure 1. Developmental stages for cherry ( © He le ne Christmann, INRA).

Table 2 .

 2 Phenology harmonisation for future flowering records of sweet and sour cherry.

	Category	Sweet cherry flowering	Sour cherry flowering
	Early flowering	Burlat, Christobalina, Rita	Meteor korai, Anglaise Hative
	Medium	Stella, Van	Chrisana Pandy, Erdibotermo
	flowering
	Late flowering	Sweetheart, Regina,	Schattenmorelle, Iiva, Ujfehrtoifurtos
		Noire de
		Meched/Germersdorfer

ACKNOWLEDGEMENTS

We are grateful to EU COST Cherry FA 1104 for the travel grant to the WG 2 (Cherry phenology and climate change) workshop in Ullensvang in April 2012 to achieve these results. The authors thank He le ne Christmann (INRA) for the pictures. Baggiolini, M. (1952). Stades repe res du pe cher. Revue Romande d'Agriculture., Viticulture and Arboriculture 4, 29. Lichou, J., Edin, M. , Tronel, C. and Sounier, R. (1990). Le Cerisier (Paris, France : CTIFL).

Literature cited


Keywords: RAD-seq, Dicentrarchus labrax, QTL, Genomic Selection

come

European sea bass (Dicentrarchus labrax) is one of the most important farmed species for Mediterranean aquaculture. The observed sexual growth and maturity dimorphism in favour of females adds additional value to deciphering its sex-determining system. Current knowledge indicates the existence of a polygenic sex-determining system interacting with temperature. Restriction-site Associated DNA (RAD) sequencing was used in a test panel of 175 offspring originating from a F2 factorial cross between two dams and four sires from one F1 family between parents with contrasted sex tendencies. In total, 1,156,659,542 raw reads (100 bases long) were produced (578,329,771 paired-end reads). After removing low quality sequences (quality score under 30), ambiguous barcodes and orphaned paired-end reads, 76.7% of the raw reads were retained (886,927,866 reads) from which 56,696 unique RADtags were retrieved. In order to maximise the number of informative markers and minimise the amount of missing or erroneous data, we used RAD-tags retrieved in at least 75% of the samples, and carrying one or two single nucleotide polymorphisms (SNPs). The family structure was identified a posteriori using R/hsphase and Vitassign with the SNP markers from the offspring and parents. A consensus SNP-based linkage map was constructed, consisting of 5,097 SNPs grouped into 24 linkage groups. Indications for putative sex-determining QTL are provided in linkage groups 13, 19 and 21. A preliminary study was conducted testing whether breeding values, estimated by including the additive SNP effects, could be used as a predictive measure of phenotypic sex. Prediction of phenotypic sex within each of the eight full-sib families, using genomic models (MCMC-BLUP, BayesCPi), consistently outperformed pedigree BLUP. The advantage of using genomic models was particularly evident in dam half-sib family 2, where one of the sex-determining QTLs had been detected. Altogether, our results further support the polygenic hypothesis for sex determination in sea bass, while they also highlight the potential use of within-family genomic selection to rank individuals within families for their estimated breeding value for sex tendency.

  Investigating the genetics of sex determination in European sea bass (Dicentrarchus labrax) using RAD-seq Christos Palaiokostas, Michaël Bekaert, John B. Taggart, Karim Gharbi, Brendan J. Mcandrew, Béatrice Chatain, David Penman, Marc Vandeputte To cite this version: Christos Palaiokostas, Michaël Bekaert, John B. Taggart, Karim Gharbi, Brendan J. Mcandrew, et al.. Investigating the genetics of sex determination in European sea bass (Dicentrarchus labrax) using RAD-seq. 12. International Symposium on Genetics in Aquaculture, Jun 2015, Santiago de Compostela, Spain. hal-02737758 INVESTIGATING THE GENETICS OF SEX DETERMINATION IN EUROPEAN SEA BASS (Dicentrarchus labrax) USING RAD-SEQ Palaiokostas C.

  1 , Bekaert M.B. 1 , Taggart J.B. 1 , Gharbi K. 2 , McAndrew B.J. 1 , Chatain B. 3 , Penman D.J. 1# , Vandeputte M. 3,4 1 Institute of Aquaculture, School of Natural Sciences, University of Stirling, Stirling FK9 4LA, Scotland, United Kingdom 2 Edinburgh Genomics, Ashworth Laboratories, King's Buildings, University of Edinburgh, Edinburgh EH9 3JT, Scotland, United Kingdom 3 Ifremer, Chemin de Maguelone, 34250 Palavas-les-Flots, France 4 INRA, UMR1313 Génétique animale et Biologie intégrative, 78350 Jouy-en-Josas, France


Investigation of the production of wet agglomerates of durum wheat semolina using a fluidized bed granulator

Results

Water content impact:

-Larger diameters -Higher values of water contents -Lower values of compactness -No impact on the sphericity Agglomeration is a key operation in the elaboration of couscous grains of durum wheat semolina (made by successive operations: wet agglomeration, cooking and drying). It is usually conducted in mechanical mixers and determines the production yield, the final size, shape and functionalities of couscous. The main objective of this work is to explore the capability of the fluidised bed technology to produce agglomerates of durum wheat semolina.

Introduction & Objectives

Characteristics of the dried agglomerates in fluidized bed

Dried agglomerates of durum wheat semolina are characterized by irregular shape, with very rough surfaces. They appeared as semi discrete particles glued to each other. The strength to compress the granules was significantly lower than those measured for the classical couscous grains.

On the hydrotextural diagram, the structures become closer to the saturation curve with an increase in their size. Mechanisms of agglomeration in a fluidized bed show a first stage of expansion, which is induced by an addition of water. The saturation degree slightly increases. During the second stage, the structures continue to grow up but with densification mechanisms due to the capillarity influence.

Conclusion

Agglomeration in fluidized bed of durum wheat semolina is characterized by a higher amount of agglomerates with diameter between 1 and 2 mm (considering as a better yield for the production of couscous). Agglomerates of durum wheat are less compact than the classical couscous grains. Agglomeration is driven by a fractal formation process followed and by densification mechanisms for the large agglomerates. Specific model is used to describe the interparticle arrangements.

! ! ! = 1 1 + ! ! * !(1 -! ! !!! ! ! )
Semolina size impact:

-Larger diameters -Lower values of water contents -Lower values of compactness -Slight decrease of the sphericity

Characteristics of the traditional couscous grains in low shear mixer

The classical couscous grains were characterized by spherical shape with regular surfaces. The native semolina particles were poorly visible in the grains. The classical couscous grains were more compact and consolidated by complete starch gelatinisation during the cooking stage.

  in the diameter values is observed on the size distribution curve (d 50 = 650 µm). The water content of the wet agglomerates is positively correlated with their diameter. The increase of the d 50 of the wet agglomerates is also associated to a decrease of their compactness. Large agglomerates are more spherical. Mechanism for the wet agglomeration can be considered as fractal morphogenesis for the smaller structures (d 50 ≤ 1600 µm), followed by some phenomena of densification for the larger structures (d 50 ≥ 1600 µm) .

  Microstructure of the agglomerates based on of durum wheat semolina prepared with the fluidized bed (a) and commercial couscous grains (b) at different magnifica<ons. authors would like to thank the Agence Nationale de la Recherche (ANR ALID 2013) for its financial support through the program "Dur Dur" the variation of the compactness (link to the size of the structure) by using only the water content. Where d * s is the real density of the semolina particle, w is the water content, w m is the water content at S=50% and d is the reverse slope of the central quasi-linear section of the sigmoidal curve

Bettina Bellocq 1 , Bernard Cuq 1 , Thierry Ruiz 1 , Agnès Duri 1 , Kevin Cronin 2 , Denis Ring 2

1 UMR IATE 1208 CIRAD/INRA/Montpellier SupAgro/Université Montpellier. 2 Department of Process & Chemical Engineering, University College, Cork, Ireland. !