This is a test version of Biostars. For the public version, visit https://www.biostars.org.
boxplot with a for loop

Hi all,

I tried to make a boxplot for each classes of (data) based on repeat names (rep_name) with a for loop but it didn't give me the desired result. What it gave me was the last boxplot of the last class

Any idea on how I can have a boxplot for each class with a loop instead of doing it manually?

Thanks for any answer

for (i in length(data)){
    boxplot(data[[i]][5:8])
}

data

[[1]]
    rep_name   rep_family    gene       distance Hepatocytes_B1 Hepatocytes_B3  Huh.7_B1  Huh.7_B2
  LTR18B        LTR           NDUFAB1      572         9.3038          9.287                     10.9524    10.8579

[[2]]
     rep_name       rep_family   gene distance Hepatocytes_B1 Hepatocytes_B3 Huh.7_B1 Huh.7_B2
        LTR19B        LTR          CASP10        0         6.5334         5.3218             7.1686              5.9633
        LTR19B        LTR               CD38        0         5.8011         6.9873              5.7066             2.5126

[[3]]
      rep_name     rep_family  gene      distance Hepatocytes_B1 Hepatocytes_B3  Huh.7_B1     Huh.7_B2
     LTR4                 LTR         ACSM3        0          8.744         6.9615                      5.9441            3.1032

[[4]]
      rep_name     rep_family    gene    distance Hepatocytes_B1 Hepatocytes_B3    Huh.7_B1    Huh.7_B2
     LTR53                  LTR       TNMD    66637         0.5905         0.9610                     0.6130        0.7612
    LTR53                    LTR      MTMR7        0         3.9547         3.7264                        3.5522        4.2649
    LTR53                    LTR        VPS41     2226         7.1890         7.2522                    7.1952        7.0576
r
• 16,331 views • link
updated by Ram Ram • written by yasjas • yasjas

2 answers

Rechecking the code I think that the correct version is (using pdf export as suggested by Sam):

pdf("Image.pdf")
for(i in 1:length(data)){
  set = data[[i]]
  boxplot(as.numeric(unlist(set[5:8])))
}
dev.off()
• 0 views • link
updated by Ram Ram • written by Nicola Casiraghi Nicola Casiraghi

You can either use par, see tutorial here

Or I do prefer just open a pdf and output images there

pdf("Image.pdf")
for (i in 1:length(data)){
    boxplot(data[[i]][5:8])
}
dev.off()
• 1 views • link
updated by Ram Ram • written by Sam Sam

Further explanation on graphing in R:

When you call boxplot() (or any graphing function) in R, it draws it in a default graphic device, which it closes after you're done. Every time you call another boxplot() function, it overwrites your previous plot. What Sam suggests above is opening your own graphics device manually (in this case, a pdf), write what output you want to it, and then close it.

A couple of things to keep in mind: be sure to close the graphic device explicitly with dev.off(), and be careful with using loops to write files, since if your loop is incorrect and turns out to be an infinite loop, you will create a huge file and potentially cause issues.

• 1 views • link
updated by Ram Ram • written by Steven Lakin Steven Lakin

Thank you guys, it is more clear now and thanks for the link of the tutorial

• 0 views • link
updated by Ram Ram • written by yasjas • yasjas

The for cycle must be:

for (i in 1:length(data)){
• 1 views • link
updated by Ram Ram • written by Nicola Casiraghi Nicola Casiraghi

Oh, ya, I forgot to check his code. That also explain why he will only ever see the last boxplot.

• 0 views • link
updated by Ram Ram • written by Sam Sam

Log in to answer this question.

More posts like this
  • Remove redundant name in a list
    written by yasjas • yasjas
    ``` [[1]] rep_name rep_family gene_name Hepatocytes_B1 Hepatocytes_B3 Huh.7_B1 Huh.7_B2 HERVL40-int LTR ANKIB1 7.2268 7.3056 7.2132 7.5750 HERVL40-int LTR ANKIB1 7.2268 7.3056 7.2132 7.5750 HERVL40-int LTR …
  • ttest on a list in R
    written by yasjas • yasjas
    Hi everyone, I have a list (data) that I want to apply a t-test on each of the component ```r [[1]] rep_name rep_family gene distance …