Question: soGGi package: Error in .validate_names(rownames, ans_rownames, "assay rownames()"
1
gravatar for nicolas.descostes
7 weeks ago by
nicolas.descostes20 wrote:

Hi,

I am following this tutorial: https://rockefelleruniversity.github.io/RUATACWorkshop.html

I tried the most recent version of soGGi but also the one of the tutorial (1.10.4).

When calculating:

nucFree <- regionPlot(bamFile = sortedBAM, testRanges = TSSs, style = "point", 
    format = "bam", paired = TRUE, minFragmentLength = 0, maxFragmentLength = 100, 
    forceFragment = 50)

I get the error:

Error in .validate_names(rownames, ans_rownames, "assay rownames()", "rowData rownames() / rowRanges names()") : 
  assay rownames() must be NULL or identical to rowData rownames() /
  rowRanges names()

I found that it has to do with "SummarizedExperiment": https://support.bioconductor.org/p/83850/

Can that be corrected in the package?

Thx

soggi • 95 views
ADD COMMENTlink modified 7 weeks ago • written 7 weeks ago by nicolas.descostes20
Answer: soGGi package: Error in .validate_names(rownames, ans_rownames, "assay rownames(
1
gravatar for nicolas.descostes
7 weeks ago by
nicolas.descostes20 wrote:

Adding this line to the code solved the issue:

profileMat <- RegionsMat
    colnames(profileMat) <- c(paste0("Point_Centre",seq(0-distanceUpStart,-1)),"Point_Centre",paste0("Point_Centre",seq(1,distanceDownEnd)))
    filteredRanges <- c(RangesPos,RangesNeg)

  ------->  ****names(filteredRanges) <- c(as.character(RangesPos$giID), as.character(RangesNeg$giID))****

    profileSample <- SummarizedExperiment(profileMat,rowRanges=filteredRanges[match(rownames(profileMat),filteredRanges$giID)])
ADD COMMENTlink written 7 weeks ago by nicolas.descostes20

hi Nicolas,

Thank you for bug report (and adding a fix!)

I will update this in release and development tomorrow.

Thanks again,

tom

ADD REPLYlink written 7 weeks ago by Thomas Carroll400

@Nicolas, Thanks for the bug report and the fix. @Thomas, I submitted a pull request containing Nicolas' fix.

Best wishes

ADD REPLYlink written 6 weeks ago by Firas0

Hello,

I am still getting the same error when I run the regionPlot function. Is there a way to go about fixing this?

Best, Sarah

ADD REPLYlink written 5 weeks ago by minkyung.lee.gr0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 123 users visited in the last hour