Comment appliquer un code par grouped data

Bonjour ,
Je viens solliciter votre aide si possible pour un cas et j'espère que vous pourriez m'aider !
Je travaille sur un dataset composé des colonnes :
Perf = response variable
factor1 = explanatory variable 1
factor 2 = explanatory variable 1
Je cherche à appliquer un "lm" et extraire lsmeans et le résultats du test tukey mais en appliquant mon code par "group".
J'ai pu créer un code qui répond parfaitement à mon besoin, mais pour un seul "group", et mon besoin est d'essayer obtenir un output renvoyant en une seule table ma dataset originale + les colonnes résultant de mon code : lsmans + tukey 'letters' colonne by "group" .
Voici le code qui fonctionne pour un seul group :
''''
library(Rmisc)
library(multcomp)
library(lsmeans)
library(dplyr)
library(broom)

output <- as.data.frame(dataset)
output$perf <- as.numeric(output$perf)
output$factor1 <- as.factor(output$factor1)
output$factor2 <- as.factor(output$factor2)

model <- lm(perf ~ factor1 + factor2, data = output)
lsmeans <- lsmeans(model, ~ factor1)
CLD = cld(lsmeans(lm(perf ~ factor1 + factor2, data = output), ~ factor1),
alpha = 0.05,
Letters = letters,
adjust = "tukey")

output$stud_resid <- studres(model)
output <- left_join(output, CLD, by = "factor1")
''''
Voici un extrait de ma dataset :

group perf factor1 factor2
1 32 1 1
1 44 1 2
1 58 1 3
1 76 2 1
1 73 2 2
1 37 2 3
1 52 3 1
1 78 3 2
1 60 3 3
2 93 1 1
2 78 1 2
2 25 1 3
2 97 2 1
2 85 2 2
2 60 2 3
2 70 3 1
2 62 3 2
2 95 3 3

Pourriez-vous m'aider à appliquer ce code par groupe et combiner les résultats en une seule table?
En vous remerciant d'avanace 🙏