df$study <- factor(df$study)

df$type_degree <- ifelse(df$type_degree=="1",1,0)

eq <- lm(pip ~ sex + study + area + type_degree, data=df, na.action=na.exclude)