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

Share:

Description

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

Usage

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
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")