From IQ signal to FM modulated carrier, how its done?
The solution to this is realizing the relationship between frequency and phase:
$$f(t) = \frac{d\phi(t)}{dt}$$
Where $f(t)$ is the instantaneous frequency as a function of time, and $\phi(t)$ is the instantaneous phase as a function of time. So by such definition, a linear change in phase represents a step change in frequency. This is exactly how we can use vector modulators as frequency translators (and equivalently single-sideband modulators). A rotating phasor on the IQ plane, if rotating at a constant rate represents either a positive or negative frequency depending on the direction of rotation. So with FSK using two frequency symbols we get positive and negative rotating phasors. The frequency realized is the rate of that rotation (cycles/sec).
Starting with MSK, which is a subset of FSK, and specifically the minimum frequency shift between the two FSK tones such that the two signals are orthogonal, for this we rotate the phase 90° over one symbol duration, such that the direction of rotation determines if a one or zero was transmitted. The reason for 90° is this results in the minimum spacing of the two frequency tones as first described (hence MSK). If we rotated further within the same time of a symbol period this would result in a higher frequency and thus increase the spacing between the tones (as in generalized FSK). Below shows the phase diagram of rotation versus symbol &quot;1&quot; or &quot;0&quot;. If we continued to send a long series of either symbol, it is quite clear that we are transmitting either a positive or negative frequency (relative to a central carrier frequency, so basically one of two possible tones: FSK).

When we transmit a series of symbols the frequency shift keying that is occurring is clear from the relationship between phase and frequency:

The sharp transition in frequency as a step is what creates a high spectral occupancy (sharp edges create high frequencies, consider all the spectral components of a square wave). By rounding these edges we minimize the spectral occupancy, and Gaussian Minimum Shift Keying specifically replaces the straight phase slopes with a Gaussian filtered transition; the phase specifically transitions as an integrated Gaussian pulse (thus the frequency follows a Gaussian transition- see diagram below). So the phase starts out slowly at the beginning of the symbol, reaches maximum rate mid way and then terminates slowly at the end where a possible otherwise abrupt direction change could occur. You can make a very simple (no multipliers so can be done in a micro-controller) GMSK modulator using a Gaussian Filter on the Frequency Control Word input to an NCO which is otherwise connected to the data stream (and similarly with a VCO in the analog world). (This is actually what I was getting at with this question: Simplest All Digital GMSK Modulator which can be combined with this solution for the simple no-multiplier Gaussian filter: Gaussian FIR filter with no multipliers?).  The GMSK phase and frequency diagrams below further and most immediately shows how this is a form of FM modulation.

GMSK is also commonly implemented as partial response signals where we overlap consecutive pulses with intentional inter-symbol interference which can later be completely undone (with Viterbi decoding for example), and thus further improves spectral efficiency since we can send more data in the same amount of time (at the expense of receiver complexity).  The strategic use of the wonderful Laurent Decompositions significantly simplifies those receivers— but that is another topic for another day.
Further FSK, MSK, GMSK are all constant envelope modulations so have the noted advantage that we can run the transmitter PA into saturation which offers significant power efficiency and without regard to immunity against intention jammers can hard limit the input instead of using an AGC (both of particular interest to low cost battery operated devices)