R/whole_ssmdzprime_output.R

In JARS3N/preciseR: data processing and calculations

whole_ssmdzprime_output <-
function(xl,rate="PER",m1=3,m2=8){
  # assumes the data is from a GRA assay,using PER Rate and measuremnet 3 & 8
  filter(xl,Measurement %in% c(m1,m2),Group!="Background") %>% 
    select(Measurement,Rate=matches(rate)) %>% 
    mutate(Measurement=as.numeric(factor(Measurement))) %>% 
    mutate(condition=c("Pos","Neg")[Measurement]) %>% 
    group_by(condition) %>% 
    summarise(mu=mean(Rate),
              sd=sd(Rate),
              n=n()
    ) %>% 
    ssmd_zprime_summary(.,rate)
}