I would like compare the job satisfaction between two industries. I've conducted a t-test and came to the conclusion that they are statistically different at the 5% significance level.
Now I would like to analyze how professional and demographic variables influence job satisfaction in the two industries.
lm_industry1 <- lm(JS~AGE+EDU+SIZE+SCOPE+MGMT+SCOPE+TENURE+COLLAR+SEX, data=data_industry1) lm_industry2 <- lm(JS~AGE+EDU+SIZE+SCOPE+MGMT+SCOPE+TENURE+COLLAR+SEX, data=data_industry2)
(COLLAR and SEX are binary, all the other variables are continuous)
Could this work with a multiple regression? What about non-linear relationships and interaction terms?