首页 > 解决方案 > How to make inactive multiple line R code in a script

问题描述

I do have a long R script and I don't want to delete a part of this R script. Is there any way to avoid running this multiple line R code in the script or should I need to insert # at the beginning of each of these lines?

I don't want to run this whole below code but want to run the code which follow it?

   alpha = 0.01
    sigtabresTreatment6Genus_T1_vs_T3 = resTreatment6Genus_T1_vs_T3[which(resTreatment6Genus_T1_vs_T3$padj < alpha), ]
    sigtabresTreatment6Genus_T1_vs_T3 = cbind(as(sigtabresTreatment6Genus_T1_vs_T3, "data.frame"), as(tax_table(physeq6Genus)[rownames(sigtabresTreatment6Genus_T1_vs_T3), ], "matrix"))
    head(sigtabresTreatment6Genus_T1_vs_T3)
    dim(sigtabresTreatment6Genus_T1_vs_T3)

    library("ggplot2")
    theme_set(theme_bw())
    scale_fill_discrete <- function(palname = "Set1", ...) {
        scale_fill_brewer(palette = palname, ...)
    }
    # Phylum order
    x = tapply(sigtabresTreatment6Genus_T1_vs_T3$log2FoldChange, sigtabresTreatment6Genus_T1_vs_T3$Phylum, function(x) max(x))
    x = sort(x, TRUE)
    sigtabresTreatment6Genus_T1_vs_T3$Phylum = factor(as.character(sigtabresTreatment6Genus_T1_vs_T3$Phylum), levels=names(x))
    # Genus order
    x = tapply(sigtabresTreatment6Genus_T1_vs_T3$log2FoldChange, sigtabresTreatment6Genus_T1_vs_T3$Genus, function(x) max(x))
    x = sort(x, TRUE)
    sigtabresTreatment6Genus_T1_vs_T3$Genus = factor(as.character(sigtabresTreatment6Genus_T1_vs_T3$Genus), levels=names(x))

    pdf("DESeq2-Soil-batch1-DEG-T1_vs_T3_pval.01.pdf", width = 8, height = 6);
    ggplot(sigtabresTreatment6Genus_T1_vs_T3, aes(x=Genus, y=log2FoldChange, color=Phylum)) + geom_point(size=6)+theme(axis.text.x = element_text(angle = -90, hjust = 0, vjust=0.5))
    dev.off();

    write.csv(sigtabresTreatment6Genus_T1_vs_T3,"sigResult-Treatment6Genus_T1_vs_T3")

###############
Further R code that needs to run

###########

Many thanks Yogesh

标签: r

解决方案


推荐阅读