読者です 読者をやめる 読者になる 読者になる

RでDeep Learning(12)

CNNの実装 library(dplyr) library(R6) #パラメータの初期化 #Conv層 W1 <- matrix(rnorm(750),nrow=25)*0.01 b1 <- matrix(numeric(30),nrow=1) #Affin層1 W2 <- matrix(rnorm(432000),nrow=4320)*0.01 b2 <- matrix(numeric(100),nrow=1) #Affine層2 W3 <-…

Rで行列計算(行列の足し算)

R

sweepは速いと思っていたが違った。 (環境:Windows 10 64bit、R3.3.3) > system.time(test1 <- out + matrix(rep(b1,576000),nrow=576000,byrow=T)) ユーザ システム 経過 0.13 0.03 0.16 > system.time(test2 <- sweep(out,2,b1,FUN="+",check.margin=F)…

RでDeep Learning(11)

Convolutin層、Pooling層の実装 Convolution <- R6Class("Convolution", #stride=1 #pad=0 #filterサイズ=5×5 #filter数=30 public = list( W = NULL, b = NULL, x = NULL, col = NULL, dW = NULL, db = NULL, initialize = function(W,b){ self$W <- W self…

Rで行列計算(行列の積)

R

RjpWikiによると「crossprod」と「%*%」では前者のほうが計算が速いとのこと。 実際にやってみた。 (環境:Windows 7 64bit、R3.4.0) > x <- matrix(rnorm(25000000),c(5000,5000)) > y <- matrix(rnorm(25000000),c(5000,5000)) > system.time(a <- t(x) …

RでDeep Learning(10)

Affineレイヤの書き換え Affine <- R6Class("Affine", public = list( W = NULL, b = NULL, x = NULL, x_shape = NULL, dW = NULL, db = NULL, initialize = function(W,b){ self$W <- W self$b <- b }, forward = function(x){ dummy_shape <- dim(x) dummy…

RでDeep Learning(9)

Rでim2colの実装(プーリング層はim2col_poolingとして別に実装) Rのdrop機能にはまった。 im2col <- function(input_data,filter_h,filter_w,stride=1,pad=0){ #imput_data:4次元配列からなるデータ #filter_h:フィルターの高さ #filter_w:フィルターの…

RでDeep Learning(8)

前回の続き。正解率はこの通り。 > test_acc [1] 0.9713 > train_acc [1] 0.9892667正解できなかった問題の一部。 y <- predict(x_test) y_1 <- apply(y,1,which.max) t <- apply(t_test,1,which.max) error <- t!=y_1 num_image <- x_test[error,] num_t <-…

RでDeep Learning(7)

前回の続き(テストデータで評価) library(dplyr) library(R6) #パラメータの初期化 W1 <- matrix(rnorm(784*50),nrow=784)*0.01 W2 <- matrix(rnorm(50*10),nrow=50)*0.01 b1 <- matrix(numeric(50),nrow=1) b2 <- matrix(numeric(10),nrow=1) #データの読…

RでDeep Learning(6)

環境はWindows10 & R-3.3.3ゼロから作るDeep Learning ―Pythonで学ぶディープラーニングの理論と実装作者: 斎藤康毅出版社/メーカー: オライリージャパン発売日: 2016/09/24メディア: 単行本(ソフトカバー)この商品を含むブログ (13件) を見る誤差逆伝播法…

RでDeep Learning(5)

環境はWindows10 & R-3.3.3ゼロから作るDeep Learning ―Pythonで学ぶディープラーニングの理論と実装作者: 斎藤康毅出版社/メーカー: オライリージャパン発売日: 2016/09/24メディア: 単行本(ソフトカバー)この商品を含むブログ (13件) を見る2層ニューラ…

RでDeep Learning(4)

環境はWindows10 & R-3.3.3ゼロから作るDeep Learning ―Pythonで学ぶディープラーニングの理論と実装作者: 斎藤康毅出版社/メーカー: オライリージャパン発売日: 2016/09/24メディア: 単行本(ソフトカバー)この商品を含むブログ (13件) を見る非常に良い本…

RでDeep Learning(3)

環境はWindows10 & R-3.3.3ゼロから作るDeep Learning ―Pythonで学ぶディープラーニングの理論と実装作者: 斎藤康毅出版社/メーカー: オライリージャパン発売日: 2016/09/24メディア: 単行本(ソフトカバー)この商品を含むブログ (13件) を見る非常に良い本…

RでDeep Learning(2)

環境はWindows10 & R-3.3.3ゼロから作るDeep Learning ―Pythonで学ぶディープラーニングの理論と実装作者: 斎藤康毅出版社/メーカー: オライリージャパン発売日: 2016/09/24メディア: 単行本(ソフトカバー)この商品を含むブログ (13件) を見る非常に良い本…

RでDeep Learning(1)

環境はWindows10 & R-3.3.3ゼロから作るDeep Learning ―Pythonで学ぶディープラーニングの理論と実装作者: 斎藤康毅出版社/メーカー: オライリージャパン発売日: 2016/09/24メディア: 単行本(ソフトカバー)この商品を含むブログ (13件) を見る非常に良い本…

RでEutilsを使ってみる

library(dplyr) library(rvest) pubmedID <- "25517282" address <- "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=pubmed&mode=XML&id=" myXML <- read_xml(paste(address,pubmedID,sep="")) #雑誌名 myXML %>% xml_nodes(xpath="//Journa…

Rのパイプ処理が使えない?

R

環境はWindows10 & R-3.3.3 原因は Help on topic 'filter' was found in the following packages: Return rows with matching conditions. (in package dplyr in library C:/R-3.3.3/library) Linear Filtering on a Time Series (in package stats in libr…

ExcelからWord内の単語を置換する

環境はWindows10 & Office 2013 Dim wdObj As Object Dim wdDoc As Object Set wdObj = CreateObject("Word.Application") wdObj.Visible = True wdObj.Activate Set wdDoc = wdObj.Documents.Open(ThisWorkbook.Path & "\test.docx") With wdDoc.Content.Fi…

RでTensorFlow

data(iris) mydata <- subset(iris,Species=="setosa") x_data <- mydata[["Sepal.Length"]] y_data <- mydata[["Sepal.Width"]] glm(y_data~x_data) Call: glm(formula = y_data ~ x_data) Coefficients: (Intercept) x_data -0.5694 0.7985 Degrees of Fre…

Windows7にTensorFlowを入れてみた

Python 3.5.3をインストールその後Windowsのコマンドプロンプトに下記のコマンドを入力 pip3 install --upgrade tensorflow これだけでPythonでTensorFlowが使えるようになった。うまくいかない人はこちら windowsにネイティブサポートしたTensorflowを動か…

Pythonでフォルダ内のファイル一覧をパス無しで取得する

import glob import os files = glob.glob('D:/CEN_list/CEN_list/bin/Debug/file/*.txt') for f in files: (dirname,filename)=os.path.split(f) print(filename)

Journalのタイトルをリストアップする

c#

using System.Text; using System.IO; namespace html { class Program { static void Main(string[] args) { StreamWriter writer = new StreamWriter("article.html", true, Encoding.UTF8); writer.WriteLine("<ol>"); for (int volume = 1; volume < 6; volu</ol>…

Journalのタイトルを一括ダウンロード

c#

namespace download { class Program { static void Main(string[] args) { for (int volume = 1; volume < 6; volume++) { for (int issue = 1; issue < 3; issue++) { for (int page = 1; page < 3; page++) { string url = "http://link.springer.com/jou…

英語学習のためのアプリを作る

c#

MainWindow.xaml <Window x:Class="WpfApplication1.MainWindow" xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation" xmlns:x="http://schemas.microsoft.com/winfx/2006/xaml" Title="英語問題" Height="350" Width="816" Initialized="Window_Initialized"> <Grid> </grid></window>

Rの結果をPDFに出力する(文字化け対応)

R

Rの結果をテキストファイルとpngファイルに出力する。 それらをWordで受け取ってPDFファイルに出力する。 Rのコード data(iris) mydata <- iris 比較 <- c("Sepal.Length","Sepal.Width","Petal.Length","Petal.Width") counter=1 for (hikaku in 1:length(…

Excelのコマンドボタンの右クリックを有効にする方法

以下の二つは共存可能。 Private Sub CommandButton1_Click() MsgBox ("左クリック") End Sub Private Sub CommandButton1_MouseUp(ByVal Button As Integer, ByVal Shift As Integer, ByVal X As Single, ByVal Y As Single) If Button = 2 Then MsgBox ("…

VBAでカレンダーコントロールを自作する

簡単にコピペで再現できるようにボタンなどのコントロールの配置やそのプロパティはすべてコードに記述。フォームやモジュールを追加して下記をコピペするだけで完成。①カレンダーフォーム「Calender」を追加する Option Explicit Private Sub UserForm_Init…

Excelユーザーフォームに動的にコマンドボタンを追加する

クラスモジュール「Button_Setting」をあらかじめ追加しておく Option Explicit Public WithEvents original_button As CommandButton Private Sub original_button_Click() 'クリックイベントをここに記載 End Sub Option Explicit Sub ユーザーフォームを…

数独を解く(別解:バックトラック法)

c#

Form1.cs using System; using System.Linq; using System.Windows.Forms; namespace 数独バックトラック法 { public partial class Form1 : Form { Button[] IDbutton = new Button[81]; Button kensyo = new Button(); Button zikko = new Button(); Butto…

数独を解く(6日目)

c#

実行ボタンのClickイベントハンドラを記述する。 private void zikko_Click(object sender, EventArgs e) { Stopwatch sw = new Stopwatch(); sw.Start(); //////////////////////////////////////////////////////////////////////////////////////////////…

数独を解く(5日目)

c#

新規ボタンのClickイベントハンドラを記述する。 private void newGame_Click(object sender, EventArgs e) { foreach (Button a in IDbutton) { a.Text = ""; a.Enabled = true; } kensyo.Enabled = true; zikko.Enabled = false; }

数独を解く(4日目)

c#

検証ボタンのClickイベントハンドラを記述。 private void kensyo_Click(object sender, EventArgs e) { //数字の拾い出し int d; for (int i = 0; i < 81; i++) { int.TryParse(IDbutton[i].Text, out d); kazu[i] = d; } //検証を実行 if (kensyo_do()) { …

数独を解く(3日目)

c#

以下のプログラムをコンソールアプリケーションとして実行する。 「resource.txt」というファイルが作られる。 出来たファイルをおおもとのプログラムにリソースとして取り込む。 using System.Collections.Generic; namespace 下準備 { class Program { sta…

数独を解く(2日目)

c#

ユーザーフォームを追加する (input.cs) Loadのイベントハンドラを作成する必要あり。 using System; using System.Drawing; using System.Windows.Forms; namespace 数独をとく { public partial class input : Form { Button[] button = new Button[9]; pu…

数独を解く(1日目)

c#

まずはGUIまわりから(GUI.cs) メインフォーム(Form1)のLoadのイベントハンドラーを作成する必要あり。 using System; using System.Drawing; using System.Windows.Forms; namespace 数独ソルバー { public partial class Form1 : Form { private void Form1…

例外処理

c#

try { } catch(Exception ex) { MessageBox.Show(ex.Message, Application.ProductName , MessageBoxButtons.OK, MessageBoxIcon.Error); }

フォームを閉じるときの終了確認

c#

private void MainForm_FormClosing(object sender, FormClosingEventArgs e) { DialogResult reply; reply = MessageBox.Show("このアプリケーションを終了します" , Application.ProductName, MessageBoxButtons.OKCancel , MessageBoxIcon.Information); …

Rの計算結果を新たなwindowに出力する方法

コンソールに文字として結果が出力される統計解析結果などを他のwindowに出力する方法。 以下の関数を定義しておく。 original_output <- function(){ x <- winDialogString("output to a new window","") filename <- tempfile() capture.output(eval(parse…

ExcelファイルをRに読み込む

R

クリップボード経由でファイルを読み込む方法 read.table(file="clipboard",header=T,sep="\t")

プロットのマーカーをきれいに表示する方法

R

plot(0, 0, xlim=c(0,3), ylim = c(0, 8), type = "n", axes = F, xlab = "", ylab = "") a<-rep(0:2,c(9,9,8)) b<-c(8:0,8:0,8:1) points(a,b,pch=0:25,cex=3) text(a+0.45,b,0:25)

biocGraph パッケージを使用して Rgraphviz で描いた図にリンクをはる方法

R

library(RCurl) library(XML) library(geneplotter) library(Rgraphviz) library(biocGraph) library(org.Hs.eg.db) gene_article <- function(geneID){ ID2symbol <- as.list(org.Hs.egSYMBOL) gene_name <- unlist(ID2symbol[geneID]) url <- paste("http:…

きれいな棒グラフを描く

ピンク (255,153,255) 黄色 (255,255,102) 水色 (204,236,255) グリーン (102,255,102)

きれいな円グラフを描く

薄い青 (134,202,232) グレー (153,165,204) ピンク (246,170,200) 薄い黄 (254,210,122) 薄いオレンジ (245,162,108) ライム (163,202,112)

C#でXML操作

using System; using System.Linq; using System.Windows.Forms; using System.Xml.Linq; 結合 XDocument xdoc1 = XDocument.Load("http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=pubmed&mode=XML&id=14993493"); XDocument xdoc2 = XDocume…

画像ファイルを複製する

using System; using System.Drawing; using System.Windows.Forms; namespace bitmap { public partial class Form1 : Form { public Form1() { InitializeComponent(); } int file_num; string file_name; string input_name; private void button1_Click(…

ImageMagickを使わずにWindowsでpng画像からgif画像を作成する

R

必要なパッケージは「png」と「caTools」 図形描画関数群を「draw.R」の名前であらかじめ保存 source("draw.R") library(png) library(caTools) for(i in 1:30){ file.name <- sprintf("%02d.png",i) png(file.name, width=960,height=720,bg="black") par(m…

PDFを操作するプログラム

using System; using System.Windows.Forms; namespace PDF { public partial class Form1 : Form { public Form1() { InitializeComponent(); } MainMenu main = new MainMenu(); MenuItem top = new MenuItem(); readonly string[] _kugiri = { ".pdf", ".…

6日間で文献管理アプリを作る(おまけ)

Form1.cs using System; using System.Windows.Forms; using System.IO; namespace 文献管理 { public partial class Form1 : Form { string choose_file; //選択したファイルの名前 int article_number; //文献数 string[] article_array; //PMIDの配列 con…

6日間で文献管理アプリを作る(6日目)

article_formに関して まずは以下の変数を定義する public FileInfo file_info; //選択しているXMLファイルの情報 public int article_number2; //文献数(Form1から受け取る) public string[] article_array2; //PMIDの配列(Form1から受け取る) int curr…

6日間で文献管理アプリを作る(5日目)

メニューバー「ファイル」/「HTML出力」を定義 まずは以下の「output_html.cs」を追加 using System; using System.Windows.Forms; using System.Xml.Xsl; namespace 文献管理 { public partial class Form1 { void output_html() { //OpenFileDialogクラス…

6日間で文献管理アプリを作る(4日目)

メニューバー「データ」/「追加」「削除」を定義 まずは以下の「add_remove.cs」を追加 using System.Windows.Forms; using System.Xml.Linq; using System.Linq; namespace 文献管理 { public partial class Form1 { void article_add(string a) { try { b…