load("/Users/oroce/Downloads/ora3.RData") Employee_data ls() getwd() library("foreign") attach(Employee_data) ls() table(GENDER,JOBCAT) chisq.test(table(GENDER,JOBCAT)) chisq.test(table(GENDER,JOBCAT)[1:7]) chisq.test(table(GENDER,JOBCAT))[1:7] install.packages("vcd") library("vcd")as assocstats(table(GENDER,JOBCAT)) table(GENDER,JOBCAT) table(EDUC,JOBCAT) table(EDUC,JOBCAT)->lamtab apply(lamtab, 2, sum) sum(apply(lamtab, 2, sim)) sum(apply(lamtab, 2, sum)) max(apply(lamtab, 2, sum)) sum(apply(lamtab, 2, sum)) - max(apply(lamtab, 2, sum)) sum(apply(lamtab, 2, sum)) - max(apply(lamtab, 2, sum))–>A (sum(apply(lamtab, 2, sum)) - max(apply(lamtab, 2, sum)))–>A (sum(apply(lamtab, 2, sum)) - max(apply(lamtab, 2, sum)))->FOO FOO apply(lamtab, 1, sum) apply(lamtab, 1, max) apply(lamtab, 1, sum) - apply(lamtab, 1, max) sum(apply(lamtab, 1, sum) - apply(lamtab, 1, max)) sum(apply(lamtab, 1, sum) - apply(lamtab, 1, max)) -> BAR 1-FOO/BAR 1-BAR/FOO table(jobcat,educ)->lamtab table(JOBCAT,EDUC)->lamtab (sum(apply(lamtab, 2, sum)) - max(apply(lamtab, 2, sum)))->FOO (sum(apply(lamtab, 2, sum)) - max(apply(lamtab, 2, sum)))->FOO sum(apply(lamtab, 1, sum) - apply(lamtab, 1, max))->BAR 1-BAR/FOO->LAMBDA 1-BAR/FOO->LAMBDA LAMBDA table(GENDER, JOBCAT, MINORITY) install.packages("MASS") 39 library("MASS") table(GENDER,JOBCAT,deparse.level=2)->K2 loglm(~GENDER+JOBCAT,data=K2) library("MASS") table(GENDER,JOBCAT,MINORITY,deparse.level=2)->K3 loglm(~GENDER+JOBCAT+MINORITY,data=K2) loglm(~GENDER+JOBCAT+MINORITY,data=K3) loglm(~GENDER+JOBCAT+MINORITY+GENDER:JOBCAT,data=K3) loglm(~GENDER+JOBCAT+MINORITY+GENDER:JOBCAT+JOBCAT:MINORITY,data=K3) loglm(~GENDER+JOBCAT+MINORITY+GENDER:MINORITY+JOBCAT:MINORITY,data=K3) SALARY[which(JOBCAT=="Clerical")] histogram(SALARY[which(JOBCAT=="Clerical")]) hist(SALARY[which(JOBCAT=="Clerical")]) shapiro.test(SALARY[which(JOBCAT=="Clerical")]) install.packages("car") library("car") levene.test(SALARY,GENDER) quit()