need step by step solution in R$-$studio with code,

also need code for cleaning up data in the first step

Study Background

The data come from a study investigating the effectiveness of Cognitive Behavioural Therapy $($CBT$)$ across different modes of delivery.

$15$ clinics took part in the study, each aimed to recruit a minimum of $20$ patients. All patients had weekly hour$-$long therapy sessions, but patients opted for a preferred mode of therapy: online chat; telephone call; or in$-$person sessions. Along with a baseline assessment, patients were followed up every $3$ months for up to $18$ months $($minimum $12$ months of follow$-$up$).$ At each assessment, participants completed the Rosenberg Self$-$Esteem Scale $($R$-$SES$),$ a unidimensional scale assessing global self$-$worth by measuring both positive and negative feelings about the self $($scores range from $0$ to $30).$ Data was also collected on two additional known predictors of self$-$esteem: age $($measured at baseline$)$ and activity$-$level, which have previously been found to be positively correlated with self$-$esteem. Table $1$ provides a description of the variables that can be found in the cbtmodes data.

TASK

Assess the extent to which improvements in patient wellbeing $($measured via a scale of self$-$esteem$)$ with CBT depends on whether it is delivered in$-$person, via telephone, or via online chat.

Table $1$:

Data Dictionary for cbtmodes

variable description

clinic Clinic Name

visit Visit number

months Duration $($months$)$ since baseline

age Age $($years$)$ of patient at baseline

activity Activity Level $($measured on a scale of $0$ to $10)$

tmode Therapy Mode $(1=$ in$-$person, $2=$ online chat, $3=$ telephone$)$

patient Patient Name

RSES Rosenberg Self$-$Esteem Scale $($R$-$SES$),$ a unidimensional scale assessing global self$-$worth $($scores range from $0$ to $30)$