Plot correlation among DNAm clockx
Arguments
- x
a tible or data.frame with the different DNAm clocks
- ...
other arguments to be passs through function 'chart.Correlation' from 'PerformanceAnalytics' package
Examples
dd <- GEOquery::getGEO("GSE109446")
#> Found 1 file(s)
#> GSE109446_series_matrix.txt.gz
#> Error in curl::curl_download(url, destfile, mode = mode, quiet = quiet, handle = h): Timeout was reached: [] Operation timed out after 120007 milliseconds with 10944512 out of 142239745 bytes received
gse109446 <- dd[[1]]
#> Error in eval(expr, envir, enclos): object 'dd' not found
controls <- pData(gse109446)$`diagnosis:ch1` == "control"
#> Error in pData(gse109446): could not find function "pData"
gse <- gse109446[, controls]
#> Error in eval(expr, envir, enclos): object 'gse109446' not found
age <- as.numeric(pData(gse)$`age:ch1`)
#> Error in pData(gse): could not find function "pData"
age.gse <- DNAmAge(gse, age = age)
#> Error in DNAmAge(gse, age = age): object 'gse' not found
plotCorClocks(age.gse)
#> Error in is.data.frame(x): object 'age.gse' not found