I have an issue with a result I just found.
I ran a sub-population analysis (for 2 age groups) with a complex survey design.
For this purpose I fitted a generalized linear model (family:normal,link:identity), using the function svyglm
from the R package survey
.
Here is my code for the model:
model=svyglm(score ~ DLA +relevel(sex,2)+raeth+pir+age+educa+(DLA:relevel(sex,2)),
design,data=imp1)
DLA is my exposure of interest and was included as categorical (tertiles)
For one of my sub-population analyses (for older age group), my results suggest a significant difference of (-7.50 (95% CI:-14.54, -0.46)) between the third tertile and the first tertile of exposure), but the p-value=0.058 suggested a non significant result.
There is not such a big difference, but I wonder how it could be possible that the confidence interval does not include 0, whereas the p-value suggests a non-significant result.
Thank you very much for your answers
svyglm
function. What package is that in? $\endgroup$