Rでコロモゴロフ-スミルノフ検定 (Kolmogrov-Smirnov test)の実装する。 比較対象がべき分布fの場合 組み込み関数ks.testを使うと #べき分布乱数の発生 h1c#べき分布の確率関数の定義 alpha xminf 1-xmin^alpha/(x)^(alpha) }test#表示 test test$statistic …

ストリーミング配信のしくみ音声や映像をすべてダウンロードしなくても、再生できる。メタファイル（たとえば、asx拡張子で示される）にストリーミングされる音声や映像 ファイル(mp3やwmvなどの拡張子)などのリンクを 記述する。ストリーミングされる音声等…

866 184 sort 列 72

To find an estimate for the standard error on ˆα in the discrete case, The MLE for the case where x is a discrete integer variable is less straightforward.

5分ニュース NPR(アメリカ国営放送) Five minituire of NPR news.http://www.npr.org/templates/player/mediaPlayer.html?action=1&t=4&islist=false CNN News update http://edition.cnn.com/services/podcasting/popups/cnn.news.update.html 参考サイト h…

べき分布 の指数の最尤推定量 (Hill estimator) } たぶん、こんな感じ。最尤推定量は漸近正規性をもつ。 その標準偏差は、となる。

NHKスペシャル 海軍反省会 http://www.nhk.or.jp/special/onair/090809.html責任を下に押しつける制度。 (ETV特集 神聖喜劇) ETVの神聖喜劇の映像 もうみれないのかなぁ。 週刊ニュース新書 信州 上田の無言館 http://www.city.ueda.nagano.jp/hp/kanko/muse…

670 132 sort 列 72

579 約60 sort 列 12

gnuplotべき乗は"**"で表現する。 例えば、xの5乗は、 plot x**5 とする。

sort並び替えの応用編2列目でまず並び替え3列目を数字で逆順に並び替える。gオプション、またはnオプション数字並び替え。rオプション、逆順に並び替え。 sort -k2,2 -k3,3gr [Filename]

gnuplotで凡例などに記号を使う。一つの方法は、 まず、 set term postscript enhanced color でpostscriptをターミナルに指定。次に、 式を書く。 上付きは"^"の記号は10^{32}。 下付きは"_"の記号は10_{32}。 記号は、{/text=\264}等記号で指定する（これ…

上場企業売上ランキング http://www.ullet.com/search.htmlこれによると、 日本の大企業の売上は約10兆円のオーダー。ちなみに、1兆円は10^12 10兆円は10^13 つまり、日本超大企業の 売上のオーダーは、10^13円。 英語では、 1兆はtrillion だそうです。参考…

gnuplotのTIPShttp://wwwnucl.ph.tsukuba.ac.jp/contents/member/inakura/gnuplot/gnuplot.html

実践女子大学図書館の新聞のレファレンス とインターネット上の情報源。図書館員の方はいい仕事をされています。 実践女子大学図書館HP http://www.jissen.ac.jp/library/frame/index.htm 外国の新聞のレファレンス http://www.jissen.ac.jp/library/frame/c…

2変量の関係がべき関係の ときの説明 (a) The relationship between the average product size and the number of products of the firm. The log-log plot of h»(K)i vs. K shows power law dependence h»(K)i » K0.38. (b) The relationship between the m…

gnuplotでの凡例の場所変更と ポイントサイズの変更方法。 凡例の移動は set key left や set key 10,20 などを使う。 ポイントサイズの移動は、 set pointsize 2 などを使う。例 set term postscript eps enhanced color "Helvetica" 30 set format y "10^{…

Rでグラフに凡例をつける例。 凡例はlegend関数を用いる。 legend(x座標,y座標,pch=c(1,2,-1,-1),col=c(1,2,3,4),lty=c(-1,-1,1,2),c("凡例1","凡例2","凡例3","凡例4")) 凡例の周りの枠を消すときは、 box.lty=0 オプションを用いる。 cexは文字の大きさを…

NHKスペシャル NHK-SP ドラマ 「気骨の判決」 裁判を公開にして、 周りから脅す。 →非公開 ETV特集 シリーズ 戦争とラジオ 第1回 放送は国民に何を伝えたのか http://www.nhk.or.jp/etv21c/update/2009/0816.html放送は、政治機関。 職員教育も。アナウンス…

8/14 368

tmp1 1 a 2 b 3 c 4 dtmp2 1 x 2 y 2 z 3 z 3 y 4 xout.txt 1 x a 2 y b 2 z b 3 z c 3 y c 4 x dを一列目のkeyで上のように併合する。 awk '{print $0,"a"}' tmp1 > tmp1b awk '{print $0,"b"}' tmp2 > tmp2bcat tmp1b tmp2b > tmp3 sort -k -k1.1n -k3 tmp…

hashクラスを使用します。 a=hash.new()

3行目の数字はほかのデータと一致する。 The figure in line three corresponds with the other data.

The researchers couldn't replicate the results of their first experiment. その研究者たちは、最初の実験の結果を再現できなかった。

rubyのbreakは内側のループを 抜けられる。

8/13 321

rubyのStringクラスのsplit関数 (String#splitを書くそうだ..。) セパレートした文字列を区切る関数カンマ区切りでは a="1,2,3,4,5" b=a.split(/\s*,\s*/) puts b[3] 関数内の正規表現の部分は、\sが空白文字。[ \t\n\r\f] と同じ。 「コメジルシ」が0回以上…

uniqは共通する要素数を数えられて 便利。uniqは連続する共通する一行を まとめるコマンド。共通する行の数を数えることもできる。データを並び替えた後 uniq -c をすれば共通する要素の数を 数えてくれる。一列のデータの場合は次のコマンド でOK sort [fil…

一列目の共通する要素数を数える awk。 awk '{a[$1]=a[$1]+1}END{for(i in a){print i,a[i];}}' netkane.csv 同じようなことは、 sort notkane.csv | uniq -c でもできる。

awkでヒストグラム暫定版。 まだこれではうまくいきません。awk '{i=int(log($2)/0.2); a[i]=a[i]+1;}END{for(i in a){print 0.1*i,a[i]}}' hi .txt > tmp