I want to use multivariate logistic regression for disease association (binary outcome) of HLA genes. I want to control effect of coexisting alleles of ethnicity haplotypes. Imagine that I have two genes (#layer1) and each gene has 4 alleles (#layer2) as bellow.
HLA-A
-A1
-A2
-A3
-A4
HLA-B
-B1
-B2
-B3
-B4
The problem is that since each person has 2 copies of allele, #layer2 items are not mutually exclusive and therefore this layer should be multivariate in turn (of course if use genotypes, only #layer1 would be multivariate and therefore #layer2 items are numbered as members of #layer1 variables. I have previously done this work Article Association of FOXP3 gene polymorphisms with risk of preecla...
table 4).I can consider #layer2 as the only layer and number the members of each variable as 0 (neg) and 1 (pos) like the above paper as bellow.
A1
0
1
A2
0
1
...
and separately for HLA-B
...
However, both layers are important to me. Because I want to remove confounding effect of both couple alleles and ethnicity associated haplotypes.
What should I do?
I am not familiar with stata commands. Please answer me via stata menus.
Kind regards!
Suggest you look at our paper which is attached. Good luck, David Booth
Thank you dear Prof. David Eugene Booth
I found this reference in your paper Article Performance of using multiple stepwise algorithms for variab...
I understood that you used "multiple logistic regression" is it the very "multivariate logistic regression"?
Hi. Lets start with the simplest things: Tell me about your Dependent variables. How many, their properties: binary, continuous, etc, how many DVs do you want to study in one regression model, and so on. That will help us get started. Looking forward to your reply. David Booth
Dear Dr, thanks a lot
The dependent variable is severity of covid-19 which is binary (sever/no). There is no further DV.
However I have two step independent variables. HLA genes and their alleles.
I have 4 genes and each gene has more than 10 alleles.
All genes are present in all individuals. But each person has 2 alleles (different or repetitive) for each gene.
I want to adjust the effect of the other allele of each gene in person, and also adjust the effect of the parallel genes. Because both of them are ethnicity dependent.
Is this what I am looking for?
https://www.youtube.com/watch?v=SbwApki_BnI
- Badges
- Science topic
- Similar topics
- Medicine
- Epidemiology