5

我正在将我的情节移动到 ggplot 中。除了这个（代码来自上一个问题）之外几乎都存在：

ggplot 图应该是这样的

#Set the bet sequence and the % lines
betseq <- 0:700 #0 to 700 bets
perlin <- 0.05 #Show the +/- 5% lines on the graph

#Define a function that plots the upper and lower % limit lines
dralim <- function(stax, endx, perlin) {
  lines(stax:endx, qnorm(1-perlin)*sqrt((stax:endx)-stax))
  lines(stax:endx, qnorm(perlin)*sqrt((stax:endx)-stax))
}

#Build the plot area and draw the vertical dashed lines
plot(betseq, rep(0, length(betseq)), type="l", ylim=c(-50, 50), main="", xlab="Trial Number", ylab="Cumulative Hits")
abline(h=0)
abline(v=35, lty="dashed") #Seg 1
abline(v=185, lty="dashed") #Seg 2
abline(v=385, lty="dashed") #Seg 3
abline(v=485, lty="dashed") #Seg 4
abline(v=585, lty="dashed") #Seg 5

#Draw the % limit lines that correspond to the vertical dashed lines by calling the
#new function dralim.
dralim(0, 35, perlin) #Seg 1
dralim(36, 185, perlin) #Seg 2
dralim(186, 385, perlin) #Seg 3
dralim(386, 485, perlin) #Seg 4
dralim(486, 585, perlin) #Seg 5
dralim(586, 701, perlin) #Seg 6

我可以显示我已经走了多远（不远）：

ggplot(a, aes(x=num,y=s, colour=ss)) +geom_line() +stat_smooth(method="lm", formula="y~poly(x,2)")

我的尝试

要清楚。我正在将我的数据绘制在参考线上（上图）。底部图像显示了我的数据以及我在获取参考线方面的糟糕尝试（显然没有奏效）。

4

1 回答 1

3

您所做的是为您的数据拟合抛物线，而不是绘制先前定义的抛物线。适应你必须适应的东西并不难ggplot。

与您相同的开始（尽管betseq实际上并未在任何地方使用）

#Set the bet sequence and the % lines
betseq <- 0:700 #0 to 700 bets
perlin <- 0.05 #Show the +/- 5% lines on the graph

不要代替绘制线条的函数，而是创建一个返回geom_line您想要的 s （在列表中）的函数。aes(x=x, y=y)稍后将在声明中给出一个暗示ggplot，但这定义了构成抛物线的数据点。

#Define a function that plots the upper and lower % limit lines
dralim <- function(stax, endx, perlin) {
  c(geom_line(data = data.frame(x=stax:endx, 
                                y=qnorm(1-perlin)*sqrt((stax:endx)-stax))),
    geom_line(data = data.frame(x=stax:endx, 
                                y=qnorm(perlin)*sqrt((stax:endx)-stax))))
}

为了避免重复，定义垂直线的位置（edges），也可以用来定义抛物线的左右端点（ranges）。

edges <- data.frame(x=c(0, 35, 185, 285, 485, 585, 700))
ranges <- data.frame(left = edges$x[-nrow(edges)],
                     right = edges$x[-1] + 1)

现在构建ggplot. 有一个geom_vline可以绘制所有垂直线（因为我们在单个数据集中定义了位置）。不寻常的步骤是遍历 of 的行（索引）并使用相应的左值和右值（和）ranges调用。这将返回 的列表列表，但可以以正常方式将其添加到绘图中，并添加所有行。最后两个 scale 调用只是设置标签，如果是 y 轴，则设置范围。dralimperlingeom_lines

ggplot(mapping=aes(x=x,  y=y)) +
  geom_vline(data=edges, aes(xintercept = x), linetype="dashed") +
  lapply(seq_len(nrow(ranges)), 
         function(r) {dralim(ranges$left[r], ranges$right[r], perlin)}) +
  scale_y_continuous("Cumulative Hits", lim=c(-50,50)) +
  scale_x_continuous("Trial Number")

在此处输入图像描述

于 2012-11-20T22:10:23.067 回答