将数据点添加到现有的散点图

Question

我有一个现有的ggplot2散点图，它显示了一个参数与普通数据库的结果。然后我想添加两个额外的点到这个图，我将作为命令行参数传递给我的脚本script age value1 value2。我想在每个点上面用r和l geom_text显示这些点为红色。我有以下的代码，到目前为止，但不知道如何添加临门一脚

pkgLoad <- function(x) 
    { 
    if (!require(x,character.only = TRUE)) 
    { 
     install.packages(x,dep=TRUE, repos='http://star-www.st-andrews.ac.uk/cran/') 
     if(!require(x,character.only = TRUE)) stop("Package not found") 
    } 

    } 

pkgLoad("ggplot2") 
#load current normals database 
df<-data.frame(read.csv("dat_normals.txt", sep='\t', header=T)) 

args<-commandArgs(TRUE) 

#specify what each argument is 
age <- args[1] 
rSBR <- args[2] 
lSBR <- args[3] 

# RUN REGRESSION AND APPEND PREDICTION INTERVALS 
lm_fit = lm(SBR ~ Age, data = df) 
sbr_with_pred = data.frame(df, predict(lm_fit, interval='prediction')) 




p <- ggplot(sbr_with_pred, aes(x=Age, y=SBR)) + 
      geom_point(shape=19, alpha=1/4) + 
      geom_smooth(method = 'lm', aes(fill = 'confidence'), alpha = 0.5) + 
      geom_ribbon(aes(y = fit, ymin = lwr, ymax = upr, 
         fill = 'prediction'), alpha = 0.2) + 
      scale_fill_manual('Interval', values = c('green', 'blue')) + 
      theme_bw() + 
      theme(legend.position = "none") 


ggsave(filename=paste("/home/data/wolf/FV_DAT/dat_results.png",sep="")) 
browseURL(paste("/home/data/wolf/FV_DAT/dat_results.png",sep"")) 

从本质上讲，我想看看2个新的点落在95％的置信区间内，从正常的数据库（蓝带）

Answer 1

您的示例不可重现。创建数据和可重复的示例非常有建设性。这不是浪费时间。对于解决方案，我写了评论中所说的内容。用新数据添加一个新图层。

newdata <- data.frame(Age = args[1], 
         SBR = c(args[2],args[3])) 
    p + geom_point(data=newdata,colour="red",size=10) 

例如：

sbr_with_pred < - data.frame（年龄=样品（15：36,50，代表= T）， SBR = RNORM（50））

p <- ggplot(sbr_with_pred, aes(x=Age, y=SBR)) + 
    geom_point(shape=19, alpha=1/4) + 
    geom_smooth(method = 'lm', aes(fill = 'confidence'), alpha = 0.5) 

    args <- c(20,rnorm(1),rnorm(2)) 
    newdata <- data.frame(Age = args[1], 
         SBR = c(args[2],args[3])) 
    p + geom_point(data=newdata,colour="red",size=10)