0
$\begingroup$

I am trying to do linear regression using retinopathy as response and 342 proteins as predictor. the model should be adjusted for age, sex and BMI. There are no missing values for any variables.

The model fitting was as follows...

mod_ret <- lm(ret~x+sex+age+bmi, data = df)

while running summary I have got following results..

summary(mod_ret)

Can you please provide insight into why I am getting no statistics and how to resolve that? Thanks!

$\endgroup$
3
  • $\begingroup$ Please don't put screenshots of text into StackExchange posts. Copy the text and paste it into your post. Don't worry about how it looks; other people can fix that up if it looks weird and you don't know how to do that. $\endgroup$
    – gringer
    Commented Apr 7 at 12:25
  • 1
    $\begingroup$ The intercept labels look odd; I'd usually expect a category name rather than a sample name. Can you show the first few lines of df? (e.g. head(df)) $\endgroup$
    – gringer
    Commented Apr 7 at 21:26
  • $\begingroup$ I suspect that you are treating age as a categorical variable rather than a numerical covariate. That would lead to the kinds of singularities due to overfitting on your categorical data- basically, you only have <=1 person in any given category due to the extreme sparsity when you treat age as categorical. $\endgroup$ Commented May 14 at 18:07

1 Answer 1

0
$\begingroup$

Are you sure this is linear regression? Gender is logistic regression, i.e. it's a boolean value. Really linear regression is y = mx + c, here this is multi-linear regression for age and bmi (GLM) and logistic regression for gender.

I would simply choose one variable either age or bmi repeat the calculation and remove everything else. I minimise my R BTW favouring scipy and Pythons statsmodels.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.