write a code in r that runs model diagnostics for the following gjrm model

#loading libraries

library$($haven$)$

library$($tidyverse$)$

library$($GJRM$)$

library$($gt$)$

library$($gtsummary$)$

library$($labelled$)$

library$($mice$)$

# loading main dhs dataset

MWIR$7$AFL $<-$ read$\_$dta$("$MW$\_2015-16\_$DHS$\_12242024\_945\_208869/$MWIR$7$ADT$/$MWIR$7$AFL.DTA"$)$

#attaching the dataset for variable easy access

attach$($MWIR$7$AFL$)$

MWIR$7$AFL$v$501$

# Creating dataset for the analysis by extracting necessary variable from dhs dataset

mydata $<-$ data.frame$($age$\_$at$\_$marrige$=$v$511,$ years$\_$since$\_$marriage$=$v$512,$

respondent$\_$age$=$v$012,$ marital$\_$status$=$ v$501$

$)$

# for the variable age at marriage, we replace with respondents age for unmarried

# the difference will be captured by the censoring indicator

mydata$age$\_$at$\_$marrige $<-$ ifelse$($is$.$na$($mydata$age$\_$at$\_$marrige$),$

mydata$respondent$\_$age,

mydata$age$\_$at$\_$marrige$)$

#generating time to first marriage censoring indicator using marital status

# $0$ for unmarried and $1$ for others

mydata$cenc$1<-$ ifelse$($mydata$marital$\_$status$==0,0,1)$

# generating divorce cencoring indicator using marital status

mydata$cenc$2<-$ ifelse$($mydata$marital$\_$status $\%$in$\%$ c$(4,5),1,0)$

# for time since first marriage, we replace $0$ for unmarried since they have never married

mydata$years$\_$since$\_$marriage $<-$ ifelse$($is$.$na$($mydata$years$\_$since$\_$marriage$),0,$ mydata$years$\_$since$\_$marriage$)$

# adding independent variables

# Assuming 'mydata' and 'MWIR$7$AFL' are already loaded data frames

mydata$place$\_$of$\_$residence $<-$ MWIR$7$AFL$v$025$

mydata$region $<-$ MWIR$7$AFL$v$024$

mydata$education $<-$ MWIR$7$AFL$v$106$

mydata$religion $<-$ MWIR$7$AFL$v$130$

mydata$ethnicity $<-$ MWIR$7$AFL$v$131$

mydata$wealth$\_$status $<-$ MWIR$7$AFL$v$190$

# unlabelling the dataset

mydata $<-$ unlabelled$($mydata$)$

# generating equations for the two models

# time to first marriage and time to first marriage dissolution

eq$1<-$ age$\_$at$\_$marrige ~ s$($log$($age$\_$at$\_$marrige$),$ bs $=$ "mpi"$)+$ respondent$\_$age $+$

place$\_$of$\_$residence$+$region$+$education$+$religion$+$ethnicity$+$wealth$\_$status

eq$2<-$ years$\_$since$\_$marriage ~ s$($log$($years$\_$since$\_$marriage$),$ bs $=$ "mpi"$)+$ respondent$\_$age $+$

place$\_$of$\_$residence$+$region$+$education$+$religion$+$ethnicity$+$wealth$\_$status

#fitting bivariate copula model

out $<-$ gjrm$($list$($eq$1,$ eq$2),$ data $=$ mydata,

copula $="$C$0",$

margins $=$ c$("-$cloglog", $"-$probit"$),$

cens$1=$ cenc$1,$ cens$2=$ cenc$2,$ model $="$B$")$