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(~drug + drug:gender, data=data)
colnames(design)
##################################################
[1] "(Intercept)" "drugpost" "drugpre:genderM" "drugpost:genderM"
###################################################
lrt1 <- glmLRT(fabry_fit, coef=2) # drugpost
returns the genes responsing to the drug.. (Q1)
======================================
lrt2 <- glmLRT(fabry_fit, coef=4) # drugpost: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 drugpre:genderF???
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,
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(~drug + drug:gender, data=data)
colnames(design)
##################################################
[1] "(Intercept)" "drugpost" "drugpre:genderM" "drugpost:genderM"
###################################################
lrt1 <- glmLRT(fabry_fit, coef=2) # drugpost
returns the genes responsing to the drug.. (Q1)
======================================
lrt2 <- glmLRT(fabry_fit, coef=4) # drugpost: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 drugpre:genderF???
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,
Comment