SNPs - dfdf4 005 nrowdf 1genes - pasteGENE 1lengthSNPs sep=_CMplotdfc134 plottype=m LOG10=TRUE col= c#3E0A52 #423D77#3F678B #468C8D #5FB47F #9FD55C#F9E956 chrden
这段代码的目的是使用CMplot函数绘制染色体图谱,并在图上突出显示具有显著性的单核苷酸多态性(SNP)以及它们所对应的基因。
代码的第一行定义了一个变量SNPs,用于存储在数据框df中第四列值小于(0.05 / nrow(df))的SNP的标识符。
第二行定义了一个变量genes,它是一个字符向量,包含了与SNPs对应的基因的名称。
接下来,代码调用了CMplot函数,并传入了以下参数:
- df[,c(1:3,4)]:表示用于绘图的数据框的前四列。
- plot.type="m":表示绘制染色体图谱的类型为密度图。
- LOG10=TRUE:表示在绘图时使用对数刻度。
- col= c("#3E0A52", "#423D77","#3F678B", "#468C8D", "#5FB47F", "#9FD55C","#F9E956"):表示绘图时使用的颜色向量。
- chr.den.col=c("darkgreen", "yellow", "red"):表示染色体密度图的颜色。
- highlight = SNPs:表示要突出显示的SNPs。
- highlight.col = NULL:表示突出显示的SNPs的颜色。
- highlight.cex = 1:表示突出显示的SNPs的大小。
- highlight.pch = c(15:17):表示突出显示的SNPs的形状。
- highlight.text = genes:表示突出显示的SNPs对应的基因名称。
- highlight.text.col = "black":表示突出显示的SNPs的基因名称的颜色。
- threshold = 0.05/nrow(df):表示显著性阈值。
- threshold.lty = 2:表示显著性阈值的线条类型。
- threshold.col = "red":表示显著性阈值的颜色。
- threshold.lwd = 3:表示显著性阈值的线条宽度。
- amplify = FALSE:表示是否放大图形以适应突出显示的SNPs。
- file = "jpg":表示输出图形的文件格式为jpg。
- dpi = 300:表示输出图形的分辨率为300 dpi。
- file.output = F:表示图形不直接输出到文件中。
- verbose = F:表示不显示绘图过程的详细信息。
- width = 14,height = 6:表示输出图形的宽度和高度。
最后一行的括号结束了CMplot函数的调用
原文地址: https://www.cveoy.top/t/topic/hMtj 著作权归作者所有。请勿转载和采集!