Hi everyone,

For my master thesis is am looking at mediation btw greenspaces/Physical activity and wellbeing, this is a bit wordy but this is just to give you an idea.

In my model i have 1 predicator separated into 2 dummy variables(coded with 1 and 0), one mediator (physical activity - numerical variable) and one dependent variable for the wellbeing (numerical).

In addition i have 4 covariates, age/gender/health condition/ employement (this last one is also seprated into 2 dummy variable)

i tested my model first with the lm() function from the simpleboot package.

But i wanted to confirm and do the sobel test using the function bmem.sobel() from the bmem package.

i followed these steps :

to precise pa3 is the mediator ;

w3.subjective.wellbeing is my dependent variable (outcome) ;

and GS_5_to_10, GS_10_more are my dummy variables for the independent variable.

the rest are my covariates.

working dataset is my dataset containing my variables.

```
library(bmem)
library(sem)
med.model <- specifyEquations(exog.variances = T) ## to precise the model to work with
1: w3.subjective.wellbeing = b*pa3 + cp1*GS_5_to_10 + cp2*GS_10_more + Age_18to60 + Condition + Employed + part.time+ Gender
2: pa3 = a1*GS_5_to_10 + a2*GS_10_more + Age_18to60 + Condition + Employed + part.time+ Gender
effects <- c('a1*b', 'cp1 + a1*b', 'a2*b', 'cp2+a2*b')
med.res<- bmem.sobel(working_dataset,med.model, effects)
```

When i try to use the bmem.sobel i get the following message:

Error in summary.objectiveML(sem.object) :

coefficient covariances cannot be computed

My guess is it doesn't like the 2 dummy variables for my independent variable but i don't know how i should introduce it for it to be ok

if anyone as an idea ?