<div dir="ltr"><div class="gmail_extra">André, mão está otimizado mas já te dá um norte.</div><div class="gmail_extra"><br></div><div class="gmail_extra">Estou assumindo que o resp do plot é o rep do objeto dados</div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_extra">b <- sort(unique(dados$B)) # ou levels(dados$B)</div><div class="gmail_extra">modelos.summary <- vector(mode = 'list', length = length(b))</div><div class="gmail_extra">for (i in 1:length(b)) {</div><div class="gmail_extra">  dados.modelo <- subset(x = dados, subset = B == b[i])</div><div class="gmail_extra">  modelos.summary[[i]] <- summary(lm(y ~ rep, data = dados.modelo))</div><div class="gmail_extra">  png(paste0("Gráfico ", b[i], ".png")) </div><div class="gmail_extra">  plot(dados.modelo$y ~ dados.modelo$rep)</div><div class="gmail_extra">  dev.off()</div><div class="gmail_extra">}</div><div class="gmail_extra">names(modelos.summary) <- b</div></div></div>