Study the following code example and add comments to describe what it does. Ensure you comment on each line of the code

library$($tidyverse$)$

library$($MASS$)$ l

oad$("./$data$/$assess$\_$data$\_0824-$resit.Rdata"$)$

idx $<-20$

c$\_$cl $<-10$:$30$

boxplot$($log$2($Y$[,$ c$\_$cl$]+1))$

x $<-$ patient$\_$data$tissue$[$c$\_$cl$]$

z $<-$ patient$\_$data$patient$[$c$\_$cl$]$

tmp $<-$ data.frame$($y $=$ Y$[$idx$,$ c$\_$cl$],$ x $=$ x$,$ z $=$ z$,$ lib$\_$size $=$ colSums$($Y$[,$ c$\_$cl$]))$

out $<-$ glm$($y ~ x $+$ z $+$ lib$\_$size, data $=$ tmp$,$ family $=$ "poisson"$)$

p$\_$val $<-$ summary$($out$)$$coefficients$[2,4]$