タグ

関連タグで絞り込む (2)

タグの絞り込みを解除

rと線形代数に関するrishidaのブックマーク (2)

  • アルマジロを可愛がることにした - サボタージュ禁止のおさぼり日記

    Armadillo: C++ linear algebra library あっ, どーも僕です. 数あるRcpp系パッケージのなかでRcppArmadilloを少し勉強したのでまとめます. RcppArmadilloによりRcpp経由で, C++ライブラリのArmadilloを使用できます. 以下にRcppArmadilloの機能をまとめます. ただし, 内容はRcppArmadilloの機能に絞ったものです. また, 複素や疎行列に関しては少し省いています. アルマジロその1. 行列操作編 - サボタージュ禁止のおさぼり日記 アルマジロその2. 演算編 - サボタージュ禁止のおさぼり日記 アルマジロその3. 線形代数編 - サボタージュ禁止のおさぼり日記 Rcppそのものの使い方は下に載せた参考文献にあたってください. なお, 私も勉強中ですので間違いに気づいていない場合があります.

    アルマジロを可愛がることにした - サボタージュ禁止のおさぼり日記
  • Rcppのすすめ

    6. gibbsR <- function(N,thin){ mat<-matrix(0,nrow=N,ncol=2) x <- 0 y <- 0 for(i in 1:N){ for(j in 1:thin){ x <- rgamma(1, 3, 1/(y*y+4)) y <- rnorm(1, 1/(x+1), 1/sqrt(2*x+2)) } mat[i,] <- c(x,y) } return(mat) } Rのコード例 ギブスサンプラー ついでに cmpfun でコンパイルしてみる gibbsC <- compiler::cmpfun(gibbsR) 7. Rcppのコード例 #include <Rcpp.h> using namespace Rcpp; ←これを書くとRcpp::は書かなくて良い // [[Rcpp::export]] ←この直下の関数がRに読み込まれる Rcp

    Rcppのすすめ
  • 1