The study of the preparation phase of large earthquakes is essential to understand the physical processes involved, and potentially useful also to develop a future reliable short-term warning system. Here we analyse electron density and magnetic field data measured by Swarm three-satellite constellation for 4.7 years, to look for possible in-situ ionospheric precursors of large earthquakes to study the interactions between the lithosphere and the above atmosphere and ionosphere, in what is called the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC). We define these anomalies statistically in the whole space-time interval of interest and use a Worldwide Statistical Correlation (WSC) analysis through a superposed epoch approach to study the possible relation with the earthquakes. We find some clear concentrations of electron density and magnetic anomalies from more than two months to some days before the earthquake occurrences. Such anomaly clustering is, in general, statistically significant with respect to homogeneous random simulations, supporting a LAIC during the preparation phase of earthquakes. By investigating different earthquake magnitude ranges, not only do we confirm the well-known Rikitake empirical law between ionospheric anomaly precursor time and earthquake magnitude, but we also give more reliability to the seismic source origin for many of the identified anomalies.