Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error in rhdf5::h5ls(fname) : HDF5. Object header. Can't open object. #253

Open
mbartl13 opened this issue Oct 23, 2020 · 0 comments
Open

Comments

@mbartl13
Copy link

Hi all- having an issue getting the sleuth object set:

Upload metadata

metadata <- read.delim("/dcl02/leased/giffin/SINV_RNAseq/fastq/metadata.txt")
#metadata <- read.table(dcl02/leased/giffin/SINV_RNAseq/fastq/'metadata.txt', sep='\t', header=TRUE, stringsAsFactors = FALSE)
metadata <- dplyr::select(metadata, c('SampleName', 'Time', 'Condition', 'Replicate'))

#rename SampleName to sample
metadata <- dplyr::rename(metadata, sample = SampleName)
head(metadata)

#set path to kallisto outputs
metadata <- dplyr::mutate(metadata,
path = file.path('.', 'trim', 'trimmed', 'kallisto_out', sample, 'abundance.h5'))

Associating transcripts to genes

mart <- useMart(biomart = "ENSEMBL_MART_ENSEMBL",
dataset = "rnorvegicus_gene_ensembl",
host = "dec2015.archive.ensembl.org")
ttg <- getBM(
attributes = c("ensembl_transcript_id", "transcript_version",
"ensembl_gene_id", "external_gene_name", "description",
"transcript_biotype"),
mart = mart)
ttg <- dplyr::rename(ttg, target_id = ensembl_transcript_id,
ens_gene = ensembl_gene_id, ext_gene = external_gene_name)
ttg <- dplyr::select(ttg, c('target_id', 'ens_gene', 'ext_gene'))
head(ttg)

#create sleuth object
so <- sleuth_prep(metadata, target_mapping = ttg,
aggregation_column = 'ens_gene', extra_bootstrap_summary = TRUE)
dropping unused factor levels
....................Error in rhdf5::h5ls(fname) : HDF5. Object header. Can't open object.

Thanks in advance!
Maggie

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant