Question: T-Test Filtering Problem In R (Simpleaffy)
1
gravatar for Andrii
7.5 years ago by
Andrii10
Andrii10 wrote:

Dear all! I need your help!!! Situation is next: I tried to perform t-test in R. I applied function

result.0vs4<-pairwise.comparison(eset,"time",c("control","4 h"))

to the expression set object (Affymetrix)

result.0vs4
........
AFFX-r2-P1-cre-5_at              AFFX-ThrX-3_at 
           1.571096e-01                2.742916e-01 
         AFFX-ThrX-5_at              AFFX-ThrX-M_at 
           2.940117e-01                8.981617e-01 
        AFFX-TrpnX-3_at             AFFX-TrpnX-5_at 
           6.314353e-01                4.371810e-01 
        AFFX-TrpnX-M_at 
           4.755571e-01 
attr(,"Csingle")
[1] TRUE

Slot "calls":

<0 x 0 matrix>

Slot "group":

[1] "time"

Slot "members":

[1] "control" "4 h"

Slot "pData":

          sample     agent    time
GSM442026 GSM442026 untreated control
GSM442027 GSM442027 untreated control
GSM442028 GSM442028 untreated control
GSM442037 GSM442037  TGF-beta     4 h
GSM442038 GSM442038  TGF-beta     4 h
GSM442039 GSM442039  TGF-beta     4 h

...............................

After that everything was fine, every slot of PairComp object (resulted object after pairwise.comparison function application) had information about fold change, means etc. Till this point no problem, but after applying filtering each slot was empty:-( I used filter contritions as follows: p> 0.001, fold change>2.

result.0vs4.f=pairwise.filter(result.0vs4,tt=0.001,fc=2)
result.0vs4.f
An object of class "PairComp"

Slot "means":

control 4 h

Slot "fc":

named numeric(0)

Slot "tt":

named numeric(0)

Slot "calls":

<0 x 0 matrix>

Slot "group":

[1] "time"

Slot "members":

[1] "control" "4 h"

Slot "pData":

          sample     agent    time
GSM442026 GSM442026 untreated control
GSM442027 GSM442027 untreated control
GSM442028 GSM442028 untreated control
GSM442037 GSM442037  TGF-beta     4 h
GSM442038 GSM442038  TGF-beta     4 h
GSM442039 GSM442039  TGF-beta     4 h

..................

I set less strict parameters and everything was the same:-( If somebody knows the reason of that don't stand aside :-) Thank you in advance! Andrii

• 1.6k views
ADD COMMENTlink modified 7.5 years ago by Neilfws48k • written 7.5 years ago by Andrii10
2
gravatar for Neilfws
7.5 years ago by
Neilfws48k
Sydney, Australia
Neilfws48k wrote:

I think you mean p < 0.001, not p > 0.001. Anyway, the most likely reason for the empty slots - assuming the code ran without errors - is simply that no probeset matched your filter conditions.

ADD COMMENTlink written 7.5 years ago by Neilfws48k
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 2.3.0
Traffic: 953 users visited in the last hour