mbic.mp: Computes the modified BIC for a given multi-platform...

In mpcbs: Multi-sample and Multiplatform CNV.

Description

Computes the modified BIC for a given multi-platform segmentation.

Usage

mbic.mp(S, imap, sigma, tau, rratio)

Arguments

S	cumulative sums.
imap	imap object from an anchor set.
sigma	platform specific variances.
tau	a vector of ascending change-points.
rratio	response ratio for each platform.

Value

mbic	mbic
term1	term 1 of the modified BIC
term2	term 2 of the modified BIC
term3	term 3 of the modified BIC

Author(s)

Nancy R. Zhang

References

Zhang, NR, Senbabaoglu, Y. and Li, J.Z. (2009) Joint Estimation of DNA Copy Number from Multiple Platforms. Under review, download manuscript from http://www-stat.stanford.edu/~nzhang/web_multiplatform/

See Also

mpcbs, plot.crossplatform, merge.pos

Examples

data(mpcbs.example)

# there are 3 platforms represented in this data example: Illumina, Affymetrix, and Agilent.
names(mpcbs.example)

# K is the number of platforms.
K=3

# Store the chromosome positions in vector pos,
# the intensities in vector y:

pos=vector("list",K)
pos[[1]] = mpcbs.example$illu[,1]
pos[[2]] = mpcbs.example$affy[,1]
pos[[3]] = mpcbs.example$agil[,1]

y = vector("list",K)
y[[1]] = mpcbs.example$illu[,2]
y[[2]] = mpcbs.example$affy[,2]
y[[3]] = mpcbs.example$agil[,2]

# Names of the platforms:
platform.names=c("Illumina","Affymetrix","Agilent")

mpcbs documentation built on May 2, 2019, 4:49 p.m.