用R语言画火山图时报错
Error in .standalone_types_check_dot_call(ffi_standalone_check_number_1.0.7, : object 'ffi_standalone_check_number_1.0.7' not found
按照2023-03-21 R语言GGPLOT2报错 ffi_standalone_check_number_1.0.7 not found - 简书 (jianshu.com)上的方法重装tidyverse后成功
原代码:
####TCGA差异分析火山图####
setwd("C:\\Users\\86189\\Desktop\\Xena")
library(tidyverse)
exp <- read.table("LIHC_fpkm_mRNA_all.txt",sep = "\t",row.names = 1,check.names = F,stringsAsFactors = F,header = T)
DEG <- as.data.frame(res)%>%
arrange(padj) %>%
dplyr::filter(abs(log2FoldChange) > 0, padj < 0.05)
logFC_cutoff <- 1
type1 = (DEG$padj < 0.05)&(DEG$log2FoldChange < -logFC_cutoff)
type2 = (DEG$padj < 0.05)&(DEG$log2FoldChange > logFC_cutoff)
DEG$change = ifelse(type1,"DOWN",ifelse(type2,"UP","NOT"))
table(DEG$change)
#remove.packages("rlang")
#remove.packages("tidyverse")
#remove.packages("ggplot2")
#install.packages("tidyverse")
# tidyverse自动安装依赖的rlang、ggplot2
#install.packages("ggpubr")
#install.packages("ggthemes")
library(ggpubr)
library(ggthemes)
DEG$logP <- -log10(DEG$padj)
ggscatter(DEG,
x = "log2FoldChange", y = "logP") +
theme_base()
Error in .standalone_types_check_dot_call(ffi_standalone_check_number_1.0.7, :
object 'ffi_standalone_check_number_1.0.7' not found
#增加基因上下调信息
ggscatter(DEG, x = "log2FoldChange", y = "logP",
color = "change",
palette = c("blue", "black", "red"),
size = 1) +
theme_base()
#添加分界线
ggscatter(DEG, x = "log2FoldChange", y = "logP", xlab = "log2FoldChange",
ylab = "-log10(Adjust P-value)",
color = "change",
palette = c("blue", "black", "red"),
size = 1) +
theme_base() +
geom_hline(yintercept = -log10(0.05), linetype = "dashed") +
geom_vline(xintercept = c(-1, 1), linetype = "dashed")
dev.off()