Simple function that counts all values in a matrix which are NA
any kind of matrix
Returns a data frame with the number of NA entries per column. Since the function is mainly written for the usage in batch effect correction of DNA methylation data, the column names of the data frame are set to "sample" and "num_pred_values". Nevertheless, the function can be used with any other matrix containing anything but beta values.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
## Shortly running example. For a more realistic example that takes ## some more time, run the same procedure with the full BEclearData ## dataset. ## Whole procedure that has to be done to use this function. data(BEclearData) ex.data <- ex.data[31:90,7:26] ex.samples <- ex.samples[7:26,] # Calculate median difference values and p-values meds <- calcMedians(data=ex.data, samples=ex.samples) pvals <- calcPvalues(data=ex.data, samples=ex.samples) # Summarize p-values and median differences for batch affected genes sum <- calcSummary(medians=meds, pvalues=pvals) clearedMatrix <- clearBEgenes(data=ex.data, samples=ex.samples, summary=sum) numberOfEntries <- countValuesToPredict(data=clearedMatrix)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.