View source: R/filter_params.R

filter_params | R Documentation |

Determine the A quantile cutoff to be used

filter_params( hic.table, SD = 2, numChanges = 300, FC = 3, alpha = 0.05, Plot = FALSE )

`hic.table` |
A hic.table object |

`SD` |
The standard deviation of the fuzzing used to produce a Hi-C matrix from your data with few true differences. |

`numChanges` |
The number of changes to add into the Hi-C matrix created. This should be proportional to the resolution of the data. High resolution data should use more changes i.e. 1MB resolution - 300 changes, 100KB resolution - 1000 changes, etc. |

`FC` |
The fold change of the changes added to the Hi-C matrix. |

`alpha` |
The alpha level for hypothesis testing. |

`Plot` |
logical, should MD plots for the normalization and difference detection be plotted? |

This function will take your data and produce an additional Hi-C matrix using the IF1 vector. Random normal noise will be added to the vector to create a "fuzzed" matrix with few true differences. Then the specified number of true changes will be added at the specified fold change level to the matrices. The HiCcompare procedure is run on the data and a plot of the MCC, TPR, and FPR based on the A minimum value filtered out will be produced. This is to aid you in determining what value you should use when analyzing your data with the hic_compare() function.

A plot of the Mathews Correlation Coefficient (MCC), true positive rate (TPR), and false positive rate (FPR) over the A minimum value filtered.

data('HMEC.chr22') data('NHEK.chr22') hic.table <- create.hic.table(HMEC.chr22, NHEK.chr22, chr = 'chr22') filter_params(hic.table)

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.