4.6.2 - same problems with sam-r and limma-r

Help
slava
2010-11-24
2013-05-20
  • slava
    slava
    2010-11-24

    Dear TMEV team!

    Downloaded, installed and tested TMEV_4_6_2_r2536_win for SAM-R and LIMMA-R.

    Test file is .tdms, 4 samples, 3-4 replicates each, 9553 genes, normalized log2 ratios, no missing values.

    SAM-R settings: use R; Number of permutations=calculated suggested; Select S0=using Tusher et al. method; Calculate q-values?=Yes.

    One-class test (one sample with 4 replicates selected from the list):
    In both Significant Genes and Non-significant Genes Tables:
    column Numerator (r)=0.0 for all genes
    column Denominator (s+s0)=0.0 for all genes
    column q-value (%)=0.0 for all genes
    First two columns, apparently, should have some meaningful values, and as genes were arranged into Significant Genes and Non-significant Genes, there must be some meaningful q-values too.
    Two-class unpaired test (two samples with 4 replicates each selected from the list) – exactly the same problems with the same columns.
    Also, for both tests, in General Information-Computed Quantities, Computed Exchangeability Factor s0=0.0, s0 Percentile=0.0, Pi0Hat=0.0, which is unlikely.
    Is there a problem with reading the output data from samr to the tables? Interestingly, when I run exactly the same analysis, but NOT using the R, all output columns in all gene tables has meaningful values, as well as in General Information-Computed Quantities.

    LIMMA-R settings : Significance Level: Alpha=.05
    One Class test (one sample with 4 replicates selected from the list):
    Gene Tables output has no columns for significance values, adj-p-values, log fold change, t-statistic, log-odds, but has irrelevant column F-values, containing only 0.0 (applicable to 3 or more samples).
    Two Class test ((two samples with 4 replicates each selected from the list)-all columns are right, but again has this irrelevant column F-values, containing only 0.0 (applicable to 3 and more samples).
    Is there a problem with reading the output data from limma-r to the tables for one-class experiment? And why the column F-value is popping up where it shouldn't?

    For me it looks like these two very useful additions to the data analysis process are again not working properly, the same as I reported for 4_6_1. Or is it only me who's having these problems?

    Regards,

    Slava.

     
  • raktim
    raktim
    2010-11-24

    Hi, We will definitely look into this and get back to you.
    Raktim