重现频率矩阵图

Question

我想在R中重新创建一个情节。这是情节：

From：Boring，E。G.（1941）。统计频率为动态均衡。心理学评论，48（4），279。

这略高于我的薪水（能力）因此在这里问。无聊状态：

  

在第一种情况下，A只能出现'从'（0）或'永远'（1）。上   第二种情况是频率   是0,1 / 2，或1;在第三个0,1 / 3，   2/3，或1等等。

显然，您不必担心标签等。只是提示生成数据以及如何绘制图片会很棒。 ;）我不知道如何开始......

Answer 1

这是一个例子：

library(plyr)
ps <- ldply(1:36, function(i)data.frame(s=0:i, n=i))
plot.new()
plot.window(c(1,36), c(0,1))
apply(ps, 1, function(x){
  s<-x[1]; n<-x[2];
  lines(c(n, n+1, n, n+1), c(s/n, s/(n+1), s/n, (s+1)/(n+1)), type="o")})
axis(1)
axis(2)

ps代表所有点。每个点都有两个孩子。 因此，从每个点到孩子们画线。

Answer 2

使用基本图形的解决方案：

x <- 1:36
boring <- function(x, n=1)n/(x+n-1)

plot(x, boring(x), type="l", usr=c(0, 36, 0, 1))
for(i in 1:36){
  lines(tail(x, 36-i+1), head(boring(x, i), 36-i+1), type="o", cex=0.5)
  lines(tail(x, 36-i+1), 1-head(boring(x, i), 36-i+1, type="o", cex=0.5))
}