Как создать цикл с openxslx и R?

У меня есть два фрейма данных. один - необработанные данные, другой - агрегированные.

Мои два фрейма данных - это варианты наборов данных радужной оболочки глаза (iris и iris_summary)

iris$ID <- 1:150

iris_summary <- iris %>% 
  group_by(Species) %>% 
  summarize(
    "Count of Species" = n(),
    "Average Sepal.Length" = mean(Sepal.Length),
    "Average Sepal.Width" = mean(Sepal.Width),
    "Average Petal.Length" = mean(Petal.Length),
    "Average Petal.Width" = mean(Petal.Width)
    )

Теперь у меня есть книга Excel, которую я хотел бы загрузить в свою среду, с предустановленным форматированием и 4 листами.

Лист 1: Лист с инструкциями 2: Сводная таблица Лист 3: Лист исходных данных 4: Контактная информация

Мне нужно создать 3 отдельных файла Excel для каждого вида (setosa, veriscolor и virginica).

В листе 2: мне нужно включить сводную статистику по каждому виду. Это пример статистики setosa в сводной таблице

И в листе 3: мне нужно включить необработанные данные только сетоса.

Листы 1 и 4 можно не трогать.

Как бы выглядела петля? Мне понадобится файл Excel - используя эту отформатированную рабочую книгу для всех 3 видов в наборе данных Iris. Мои фактические данные содержат около 30 категориальных переменных (по сравнению с 3 «видами» в наборах данных по радужной оболочке), и копирование и вставка их вручную займет часы.

Код обновления пробовал

wb <- loadWorkbook("Sample Workbook.xlsx")

for (i in iris$Species){
  print(i)
  FilterData <- iris %>% filter(Species == i) 
  filename <- paste0(i, "-data.xlsx")
  
  writeData(wb, sheet="Summary Table", i)
}