THYMmethy.m<-as.matrix(THYMmethy) #
samlename<-THYMmethy.m[,1] ##添加行名
THYMmethy.mn<-THYMmethy.m[,-1] ##添加行名
#查看数据结构
head(THYMmethy.mn)
class(THYMmethy.mn[,1])
THYMmethy.mn[1:6,1:5]
THYMmethy.mnn=apply(THYMmethy.mn,2,as.numeric)
THYMmethy.mnn[1:6,1:5]
rownames(THYMmethy.mnn) <- samlename
THYMmethyout.l <- epidish(THYMmethy.mnn, centEpiFibIC.m, method = 'RPC')
frac.m <- THYMmethyout.l$estF
frac.m
dim(frac.m)
pheno.v <- rep(c(0, 1), each = 63)
pheno.v
celldmc.o <- CellDMC(THYMmethy.mnn, pheno.v, frac.m)
dotchart(celldmc.o)
