我似乎无法重新生成ggplot2我之前使用下面的代码生成的图(使用).我现在收到错误消息"No stat called StatHline".有替代品吗？
data <- data.frame(
CorrectedIntensity=c(0, -0.66, -0.37, 0, -1.04, -0.38, 0, -1.89, -1.50),
Day=c("Day 1", "Day 9", "Day 5", "Day 1", "Day 9", "Day 5", "Day 1", "Day 9", "Day 5"))
library(ggplot2)
plot_data <- ggplot() +
ylim(-2.5, 0.5) +
geom_point(data=data, aes(x=Day, y=CorrectedIntensity), size=7,
colour="royalblue3", alpha=0.30) +
geom_errorbar(data=data, aes(x=Day, y=CorrectedIntensity, ymax=..y.., ymin=..y..),
stat = "hline", yintercept = "mean", width=0.3, colour="royalblue3",
size=1.25) +
stat_summary(data=data, aes(x=Day, y=CorrectedIntensity),
fun.ymax=function(i) mean(i) + sd(i),
fun.ymin=function(i) mean(i) - sd(i),
geom="errorbar", …我想从数据框的列中提取特定信息,并将其添加到同一数据框中的新列.复杂性在于某些行根本没有我要提取的信息("UniProt:"之后的6个字符),而其他行有多次出现 - 我希望这些行相应地显示,因为此列包含标识符在我的数据框中.
这是一个例子; 我从我的数据框中复制了几行Fasta.headers:
第1行:
H05C05.1c; CE43771; WBGene00019157;状态:Partially_confirmed; UniProt的:H2L0A8; protein_id:CCD72193.1;> H05C05.1a; CE37385; WBGene00019157;状态:Partially_confirmed; UniProt的:Q9TXU2; protein_id:CCD72188.1
第2行:
C02B10.5; CE16802; WBGene00015330;状态:Partially_confirmed; UniProt的:O44447; protein_id:CCD61167.1
第3行:
ZK1127.4; CE07643; WBGene00022851;状态:确认; protein_id:CCD73716.1
第4行:
T27C4.4a; CE21211; WBGene00003025;轨迹:LIN-40;状态:成熟; UniProt的:O61907; protein_id:CCD74255.1;> T27C4.4b; CE21212; WBGene00003025;轨迹:LIN-40;状态:成熟; UniProt的:Q76NP4 ; protein_id:CCD74256.1;> T27C4.4d; CE33331;> F54F2.9; CE39158; WBGene00018836;状态:成熟; UniProt的:P34454; protein_id:CCD71243.1
我希望输出为:
H2L0A8;Q9TXU2
O44447
O61907;Q76NP4;P34454