4

所以我正在尝试编写一个创建散点图的函数。

根据数据,我可能希望 V1 或 V2 的值来确定点的形状/颜色。V1 中的数据只取特定值,V2 也是如此,但两者对于 R/ggplot 看起来都是连续的,所以我对待它们as.factor。我也可能会退步购买 V1 或 V2。

这工作正常:

p1<-ggplot(sigmas.table,
           aes(x=V4,y=V5, color=as.factor(V2),shape=as.factor(V1)))+  
    geom_point(size=4)+
    stat_smooth(aes(group=as.factor(V2)),method=lm,formula=y~poly(x,1),size=1)

但是要将它包装在一个函数中,首先我必须切换到aes_string以便我可以读取列标题(我也改变 x 轴值),但现在我无法说服 ggplot V1/V2 只需要有限数量的分配不同点的值。颜色排序的作品,但分配一个连续的阴影而不是不同的颜色,所以它也将数据视为连续的:

scatter.plot<- function(s.table,x.var, col.var, shape.var, line.var){
              ploti <- ggplot(s.table,
                              aes_string(x=x.var,y="V5", color=as.factor(col.var),shape=as.factor(shape.var)))+
                       geom_point(size=4)+
                       stat_smooth(aes_string(group=as.factor(line.var)),method=lm,formula=y~poly(x,1),size=1)+
               return(ploti)}

p1 <- scatter.plot(sigmas.table,"V4","V2","V1","V2")

返回

Error: A continuous variable can not be mapped to shape

我确信有很多可用的工具,但在我看来(也许在我的天真中)好像有一个真正的解决方案。任何和所有的帮助表示赞赏。

Ps 这是 sigmas 表的示例:

     V1  V2         V3        V4         V5         V6
 1: 0.3  16 0.12791584 0.3454941 0.19463432 0.04231422
 2: 0.5  16 0.09908318 0.4460310 0.06286376 0.05462742
 3: 0.7  16 0.08374057 0.5277510 0.03820782 0.06463604
 4: 0.9  16 0.07385224 0.5984134 0.03026121 0.07329038
 5: 0.3  32 0.09045016 0.2443013 0.17003319 0.02992067
 6: 0.5  32 0.07006239 0.3153916 0.04250670 0.03862742
 7: 0.7  32 0.05921353 0.3731763 0.02563209 0.04570458
 8: 0.9  32 0.05222142 0.4231422 0.02037799 0.05182412
 9: 0.3  64 0.06395792 0.1727471 0.14072683 0.02115711
10: 0.5  64 0.04954159 0.2230155 0.03286223 0.02731371
11: 0.7  64 0.04187029 0.2638755 0.01946985 0.03231802
12: 0.9  64 0.03692612 0.2992067 0.01475039 0.03664519
13: 0.3 128 0.04522508 0.1221506 0.11012266 0.01496034
14: 0.5 128 0.03503120 0.1576958 0.02260296 0.01931371
15: 0.7 128 0.02960676 0.1865882 0.01317059 0.02285229
16: 0.9 128 0.02611071 0.2115711 0.01036970 0.02591206
4

1 回答 1

3

正如我的评论中所述,如果原始数据实际上包含因子,则应该转换它。但是,您也可以在您的函数中执行此操作:

scatter.plot<- function(s.table,x.var, col.var, shape.var, line.var){
  s.table[,col.var] <- as.factor(s.table[,col.var])
  s.table[,shape.var] <- as.factor(s.table[,shape.var])
  s.table[,line.var] <- as.factor(s.table[,line.var])

  ggplot(s.table, aes_string(x=x.var, y="V5", color=col.var, shape=shape.var)) +
    geom_point(size=4) +
    stat_smooth(aes_string(group=line.var), 
                method=lm, formula=y~poly(x,1), size=1)
    }
于 2013-09-17T11:03:28.267 回答