cummeRbund package function diffData returns no rows for TSS features

0

Entering edit mode

Asta Laiho ▴ 120

@asta-laiho-4271

Last seen 8.2 years ago

Finland

Many thanks for a quick response! I installed the version 1.1.3 and tried rebuilding the database but got the following error: cuff <- readCufflinks(rebuild=T) Creating database /Users/alaiho/Documents/Data-analysis/Projects/NGS/1 10007_Liisa_Nissinen/cufflinks/cuffdiff_output/cuffData.db Reading /Users/alaiho/Documents/Data-analysis/Projects/NGS/110007_Liis a_Nissinen/cufflinks/cuffdiff_output/genes.fpkm_tracking Checking samples table... Populating samples table... Writing genes table Reshaping geneData table Recasting Error in as.data.frame(cast(genemelt, ... ~ measurement)) : could not find function "cast" > sessionInfo() R version 2.14.2 (2012-02-29) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] C/UTF-8/C/C/C/C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] cummeRbund_1.1.3 reshape2_1.2.1 ggplot2_0.9.0 RSQLite_0.11.1 [5] DBI_0.2-5 loaded via a namespace (and not attached): [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 [9] plyr_1.7.1 proto_0.3-9.2 scales_0.2.0 stringr_0.6 On Mar 15, 2012, at 5:28 AM, Loyal Goff wrote: > Hi Asta, > I believe that I had already fixed this bug in the newer releases of cummeRbund. Can you try to upgrade to the development cummeRbund version (1.1.3)? You will have to rebuild your cuffData.db database with this new version. Please let me know if this does not resolve your issue and I will try to work with you to get it resolved. > > Thanks! > > Cheers, > Loyal > > On Mar 14, 2012, at 9:05 AM, Asta Laiho [guest] wrote: > >> >> Hi, >> >> I am using cummeRbund package to analyze human RNAseq data. Everything else is working very smoothly, except that the function diffData(TSS(cuff)) returns zero rows which to me seems like an error. Other features; genes, isoforms and CDSs are returned correctly. What could be wrong? >> >> This is what I'm doing: >> cuff <- readCufflinks() >> cuff >> CuffSet instance with: >> 4 samples >> 22260 genes >> 37557 isoforms >> 27776 TSS >> 25411 CDS >> 133002 promoters >> 166656 splicing >> 101436 relCDS >> >> As you can see, the resulting data frame has no rows: >> diffData(TSS(cuff)) >> [1] TSS_group_id sample_1 sample_2 status value_1 >> [6] value_2 ln_fold_change test_stat p_value q_value >> [11] significant >> <0 rows> (or 0-length row.names) >> >> The following works correctly: >> head(> head(diffData(CDS(cuff)) ) >> CDS_id sample_1 sample_2 status value_1 value_2 ln_fold_change >> 1 P1 S1 S2 OK 11.30910 19.39660 7.78320e-01 >> 2 P10 S1 S2 NOTEST 3.25391 0.00000 -1.79769e+308 >> 3 P100 S1 S2 NOTEST 2.26812 2.28999 1.38405e-02 >> 4 P1000 S1 S2 NOTEST 0.00000 0.00000 0.00000e+00 >> 5 P10000 S1 S2 LOWDATA 12.03700 17.89210 5.71841e-01 >> 6 P10001 S1 S2 NOTEST 6.83010 4.39573 -6.35805e-01 >> test_stat p_value q_value significant >> 1 -1.27958e+00 0.200693000 0.607328 no >> 2 -1.79769e+308 0.000624649 1.000000 no >> 3 -2.06580e-02 0.983519000 1.000000 no >> 4 0.00000e+00 1.000000000 1.000000 no >> 5 0.00000e+00 1.000000000 1.000000 no >> 6 9.72050e-01 0.331026000 1.000000 no >> >> >> >> >> >> -- output of sessionInfo(): >> >> R version 2.14.2 (2012-02-29) >> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >> >> locale: >> [1] C/UTF-8/C/C/C/C >> >> attached base packages: >> [1] stats graphics grDevices utils datasets methods base >> >> other attached packages: >> [1] cummeRbund_1.0.0 ggplot2_0.9.0 reshape_0.8.4 plyr_1.7.1 >> [5] RSQLite_0.11.1 DBI_0.2-5 >> >> loaded via a namespace (and not attached): >> [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 >> [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 >> [9] proto_0.3-9.2 reshape2_1.2.1 scales_0.2.0 stringr_0.6 >> >> -- >> Sent via the guest posting facility at bioconductor.org. >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >

RNASeq cummeRbund RNASeq cummeRbund • 1.7k views

ADD COMMENT • link updated 12.1 years ago by Loyal Goff ▴ 140 • written 12.1 years ago by Asta Laiho ▴ 120

0

Entering edit mode

Loyal Goff ▴ 140

@loyal-goff-4989

Last seen 9.6 years ago

Hi Asta, This is a bug that has been fixed in future versions, the current workaround is to directly import the reshape library before running readCufflinks, as follows: >library(cummeRbund) >library(reshape) >cuff<-readCufflinks(rebuild=T) Cheers, Loyal ------------------ Loyal A. Goff Postdoctoral Fellow CSAIL, MIT & SCRB, Harvard On Mar 15, 2012, at 4:12 AM, Asta Laiho <asta.laiho at="" btk.fi=""> wrote: > Many thanks for a quick response! > > I installed the version 1.1.3 and tried rebuilding the database but got the following error: > > cuff <- readCufflinks(rebuild=T) > Creating database /Users/alaiho/Documents/Data-analysis/Projects/NGS /110007_Liisa_Nissinen/cufflinks/cuffdiff_output/cuffData.db > Reading /Users/alaiho/Documents/Data-analysis/Projects/NGS/110007_Li isa_Nissinen/cufflinks/cuffdiff_output/genes.fpkm_tracking > Checking samples table... > Populating samples table... > Writing genes table > Reshaping geneData table > Recasting > Error in as.data.frame(cast(genemelt, ... ~ measurement)) : > could not find function "cast" > >> sessionInfo() > R version 2.14.2 (2012-02-29) > Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) > > locale: > [1] C/UTF-8/C/C/C/C > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] cummeRbund_1.1.3 reshape2_1.2.1 ggplot2_0.9.0 RSQLite_0.11.1 > [5] DBI_0.2-5 > > loaded via a namespace (and not attached): > [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 > [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 > [9] plyr_1.7.1 proto_0.3-9.2 scales_0.2.0 stringr_0.6 > > On Mar 15, 2012, at 5:28 AM, Loyal Goff wrote: > >> Hi Asta, >> I believe that I had already fixed this bug in the newer releases of cummeRbund. Can you try to upgrade to the development cummeRbund version (1.1.3)? You will have to rebuild your cuffData.db database with this new version. Please let me know if this does not resolve your issue and I will try to work with you to get it resolved. >> >> Thanks! >> >> Cheers, >> Loyal >> >> On Mar 14, 2012, at 9:05 AM, Asta Laiho [guest] wrote: >> >>> >>> Hi, >>> >>> I am using cummeRbund package to analyze human RNAseq data. Everything else is working very smoothly, except that the function diffData(TSS(cuff)) returns zero rows which to me seems like an error. Other features; genes, isoforms and CDSs are returned correctly. What could be wrong? >>> >>> This is what I'm doing: >>> cuff <- readCufflinks() >>> cuff >>> CuffSet instance with: >>> 4 samples >>> 22260 genes >>> 37557 isoforms >>> 27776 TSS >>> 25411 CDS >>> 133002 promoters >>> 166656 splicing >>> 101436 relCDS >>> >>> As you can see, the resulting data frame has no rows: >>> diffData(TSS(cuff)) >>> [1] TSS_group_id sample_1 sample_2 status value_1 >>> [6] value_2 ln_fold_change test_stat p_value q_value >>> [11] significant >>> <0 rows> (or 0-length row.names) >>> >>> The following works correctly: >>> head(> head(diffData(CDS(cuff)) ) >>> CDS_id sample_1 sample_2 status value_1 value_2 ln_fold_change >>> 1 P1 S1 S2 OK 11.30910 19.39660 7.78320e-01 >>> 2 P10 S1 S2 NOTEST 3.25391 0.00000 -1.79769e+308 >>> 3 P100 S1 S2 NOTEST 2.26812 2.28999 1.38405e-02 >>> 4 P1000 S1 S2 NOTEST 0.00000 0.00000 0.00000e+00 >>> 5 P10000 S1 S2 LOWDATA 12.03700 17.89210 5.71841e-01 >>> 6 P10001 S1 S2 NOTEST 6.83010 4.39573 -6.35805e-01 >>> test_stat p_value q_value significant >>> 1 -1.27958e+00 0.200693000 0.607328 no >>> 2 -1.79769e+308 0.000624649 1.000000 no >>> 3 -2.06580e-02 0.983519000 1.000000 no >>> 4 0.00000e+00 1.000000000 1.000000 no >>> 5 0.00000e+00 1.000000000 1.000000 no >>> 6 9.72050e-01 0.331026000 1.000000 no >>> >>> >>> >>> >>> >>> -- output of sessionInfo(): >>> >>> R version 2.14.2 (2012-02-29) >>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >>> >>> locale: >>> [1] C/UTF-8/C/C/C/C >>> >>> attached base packages: >>> [1] stats graphics grDevices utils datasets methods base >>> >>> other attached packages: >>> [1] cummeRbund_1.0.0 ggplot2_0.9.0 reshape_0.8.4 plyr_1.7.1 >>> [5] RSQLite_0.11.1 DBI_0.2-5 >>> >>> loaded via a namespace (and not attached): >>> [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 >>> [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 >>> [9] proto_0.3-9.2 reshape2_1.2.1 scales_0.2.0 stringr_0.6 >>> >>> -- >>> Sent via the guest posting facility at bioconductor.org. >>> >>> _______________________________________________ >>> Bioconductor mailing list >>> Bioconductor at r-project.org >>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >> >

ADD COMMENT • link 12.1 years ago Loyal Goff ▴ 140

0

Entering edit mode

Hi, After importing the reshape the database is succesfully rebuilt: cuff CuffSet instance with: 4 samples 22260 genes 37557 isoforms 27776 TSS 25411 CDS 133002 promoters 166656 splicing 101436 relCDS However, the diffData for TSS is still empty: head(diffData(TSS(cuff))) [1] TSS_group_id sample_1 sample_2 status [5] value_1 value_2 log2_fold_change test_stat [9] p_value q_value significant <0 rows> (or 0-length row.names) other features still work: head(diffData(CDS(cuff))) CDS_id sample_1 sample_2 status value_1 value_2 log2_fold_change 1 P1 HaCat A5 OK 11.30910 19.39660 7.78320e-01 2 P10 HaCat A5 NOTEST 3.25391 0.00000 -1.79769e+308 3 P100 HaCat A5 NOTEST 2.26812 2.28999 1.38405e-02 4 P1000 HaCat A5 NOTEST 0.00000 0.00000 0.00000e+00 5 P10000 HaCat A5 LOWDATA 12.03700 17.89210 5.71841e-01 6 P10001 HaCat A5 NOTEST 6.83010 4.39573 -6.35805e-01 test_stat p_value q_value significant 1 -1.27958e+00 0.200693000 0.607328 no 2 -1.79769e+308 0.000624649 1.000000 no 3 -2.06580e-02 0.983519000 1.000000 no 4 0.00000e+00 1.000000000 1.000000 no 5 0.00000e+00 1.000000000 1.000000 no 6 9.72050e-01 0.331026000 1.000000 no On Mar 15, 2012, at 1:03 PM, Loyal A. Goff wrote: > Hi Asta, > This is a bug that has been fixed in future versions, the current workaround is to directly import the reshape library before running readCufflinks, as follows: > >> library(cummeRbund) >> library(reshape) >> cuff<-readCufflinks(rebuild=T) > > Cheers, > Loyal > > > ------------------ > Loyal A. Goff > Postdoctoral Fellow > CSAIL, MIT & > SCRB, Harvard > > On Mar 15, 2012, at 4:12 AM, Asta Laiho <asta.laiho at="" btk.fi=""> wrote: > >> Many thanks for a quick response! >> >> I installed the version 1.1.3 and tried rebuilding the database but got the following error: >> >> cuff <- readCufflinks(rebuild=T) >> Creating database /Users/alaiho/Documents/Data-analysis/Projects/NG S/110007_Liisa_Nissinen/cufflinks/cuffdiff_output/cuffData.db >> Reading /Users/alaiho/Documents/Data-analysis/Projects/NGS/110007_L iisa_Nissinen/cufflinks/cuffdiff_output/genes.fpkm_tracking >> Checking samples table... >> Populating samples table... >> Writing genes table >> Reshaping geneData table >> Recasting >> Error in as.data.frame(cast(genemelt, ... ~ measurement)) : >> could not find function "cast" >> >>> sessionInfo() >> R version 2.14.2 (2012-02-29) >> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >> >> locale: >> [1] C/UTF-8/C/C/C/C >> >> attached base packages: >> [1] stats graphics grDevices utils datasets methods base >> >> other attached packages: >> [1] cummeRbund_1.1.3 reshape2_1.2.1 ggplot2_0.9.0 RSQLite_0.11.1 >> [5] DBI_0.2-5 >> >> loaded via a namespace (and not attached): >> [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 >> [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 >> [9] plyr_1.7.1 proto_0.3-9.2 scales_0.2.0 stringr_0.6 >> >> On Mar 15, 2012, at 5:28 AM, Loyal Goff wrote: >> >>> Hi Asta, >>> I believe that I had already fixed this bug in the newer releases of cummeRbund. Can you try to upgrade to the development cummeRbund version (1.1.3)? You will have to rebuild your cuffData.db database with this new version. Please let me know if this does not resolve your issue and I will try to work with you to get it resolved. >>> >>> Thanks! >>> >>> Cheers, >>> Loyal >>> >>> On Mar 14, 2012, at 9:05 AM, Asta Laiho [guest] wrote: >>> >>>> >>>> Hi, >>>> >>>> I am using cummeRbund package to analyze human RNAseq data. Everything else is working very smoothly, except that the function diffData(TSS(cuff)) returns zero rows which to me seems like an error. Other features; genes, isoforms and CDSs are returned correctly. What could be wrong? >>>> >>>> This is what I'm doing: >>>> cuff <- readCufflinks() >>>> cuff >>>> CuffSet instance with: >>>> 4 samples >>>> 22260 genes >>>> 37557 isoforms >>>> 27776 TSS >>>> 25411 CDS >>>> 133002 promoters >>>> 166656 splicing >>>> 101436 relCDS >>>> >>>> As you can see, the resulting data frame has no rows: >>>> diffData(TSS(cuff)) >>>> [1] TSS_group_id sample_1 sample_2 status value_1 >>>> [6] value_2 ln_fold_change test_stat p_value q_value >>>> [11] significant >>>> <0 rows> (or 0-length row.names) >>>> >>>> The following works correctly: >>>> head(> head(diffData(CDS(cuff)) ) >>>> CDS_id sample_1 sample_2 status value_1 value_2 ln_fold_change >>>> 1 P1 S1 S2 OK 11.30910 19.39660 7.78320e-01 >>>> 2 P10 S1 S2 NOTEST 3.25391 0.00000 -1.79769e+308 >>>> 3 P100 S1 S2 NOTEST 2.26812 2.28999 1.38405e-02 >>>> 4 P1000 S1 S2 NOTEST 0.00000 0.00000 0.00000e+00 >>>> 5 P10000 S1 S2 LOWDATA 12.03700 17.89210 5.71841e-01 >>>> 6 P10001 S1 S2 NOTEST 6.83010 4.39573 -6.35805e-01 >>>> test_stat p_value q_value significant >>>> 1 -1.27958e+00 0.200693000 0.607328 no >>>> 2 -1.79769e+308 0.000624649 1.000000 no >>>> 3 -2.06580e-02 0.983519000 1.000000 no >>>> 4 0.00000e+00 1.000000000 1.000000 no >>>> 5 0.00000e+00 1.000000000 1.000000 no >>>> 6 9.72050e-01 0.331026000 1.000000 no >>>> >>>> >>>> >>>> >>>> >>>> -- output of sessionInfo(): >>>> >>>> R version 2.14.2 (2012-02-29) >>>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >>>> >>>> locale: >>>> [1] C/UTF-8/C/C/C/C >>>> >>>> attached base packages: >>>> [1] stats graphics grDevices utils datasets methods base >>>> >>>> other attached packages: >>>> [1] cummeRbund_1.0.0 ggplot2_0.9.0 reshape_0.8.4 plyr_1.7.1 >>>> [5] RSQLite_0.11.1 DBI_0.2-5 >>>> >>>> loaded via a namespace (and not attached): >>>> [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 >>>> [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 >>>> [9] proto_0.3-9.2 reshape2_1.2.1 scales_0.2.0 stringr_0.6 >>>> >>>> -- >>>> Sent via the guest posting facility at bioconductor.org. >>>> >>>> _______________________________________________ >>>> Bioconductor mailing list >>>> Bioconductor at r-project.org >>>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>>> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >>> >>

ADD REPLY • link 12.1 years ago Asta Laiho ▴ 120

0

Entering edit mode

Hmm, can you send me a copy of your cuffData.db file so that I can try to figure out exactly what's going on? Thanks, loyal ------------------ Loyal A. Goff Postdoctoral Fellow CSAIL, MIT & SCRB, Harvard On Mar 15, 2012, at 7:27 AM, Asta Laiho <asta.laiho at="" btk.fi=""> wrote: > Hi, > > After importing the reshape the database is succesfully rebuilt: > cuff > CuffSet instance with: > 4 samples > 22260 genes > 37557 isoforms > 27776 TSS > 25411 CDS > 133002 promoters > 166656 splicing > 101436 relCDS > > However, the diffData for TSS is still empty: > > head(diffData(TSS(cuff))) > [1] TSS_group_id sample_1 sample_2 status > [5] value_1 value_2 log2_fold_change test_stat > [9] p_value q_value significant > <0 rows> (or 0-length row.names) > > other features still work: > > head(diffData(CDS(cuff))) > CDS_id sample_1 sample_2 status value_1 value_2 log2_fold_change > 1 P1 HaCat A5 OK 11.30910 19.39660 7.78320e-01 > 2 P10 HaCat A5 NOTEST 3.25391 0.00000 -1.79769e+308 > 3 P100 HaCat A5 NOTEST 2.26812 2.28999 1.38405e-02 > 4 P1000 HaCat A5 NOTEST 0.00000 0.00000 0.00000e+00 > 5 P10000 HaCat A5 LOWDATA 12.03700 17.89210 5.71841e-01 > 6 P10001 HaCat A5 NOTEST 6.83010 4.39573 -6.35805e-01 > test_stat p_value q_value significant > 1 -1.27958e+00 0.200693000 0.607328 no > 2 -1.79769e+308 0.000624649 1.000000 no > 3 -2.06580e-02 0.983519000 1.000000 no > 4 0.00000e+00 1.000000000 1.000000 no > 5 0.00000e+00 1.000000000 1.000000 no > 6 9.72050e-01 0.331026000 1.000000 no > > > On Mar 15, 2012, at 1:03 PM, Loyal A. Goff wrote: > >> Hi Asta, >> This is a bug that has been fixed in future versions, the current workaround is to directly import the reshape library before running readCufflinks, as follows: >> >>> library(cummeRbund) >>> library(reshape) >>> cuff<-readCufflinks(rebuild=T) >> >> Cheers, >> Loyal >> >> >> ------------------ >> Loyal A. Goff >> Postdoctoral Fellow >> CSAIL, MIT & >> SCRB, Harvard >> >> On Mar 15, 2012, at 4:12 AM, Asta Laiho <asta.laiho at="" btk.fi=""> wrote: >> >>> Many thanks for a quick response! >>> >>> I installed the version 1.1.3 and tried rebuilding the database but got the following error: >>> >>> cuff <- readCufflinks(rebuild=T) >>> Creating database /Users/alaiho/Documents/Data-analysis/Projects/N GS/110007_Liisa_Nissinen/cufflinks/cuffdiff_output/cuffData.db >>> Reading /Users/alaiho/Documents/Data-analysis/Projects/NGS/110007_ Liisa_Nissinen/cufflinks/cuffdiff_output/genes.fpkm_tracking >>> Checking samples table... >>> Populating samples table... >>> Writing genes table >>> Reshaping geneData table >>> Recasting >>> Error in as.data.frame(cast(genemelt, ... ~ measurement)) : >>> could not find function "cast" >>> >>>> sessionInfo() >>> R version 2.14.2 (2012-02-29) >>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >>> >>> locale: >>> [1] C/UTF-8/C/C/C/C >>> >>> attached base packages: >>> [1] stats graphics grDevices utils datasets methods base >>> >>> other attached packages: >>> [1] cummeRbund_1.1.3 reshape2_1.2.1 ggplot2_0.9.0 RSQLite_0.11.1 >>> [5] DBI_0.2-5 >>> >>> loaded via a namespace (and not attached): >>> [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 >>> [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 >>> [9] plyr_1.7.1 proto_0.3-9.2 scales_0.2.0 stringr_0.6 >>> >>> On Mar 15, 2012, at 5:28 AM, Loyal Goff wrote: >>> >>>> Hi Asta, >>>> I believe that I had already fixed this bug in the newer releases of cummeRbund. Can you try to upgrade to the development cummeRbund version (1.1.3)? You will have to rebuild your cuffData.db database with this new version. Please let me know if this does not resolve your issue and I will try to work with you to get it resolved. >>>> >>>> Thanks! >>>> >>>> Cheers, >>>> Loyal >>>> >>>> On Mar 14, 2012, at 9:05 AM, Asta Laiho [guest] wrote: >>>> >>>>> >>>>> Hi, >>>>> >>>>> I am using cummeRbund package to analyze human RNAseq data. Everything else is working very smoothly, except that the function diffData(TSS(cuff)) returns zero rows which to me seems like an error. Other features; genes, isoforms and CDSs are returned correctly. What could be wrong? >>>>> >>>>> This is what I'm doing: >>>>> cuff <- readCufflinks() >>>>> cuff >>>>> CuffSet instance with: >>>>> 4 samples >>>>> 22260 genes >>>>> 37557 isoforms >>>>> 27776 TSS >>>>> 25411 CDS >>>>> 133002 promoters >>>>> 166656 splicing >>>>> 101436 relCDS >>>>> >>>>> As you can see, the resulting data frame has no rows: >>>>> diffData(TSS(cuff)) >>>>> [1] TSS_group_id sample_1 sample_2 status value_1 >>>>> [6] value_2 ln_fold_change test_stat p_value q_value >>>>> [11] significant >>>>> <0 rows> (or 0-length row.names) >>>>> >>>>> The following works correctly: >>>>> head(> head(diffData(CDS(cuff)) ) >>>>> CDS_id sample_1 sample_2 status value_1 value_2 ln_fold_change >>>>> 1 P1 S1 S2 OK 11.30910 19.39660 7.78320e-01 >>>>> 2 P10 S1 S2 NOTEST 3.25391 0.00000 -1.79769e+308 >>>>> 3 P100 S1 S2 NOTEST 2.26812 2.28999 1.38405e-02 >>>>> 4 P1000 S1 S2 NOTEST 0.00000 0.00000 0.00000e+00 >>>>> 5 P10000 S1 S2 LOWDATA 12.03700 17.89210 5.71841e-01 >>>>> 6 P10001 S1 S2 NOTEST 6.83010 4.39573 -6.35805e-01 >>>>> test_stat p_value q_value significant >>>>> 1 -1.27958e+00 0.200693000 0.607328 no >>>>> 2 -1.79769e+308 0.000624649 1.000000 no >>>>> 3 -2.06580e-02 0.983519000 1.000000 no >>>>> 4 0.00000e+00 1.000000000 1.000000 no >>>>> 5 0.00000e+00 1.000000000 1.000000 no >>>>> 6 9.72050e-01 0.331026000 1.000000 no >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> -- output of sessionInfo(): >>>>> >>>>> R version 2.14.2 (2012-02-29) >>>>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) >>>>> >>>>> locale: >>>>> [1] C/UTF-8/C/C/C/C >>>>> >>>>> attached base packages: >>>>> [1] stats graphics grDevices utils datasets methods base >>>>> >>>>> other attached packages: >>>>> [1] cummeRbund_1.0.0 ggplot2_0.9.0 reshape_0.8.4 plyr_1.7.1 >>>>> [5] RSQLite_0.11.1 DBI_0.2-5 >>>>> >>>>> loaded via a namespace (and not attached): >>>>> [1] MASS_7.3-17 RColorBrewer_1.0-5 colorspace_1.1-1 dichromat_1.2-4 >>>>> [5] digest_0.5.1 grid_2.14.2 memoise_0.1 munsell_0.3 >>>>> [9] proto_0.3-9.2 reshape2_1.2.1 scales_0.2.0 stringr_0.6 >>>>> >>>>> -- >>>>> Sent via the guest posting facility at bioconductor.org. >>>>> >>>>> _______________________________________________ >>>>> Bioconductor mailing list >>>>> Bioconductor at r-project.org >>>>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>>>> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >>>> >>> >

ADD REPLY • link 12.1 years ago Loyal Goff ▴ 140

Login before adding your answer.