http://virsalus.hatenablog.com/entry/2013/11/02/125200やはりipythonのpylab=inlineモードでは使えなかった。

http://stackoverflow.com/questions/9906938/change-ipython-qtconsole-settings ipython profile createで設定ファイルが~/.ipython/profile_defaultに生成されるpython_qtconsole_config.py c.ZMQInteractiveShell.colors = 'Linux'でカラースキームを設…

sudo pip install ipython sudo pip install spydersipとPyQt4のインストールが必要だったのでインストール sipのインストール ソースから python configure.py make sudo make installPyQt4のインストール ソースから configure-ng.pyは関係無かった python…

まだうまくいってない lengthとかnrowとかがこんがらがってきた途中のコード window=20; update = 0.5; df <- data.frame(); #kalman parameter dt <- matrix(0) ct <- matrix(0) Zt <- matrix(1) Tt <- matrix(1) a0 <- 0 P0 <- matrix(0.01); alpha <- cum…

ここを参考に http://d.hatena.ne.jp/teramonagi/20110508/1304821958真の状態も合わせてプロットしてみた library(FKF) ## alpha[t+1] = alpha[t] + eta[t], eta[t] ~ N(0, 2) ## y[t] = alpha[t] + eps[t], eps[t] ~ N(0, 1) alpha <- cumsum(2.0 * rnorm(…

ポイントはArduinoから取得した値をRpy2でRにアサインした後、numericに変換する所Pythonコード from rpy2.robjects import r import serial #global parameters window = 20 update = 0.5 r.assign('window',window); r.assign('update',update) r(''' df <…

sourceコマンドを使う。irisplot.R library(ggplot2) p <- ggplot(iris,aes(x=Sepal.Length,y=Sepal.Width,color=Species))+geom_point() R > source("irisplot.R") > p

http://blog.yhathq.com/posts/ggplot-for-python.html sudo pip install pandas sudo pip install patsy sudo pip install ggplotfrom ggplot import * print ggplot(mtcars, aes('mpg', 'qsec')) + \ geom_point(colour='steelblue') + \ scale_x_continuo…

RPy2を使う。 sudo pip install rp2#/usr/bin/env python #-*- coding:utf-8 -*- from rpy2.robjects import r r.assign('a',1) p = r('a') print p

http://r.789695.n4.nabble.com/Serial-connections-td2017474.html pipeを使ってファイルを読めるみたいだけど、使い方がわからない。 f <- pipe("ttys /dev/tty.usbmodemfd121") > f description class "ttys /dev/tty.usbmodemfd121" "pipe" mode text "r…

Gaussianライブラリを使う https://github.com/ivanseidel/Gaussian #include <Gaussian.h> double val; void setup(){ Serial.begin(9600); Serial.println("\n\nStarting Random Gaussian Distribuition..."); delay(20); randomSeed(analogRead(0)); } void loop(){ G</gaussian.h>…

http://stackoverflow.com/questions/11365857/real-time-auto-updating-incremental-plot-in-r n=1000 df=data.frame(time=1:n,y=runif(n)) window=100 for(i in 1:(n-window)) { flush.console() plot(df$time,df$y,type='l',xlim=c(i,i+window)) Sys.slee…

Mavericks、Homebrewでwxmaximaがビルドできないみたいだけど、インタプリタのmaximaはビルドできた brew install maxima途中 brew install --universal libpng brew install --universal gdが必要だった $ maxima Maxima 5.32.1 http://maxima.sourceforge.…