読者です 読者をやめる 読者になる 読者になる

INFINITY -数学とかプログラミングとか-

統計とプログラムを使って役に立たせたい

TeX用コマンド入力を支援するための辞書をご利用ください。
sanctuary's blogは,適当なことが書いてあります。

ラプラス分布

#指数べき乗族のvを0.4から1までの図(phi=1)。0.4からにしたのは、0.1~0.3では上に突き抜けすぎてしまったため。

x=seq(-2,2,length=50)
y=seq(-1,2,length=50)

phi=1
for(i in 4:9){
	v=0.1*i
	density1=function(x) (1/(2*v^(1/v-1)*phi^(1/v)*gamma(1/v)))*exp(-abs(x)^v/(v*phi))
	plot(density1,xlab="x",ylab="",xlim=c(-2,2),ylim=c(0,1.5))
	par(new=T)
}