2013年6月 – IT研究所

統計ソフトRを使用して株価チャートを簡単に作成(3)

投稿者: utsubo 投稿日: 2013-06-30 in R

データをとれるようになったのでチャートを書いてみる

library(quantmod)
library(RFinanceYJ)
stockCode&lt;-&quot;6758.t&quot;
today&lt;-Sys.Date()
dataTerm&lt;-100 # 100日分取得
sinceDate&lt;-today-dataTerm
colnames(stockData)&lt;-c(&quot;Date&quot;,&quot;Open&quot;,&quot;High&quot;,&quot;Low&quot;,&quot;Close&quot;,&quot;Volume&quot;)
stockData.zoo&lt;-read.zoo(stockData)
candleChart(stockData.zoo,theme=&quot;white&quot;)

library(quantmod)

library(RFinanceYJ)

stockCode<-"6758.t"

today<-Sys.Date()

dataTerm<-100 # 100日分取得

sinceDate<-today-dataTerm

colnames(stockData)<-c("Date","Open","High","Low","Close","Volume")

stockData.zoo<-read.zoo(stockData)

candleChart(stockData.zoo,theme="white")

ファイルにする場合はこちら

library(quantmod)
library(RFinanceYJ)
stockCode&lt;-&quot;6758.t&quot;
today&lt;-Sys.Date()
dataTerm&lt;-100 # 100日分取得
sinceDate&lt;-today-dataTerm
colnames(stockData)&lt;-c(&quot;Date&quot;,&quot;Open&quot;,&quot;High&quot;,&quot;Low&quot;,&quot;Close&quot;,&quot;Volume&quot;)
stockData.zoo&lt;-read.zoo(stockData)
pngDir&lt;-&quot;png/&quot;
fileName&lt;-paste(pngDir,stockCode,&quot;.png&quot;,sep=&quot;&quot;)
png(fileName,width=600,height=400)
candleChart(stockData.zoo,theme=&quot;white&quot;)
dev.off()

library(quantmod)

library(RFinanceYJ)

stockCode<-"6758.t"

today<-Sys.Date()

dataTerm<-100 # 100日分取得

sinceDate<-today-dataTerm

colnames(stockData)<-c("Date","Open","High","Low","Close","Volume")

stockData.zoo<-read.zoo(stockData)

pngDir<-"png/"

fileName<-paste(pngDir,stockCode,".png",sep="")

png(fileName,width=600,height=400)

candleChart(stockData.zoo,theme="white")

dev.off()

これをシェルでまわせば複数銘柄とれます

chart.R

stockCode&lt;-commandArgs()[5]
library(quantmod)
library(RFinanceYJ)
today&lt;-Sys.Date()
dataTerm&lt;-100 # 100日分取得
sinceDate&lt;-today-dataTerm
colnames(stockData)&lt;-c(&quot;Date&quot;,&quot;Open&quot;,&quot;High&quot;,&quot;Low&quot;,&quot;Close&quot;,&quot;Volume&quot;)
stockData.zoo&lt;-read.zoo(stockData)
pngDir&lt;-&quot;png/&quot;
fileName&lt;-paste(pngDir,stockCode,&quot;.png&quot;,sep=&quot;&quot;)
png(fileName,width=600,height=400)
candleChart(stockData.zoo,theme=&quot;white&quot;)
dev.off()

stockCode<-commandArgs()[5]

library(quantmod)

library(RFinanceYJ)

today<-Sys.Date()

dataTerm<-100 # 100日分取得

sinceDate<-today-dataTerm

colnames(stockData)<-c("Date","Open","High","Low","Close","Volume")

stockData.zoo<-read.zoo(stockData)

pngDir<-"png/"

fileName<-paste(pngDir,stockCode,".png",sep="")

png(fileName,width=600,height=400)

candleChart(stockData.zoo,theme="white")

dev.off()

chart.sh

#!/bin/sh
codes="
6758
9501
9433
"
for code in $codes;do
R --vanilla --slave --args ${code}.t &lt;	chart.R
done

