a$).$ Explain how the function myplot$()$ works $($this is not buit$-$in R function. You

can also Learn how to write your own function$)$

b$).$ Compare this two models and write a summary of the findings.

Problem $2$

Read the following R codes and the corresponding output then answer the questions.

Model $1$ without interaction

$11$ brary$($HSAUR$3)$

data$("$womensrole$")$

womensrole cbind $($agree$,$ disagree$)-$ gender $+$ education

womensrole$\_$glm glm$($fm$,$ data $=$ vomensrole, family $=$ binomial$())$

#role. fitted predict$($womensig$\_$glm

Plot function

myplot function$($role$.$fitted$)\{$

f womensrolesgender $==$ "Female"

plot$($womensrole$education, role.fitted, type $="$n$",$

ylab $=$ "Probability of agreeing",

Iines$($womensroleseducation$[|$f$],$ role.fitted$[|$f$],$ lty $=1)$

lines$($womensrole$education$[$f$],$ role.fitted$[$f$],1$ty $=2)$

legend$("$topright$",$ lgtxt$,$ lty $=($Females$)")$

y vomensrolesagree ; $($vomensrole$agree $+"\backslash )$

size womensroleagree$+$womensroledisagree

size size $-$ min$($size$)$

text $($womensre$)$ max$($size$))*3+$

family $=$ "HersheySerif", y$[$f$],"\backslash$VE$",$

text$($womensrole$education$[$If$],$ y$[$lf$]"$ex $=$ size, col$=$"red"$)$

family $=$ "HersheySerif", y$[$lf$],"\backslash$MA$"$

$\}$

Model $2$ with interaction

fm cbind$($agree$,$ di sagree$)-$ gendor $*$ education

vomensrole$\_$glm fm $2,$ data $=$ romensrole, family $=$ binomial $())$

role.fitted fitted$($womensrole$\_$glm$\_2)$

Call myplot fuction

par$(mfrowec(2,1))$

myplot$($role$.$fitted$1)$

myplot$($role$.$fitted$2)$