PPT図をTEXに張り込むのにはまったので記録。まず、少し前までつかえた、PSプリンタをパワポにとりこむ、 EPSにして出力、バウンダリーボックスを調整して入れ込むという 方法はもう使えないみたいです。これらの方法はもうだめ http://www4.cs.miyazaki-u.a…

scalaでファイルの内容一括読み込みと表示 var text1=Source.fromFile("text1.txt").getLines.toArray Source.fromFile("text2.txt").getLines.foreach{println} Source.fromFile("text3.txt","utf-8").getLines.toArray.map{b=>| b.split("\t")(4).toDouble…

処理結果が-1の場合は0を返す。それ以外は値を２倍して５をたす． http://xawa99.blogspot.jp/2012/09/Scala-MatchBasic.html var a=-1 var b=(2*a+5).match{case -1=>0; case other => other}

scalaでrubyのjoin的なものを書く場合はmkString Array(1,2,3,4).mkString(",")

Scalaで単位根検定 unit root test (Phillips-Perron Tests in scala) RのPP.testをscalaに翻訳 import scala.math._ import org.apache.commons.math3.stat.correlation.SpearmansCorrelation; import org.apache.commons.math3.stat.regression.OLSMultipl…

英作文の仕方（１）主体を考える （２）主体の修飾を考える （３）主体の動作＋時間を考える （４）客体を考える （５）客体の修飾を考える （６）場所や状態を考える。 （７）理由等を考えるWhy は SVをした理由の省略形 前置詞のin は接続詞の代わり When …

重回帰分析 > x [1] 1 2 3 4 5 6 > q1<-cbind(rep(1,6),x) > q1 x [1,] 1 1 [2,] 1 2 [3,] 1 3 [4,] 1 4 [5,] 1 5 [6,] 1 6 > solve(t(q1)%*%q1) x 0.8666667 -0.20000000 x -0.2000000 0.05714286 >|sh| [,1] [1,] 0.0616738 > > q2 > 0.258^2*solve(t(q2)%…

単回帰の係数の標準偏差は， y=bx+f 平均値は， bcov(x,y)/var(x) で推定できる． 誤差：y-bxは， sd な感じ（本当は不偏推定量を使うべき）． sd_b で与えられる sdは，誤差の推定値.チェック summary(lm(y~x)) http://www.ic.daito.ac.jp/~tkadoda/2005/ec…

http://bitingcode.blogspot.jp/2012/01/simplest-olsmultiplelinearregression.html 参考

参考： http://www.jaceksuda.com/teach/pse2011/class9_slides.pdf単位根検定→定常性の検定→RW性の検定70年代後半まで： 線形トレンドでlog(GDP)をフィッティング、確率部分はここからのゆらぎ 70年代問題 log(GDP)が線形増加がスローダウン→予見できなくな…

組み合わせの差 ax-yb 例 x=(1,2,3); y=(4,5,6); a=1 ;b=2 a<-1;b<-2;log(outer(exp(1:3)^a,exp(4:6)^-b)) 出力 [,1] [,2] [,3] [1,] -7 -9 -11 [2,] -6 -8 -10 [3,] -5 -7 -9 応用例：２つのポアソン分布の差の確率 Z=0.4X-0.3Y mul_pdf2<-function(x,c,sds…

scalaで移動ボックス平均 object test{ def box_mean(x:List[Double],n:Int=2):List[Double]={ var y=List[Double](); var k:Int=0; var nn=n-1; var sum=0.0; x.foreach{ a => k match{ case k if k%n==n-1 => sum=sum+a/n;y=y::sum; sum=0; case _ => sum=…

入力 $ cat test1.txt abc xxx 11 efg yyy 23 hij zzz 33 lmk qqq 41.3 $ cat test2.txt abc xxx 11 efg yyy 39 hij zzz 22 lmk qqq 17 プログラム import scala.io.Source import org.apache.commons.math3.stat.correlation.SpearmansCorrelation; object …

例１） 基本よみとって出力 import scala.io.Source object test{ def main(args: Array[String]){ var source =Source.fromFile("test.txt","utf-8") val lines=source.getLines var a1: String =""; var k=0; lines.foreach{println(_);} } } foreachの- …

http://sourceforge.jp/projects/milm-search/lists/archive/public/2012-April/000174.html がわかりやすい。 mkdir -p ~/.vim/{ftdetect,indent,syntax} で、できた各フォルダに それぞれのフォルダの http://lampsvn.epfl.ch/svn-repos/scala/scala-tool-…

関数などはjavaで順位相関を参照 http://d.hatena.ne.jp/arupaka-_-arupaka/20130218/1361176768 import org.apache.commons.math3.stat.correlation.SpearmansCorrelation; object test6{ def main(args: Array[String]){ val data1=Array(20,17,30,42.3,17…

Javaで順位相関 (Spearman order correlation on Java)： Apache Commons のCommons Mathの利用（↓コンパイル、導入法などはこれを参考） http://d.hatena.ne.jp/arupaka-_-arupaka/20130218/1361169072 本家を参考に http://commons.apache.org/math/apidoc…

Javaで順位統計： ランキングの計算 Apache Commons のCommons Mathの利用 http://d.hatena.ne.jp/arupaka-_-arupaka/20130218/1361169072 本家を参考に http://commons.apache.org/math/userguide/stat.html#a1.6_Rank_transformations import org.apache.c…

Commons Mathはjavaの準正規ライブラリ的なのもの数学ライブラリ 統計計算などいろいろできて便利。ライブラリできること http://ja.wikipedia.org/wiki/Apache_Commons_Math 公式 http://commons.apache.org/導入は http://wiki.unfindable.net/webbook2/in…

(1)まず、ダウンロード http://www.scala-lang.org/downloads (2)次に、解凍 tar zxvf scala-2.9.0.1.tgz .basrcに以下を追加。 （ホームに展開した場合） export SCALA_HOME=~/scala-2.10.0 export PATH=$PATH:$SCALA_HOME/bin https://sites.google.com/si…

C++からRの関数を呼び出す. cygwin，gcc,g++版： 時間計測： これ http://d.hatena.ne.jp/teramonagi/20101213/1292253192 を参考にtest.cpp #include <windows.h> #include <iostream> //doubleの引数一個取ってdoubleの値返す関数ポインタをRFuncと定義 typedef double(*RFunc)(</iostream></windows.h>…

kugiriapply(matrix(x[1:(kugiri*((length(x)%/%kugiri)))],kugiri),2,"mean")