#!/bin/sh

codes="

6758

9501

9433

for code in $codes;do

R --vanilla --slave --args ${code}.t < chart.R

done

PlayFramework2.1で既存のMySQLにつないでみる

投稿者: utsubo 投稿日: 2013-06-25 in PlayFramework

ちょっと悩んだのでメモ

PlayFrameworkで既存のMySQLにつなぐためには

テーブルはこれ

use sampled
;
create table sample(
id integer,
name varchar(32)
)
;
create unique index idx_sampleTable on sampleTable(id)
;

use sampled

;

create table sample(

id integer,

name varchar(32)

)

;

create unique index idx_sampleTable on sampleTable(id)

;

まずモデル

package models;

import javax.persistence.Column;
import javax.persistence.Entity;
import javax.persistence.Id;
import javax.persistence.Table;

import play.db.ebean.model;

@Entity
@Table(name="sample")
public class Sample extends Model
{
	@Column(name="id")
	@Id
	public integer id;

	@Column(name="name")
	public String name;

	public String toString(){
		return id+","+name;
	}
}

package models;

import javax.persistence.Column;

import javax.persistence.Entity;

import javax.persistence.Id;

import javax.persistence.Table;

import play.db.ebean.model;

@Entity

@Table(name="sample")

public class Sample extends Model

