时间流逝密谋

Question

我是R环境的新手。我使用以下代码生成了模拟时间流逝图。

seq_x<-seq(1,10)
seq_y<-function(y)
{
z<-y^2+y+1
return (c(z))
}
yrange<-seq_y(1)
yrange[2]<-seq_y(length(seq_x))
for(i in 1:length(seq_x) )
{
xdata<-seq_x[1:i]
ydata<-seq_y(xdata)
plot(xdata,ydata,xlim=range(seq_x),ylim=range(yrange),type="o",col="royalblue",plot.first=grid())
Sys.sleep(1)
}

我得到以下图（在最后一次迭代之后）。

现在我需要在图中绘制一条直线和一个圆，如下所示。

直线应随数据序列增长。圆应位于用户指定半径的数据序列的中心。任何有关这方面的建议都将受到高度赞赏。

Answer 1

试试这个：

#draw circle
require(plotrix)
userRadius <- 1
draw.circle(median(xdata), median(ydata), userRadius)

#draw line
segments(x0=xdata[1],y0=ydata[1],
         x1=xdata[length(xdata)],y1=ydata[length(ydata)])

Answer 2

您可以使用lines为地块添加其他线条。然后你只需要计算圆（或省略号）上的点，你可以绘制额外的线和圆。

# data
x <- seq(1,10)
y <- x^2 + x + 1
# function to calculate points on the ellipsis
ellipsis_fct <- function(mx, my, rx, ry){
  phi <- seq(0, 2*pi, length = 100) # change length if you need better resolution
  data.frame(x = mx + rx*sin(phi), 
             y = my + ry*cos(phi))
}
# actually calculate the points. 
circ <- ellipsis_fct(mean(range(x)), mean(range(y)), diff(range(x))/5, diff(range(y))/5)
# plotting commands
plot(x, y, xlim=range(x), ylim=range(y), type="o", col="royalblue", plot.first=grid())
lines(range(x), range(y), col = "darkred", lty = "dashed")
lines(circ, col = "orange")