请教各位哥问题，只有1个论坛币了~谢谢大家

各位老师，请教下如何按月份求以下数据的sd(p)/mean(p)
date                        ym           a         b         c         d
2017/1/1        2017年1月        200        33        44        160
2017/1/2        2017年1月        201        34        44        160
2017/1/3        2017年1月        202        35        44        160
2017/1/4        2017年1月        203        36        44        160
2017/2/1        2017年2月        204        37        44        160
2017/2/2        2017年2月        205        38        44        160
2017/2/3        2017年2月        206        39        44        160
2017/2/4        2017年2月        207        40        44        160
2017/3/1        2017年3月        208        41        44        160
2017/3/2        2017年3月        209        42        44        160
2017/3/3        2017年3月        210        43        44        160
2017/3/4        2017年3月        211        44        44        160
2017/4/1        2017年4月        212        45        44        160
2017/4/2        2017年4月        213        46        44        160
2017/4/3        2017年4月        214        47        44        160
2017/4/4        2017年4月        215        48        44        160


结果是：
    ym        variable                   p_rate
2017年1月        a             sd（当月p值)/mean(当月p值)
2017年2月        a             sd（当月p值)/mean(当月p值)
2017年3月        a             sd（当月p值)/mean(当月p值)
2017年4月        a             sd（当月p值)/mean(当月p值)
2017年1月        b             sd（当月p值)/mean(当月p值)
2017年2月        b             sd（当月p值)/mean(当月p值)
2017年3月        b             sd（当月p值)/mean(当月p值)
2017年4月        b             sd（当月p值)/mean(当月p值)
2017年1月        c             sd（当月p值)/mean(当月p值)
2017年2月        c             sd（当月p值)/mean(当月p值)
2017年3月        c             sd（当月p值)/mean(当月p值)
2017年4月        c             sd（当月p值)/mean(当月p值)




恳求各位老师帮助，谢谢了。

library(tidyverse)
data %>%
    gather(variable, p_value, a:d) %>%
    group_by(ym, variable) %>%
    summarise(p_rate = sd(p_value) / mean(p_value))