blog：http://d.hatena.ne.jp/infinity_th4/ skype：ryu_th4_dai mail：infinity_ar_net_08_thアットyahoo.co.jp（アットの部分は＠） mixi：http://mixi.jp/show_profile.pl?id=35636856&from=navi twitter：http://twitter.com/#!/tashirohiro4 本棚：http…

指数べき分布,(別名、誤差分布、他にもある） In[16]:= f[x_] := Exp[-Abs[x]^v/(v)] In[20]:= g[v_] := Integrate[f[x_], {x, -Infinity, Infinity}] In[23]:= Solve[g[v_] == 1, v] 無理やな。

が一般化ガウス族または指数べき乗族と言われる族に属する密度関数というらしい。(まだ、この正体が不明なので調査中) の辺が理解不能なので、試しに、 の図を描いてみた。vを0.1~1まで0.1刻みで。 x=seq(-2,2,length=50) y=seq(-1,2,length=50) for(i in 1:…

分布が対称なら歪度が0であることの証明を長い時間考えていたが、簡単だったな。尖度の最小値が−2であることの証明は、多分出来ると思うけど。pf.よりおｋ。シュワルツの不等式をつかうと、 とおいてシュワルツの不等式を適用する。

plotが変。あとで書こう。 x=seq(-15,15,length=500) lambda=0.8 laplace=function(x){ (lambda/2)*exp(-lambda*abs(x)) } plot(x,laplace(x),type='l',ylim=c(0,0.4)) par(new=T) lambda=0.6 plot(x,laplace(x),type='l',ylim=c(0,0.4)) par(new=T) lambda=…

そういえば、kurtosisの加法性の計算がおかしいと思ったら、平均が0だということを忘れてました。 そして、白色化されていたら、2次モーメントは1ですね。やってみましょう。xとyが独立なとき、(今、白色性は仮定しません) が成立する。ってやつなんですけど…

#共分散が0 x=rnorm(1000) y=rnorm(1000) plot(x,y) cov(x,y) #-0.004686349 共分散というのは、線形関係が見られるかどうかを見ているようです。例えば、原点を通る2次曲線上のある点と対称となる点も考え、共分散を計算すると0になります。これは、曲線関…

pysqlite2をダウンロード： http://code.google.com/p/pysqlite/downloads/listからpysqlite-2.6.3.tar.gzをダウンロードし、 tar xvf pysqlite-2.6.3.tar.gz cd pysqlite-2.6.3 python setup.py build_static install としたら、以下のエラーが起こった。 ~…

x=seq(-15,15,length=50) y=seq(-15,15,length=50) pdf=function(x,y){ return((1/(2*pi*70))*exp(-(1/2)*(((x-(-1))^2)/49))*exp((-1/2)*(((y-3)^2)/10))) } z=matrix(0,50,50) #z for(i in 1:50){ for(j in 1:50){ z[i,j]=pdf(x[i],y[j]) } } persp(x,y,z,…

アンドロイド タブレット PC MOMO9 加強版 Android 4.0 AllWinner A10 タブレット 7インチ CPU 1.5GHz ANDROID TABLET 静電式出版社/メーカー: MOMO9メディア: エレクトロニクス購入: 9人 クリック: 366回この商品を含むブログ (12件) を見るを購入した(だい…

#pdf x=seq(-0.5,2.5,length=50) y=seq(-0.5,1.5,length=50) pdf=function(x,y){ if(((0<=x)&&(x<=2))&&((0<=y)&&(y<=1))){ (3/7)*(2-x)*(x+y) }else{ return(0) } } z=matrix(0,50,50) z for(i in 1:50){ for(j in 1:50){ z[i,j]=pdf(x[i],y[j]) } } persp(…

theta=3.14/6 mm <- matrix(c(cos(theta),-sin(theta), sin(theta),cos(theta)),2,2) mm=t(mm) n=200 x=runif(n) y=runif(n) xt=matrix(0,ncol=n,nrow=2) for(i in 1:n){ for(p in 1:2){ for(q in 1:2){ xt[,i]=mm%*%c(x[i],y[i]) } } } plot(x,y,xlim=c(-1.…

#平均0,分散2のLaplace distribution phi=1 mu=0 x=seq(-4,4,,100) #x curve((1/(2*phi))*exp(-abs(x-mu)/phi),xlim=c(-4,4),ylim=c(0,0.5),ylab="Laplace distribution") par(new=T) plot(dnorm, -4, 4, type="l",xlim=c(-4,4),ylim=c(0,0.5),lty=2,ylab=""…

とは、との商空間である、つまりである。とって(;・∀・) ﾅﾝ！ (; ∀・)・ ﾃﾞｽ!! (; ∀ )・・ ｶｰ!!!まず、について。上の複素数値関数がp乗可積分であるとは、 が可測でかつ、を満たすことをいう。上のp乗可積分関数の全体の集合をと表す。 次に、について。こ…

:Hilbert 空間 線形汎関数とは,からへの線形作用素のこと。:有界線形汎関数からなるベクトル空間(双対空間)このとき, s.t. 関数を内積で計算できることが,嬉しい(´∀｀) ？記憶が曖昧だ。

grub rescue で解決できそうにないので、また、HDDをコピーすることにした。 新HDDを1日しか使ってないので、めんどくさいことにはならないでしょう。多分。コピーしようと思ったが、EASEUS Partition Masterがdid not find pending operation partitionと怒…

#Copyright(c) Hiroshi Tashiro data<-source("なんとか.なんとか")$value data data=t(data) data # #相関行列 data.cor=cor(t(data)) data.cor X=data X Xvar=function(X){ return((1/n)*rowSums(X))#各行の和 } xvar=Xvar(X) xvar#標本平均ベクトル n=nco…

dat[,1]が名前として与えられているときに、datの行の名前を書き換えるときは、以下のようにする。 for(i in 1:392){ rownames(dat)[i]=sprintf("%03d%s",i,dat[i,1]) } rownames(dat) 他に方法があるのかもしれないが、一例。

\documentclass[a4paper]{jarticle} \setlength{\topmargin}{-0.3in} \setlength{\oddsidemargin}{0pt} \setlength{\evensidemargin}{0pt} \setlength{\textheight}{46\baselineskip} \setlength{\textwidth}{47zw}\def\vect#1{\mbox{\boldmath $#1$}}\usepa…

\documentclass[a4paper]{jarticle} \setlength{\topmargin}{-0.3in} \setlength{\oddsidemargin}{0pt} \setlength{\evensidemargin}{0pt} \setlength{\textheight}{46\baselineskip} \setlength{\textwidth}{47zw}\def\vect#1{\mbox{\boldmath $#1$}}\usepa…

ノルム： norm=function(x){ return(sqrt(t(x)%*%x)) } 線形カーネル(標準内積) linearKernel=function(x1,x2){ return (t(x1)%*%x2) } x1=c(1,2) x2=c(3,4) linearKernel(x1,x2) >11 正規化線形カーネル (normalized linear kernel) normalizedLinearKernel…

テスト(´∀｀)！！

//再帰関数( recursive function) #include <stdio.h> int cnt=0; void recursive(int); void main(){ int i=0; recursive(i); } void recursive(int i){ printf("i:%d\n",i); i++; if(cnt<100){ cnt++;//recursive(i)の前じゃないと止まらない。 recursive(i); } }</stdio.h>

状態空間グラフの隣接行列による表現。今の場合は5×5のmatrixの場合だけど、 10×10とかの行列を手打ちするのは、ｸｿ ﾏﾝﾄﾞｸｾ('A`)ってなるので、描いてみた。(持ち前のｸｿみたいな忍耐力で手打ちしようと思ったが、ｸｿみたいな面倒臭さが勝ってやめた)この行列は…

Google アラート - クローラ 受信トレイ x Google アラート googlealerts-noreply@google.com 15:47 (1時間前)To 自分 ウェブ 「クローラ」に関する 3 件の新しい検索結果 かわいいアクアガール クローラ グリッターパンプス36 かわいい☆アクアガール クロー…

inprd=function(x,y){ return(t(x)%*%y)#x,y must be col vectors. } inprd(x,y) norm=function(x){ return(sqrt(t(x)%*%x)) } norm(x)