R语言中，用多组散点图加折线，展现不同地区美女数量变化趋势？

这篇博客是想向各位大佬请教问题，请看到的大佬不吝赐教！万分感谢！
有点标题党，为吸引大家进来看，请大家谅解！

我的目标是：在一个图中，将多组数据画成散点加折线图。

我用的语句和现在已经画出来的散点图：

ggplot(ceshia, aes(ceshia[,1],ceshia[,2],col=ceshia[,3])) + geom_point(size = 1)

ceshia是我的数据表，格式是这样的（完整的数据我放到文章最后面）：

    统计时间 美女数量   地区
1    2019-07           37 广州市
2    2019-08           37 广州市
3    2019-09           38 广州市
16   2019-07           10 珠海市
17   2019-08           10 珠海市
18   2019-09           10 珠海市
31   2019-07           33 汕头市
32   2019-08           33 汕头市
33   2019-09           33 汕头市

我想要画成下面这种类型（画得不好看，就是那种意思，把同一类的点用折线连起来）：
用折线表示同一个市的美女数量变化。

但是，实际上我画出来的图是这样的：

ggplot(ceshia, aes(ceshia[,1],ceshia[,2],col=ceshia[,3])) + geom_point(size = 1)+
  geom_line(group=1)

不知道怎么回事。在网上搜了很久，用了很多方法，都没成功。
有没有大牛能帮帮我，看看要怎么实现我想要的那个效果？
万分感谢！

完整数据如下：

    统计时间 美女数量   地区
1    2019-07           37 广州市
2    2019-08           37 广州市
3    2019-09           38 广州市
4    2019-10           37 广州市
5    2019-11           38 广州市
6    2019-12           37 广州市
7    2020-01           35 广州市
8    2020-02           19 广州市
9    2020-03           19 广州市
10   2020-04           19 广州市
11   2020-05           20 广州市
12   2020-06           24 广州市
13   2020-07           26 广州市
14   2020-08           28 广州市
15   2020-09           30 广州市
16   2019-07           10 珠海市
17   2019-08           10 珠海市
18   2019-09           10 珠海市
19   2019-10           10 珠海市
20   2019-11            8 珠海市
21   2019-12            8 珠海市
22   2020-01            8 珠海市
23   2020-02            8 珠海市
24   2020-03            8 珠海市
25   2020-04            8 珠海市
26   2020-05            8 珠海市
27   2020-06            8 珠海市
28   2020-07            8 珠海市
29   2020-08            8 珠海市
30   2020-09            8 珠海市
31   2019-07           33 汕头市
32   2019-08           33 汕头市
33   2019-09           33 汕头市
34   2019-10           33 汕头市
35   2019-11           33 汕头市
36   2019-12           33 汕头市
37   2020-01           33 汕头市
38   2020-02           33 汕头市
39   2020-03