How to plot coordinates on tmap output?
<p>So far, I have managed to visually present data for different socio-economic factors at a district level in Pakistan with the help of the <code>tmap</code> package to create the following output:</p>
<p><a href="https://i.stack.imgur.com/nVPdX.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/nVPdX.png" alt="enter image description here" /></a></p>
<pre class="lang-r prettyprint-override"><code>accesstowater &lt;- tm_shape(pak) + 
    tm_fill(&quot;Condition&quot;, n = 5, palette = &quot;-Blues&quot;, colorNA = &quot;grey&quot;,
             textNA = &quot;Missing Values&quot;, title = &quot;Access to Water&quot;) + 
    tm_borders(alpha = 0.5) +
    tm_credits(&quot;Source: PBS - PSLM 2014 - 15\nDH Corp Ltd.&quot;, size = .64,
                position = c(.62,0.02))
</code></pre>
<p>In the code above, <code>pak</code> is the name of my Spatial Dataframe and <code>Condition</code> is a variable within that dataframe showing access to water.</p>
<hr />
<p>At the same time, I have used the <code>ggplot2</code> package to plot the geographic coordinates of various towers for different Mobile Network Operators in Pakistan (roughly 13000 towers). The output is as follows:</p>
<p><a href="https://i.stack.imgur.com/jDgpG.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/jDgpG.png" alt="enter image description here" /></a></p>
<pre class="lang-r prettyprint-override"><code>pakTowers &lt;- ggmap(pakmap6) + 
               geom_point(aes(Longitude, Latitude, colour = MNO), 
                          data = Towers, size = 0.7, alpha = 0.5,
                          show.legend = FALSE)
</code></pre>
<p>In the code above , <code>pakmap6</code> is the output of the <code>get-map()</code> function, <code>MNO</code> represents the Mobile Network Operator and the <code>Towers</code> dataset includes entries for each tower's coordinates.</p>
<p>My next requirement is to overlay the second plot over the first, essentially showing the towers on top of the district wise electricity scores.</p>
<p>Any ideas on which packages I could use (as well as sample code) in order to achieve this?</p>