重现频率矩阵图

Question

我有一个我想在R中重新创建的情节.这是情节:

来自:Boring,EG(1941).统计频率为动态均衡.心理学评论,48(4),279.

这略高于我的薪水(能力)因此在这里问.无聊状态:

在第一种情况下,A只能出现'从'(0)或'永远'(1).在第二种情况下,频率为0,1/2或1; 在第三个,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))
}