1

我试图绘制成对的点,因为这是它们出现的格式。该数据描述了商业部件(电子元件),其行为被描述为介于两个频率之间。这是一个数据集的小例子:

freq1  freq2  gain  
2.0    6.0    43
6.0   18.0    40
8.5   10.5    50 
8.5    9.3    52

因为数据是在一个频率范围内呈现的,所以我想将其显示为连接两个点的线。例如,第一行描述了两个点:(2, 43) 和 (6, 43)。

另一种可能性是计算中心频率,定位那里的点并计算我在每个点绘制的图像的宽度,使其范围为freq1freq2。这引发了很多问题,我无法开始编写代码。

所以,我的问题是:哇,我可以根据连接线和(最好使用ggplot2)绘制每个gain值吗?freq1freq2

4

3 回答 3

9

只是一个猜测,latticeExtra有一个有趣的segplot(它是一个扩展的点图)

在此处输入图像描述

读取一些数据:

  dat <- read.table(text ='freq1  freq2  gain  
2.0    6.0    43.
6.0   18.0    40.
8.5   10.5    50. 
8.5    9.3    52.',head=T)

dat$compo <- paste('compo',1:nrow(dat),sep='')
library(latticeExtra)
segplot(reorder(factor(compo), gain)~freq1+freq2,
        data=dat,draw.bands=FALSE,centers=gain,
        segments.fun = panel.arrows,ends = "both", 
        angle = 90, length = 1, scales=list(y=list(cex=1.5)),
        unit = "mm",
        main = ' Range frequencies'  ,sub= 'electronic components', 
        ## a serious theme here :)
        par.settings = theEconomist.theme(with.bg = TRUE))

编辑添加一些图片

library(png)

ll <- list.files(path=path_picts,patt='compo[0-9].*',full.names=T)
imgs <- lapply(ll,readPNG)
#ll <- gsub('.*(compo[0-9]).png','\\1',ll)
#names(imgs) <- ll
dat$compo <- paste('compo',1:nrow(dat),sep='')
segplot(factor(compo)~freq1+freq2,
        data=dat,draw.bands=FALSE,centers=gain,
        segments.fun = panel.arrows,ends = "both", 
        angle = 90, length = 1, scales=list(y=list(cex=1.5)),
        unit = "mm",
        main = ' Range frequencies'  ,sub= 'electronic components',
        par.settings = ggplot2like(),axis = axis.grid,

        panel = function(x,y,...){
            panel.segplot(x,y,...)
            browser()
            lapply(seq_along(ll),function(img){
                 x1 <- x[img];y1 <- y[img];
                  grid.raster(image=imgs[[img]],x=(x1+y1)*0.5,y=img,width=y1-x1, height=0.5,interpolate=F,
                        default.units = 'native')})

        })

在此处输入图像描述

于 2013-03-04T20:59:44.020 回答
8

这是您问题的答案

> dat
  freq1 freq2 gain
1   2.0   6.0   43
2   6.0  18.0   40
3   8.5  10.5   50
4   8.5   9.3   52

> attach(dat)

#Don't actually need to calculate the midpoint but since you suggested it
#that's the way I did it
midpoint = (freq1+freq2)/2
plot(midpoint,gain,xlim=c(min(freq1),max(freq2)),col="white",xlab="")

points(freq1,gain,col=1:length(gain),pch=19)
points(freq2,gain,col=1:length(gain),pch=19)

for(i in 1:length(gain)){
    lines(c(freq1[i],freq2[i]),c(gain[i],gain[i]),col=i)
}

在此处输入图像描述

于 2013-03-04T20:54:59.647 回答
5

对于它的价值,这里有一个想法ggplot2。我会让你在它上面建造任何东西,尽管它或多或少拥有一切。

ggplot(data = df, aes(x=freq1, y=gain)) + 
      geom_segment(aes(xend=freq2, yend=gain, colour=factor(1:4)))

在此处输入图像描述

于 2013-03-04T21:40:12.807 回答