2013-10-02から1日間の記事一覧
重みつきスコア法によるポアソンパラメータのロバスト推定 参考（べき分布の場合） http://d.hatena.ne.jp/arupaka-_-arupaka/20130616/1371343477 #Saiyuuhou z<-rpois(10000,lambda=5) f<-function(lambda){ -sum(log(dpois(z,lambda))) } optimize(f,uppe…
c1 #c1 #c1 #qcsv("data/sikousakugo/siokouji_s.csv")[,3] n z x #x x1 yy #yy #y1 #plot(x1,yy) knot f12 #y# plot(x,y) h len v m for(i in 1:(len-1)){ if(i==1){ m[i,] } if(i>1 & i*1/6/diff(x_knots) d c f sum for(i in 1:(length(x_knots)-1)){ sum…
c1 #c1 #c1 #qcsv("data/sikousakugo/siokouji_s.csv")[,3] n z x #x x1 yy #yy #y1 #plot(x1,yy) knot f12 #y# plot(x,y) h len v m for(i in 1:(len-1)){ if(i==1){ m[i,] } if(i>1 & i*1/6/diff(x_knots) d c f sum for(i in 1:(length(x_knots)-1)){ sum…
