getPloidy: Genomewide average ploidy

Description Usage Arguments Value Author(s) Examples

Description

This function returns the genomewide average ploidy of aneuploid fraction and overall of the sample.

Usage

1
getPloidy(p, type = 1, sam.dat, oo)

Arguments

p

AGP value. see getGrid() for details.

type

integer, ploidy indicator. 1, diploid; 2, tetraploid; 3, hexaploid

sam.dat

numeric matrix, as returned from getSeg() or getSegChr()

oo

Origin Cluster, as returned from getOrigin()

Value

Ploidy of the overall sample and ploidy of the aneuploid proportion are returned.

Author(s)

Bo Li

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
data(A0SD.BAF)
data(A0SD.LRR)
seg.dat=c()
for(CHR in c(8,9,10)){
	baf=A0SD.BAF[A0SD.BAF[,2]==CHR,]
	lrr=A0SD.LRR[A0SD.LRR[,2]==CHR,]
	x=getSegChr(baf,lrr)
	seg.dat=rbind(seg.dat,x)
}
dd.dat=seg.dat[,2:8]
rownames(dd.dat)=seg.dat[,1]
mode(dd.dat)='numeric'
para=getPara()
oo=getOrigin(dd.dat,para=para)
getPloidy(0.5,sam.dat=dd.dat,oo=oo)


Search within the CHAT package
Search all R packages, documentation and source code

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

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

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