Rcppを試してみる

R C

この前紹介したGoogle Summer of Codeのやつですが、メンターがDirk Eddelbuettelさんで、この人はRcppというパッケージを作ってる人で僕はこのパッケージを使ったことがなかったので、使ってみようじゃないかというお話。以下を参考にしつつ。 Rcpp: Seamle…

2010-03-01

uint8とかMiniseのコード読んだりとか

C C++

Miniseのコード読みの続き。uint8とかuint64とか出てきてなんじゃいなという感じだったので調べてみる。Minise::addFilehttp://developers.curlap.com/curl/docs/rte/6.0/ja/docs/ja/dguide/primitive.htmlによるとデータ型説明既定値サイズ(ビット) 有…

2009-02-22

俺俺lapply&sapplyを書いてみる【Rの中身はlisp編】

R C

というかWriting R Extensionsにある例をやってみているという感じなんですけど。C言語が(も)あんまり分かっていない僕ですが、それでもRのベクトルや行列のオブジェクトをCレベルで触るにはどうすればいいか段々分かってきました。こうなってくるとRの色々…

2009-02-21

Kmeans法のソースを見つつ

R 機械学習 C

Cの勉強もかねつつ。kmeans.cとkmeans.Rから最低限必要そうなところを引っ込抜いてきて、コメントを適度に埋めつつやってみました。Cのソース。100行くらい。短かい。 #include <R.h> #include <Rinternals.h> void mykmeans(double *x, int *pn, int *pp, double *cen, int *pk</rinternals.h></r.h>…

2009-02-21

Rの言語仕様とか色々参考にできそうな関数とかのメモ

R C

なんというS式R Language Deﬁnitionのp3のシンボルオブジェクト。 > as.list(quote(x + y)) [[1]] `+` [[2]] x [[3]] y > as.list(quote((x + y) * z)) [[1]] `*` [[2]] (x + y) [[3]] z モデル式言語オブジェクト R 言語を構成する三つのオブジェクトの型…

2009-02-19

Confidence interval for kernel density estimator with optimal bandwidth under the cross-validation criterion

R C

ここの続き。カーネル密度推定の場合においての最適なバンド幅をクロスバリデーションで求める、というのはまだやっていなかったのでやりました。カーネル密度推定の場合のクロスバリデーションは、積分の計算が入ってくるので、数値積分をやったりしました…

2009-02-05

Rのソースコードを読むために

R C

gtagsとかを入れてみた。関数の定義に飛んだりすることができるなどする。 GNU GLOBAL（gtags）ソースコードタグシステムの使い方

2009-01-19

RのコードをCに書き変えてみたら1500倍くらい早くなったけど、早すぎて将来が不安

C R

昨日の続きのところ。Cでナダラヤワトン推定量*1計算させて、Rのoptimを使って計算したら早くなったっぽい気がした。でも、一変数はRだけでまともなコードがなかった(文字列evalとかのやつしか残っていなかった)ので速度比較ができなかった。これじゃ最適化…

2009-01-18

クロスバリデーション最小化のコード関係をCで書いてみた

R C

Cで書いてどうなるかを見てみるテスト。やってみて気がついたが、サンプル数が少ないとき、epanechnikovだと分母が0の領域が出てきて微妙な感じになってしまうのでどうにかしないといけないが、面倒なので放置している。。。サンプル少ない時はガウシアンカ…

2009-01-18

Cの勉強などなど

C

Rを読む時に必要そうなところを中心にピックアップ。C言語入門 (ASCII SOFTWARE SCIENCE Language)を参考にしつつ。構造体と配列なんちゃら型の配列というのができるらしい。Javaの型付きのArrayListみたいな感じか。 #include <stdio.h> struct pencil { int hardnes</stdio.h>…

2009-01-17

Cを使ってRの拡張をする関係の参考リンク

R C

Rのソースコードを読むオフなどしたいなどする。Cに詳しい人がいるとよさそうですか、分かりません。 Introduction to Calling C and Fortran from R ppt注意暇を見つつ日本語訳などしたい R-C-interface なかなか分かりやすそう読みたいソース R-2.8.0/sr…