{

@Column(name="id")

@Id

public integer id;

@Column(name="name")

public String name;

public String toString(){

return id+","+name;

}

application.conf

..
db.default.driver=com.mysql.jdbc.Driver
db.default.url="jdbc:mysql://dbserver/sampledb"
db.default.user=root
db.default.password=password

evolutionplugin=disabled


ebean.default="models.*"
..

db.default.driver=com.mysql.jdbc.Driver

db.default.url="jdbc:mysql://dbserver/sampledb"

db.default.user=root

db.default.password=password

evolutionplugin=disabled

ebean.default="models.*"

これを使って、ViewとControllerにつなげれやればおっけー

統計ソフトRを使用して株価チャートを簡単に作成(2)

投稿者: utsubo 投稿日: 2013-06-24 in R 1件のコメント

RFinanceYJパッケージですがどうやらYahoo側の変更のせいでうまく動かないみたいです

> stockCode<-"6758.t"
> quoteStockTsData(stockCode)
 以下にエラー as.POSIXlt.character(x, tz, ...) : 
	character string is not in a standard unambiguous format

> stockCode<-"6758.t"

> quoteStockTsData(stockCode)

以下にエラー as.POSIXlt.character(x, tz, ...) :

character string is not in a standard unambiguous format

いろいろ調べてみたところ、どうやら日付部分が日本語のためうまく処理できていないようでした。

ちょっといじったソースはこれ

quoteStockTsData <- function(x, since=NULL,start.num=0,date.end=NULL,time.interval='daily')
{
	function.stock <- function(quote.table.item){
		if( xmlSize(quote.table.item) < 5) return(NULL)
		d <- convertToDate(xmlValue(quote.table.item[[1]]),time.interval)
		o <- as.number(xmlValue(quote.table.item[[2]]))
		h <- as.number(xmlValue(quote.table.item[[3]]))
		l <- as.number(xmlValue(quote.table.item[[4]]))
		c <- as.number(xmlValue(quote.table.item[[5]]))
		v <- ifelse(xmlSize(quote.table.item) >= 6,as.number(xmlValue(quote.table.item[[6]])),0)
		return(data.frame(date=d,open=o,height=h,low=l,close=c,volume=v))
	}
	return(quoteTsData(x,function.stock,since,start.num,date.end,time.interval))
}
quoteFundTsData <- function(x, since=NULL,start.num=0,date.end=NULL,time.interval='daily')
{
	function.fund <- function(quote.table.item){
		d <- convertToDate(xmlValue(quote.table.item[[1]]),time.interval)
		if(time.interval=='monthly'){
			d <- endOfMonth(d)
		}
		c <- as.number(xmlValue(quote.table.item[[2]]))
		v <- as.number(xmlValue(quote.table.item[[3]]))
		return(data.frame(date=d,constant.value=c,NAV=v))
	}
	return(quoteTsData(x,function.fund,since,start.num,date.end,time.interval))
}
quoteFXTsData <- function(x, since=NULL,start.num=0,date.end=NULL,time.interval='daily')
{
	function.fx <- function(quote.table.item){
		d <- convertToDate(xmlValue(quote.table.item[[1]]),time.interval)
		o <- as.number(xmlValue(quote.table.item[[2]]))
		h <- as.number(xmlValue(quote.table.item[[3]]))
		l <- as.number(xmlValue(quote.table.item[[4]]))
		c <- as.number(xmlValue(quote.table.item[[5]]))
		return(data.frame(date=d,open=o,height=h,low=l,close=c))
	}
	return(quoteTsData(x,function.fx,since,start.num,date.end,time.interval))
}


quoteTsData <- function(x,function.financialproduct,since,start.num,date.end,time.interval){
	r <- NULL
	result.num <- 51
	quote.table.list <- list(NULL)
	quote.table <- NULL
	financial.data <- data.frame(NULL)
	start <- (gsub("([0-9]{4,4})-([0-9]{2,2})-([0-9]{2,2})","&c=\\1&a=\\2&b=\\3",since))
	end	 <- (gsub("([0-9]{4,4})-([0-9]{2,2})-([0-9]{2,2})","&f=\\1&d=\\2&e=\\3",date.end))

	if(!any(time.interval==c('daily','weekly','monthly'))) stop("Invalid time.interval value")
	
	while( result.num >= 51 ){
		quote.url <- paste('http://table.yahoo.co.jp/t?s=',x,start,end,'&y=',start.num,'&g=',substr(time.interval,1,1),sep="")
		try( r <- xmlRoot(htmlTreeParse(quote.url,encoding="EUC-JP",error=xmlErrorCumulator(immediate=F))), TRUE)
		if( is.null(r) ) stop(paste("Can not access :", quote.url))

		try( quote.table <- r[[2]][[1]][[1]][[13]][[1]][[1]][[1]][[4]][[1]][[1]][[1]], TRUE )
		if( is.null(quote.table) ) stop(paste("Can not quote :", x))

		size <- xmlSize(quote.table) 

		for(i in 2:size){
			financial.data <- rbind(financial.data,function.financialproduct(quote.table[[i]]))
		}	 
					
		result.num <- xmlSize(quote.table)
		start.num <- start.num + 50
		Sys.sleep(1)
	}
	financial.data <- financial.data[order(financial.data$date),]
	return(financial.data)
}

#convert string formart date to POSIXct object
convertToDate <- function(date.string,time.interval)
{
	#date.string <- iconv(date.string,"EUC-JP","UTF-8","")
	date.string<-gsub("日","",gsub("[年月]","-",date.string))
	#data format is different between monthly and dialy or weekly
	if(any(time.interval==c('daily','weekly'))){
		result <- gsub("^([0-9]{4})([^0-9]+)([0-9]{1,2})([^0-9]+)([0-9]{1,2})([^0-9]+)","\\1-\\3-\\5",date.string)
	}else if(time.interval=='monthly'){
		result <- gsub("^([0-9]{4})([^0-9]+)([0-9]{1,2})([^0-9]+)","\\1-\\3-01",date.string)
	}
	return(as.POSIXct(result))
}

#convert string to number.
as.number <- function(string)
{
	return(as.double(as.character(gsub("[^0-9.]", "",string))))
}
#return end of month date.
endOfMonth <- function(date.obj)
{
	startOfMonth		 <- as.Date(format(date.obj,"%Y%m01"),"%Y%m%d")
	startOfNextMonth <- as.Date(format(startOfMonth+31,"%Y%m01"),"%Y%m%d")
	return(startOfNextMonth-1)
}

100

quoteStockTsData <- function(x, since=NULL,start.num=0,date.end=NULL,time.interval='daily')

{

function.stock <- function(quote.table.item){

if( xmlSize(quote.table.item) < 5) return(NULL)

d <- convertToDate(xmlValue(quote.table.item[[1]]),time.interval)

o <- as.number(xmlValue(quote.table.item[[2]]))

h <- as.number(xmlValue(quote.table.item[[3]]))

l <- as.number(xmlValue(quote.table.item[[4]]))

c <- as.number(xmlValue(quote.table.item[[5]]))

v <- ifelse(xmlSize(quote.table.item) >= 6,as.number(xmlValue(quote.table.item[[6]])),0)

return(data.frame(date=d,open=o,height=h,low=l,close=c,volume=v))

}

return(quoteTsData(x,function.stock,since,start.num,date.end,time.interval))

}

quoteFundTsData <- function(x, since=NULL,start.num=0,date.end=NULL,time.interval='daily')

{

function.fund <- function(quote.table.item){

d <- convertToDate(xmlValue(quote.table.item[[1]]),time.interval)

if(time.interval=='monthly'){

d <- endOfMonth(d)

}

c <- as.number(xmlValue(quote.table.item[[2]]))

v <- as.number(xmlValue(quote.table.item[[3]]))

return(data.frame(date=d,constant.value=c,NAV=v))

}

return(quoteTsData(x,function.fund,since,start.num,date.end,time.interval))

}

quoteFXTsData <- function(x, since=NULL,start.num=0,date.end=NULL,time.interval='daily')

{

function.fx <- function(quote.table.item){

d <- convertToDate(xmlValue(quote.table.item[[1]]),time.interval)

o <- as.number(xmlValue(quote.table.item[[2]]))

h <- as.number(xmlValue(quote.table.item[[3]]))

l <- as.number(xmlValue(quote.table.item[[4]]))

c <- as.number(xmlValue(quote.table.item[[5]]))

return(data.frame(date=d,open=o,height=h,low=l,close=c))

}

return(quoteTsData(x,function.fx,since,start.num,date.end,time.interval))

}

quoteTsData <- function(x,function.financialproduct,since,start.num,date.end,time.interval){

r <- NULL

result.num <- 51

quote.table.list <- list(NULL)

quote.table <- NULL

financial.data <- data.frame(NULL)

start <- (gsub("([0-9]{4,4})-([0-9]{2,2})-([0-9]{2,2})","&c=\\1&a=\\2&b=\\3",since))

end <- (gsub("([0-9]{4,4})-([0-9]{2,2})-([0-9]{2,2})","&f=\\1&d=\\2&e=\\3",date.end))

if(!any(time.interval==c('daily','weekly','monthly'))) stop("Invalid time.interval value")

while( result.num >= 51 ){

quote.url <- paste('http://table.yahoo.co.jp/t?s=',x,start,end,'&y=',start.num,'&g=',substr(time.interval,1,1),sep="")

try( r <- xmlRoot(htmlTreeParse(quote.url,encoding="EUC-JP",error=xmlErrorCumulator(immediate=F))), TRUE)

if( is.null(r) ) stop(paste("Can not access :", quote.url))

try( quote.table <- r[[2]][[1]][[1]][[13]][[1]][[1]][[1]][[4]][[1]][[1]][[1]], TRUE )

if( is.null(quote.table) ) stop(paste("Can not quote :", x))

size <- xmlSize(quote.table)

for(i in 2:size){

financial.data <- rbind(financial.data,function.financialproduct(quote.table[[i]]))

}

result.num <- xmlSize(quote.table)

start.num <- start.num + 50

Sys.sleep(1)

}

financial.data <- financial.data[order(financial.data$date),]

return(financial.data)

}

#convert string formart date to POSIXct object

convertToDate <- function(date.string,time.interval)

{

#date.string <- iconv(date.string,"EUC-JP","UTF-8","")

date.string<-gsub("日","",gsub("[年月]","-",date.string))

#data format is different between monthly and dialy or weekly

if(any(time.interval==c('daily','weekly'))){

result <- gsub("^([0-9]{4})([^0-9]+)([0-9]{1,2})([^0-9]+)([0-9]{1,2})([^0-9]+)","\\1-\\3-\\5",date.string)

}else if(time.interval=='monthly'){

result <- gsub("^([0-9]{4})([^0-9]+)([0-9]{1,2})([^0-9]+)","\\1-\\3-01",date.string)

}

return(as.POSIXct(result))

}

#convert string to number.

as.number <- function(string)

{

return(as.double(as.character(gsub("[^0-9.]", "",string))))

}

#return end of month date.

endOfMonth <- function(date.obj)

{

startOfMonth <- as.Date(format(date.obj,"%Y%m01"),"%Y%m%d")

startOfNextMonth <- as.Date(format(startOfMonth+31,"%Y%m01"),"%Y%m%d")

return(startOfNextMonth-1)

}

これをRFinanceYJパッケージに続けて呼べばデータが取得できました

> 
> quoteStockTsData("6758.t")
				 date open height	low close	 volume
62 2013-03-25 1694	 1729 1662	1712 34382500
61 2013-03-26 1682	 1688 1660	1666 26291500
60 2013-03-27 1666	 1691 1663	1675 18503100
59 2013-03-28 1663	 1666 1587	1625 33812600
...

> quoteStockTsData("6758.t")

date open height low close volume

62 2013-03-25 1694 1729 1662 1712 34382500

61 2013-03-26 1682 1688 1660 1666 26291500

60 2013-03-27 1666 1691 1663 1675 18503100

59 2013-03-28 1663 1666 1587 1625 33812600

...

統計ソフトRを使用して株価チャートを簡単に作成(1)

投稿者: utsubo 投稿日: 2013-06-18 in R 1件のコメント

統計ソフトのRには様々なパッケージがありExcelなどでは結構手間がかかることでもあっという間にできる場合があります

今回は株価チャートを作成してみました

まずはRのインストールから

環境

mac OS10.8.3
R 3.0.1

http://cran.r-project.org/

こちらからMacOX用のものをダウンロード。（以前は３２ビット、６４ビットと分かれていたが統合されたみたい）
インストールして起動
チャート用のパッケージとデータダウンロード用のパッケージを入れる

> #プロキシサーバ越えの場合はこれを入れる
> Sys.setenv(http_proxy="http://proxyserver:port")
> 
> install.packages("quantmod")
> install.packages("RFinanceYJ")
> library(RFinanceYJ)
 要求されたパッケージ XML をロード中です 
> library(quantmod)
 要求されたパッケージ Defaults をロード中です 
 要求されたパッケージ xts をロード中です 
 要求されたパッケージ zoo をロード中です 

