in R code please

[25 marks] Load the iris data (included by default in R). (a) [10 marks) Write magrittr and/or dplyr code using pipe operators which: takes only those observations with sepal width > 3, computes the petal area as 7 X (petal length/2) x (petal width/2). (Hint: see ?pi) aggregates to produce the median petal area across species. (Hint: you may use either i) the functions subset, within, multiply_by, and aggregate, or ii) a combination of filter, mutate, multiply_by, group_by, and summarise)