Description Usage Arguments Details Value Note Examples

View source: R/getQualityScore.R

This function provides different types of measurements of uncertainty after CNV calling

1 2 3 4 5 | ```
getQualityScore(x, ...)
## Default S3 method:
getQualityScore(x, sds, w, type, iter = 10000, threshold = 0.1, ...)
## S3 method for class 'cnv'
getQualityScore(x, type = "class", iter = 10000, threshold = 0.1, ...)
``` |

`x` |
and object of class |

`...` |
further arguments passed to or from |

`type` |
the type of quality score measurement computed. Possible values are "class", "CNVtools" or "CANARY" (see Details) |

`iter` |
number of iterations when |

`threshold` |
a value to compute the proportion of sample individuals with confidence score bigger than it (see Details) |

`sds` |
standard deviations vector of intensity signal for each copy number status |

`w` |
copy number status proportions vector |

The quality scores measures how well the clusters are separated.
If compares the locations of the means with the standard error for each pair of adjacent cluster.
Obviously, except for probability of good classification (`type="class")`

,
the lower quality score the highest uncertainty.
There are 3 possible types of quality score measurements:
"class": probability of good classification),
"CNVtools": the score defined in 'CNVtools' package) and
"CANARY": proportions of sample individuals with confidence score bigger than `threshold`

. The confidence
score is defined as the ratio between the second biggest copy number call probability divided by
the biggest one.

An object of class `getQualityScore`

with a single number of quality score.

For `cnv`

objects created directly from probabilities and not from fitting a univariate intensity signal,
only "class" quality score type can be calculated.

1 2 3 4 5 | ```
data(dataMLPA)
CNV<-cnv(x = dataMLPA$Gene2, threshold.0 = 0.01, mix.method = "mixdist")
getQualityScore(CNV,type="class")
getQualityScore(CNV,type="CNVtools")
getQualityScore(CNV,type="CANARY")
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.

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.