L'analisi dei dati è stata effettuata in collaborazione con il Professor Silvano Milani dell'Istituto di statistica Medica e Biometria dell'Università di Milano.( ) Le costanti di crescita individuali sono state interpolate con un modello lineare generale multivariato per stimare gli effetti sulla cinetica della crescita esercitati da età gestazionale, sesso, adeguatezza/inadeguatezza del peso per l'età gestazionale (soggetti AGA od SGA), sepsi tardiva, displasia broncopolmonare, sequele neurologiche maggiori a 2 anni.[ ] [ ] ) 1 ( ) 1 t ( 1 4 ) 1 ( ) 1 t ( 1 0 0 0 1 ) 1 t ( y(t) E + − + − + − + − 1 + + ⋅ + − = τ τ γ αIl modello è atto a descrivere la crescita sia nelle prime settimane (calo ponderale, recupero del peso neonatale, primo picco di velocità, riduzione della velocità, take-off, secondo picco di velocità, in sequenza), sia nel periodo di lento decremento della velocità che inizia dopo il terzo mese di vita. Il modello ha 4 costanti di crescita, ovvero, peso alla nascita (), calo ponderale (γ 0 ), età al minimo peso postnatale (τ 0 ), età al picco di velocità (τ 1 ) e 3 costanti di velocità per i periodi nell'intorno del peso minimo postnatale (β 0 ), del picco della velocità di crescita (β 1 ) e della fine del secondo anno di vita (γ 1 ). La parametrizzazione del modello (ad esempio e  anziché ) è stata scelta per ottenere una distribuzione delle costanti individuali pressoché gaussiana).Le curve medie ottenute con tale modello e relative ad uno dei predetti fattori sono "corrette" per gli effetti di tutti gli altri fattori presenti nel modello: ad esempio, la differenza tra la curva media delle femmine e quella dei maschi dipende solo dal sesso e non da eventuali differenze di distribuzione degli altri fattori tra maschi e femmine.