Question: GLM model design question
0
gravatar for illinois.ks
4.7 years ago by
illinois.ks140
Korea, Republic Of
illinois.ks140 wrote:

Dear, 

I have a question related with my design matrix.
My RNA-seq experimental design looks like this. (paried design)


patient gender drug
sample1 patient1 F Pre
sample2 patinet1 F Post
sample3 patient2 M Pre
sample4 patient2 M Post
sample5 patient3 F Pre
sample6 patient3 F Post
sample7 patient4 F Pre
sample8 patient4 F Post
sample9 patient5 F Pre
sample10 patient5 F Post
sample11 patient6 M Pre
sample12 patient6 M Post

=====================================

What I want to test is this.

Q1) DEG responsing to drug... (regardless of gender)
Q2) DEG responsing to drug on female
DEG responsing to drug on male.. 

More ultimately, Does genes reponse to the drug differently based on gender??? (or not really?? ) In other words, gender is the significant factor to drug response?? 

Q3) Which gender is more sensitive to the drug...?? 

==================================================
In order to test or answer Q1 , my model looks like this. 

design<- model.matrix(~enzyme + enzyme:gender, data=data)

colnames(design)
##################################################
[1] "(Intercept)" "enzymepost" "enzymepre:genderM" "enzymepost:genderM"
###################################################

lrt1 <- glmLRT(fabry_fit, coef=2) # enzymepost

returns the genes responsing to the drug.. (Q1)


======================================


lrt2 <- glmLRT(fabry_fit, coef=4) # enzymepost:genderM

does it returns that genes responding the drug depending on gender differently ?? In this case, we can answer "DEG responsing to drug on female
DEG responsing to drug on male.. " ?? Q2 ...........

=======================================================

Then, what is the meaing of enzymepre:genderF??? 
I also found that enzymepost:genderM VS enzymepre:genderM returuns relatively similar results!!! 


Could you please somebody to explain the meaning of this two model??

Again, I would like to differentally expressed genes under the drug treatment as well as gender.
(e.g. up after drug treatment for male, but not for female.. something like this... )

thanks in advance,

rna-seq glm • 1.7k views
ADD COMMENTlink modified 4.7 years ago by Devon Ryan90k • written 4.7 years ago by illinois.ks140

Are you sure it's even worth asking the question 'does a gene's response to the drug differ according to gender' when you've only two male subjects?

ADD REPLYlink written 4.7 years ago by russhh4.4k

Hello illinois.ks!

It appears that your post has been cross-posted to another site: SEQanswers, where Gordon Smyth already replied.

This is typically not recommended as it runs the risk of annoying people in both communities.

ADD REPLYlink written 4.7 years ago by Devon Ryan90k
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: 1827 users visited in the last hour