paired samples, comparing expression given a continuous variable
1
0
Entering edit mode
naz • 0
@b77d67c2
Last seen 5 weeks ago
Canada

I am comparing expression values in paired samples given a continuous variable. In theory, I should see an increase in expression given an increase in the total score between the paired samples.

example data

This is the code that I have generated to run the analysis. I am using the section "9.4.1 Paired Samples" of the limma user guide as a reference.


ID <- factor(mSetSqFlt$Collaborator.ID)

Visit <- factor(mSetSqFlt$TimePoint, levels = c("Baseline","Followup"))

Score <- factor(mSetSqFlt$Total.Score)


design <- model.matrix(~ID + Visit + Visit:Score)

fit <- lmFit(mVals, design)

fit2 <- eBayes(fit)

topTable(fit2, coef = "VisitFollowup")

I get outputs using this method and would just like some insights into if I am on the right track.

Any insights would be very much appreciated as I am new to limma.

limma • 200 views
ADD COMMENT
2
Entering edit mode
@gordon-smyth
Last seen 4 hours ago
WEHI, Melbourne, Australia

You are including too many terms in the linear model. TimePoint and Total.Score give almost the same information. You can include either one of them in the model but not both.

The model you have defined in your question is recommended in the User's Guide for time-course experiments in independent groups, rather than for paired experiments. Defining an interaction is not at all relevant in the paired comparison context.

ADD COMMENT
0
Entering edit mode

Thank you for your help, Dr. Smyth. I appreciate your insights and I will adjust my analysis as you recommend.

ADD REPLY

Login before adding your answer.

Traffic: 443 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6