If you assign your call to the qic function to a variable, it will return a plot, and a dataframe with all the relevant info you need.
E.G.
p <- qic(....) returns a list with a plot and dataframe
my_df <- p$data (extract the data with all the calculations)
my_df >%> select (x, y, sigma.signal) should give you the date, y value and any out of control points.
On my phone at the moment, so can't give precise info, but this should work..
I often take the results from the qic function and build my own plots in ggplot2 using this approach
Let me know how you get on
John