Modèle linéaire mixte et predictions

Bonjour à tous,
C'est la première fois que je poste ici donc j'espèce que tout sera dans les règles et compréhensible.
Je cherche à prédire des valeurs de pression sonore en fonction de différents paramètres.

Variable réponse : pression sonore qui est de type continue --> pression.
Variable explicative :
- le lieu avec à 2 facteurs 1 ou 2 --> lieu2
- le niveau d'activité avec 2 facteurs : haut ou bas --> activité (baja ou alta)
- le vent avec 3 facteurs :0 (calme), 1(leger), 2(fort) --> vent
- la distance qui est continue (entre 2 et 400 mètres) --> distance
- et une autre variable lieu à effet aléatoire (12 catégories de lieux) --> lieu1

J'ai crée plusieurs modèles et après comparaisons, j'ai retenu celui :
> lmer(formula =pression~lieu2*activité+ vent:activité + distance+ (1|lieu1), data = donnéesnew)

Je souhaiterais ensuite faire des prédictions de la pression sonore en fonction de ce modèle. Et pouvoir les représenter graphiquement.
C'est la première fois que je travaille sur les modèles mixtes donc il y pas mal de choses que je ne sais pas ou ne sais pas faire. Donc voici mes différentes questions :
1) Est-il possible d'effectuer des prédictions sur un tel modèle ?
2) Si oui, auriez vous des pistes pour m'indiquer comment les obtenir et les représenter graphiquement ?

J'espère que ma demande est compréhensible et a un sens, je laisse une partie de mon jeu de données ici.
Je précise que j'effectue quelques transformations avant d'obtenir mon modèle:

données$lieu2 <- factor(données$lieu2)
données$vent <- factor(données$vent)

données$activité <-relevel(données$activité ,ref="Baja")

donnéesnew <- na.omit(données)


En vous remerciant d'avance,

Oui, vous pouvez faire des prédictions avec ce genre de modèles et même construire des intervalles de confiance sur ces prédictions. Lorsqu'il s'agit de modèles mixte, une quantité importante (énorme, exhaustive) d'information se trouve ici.

Vous deviez y trouver votre bonheur..

HTH, Eric.