Hi, Great tools.
I got stuck at the merge_fasta step, where I used data from your tutorial: https://jiang-junyao.github.io/IReNA/bulk-ATAC+scRNA#integrate-bulk-atac-seq-and-scrna-seq-data-to-reconstruct-regulatory-networks
filtered_footprints <- read.table(file.path(analysis_dir, "ATAC-seq", 'footprints.bed'),sep = '\t')
fastadir <- file.path(analysis_dir, "Genome")
merged_fasta <- get_merged_fasta(filtered_footprints, fastadir)
Warning in merge_footprints(footprints, ditance = distance) :
NAs introduced by coercion
Error in pva1 + pva2 : non-numeric argument to binary operator
###merge footprints whose distance is less than 4
filtered_footprints <- read.table(file.path(analysis_dir, "ATAC-seq", 'footprints.bed'),sep = '\t')
fastadir <- file.path(analysis_dir, "Genome/hg38.fa")
merged_fasta <- get_merged_fasta(filtered_footprints, fastadir)
Error in validInput(fastadir, "fastadir", "direxists") :
Input value for 'fastadir' is not a direxists, (fastadir = character) please supply valid input!
Could you help me out with this issue?
Thank you!
Hi, Great tools.
I got stuck at the merge_fasta step, where I used data from your tutorial: https://jiang-junyao.github.io/IReNA/bulk-ATAC+scRNA#integrate-bulk-atac-seq-and-scrna-seq-data-to-reconstruct-regulatory-networks
filtered_footprints <- read.table(file.path(analysis_dir, "ATAC-seq", 'footprints.bed'),sep = '\t')
Could you help me out with this issue?
Thank you!