search

hanchenphd / GMMAT

Generalized linear Mixed Model Association Tests
Other
36 stars 22 forks source link

glmm.wald exporting all NAs #20

Closed ntilhou closed 3 years ago

ntilhou commented 3 years ago

Hello,

For some reason "glmm.wald" is exporting only NAs when I provide the inputs below. For reference, it's returning the full data frame full of NAs. I'm not sure where my error is.

Thank you, Neal

snps<-c("Chr01K_671719","Chr01K_676401","Chr01K_676418","Chr01K_676677","Chr01K_676813") glmm.wald(fixed = trt ~ 1, data = y, kins = K, id = "ID",verbose=T,snp.col=1, family = quasibinomial(link = "logit"), infile = "tmp_infile.txt", snps=snps, infile.nrow.skip = 0, infile.ncol.skip = 1, infile.ncol.print = 1,infile.sep = ",")

K.txt tmp_y.txt tmp_infile.txt

ntilhou commented 3 years ago

Oh no! Github gave me an error before so I did not think it accepted the issue I previously submitted. Sorry for submitting this three times!!