saezlab / bradiplus Goto Github PK

For selectSamplesPeaks, the casemetric="AUC" is not working. It does not throw an error, but instead returns peaks which have AUC of 0 in every channel...

The problem is that the AUC-metric takes as time points not the time points of the peak, but time points 1 and 2 from the complete sample. Those time points are very close together, so everything is multiplied by 0:

# select corresponding data and compute metric
if (metric=="AUC"){  # AUC is computed as the sum of data points multiplied by the sampling time
    metric<-function(x){sum(x)*(data$time[2]-data$time[1])}
}else{
    metric <- match.fun(metric)
}

Instead of data$time[x], it should be the start and end time of the peak.