4/25のセミナーの続き。
記事が前後したがスクリプトはこちらを使う。
3人の取引を考える。
N<- 100 M<- 100 team<- 3 cycles<- 1000*5 ylim<- c(0,N*M) Res<- lapply(mapply(rep,0,1:3),unique) for(i in 1:3){ Res[[i]]<- trademodel(N,M,team,cycles,i) } par(mfcol=c(2,3)) for(pl in 1:3){ # 3種類の取引を行う。 matplot(Res[[pl]],type="l",pch=1,col=rainbow(N),ylim=ylim,frame=FALSE, xlab="trade cycle",ylab="property",main="money trade") hist(Res[[pl]][cycles+1,],breaks=N,col=5,ylim=c(0,100), xlab="money",ylab="frequency",main="money distribution") }
左から、フリマ、詐欺師、離婚モデル。
3人で取引したら、資産を独占しないのか、と思って次の日
