I have been reading about gtsummary to make regression tables, also huxtable and jtools. I would like to make a table of the p-values for 6 different cell types. The cell types are the outcome variable and I want to see if the proportion of cells associated with Cases (model 1) is different when I adjust for Time of blood draw (model 2)
Model 1: lm(cellTypes ~ factor(Case_Control) + smoking score + Age + factor(Gender) + 15PCs, data = df
Model 2: lm(cellTypes ~ factor(Case_Control) + smoking score + Age + factor(Gender) + 15PCs + Time_of_blood_draw, data = df
I made the table by hand but would like to learn to make a publication-ready table.
Thanks for any tips. I would also like to know how to omit the 15 PCs when I print the table. I tried adding a vector of the covariates I wanted to include. This worked, but it wouldn´t accept Case_Control or Gender which are categorical.
Thanks for your help.