実験計画法と連立方程式 （１）実験計画法は以下の連立方程式をとくことに対応 y1=mu+a+b+ab+q[ab1] y2=mu-a+b+ab+q[ab2] y3=mu-a-b+ab+q[ab3]未知数 mu[1つ],a[1つ], b[1つ], ab[1つ？],q[N-3つ] を解ければよい。 （２）効率のよい連立方程式のたてかたが…

L8の実験関数 v1:第一水準 v2:第二水準 f:関数 L8<-function(x=1.0,fun=f,v1,v2){ mat<-cbind(v1,v2) e1<-c(1,1,1,1,1,1,1) e2<-c(1,1,1,2,2,2,2) e3<-c(1,2,2,1,1,2,2) e4<-c(1,2,2,2,2,1,1) e5<-c(2,1,2,1,2,1,2) e6<-c(2,1,2,2,1,2,1) e7<-c(2,2,1,1,2,2,…

２つのディレクトリからリストに含まれるファイルをコピー以下のプロセス list1.txtに含まれるファイルは 行数1_s.csv からコピー ディレクトリは../../list1dir/list2.txtに含まれるファイルは 行数1_s.csv からコピー ディレクトリは../../list2dir/list3.…

rubyでファイルの併合３ （共通するもの対応関係の出力） f=open("A.txt"); a=Hash::new(); k=1; f.readlines.each{|i| i.chomp! q=i.split(" ") a[q[0]]=q[0]+" "+k.to_s; k=k+1; } f3=open("B.txt"); b=Hash::new(); k=1; f3.readlines.each{|i| i.chomp! …

さまざまなエスケープ シェルスクリプト nkf -u xxxx.sh > tmp.sh sed "s/\`/\\\\\`/g;s/\"/\\\\\"/g;s/\\$/\\\\\$/g;" tmp.sh > tmp2.sh sed "s/\\\/\\\\\\\/g;s/\`/\\\\\`/g;s/\"/\\\\\"/g;s/\\$/\\\\\$/g;" tmp.sh > tmp3.sh sed "s/RRR/ \"/g;s/SSS/\" /…

scalac -optimise xxx.scala

メルセンヌツイスタ64bit版の導入ソースの入手 wget http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/mt19937-64.tgz gzip -d mt19937-64.tgz tar -xvf mt19937-64.tar コピー cd /mt19937-64 cp mt64.h /usr/include コンパイル gcc -c -O3 mt19937-64.c …

setの使い方 #include<stdio.h> #include<stdlib.h> #include<iostream> #include<set> #include<map> #include<vector> #include<algorithm> #include<math.h> int main(){ using namespace std; std::set<int> nodelist; int i; int nodeno=10; for(i=0;i</int></math.h></algorithm></vector></map></set></iostream></stdlib.h></stdio.h>

一元配置実験 > x<-1:4 > y<-1:4+rnorm(4,mean=0,sd=0.1) > q<-aov(y~x) > anova(q) Analysis of Variance Table Response: y Df Sum Sq Mean Sq F value Pr(>F) x 1 5.3171 5.3171 12046 8.3e-05 *** Residuals 2 0.0009 0.0004 --- Signif. codes: 0 ‘***’…

まずは普通に su do yum install R 失敗する Error: Package: R-core-2.15.2-1.el5.x86_64 (epel) Requires: libtcl8.4.so()(64bit) Error: Package: R-core-2.15.2-1.el5.x86_64 (epel) Requires: libtk8.4.so()(64bit) epel? とはEPELとは http://kb.xpres…