二次的現象なのか

  • 特定の要因がなくても発生する事象があって、ある要因があるとその事象が引き起こされるような印象がある。その印象はどんなときにどのくらいもたらされる?
  • 印象を強める要素は何だろう?

N<-20
# 単位期間にDという有害事象が起きる確率
pD<-0.1
D<-runif(N)/pD
# 単位期間にVという有害事象が起きる確率
pV<-0.4

# Vの後、probの確率で、Dが短期間に引き起こされる(有害事象が起きる)
V<-runif(N,min=0,max=1)/pV
VearlierD<-which(V<D)
# side effectのタイムラグは abs(正規分布(平均:SideEffectM,分散:SideEffectV))
SideEffectM<-0.01
SideEffectV<-0.01
# side effectとしてDが発生する割合
M <- N/2
probs<-seq(from=0,to=1,length=M)
par(ask=FALSE)
for(x in 1:length(probs)){
	prob<-probs[x]
Sideeffected<-sample(VearlierD,length(VearlierD)*prob)
D[Sideeffected]<-V[Sideeffected]+abs(rnorm(length(Sideeffected),SideEffectM,sqrt(SideEffectV)))
t<-seq(from=0,to=1,length=100)
xlim<-ylim<-c(0,1)
maintitle=x
for(i in 1:length(t)){
	if(i==length(t)){
		par(ask=TRUE)
	}else{
		par(ask=FALSE)
	}
	plot(rep(0,length(t)),t,type="l",xlim=xlim,ylim=ylim,main=maintitle)
	for(j in 1:N){
		if(t[i]<min(D[j],V[j])){
			segments(0,j/N,t[i],j/N)
		}else{
			if(D[j]<V[j]){
				segments(0,j/N,D[j],j/N)
			}else{
				segments(0,j/N,V[j],j/N)
				col<-1
				if(t[i]>D[j])col<-2
				segments(V[j],j/N,min(t[i],D[j]),j/N,col=col)
			}
		}
		if(D[j]<t[i])points(D[j],j/N,pch=15,col=2)
		if(V[j]<t[i] && V[j]<D[j])points(V[j],j/N,pch=19,col=1)
	}
}

}