Изменить print () на операторы присваивания?

Question

Я хочу использовать результат печати (a c), полученный в результате повтора l oop, для дальнейшего анализа, и я пробовал различные способы, но это не было достигнуто. Пожалуйста, помогите мне изменить print () на оператор присваивания, который можно вызвать в дальнейшем анализе.

z=2
cut_avg=cutree(hclust.avg,k=z)     #cutting dendrogram
clst <- function(dat1, dat2, k=z){
ct <- vector("list", length =k)
rp <- vector("list", length = k)
for(i in 1:k){

   ct[[i]]=dat1[rownames(dat2)[cut_avg==i],]
   rp[[i]]=replace(ct[i]$sch_no,c(1:length(ct[i]$sch_no)),i)
   ct[[i]]$sch_no=rp[[i]]
  }
    cts = function(x) Reduce(rbind, x)
    cts = cts(ct)
    return(cts)
 }
 out2=clst(school, clt_sch, k=z)[,1:8]
 csch_intercept=lmer(phySCORE2~mthSCORE+(1|sch_no),data=out2,REML=FALSE)
 ac=AIC(csch_intercept)

 repeat {
    cut_avg=cutree(hclust.avg,k=z)     #cutting dendrogram
    clst <- function(dat1, dat2, k=z){
    ct <- vector("list", length =k)
    rp <- vector("list", length = k)
    for(i in 1:k){

       ct[[i]]=dat1[rownames(dat2)[cut_avg==i],]
       rp[[i]]=replace(ct[i]$sch_no,c(1:length(ct[i]$sch_no)),i)
       ct[[i]]$sch_no=rp[[i]]
         }
    cts = function(x) Reduce(rbind, x)
    cts = cts(ct)
    return(cts)
  }

 out2=clst(school, clt_sch, k=z)[,1:8]
 csch_intercept=lmer(phySCORE2~mthSCORE+(1|sch_no),data=ou2,REML=FALSE)
 ac=AIC(csch_intercept)

 z <- z+1
  print(ac)

 if(z > 21) {
  break
  }
 }

  The output
   [1] 1468.853
   [1] 1467.783
    -----------
   [1] 1421.206