-
Notifications
You must be signed in to change notification settings - Fork 11
Description
Hello--
I'm running into issues with the logistic regression model - no error, but the program is returning all NA's for those parameters. All other aspects of PyHLA are working for me, so I don't think it's an issue with my input files.
Here is my code:
python2 PyHLA.py
--input HLA_input.txt
--exclude exclude_alleles.txt
--adjust FDR
--assoc
--test logistic
--model additive
--covar covar.txt
--out logistic.out
Here is the first part of the output (truncated):
Allele A_case B_case A_ctrl B_ctrl F_case F_ctrl Freq P_Logit OR L95 U95 P_adj
HLA-A01:01 38 240 47 185 0.1367 0.2026 0.1667 NA NA NA NA NA
HLA-A02:01 77 201 65 167 0.2770 0.2802 0.2784 NA NA NA NA NA
HLA-A03:01 34 244 21 211 0.1223 0.0905 0.1078 NA NA NA NA NA
HLA-A11:01 19 259 12 220 0.0683 0.0517 0.0608 NA NA NA NA NA
HLA-A*24:02 28 250 23 209 0.1007 0.0991 0.1000 NA NA NA NA NA
I get the same NA values for all HLA alleles.
Thanks for any advice you can provide.
Melanie