R 再现频率矩阵图

我想在R中重新创建一个绘图。以下是绘图：

来源：无聊，例如（1941年）。统计频率作为动态平衡。《心理学评论》，48（4），279

这比我的工资等级（能力）稍微高一点，因此在这里提问。令人厌烦的国家：

在第一种情况下，只能出现“从不”（0）或“始终”（1）。在…上 第二次是频率 为0,1/2或1；在第三个0，1/3， 2/3或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)

库（plyr）
ps使用基本图形的解决方案：

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))
}

xniiice。这给了我一些可以玩的东西。非常感谢。科什克，你介意在你的代码中添加一些解释性的注释来解释你在这里做什么吗？@koshke：我用你的代码作为发问者后续问题的起点，这也很好。我将调查这两个答案。谢谢