スポンサーサイト

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

一様乱数のヒストグラム

class <- 600
data <- hist(trunc(runif(344)*class),breaks=seq(0,class,1))[[2]]
data <- data[data!=0]
plot(sort(data, decreasing=T))
plot(log(seq(1,length(data)+1,1)),log(sort(data, decreasing=T)))
スポンサーサイト

てめェはパパ!!

この記事を閲覧するにはパスワードが必要です
パスワード入力

憲法 第29条

憲法 第29条

(1)財産権は、これを侵してはならない。
(2)財産権の内容は、公共の福祉に適合するやうに、法律でこれを定める。
(3)私有財産は、正当な補償の下に、これを公共のために用ひることができる。

不法行為による損害賠償

民法 第709条

故意又は過失によって他人の権利又は法律上保護される利益を侵害した者は、これによって生じた損害を賠償する責任を負う。

w-oN

・・・・やはり。確認した。○△×○△×に従う○△×から○△×○△×は○△×○△×○△×らしき○△×に○△×(○△×○△×○△×○△×○△×)。
つまり、○△×からの乱数生成法は得られているから(○△×について考えなくてよい)、パラメータαとBについて最適化すればよい。


ZMdistribution_21Sep2013-mini.png


方法:

α, B定義

Zipf分布からの乱数生成(340個)

観測データの尤度推定(Kolmogorov-Smirnov検定、あるいはAnderson-Darling検定でも使えばよい。AD検定の場合は最尤推定したZM分布を使うこと)・・・複数回繰り返し、平均を評価する

α, B変更(最初に戻る)

ファイル操作

ls機能

> data <- list.files()
> data
[1] "Chao0_18Sep2013-mini.png" "Chao0_18Sep2013.pdf"
[3] "Chao0_18Sep2013.png"
> data[1]
[1] "Chao0_18Sep2013-mini.png"

dir()も同じ効果。

list.dirs()


cd ../

setwd("../")

Chao 0

f<-function(x){
x*(1-exp(-200/x))
}
x <- c(1:110000)
pdf("Chao0_18Sep2013.pdf",height=5.5,width=6,pointsize=14)
plot(f(x),x,type="l",xlim=c(200,150),ylim=c(0,11000),lwd = 2,xlab="D",ylab="N")
abline(v=seq(200,150,-2.5),col="#0006E560",lty=2)
dev.off()

Chao0_18Sep2013-mini.png

If you couldn't judge, please do nothing.

rglのtriangles3dでポリゴンメッシュを復元する方法:

forループでは無間地獄となる。

.plyファイルについて:



segments3dの操作方法


国境なき記者団 ⇔ 法人格なき社団等  ← 似てる!!(特6条)

痔のくすり
Profile

zoa

Author:zoa
Metropolis (1927)

Calendar
08 | 2013/09 | 10
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ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。