new pat_timeAverageDiagnostics() function about airsensor HOT 2 CLOSED

I disagree. This function is all about applying statistical functions to consecutive chunks of the overall time series and returning vector that represents the value of that statistic calculated for each chunk. Currently, the statical functions you have enabled are mean, sd and count, etc.

The idea for the t-test is to have another function that would behave similar to sd in that it takes the data within a chunk of time, applies an algorithm and returns a number. So you will need to implement another internal function that looks like this:

if ( stats == "ttest_qc" ) {
func <- function(x) {
# use t.test with pm25_A and pm25_B and appropriate parameters
# return the t test statistic or p value as the result
}

The use case is to generate clean, hourly pm25 data with the following steps:

• remove outliers from raw data
• pm25 <- calculate AB mean
• calculate t_statistic
• badQCIndices <- t_statistic > someThreshold
• pm25[baQCIndices] <- NA

from airsensor.

I think that any sort of statistical analysis should exist independently of the this function. Perhaps a pat_ttest() or in a pat_statistics()?

from airsensor.