When I fit a logistic regression in R using glm() (I have 1 binary outcome and 1000 predictors). I do not want to reduce my predictors, but glm() yields beta coefficients for the first 50 predictors ONLY. However, when I use Python Sklearn, I can get all beta coefficients for all predictors. How to solve this problem? I have to use R, and I cannot do dimension reduction.
Tried in Python and it works. But R does not.