2013-04-07から1日間の記事一覧

WikiのようなProjective Varietyのような絵(こちら)が書きたかったが… my.pr <- function(xy,x0,y0,y1){ k <- 1/(xy[,2]-y1+1) new.y <- y0 * k new.x <- x0 + (xy[,1]-x0) * k return(cbind(new.x,new.y)) } x <- y <- seq(from=-10,to=10,length=1000) xyz…

Solving the Likelihood Euations

昨日の記事の続き 論文はこちら 計算機代数統計的に尤度関数を解くのはこんな感じ 尤度関数は必ずしも凸ではなく、そうすると複数の極大値があり得る 探索すると最大値にたどり着き損ねるので、すべての微分0の点を列挙して それらを比較して最大値を得るこ…

尤度関数を解く