スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

2度目のZipf-Mandelbrot分布

今日は休んで、午前は反応拡散パターン探索、午後からは、あるデータの解析を続けていた。

$$

修正ZM分布版(完成。稀なデータでズレ)
setwd("F:/Seimei")
data <- read.table("seimei-08Feb2014.csv", header=T, sep=",")
data[,2] <- as.numeric(gsub("位","",data[,2]))
x <- c(1:length(data[,1]))
y <- data[,3]/sum(data[,3])
resid <- function(par){
xx <- c(1:1000000)
C <- sum((par[3]^x)/(xx+par[2])^par[1])
yhat <- (par[3]^x)/((x+par[2])^par[1])/C
sum((y-yhat)^2) # 残差平方和を返す
}
result <- optim(c(1,1,1),resid,method="Nelder-Mead")
C <- sum((result$par[3]^x)/(x+result$par[2])^result$par[1])
plot(x, (result$par[3]^x)/((x+result$par[2])^result$par[1])/C)
points(x, y, col="red")
ey <- (result$par[3]^x)/((x+result$par[2])^result$par[1])/C
plot(log(x), log(ey))
points(log(x), log(y), col="red")

> pdf("seimei-distribution_26Feb2014.pdf", height=7, width=10, pointsize=14)
> plot(log(x), log(ey))
> points(log(x), log(y), col="red")
> dev.off()

seimei-distribution_26Feb2014.png


岸辺露伴は動かない (ジャンプコミックス)岸辺露伴は動かない (ジャンプコミックス)
(2013/11/19)
荒木 飛呂彦

スポンサーサイト

コメントの投稿


非公開コメント

コメント

Profile

zoa

Author:zoa
Metropolis (1927)

Calendar
05 | 2017/06 | 07
- - - - 1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30 -
Labels
Admin
Previous Posts
Recent Comments
Recent Trackback
Archive
このページのトップへ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。