1つのwのみ。 #確率的最急勾配法 f<-file("rnorm_data.txt","r") data=function(){ a<-readLines(con=f,1) b=unlist(strsplit(a, "\\,")) # 文字 "," で分割 return(as.numeric(b)) } #data()#ファイルの1行目が読み込まれる #data()#ファイルの2行目が読み…

getwd() setwd("作業ディレクトリ") N=1000 out <- file("rnorm_data.txt", "w") # 書き込みモードで開く for(cnt in 1:N){ x=rnorm(5) for(i in 1:length(x)){ if (i < 5) writeLines(paste(x[i]),out,sep=", ") else if(i==5) writeLines(paste(x[i]),out,…

分布を平均0,分散1にしたときのネゲントロピーを計算し,横軸が0.001刻みでを0から1まで動かし,縦軸がネゲントロピーの近似した図 (ちゃんとやったのに図が間違っているかも？) ネゲントロピーが大きいという事はガウス分布から形状が離れているという事であ…

#一般化ガウス族(GGF) p=function(xi){ (1/(2*a^(1/a-1)*gamma(1/a)))*exp(-abs(xi)^a/a) } a=2 mean=(a^(1/a))*(gamma(2/a)/gamma(1/a)) var=(a^(2/a))*(gamma(3/a)/gamma(1/a))-((a^(1/a))*(gamma(2/a)/gamma(1/a)))^2 Gauss = function(x){ (1/sqrt(2*pi*…

my $row = $c->model('RemaindDB::Contents0')->create({ content => $content }); とかやったときに、データベースにInsertされない。不思議。 原因がよくわからないが、下が気になる。 $ script/remaind_server.pl -r DBIx::Class::Schema::connection(): …