 次のパッケージを付け加えます: ‘zoo’ 

 以下のオブジェクトはマスクされています (from ‘package:base’) : 

		 as.Date, as.Date.numeric 

 要求されたパッケージ TTR をロード中です 
Version 0.4-0 included new data defaults. See ?getSymbols.
> stockData<-quoteStockTsData("6758.t")
 以下にエラー order(financial.data$date) :	引数 1 がベクトルではありません

> #プロキシサーバ越えの場合はこれを入れる

> Sys.setenv(http_proxy="http://proxyserver:port")

> install.packages("quantmod")

> install.packages("RFinanceYJ")

> library(RFinanceYJ)

要求されたパッケージ XML をロード中です

> library(quantmod)

要求されたパッケージ Defaults をロード中です

要求されたパッケージ xts をロード中です

要求されたパッケージ zoo をロード中です

次のパッケージを付け加えます: ‘zoo’

以下のオブジェクトはマスクされています (from ‘package:base’) :

as.Date, as.Date.numeric

要求されたパッケージ TTR をロード中です

Version 0.4-0 included new data defaults. See ?getSymbols.

> stockData<-quoteStockTsData("6758.t")

以下にエラー order(financial.data$date) : 引数 1 がベクトルではありません

なぜかエラー。。。

xcodeでフレームワーク作成

投稿者: utsubo 投稿日: 2013-06-17 in iOS

フレームワークの作成ですが、ちょっとググるといろいろ書いてあるので、動くものを作成してみました。

https://github.com/anagotan/ios_framework_sample

カーソルがウインドウの外に出たのを検知する

投稿者: utsubo 投稿日: 2013-06-06 in win32

win32apiではWM_MOUSEMOVEでマウスが移動しているというのは判定できるのですが、マウスがそのWindowの外に出たのがそのままではわかりません。

そこでTrackMouseEventを使用するらしい

この辺りがすごく参考になる

http://hp.vector.co.jp/authors/VA016117/hoverbtn2.html

http://support.microsoft.com/kb/183107/ja

GDIでPNGを読み込む

投稿者: utsubo 投稿日: 2013-06-05 in win32 2件のコメント

GDIでは通常PNGファイルを読み込むことはできないらしい。

どうしても読み込む必要があったのでちょっと調べてみた

http://code.logos.com/blog/2008/09/displaying_a_splash_screen_with_c_part_i.html

このやり方でうまく読み込めるのだが、どうやらマシンによっては読み込みにFailする。一応XPSP３から使えるとなっているのだが。。。

仕方がないのでGDI+を用いました。

HBITMAP LoadImage( HINSTANCE hinst, LPCTSTR pszName, LPCTSTR pszType )
{
		HRSRC hRes = FindResource(hinst, pszName, pszType);
		DWORD Size = SizeofResource(hinst, hRes);
		HGLOBAL hData = LoadResource(hinst, hRes);
		const void *pData = LockResource(hData);
		HGLOBAL hBuffer = GlobalAlloc(GMEM_MOVEABLE, Size);
		void *pBuffer = GlobalLock(hBuffer);
		CopyMemory(pBuffer, pData, Size);
		GlobalUnlock(hBuffer);

		IStream *pStream;
		CreateStreamOnHGlobal(hBuffer, TRUE, &pStream) ;
	 IStream *pStream=CreateStreamOnResource(hinst,pszName,pszType);

		Gdiplus::Bitmap* pBitmap = Gdiplus::Bitmap::FromStream(pStream);
		pStream->Release();
		HBITMAP ret;
		pBitmap->GetHBITMAP(Color(0,0,0,0),&ret);
		delete pBitmap;
	return ret;
}

HBITMAP LoadImage( HINSTANCE hinst, LPCTSTR pszName, LPCTSTR pszType )

{

HRSRC hRes = FindResource(hinst, pszName, pszType);

DWORD Size = SizeofResource(hinst, hRes);

HGLOBAL hData = LoadResource(hinst, hRes);

const void *pData = LockResource(hData);

HGLOBAL hBuffer = GlobalAlloc(GMEM_MOVEABLE, Size);

void *pBuffer = GlobalLock(hBuffer);

CopyMemory(pBuffer, pData, Size);

GlobalUnlock(hBuffer);

IStream *pStream;

CreateStreamOnHGlobal(hBuffer, TRUE, &pStream) ;

IStream *pStream=CreateStreamOnResource(hinst,pszName,pszType);

Gdiplus::Bitmap* pBitmap = Gdiplus::Bitmap::FromStream(pStream);

pStream->Release();

HBITMAP ret;

pBitmap->GetHBITMAP(Color(0,0,0,0),&ret);

delete pBitmap;

return ret;

}

Just another 靱技術士事務所 site

月別: 2013年6月

統計ソフトRを使用して株価チャートを簡単に作成(3)

PlayFramework2.1で既存のMySQLにつないでみる

統計ソフトRを使用して株価チャートを簡単に作成(2)

統計ソフトRを使用して株価チャートを簡単に作成(1)

xcodeでフレームワーク作成

カーソルがウインドウの外に出たのを検知する

GDIでPNGを読み込む