getquantilesandranks: Compute the quantiles and ranks for a given gdsn.node

Description Usage Arguments Details Value Author(s) Examples

View source: R/ranknorm.R

Description

Used inside dasenrank to generate the quantiles for both type 'I' and type 'II' probes to normalise DNA methylation data using bigmelon.

Usage

1
getquantilesandranks(gds, node, onetwo, rank.node = NULL, perc = 1)

Arguments

gds

A gds.class object

node

A gdsn.class object, or a character string that refers to a node within supplied gds.

onetwo

gdsn.class object OR character string that refers to location in gds that contains information relating to probe design OR vector of length equal to the number of rows in the array that contains 'I' and 'II' in accordance to Illumina HumanMethylation micro-array design. Can be obtained with fot(gds)

rank.node

Default = NULL. If supplied with character string, function will calculate the ranks of given node and store them in gds. Additionally, the computed quantiles will now instead be attributed to rank.node which can be accessed with get.attr.gdsn

perc

A number between 0 and 1 that relates to the given proportion of columns that are used to normalise the data. Default is set to 1, but in cases where there many of samples to normalise this number can be reduced to increase speed of code.

Details

Used in dasenrank, can be used externally for testing purposes.

Value

If rank.node is NULL. A list containing quantiles, intervals and supplied probe design will be returned. If rank.node was supplied, nothing will be returned. Instead a new node will be created in given gds that has the otherwise returned list attached as an attribute. Which can be accessed with get.attr.gdsn

Author(s)

Tyler J. Gorrie-Stone Who to contact: <[email protected]>

Examples

1
2
3
4
5
6
7
8
data(melon)
e <- es2gds(melon, "melon.gds")
output <- getquantilesandranks(gds = e, 'methylated', onetwo = fot(e), perc = 1, rank.node = NULL)
str(output)
# with-out put.
getquantilesandranks(gds = e, 'methylated', onetwo = fot(e), perc = 1, rank.node = 'mnsrank')
closefn.gds(e)
unlink("melon.gds")

Bioconductor-mirror/bigmelon documentation built on June 1, 2017, 4:26 a.m.