On a pulse-to-pulse timescale, pulsar emission is highly variable. However integrating over many pulses yields an average pulse profile that is very stable and unique to that pulsar.
The profiles of some pulsars however, change between two or more different shapes in a process known as mode-switching. Lyne et al (2010) showed that profile changes can correlate with the spin-down rate, implying that this phenomenon isn’t just rooted in the emission region but instead is a global magnetospheric process.
I’m using Gaussian Process techniques to investigate these correlations in detail using Lovell telescope data, as well as to search for emission changes near glitch events.