Description Usage Arguments Details Value Author(s) Examples
Uses maximum likelihood to estimate (C,5mC,5hmC) vectors for one specimen
1 | fitOxBS(betaBS, betaOxBS, signalBS, signalOxBS, eps=1E-5)
|
betaBS |
beta value from conventional bisulfite conversion |
betaOxBS |
beta value from oxy-bisulfite conversion |
signalBS |
total signal from conventional bisulfite conversion |
signalOxBS |
total signal from oxy-bisulfite conversion |
eps |
small positive value representing numerical zero |
Uses maximum likelihood to estimate (C,5mC,5hmC) one specimen (many CpGs).
matrix of (C,5mC,5hmC) values (each row corresponds to a separate CpG).
E. Andres Houseman
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | ## Not run:
data(OxyBSSampleData)
nSpecimens <- 30
nCpGs <- 30
# Calculate Total Signals
signalBS <- exampleMethBS+exampleUnmethBS
signalOxBS <- exampleMethOxBS+exampleUnmethOxBS
# Calculate Beta Values
betaBS <- exampleMethBS/signalBS
betaOxBS <- exampleMethOxBS/signalOxBS
# Create container for results
MethOxy <- array(NA,dim=c(nCpGs,nSpecimens,3))
dimnames(MethOxy) <- list(
rownames(exampleMethBS)[1:nCpGs],
colnames(exampleMethBS)[1:nSpecimens],
c("C","5mC","5hmC"))
# Process results (one array at a time)
for(i in 1:nSpecimens){
MethOxy[,i,] <- fitOxBS(betaBS[,i],betaOxBS[,i],signalBS[,i],signalOxBS[,i])
}
# Check that results sum to one
table(apply(MethOxy,1:2,sum))
# First specimen
MethOxy[,1,]
# Ranges
range(MethOxy[,,1])
range(MethOxy[,,2])
range(MethOxy[,,3])
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.