Hi!
I have two groups:
Healthy=26
Diabetic=22
> colnames(design)
[1] "(Intercept)" "DiseaseDiabetic"
[3] "DiseaseHealthy:Patients2" "DiseaseDiabetic:Patients2"
[5] "DiseaseHealthy:Patients3" "DiseaseDiabetic:Patients3"
[7] "DiseaseHealthy:Patients4" "DiseaseDiabetic:Patients4"
[9] "DiseaseHealthy:Patients5" "DiseaseDiabetic:Patients5"
[11] "DiseaseHealthy:Patients6" "DiseaseDiabetic:Patients6"
[13] "DiseaseHealthy:Patients7" "DiseaseDiabetic:Patients7"
[15] "DiseaseHealthy:Patients8" "DiseaseDiabetic:Patients8"
[17] "DiseaseHealthy:Patients9" "DiseaseDiabetic:Patients9"
[19] "DiseaseHealthy:Patients10" "DiseaseDiabetic:Patients10"
[21] "DiseaseHealthy:Patients11" "DiseaseDiabetic:Patients11"
[23] "DiseaseHealthy:Patients12" "DiseaseDiabetic:Patients12"
[25] "DiseaseHealthy:Patients13" "DiseaseDiabetic:Patients13"
[27] "DiseaseHealthy:TreatmentAfterTreatment" "DiseaseDiabetic:TreatmentAfterTreatment"
> fit <- glmFit(d,design)
Error in glmFit.default(y = y$counts, design = design, dispersion = dispersion, :
Design matrix not of full rank. The following coefficients not estimable:
DiseaseDiabetic:Patients12 DiseaseDiabetic:Patients13
How can I adjust for uneven groups?
Thanks!
I have two groups:
Healthy=26
Diabetic=22
> colnames(design)
[1] "(Intercept)" "DiseaseDiabetic"
[3] "DiseaseHealthy:Patients2" "DiseaseDiabetic:Patients2"
[5] "DiseaseHealthy:Patients3" "DiseaseDiabetic:Patients3"
[7] "DiseaseHealthy:Patients4" "DiseaseDiabetic:Patients4"
[9] "DiseaseHealthy:Patients5" "DiseaseDiabetic:Patients5"
[11] "DiseaseHealthy:Patients6" "DiseaseDiabetic:Patients6"
[13] "DiseaseHealthy:Patients7" "DiseaseDiabetic:Patients7"
[15] "DiseaseHealthy:Patients8" "DiseaseDiabetic:Patients8"
[17] "DiseaseHealthy:Patients9" "DiseaseDiabetic:Patients9"
[19] "DiseaseHealthy:Patients10" "DiseaseDiabetic:Patients10"
[21] "DiseaseHealthy:Patients11" "DiseaseDiabetic:Patients11"
[23] "DiseaseHealthy:Patients12" "DiseaseDiabetic:Patients12"
[25] "DiseaseHealthy:Patients13" "DiseaseDiabetic:Patients13"
[27] "DiseaseHealthy:TreatmentAfterTreatment" "DiseaseDiabetic:TreatmentAfterTreatment"
> fit <- glmFit(d,design)
Error in glmFit.default(y = y$counts, design = design, dispersion = dispersion, :
Design matrix not of full rank. The following coefficients not estimable:
DiseaseDiabetic:Patients12 DiseaseDiabetic:Patients13
How can I adjust for uneven groups?
Thanks!
Comment