在 function 中使用多個變量或在 R 中使用 for 循環

Question

我一直無法在以前的帖子中找到我需要的答案，所以我開始了一個新問題。 我正在嘗試將多個變量輸入到 for 循環或 function 中，以便通過變量的每個步驟獲得一個 output。 我不需要變量的所有排列，我只需要循環或 function 對向量的每個元素並行執行一次（以相同的順序）。 這是我所做的...

a=c("BP GOterms","CC GOterms","MF GOterms")
i = c(0,2,4)
j = c(2,4,6)
n = rep(n,3)

func<-function(a,i,j,n) {
  temp<-arrange(df[df$ID %in% hit[(i*n+1):(n*j)],], -Cluster,Count)
  temp$Description<-factor(temp$Description,levels=unique(temp$Description,fromLast=T))
  plot<-ggplot(temp,aes(Cluster,Description,size=Count,color=p.adjust))+geom_point()+ggtitle(a)
  print(plot)
}

mapply(func,i=i,j=j,n=n,a=a)

這給了我正確的 output 3 次，我只想要一次

或者，我一直在嘗試這樣的事情......

for (i in c(0,2,4) & j in c(2,4,6) & title in c("BP GOterms","CC GOterms","MF GOterms")) {
  temp<-arrange(df[df$ID %in% hit[(i*n+1):(n*j)],], -Cluster,Count)
  temp$Description<-factor(temp$Description,levels=unique(temp$Description,fromLast=T))
  plot<-ggplot(temp,aes(Cluster,Description,size=Count,color=p.adjust))+geom_point()+ggtitle(title)
  print(plot)
}

我收到以下錯誤：

Warning messages:
1: In (i * n + 1):(n * j) :
  numerical expression has 81 elements: only the first used
2: In (i * n + 1):(n * j) :
  numerical expression has 81 elements: only the first used
3: In Ops.factor(Cluster) : ‘-’ not meaningful for factors

要運行我的代碼，您需要非常大的 df。 但這里有一個例子說明它的樣子......

structure(list(Cluster = structure(c(2L, 2L, 2L, 2L, 2L, 2L), .Label = c("M1", 
"M2"), class = "factor"), ONTOLOGY = structure(c(1L, 1L, 1L, 
1L, 1L, 1L), .Label = c("BP", "CC", "MF"), class = "factor"), 
    ID = c("GO:1901673", "GO:0034314", "GO:0006892", "GO:0036503", 
    "GO:0035967", "GO:0006986"), Description = structure(c(1L, 
    2L, 3L, 6L, 7L, 8L), .Label = c("regulation of mitotic spindle assembly", 
    "Arp2/3 complex-mediated actin nucleation", "post-Golgi vesicle-mediated transport", 
    "membrane fusion", "establishment of protein localization to membrane", 
    "ERAD pathway", "cellular response to topologically incorrect protein", 
    "response to unfolded protein", "vacuolar transport", "response to topologically incorrect protein", 
    "cytosolic transport", "endosomal transport", "vacuole organization", 
    "DNA-dependent DNA replication", "macroautophagy", "response to endoplasmic reticulum stress", 
    "Golgi vesicle transport", "rRNA processing", "myeloid leukocyte activation", 
    "rRNA metabolic process", "mitotic nuclear division", "DNA replication", 
    "regulation of mitotic cell cycle phase transition", "regulation of cell cycle phase transition", 
    "ribosome biogenesis", "proteasome-mediated ubiquitin-dependent protein catabolic process", 
    "chromosome segregation", "regulation of chromosome organization", 
    "autophagy", "process utilizing autophagic mechanism", "proteasomal protein catabolic process", 
    "ncRNA processing", "mitotic cell cycle phase transition", 
    "cell cycle phase transition", "regulation of DNA metabolic process", 
    "negative regulation of cell cycle", "ribonucleoprotein complex biogenesis", 
    "DNA repair", "ncRNA metabolic process", "regulation of mitotic cell cycle"
    ), class = "factor"), GeneRatio = c("8/1090", "10/1090", 
    "18/1090", "18/1090", "20/1090", "21/1090"), BgRatio = c("21/19527", 
    "32/19527", "95/19527", "95/19527", "97/19527", "104/19527"
    ), pvalue = c(9.7302693101372e-06, 5.88348341881873e-06, 
    4.47976958483539e-06, 4.47976958483539e-06, 3.24568370146555e-07, 
    2.40596847523017e-07), p.adjust = c(0.00248177711775767, 
    0.00171451628569576, 0.0014795185682183, 0.0014795185682183, 
    0.000229701672243719, 0.000198652797104838), qvalue = c(0.00230759974232573, 
    0.00159418720995361, 0.00137568222478875, 0.00137568222478875, 
    0.000213580629738545, 0.000184710842940477), Count = c(8L, 
    10L, 18L, 18L, 20L, 21L)), row.names = c(NA, 6L), class = "data.frame")

謝謝你盡你所能的幫助！