我正在尝试使用getAgeR
中的cgageR
函数。我首先尝试重现描述中提供的示例,但是我一直都在收到相同的错误消息。有没有人对getAgeR
有过经验并且能够运行它?
library("cgageR")
library("wateRmelon")
data(melon)
bval <- betas(melon)
est.ages <- getAgeR(bval,epitoc=TRUE,horvath=TRUE,hannum=TRUE,drift=FALSE,showStatusHannum=TRUE,
keepcpgs.epitoc=TRUE,keepcpgs.hannum=TRUE,keepres=FALSE,chrage=NULL)
我得到:
获取epiTOC年龄估算...
可用的epiTOC CpG数量:2/385。
做完了!继续...
获取Horvath年龄估算值...
使用Horvaths时钟未检测到系数
做完了!继续...
正在初始化Hannum Clock年龄估算...
df [cgi,i]中的错误:下标超出范围