Question: mageckflute can't output my own data but the model data
0
4 weeks ago by
pamelalee0 wrote:

When I run the MAGeCKFlute, the output graph and data are not based on my own data but the model data. For Example:

> library('MAGeCKFlute')

Registered S3 method overwritten by 'enrichplot':
method               from
fortify.enrichResult DOSE
##############################################################################
Public License version 3 (GPLv3). Details of GPLv3 is available at
http://www.gnu.org/licenses/gpl-3.0.html. Particullary, users are required to
formally cite the original Pathview paper (not just mention it) in publications
or products. For details, do citation("pathview") within R.

##############################################################################
MAGeCKFlute is an open source software package. Particullary, users are
required to formally cite the original MAGeCKFlute paper in publications or
products. For details, do citation("MAGeCKFlute") within R.
> View(countsummary)
1   File    Label   Reads   Mapped  Percentage  TotalsgRNAs Zerocounts  GiniIndex   NegSelQC    NegSelQCPval    NegSelQCPvalPermutation NegSelQCPvalPermutationFDR  NegSelQCGene
2   G20.bam G20 5214605 3602121 0.6908  11107   56  0.07057 0   1   1   1   0.0
3   G1.bam  G1  4272051 3232724 0.7567  11107   5   0.05494 0   1   1   1   0.0
Showing 1 to 3 of 3 entries, 13 total columns
##Data above are my Data created from mageck
##then I started to run MAGeckFlute
> data("countsummary")
1 ../data/GSC_0131_Day23_Rep1.fastq.gz day23_r1 62818064 39992777     0.6366
2  ../data/GSC_0131_Day0_Rep2.fastq.gz  day0_r2 47289074 31709075     0.6705
3  ../data/GSC_0131_Day0_Rep1.fastq.gz  day0_r1 51190401 34729858     0.6784
4 ../data/GSC_0131_Day23_Rep2.fastq.gz day23_r2 58686580 37836392     0.6447
TotalsgRNAs Zerocounts GiniIndex NegSelQC NegSelQCPval NegSelQCPvalPermutation
1       64076         57   0.08510        0            1                       1
2       64076         17   0.07496        0            1                       1
3       64076         14   0.07335        0            1                       1
4       64076         51   0.08587        0            1                       1
NegSelQCPvalPermutationFDR NegSelQCGene
1                          1            0
2                          1            0
3                          1            0
4                          1            0


the table is not my table any more and the same problem occurred when I run FluteRRA, the output graphs are not based on my data.

mageckflute • 84 views
written 4 weeks ago by pamelalee0
Answer: C: mageckflute can't output my own data but the model data
1
4 weeks ago by
Mike Smith4.0k
EMBL Heidelberg / de.NBI
Mike Smith4.0k wrote:

When you run the command data("countsummary") it loads the example data from the MAGeCKFlute package, and this overwrites your existing countsummary table.

It's common to see this used in package vignettes as it allows the code to be run on some example dataset, but you don't have to run that step if you have your own data to provide instead.