Yusuke Saito ysaito8015

ysaito8015 / calcTaronesTest.R

Last active August 29, 2015 14:01

層によるオッズ比が異なるかどうかを検定するには,Tarone の検定を利用する. Taroneの検定スクリプト

	##metaforパッケージをインストールする
	#install.packages("metafor")
	#library("metafor")

	calcTaronesTest <- function(mylist,referencerow=2)
	{
	numstrata <- length(mylist)
	# 各層における曝露グループの人数のアレイ"ntrt"を作る
	# 各層における非曝露グループの人数のアレイ"nctrl"を作る
	# 各層における疾病を持つ曝露グループの人数のアレイ"ptrt"を作る

ysaito8015 / calcMHRatio.R

Last active August 29, 2015 14:01

対症例対照研究(Matched Case-Control Study) Mantel-Haenszel のオッズ比計算スクリプト

	calcMHRatio <- function(mymatrix, alpha=0.05)
	{
	caseExposedControlUnexposed <- mymatrix[1, 2]
	caseUnexposedControlExposed <- mymatrix[2, 1]
	MHRatio <- caseExposedControlUnexposed / caseUnexposedControlExposed
	print(paste("Mantel-Haenszelのオッズ比 =", MHRatio))
	# 信頼区間の計算
	confidenceLevel <- (1 - alpha) * 100
	sigma <- sqrt((1 / caseExposedControlUnexposed) + (1 / caseUnexposedControlExposed))
	# sigmaは,オッズ比の対数の推定値の標準誤差

ysaito8015 / calcRerativeRist.R

Created May 11, 2014 01:21

相対危険度スクリプト

	calcRelativeRisk <- function(mymatrix,alpha=0.05,referencerow=2)
	{
	numrow <- nrow(mymatrix)
	myrownames <- rownames(mymatrix)
	for (i in 1:numrow)
	{
	rowname <- myrownames[i]
	DiseaseUnexposed <- mymatrix[referencerow,1]
	ControlUnexposed <- mymatrix[referencerow,2]
	if (i != referencerow)

ysaito8015 / calcOddsRatio.R

Created May 11, 2014 01:19

オッズレシオ計算スクリプト

	calcOddsRatio <- function(mymatrix,alpha=0.05,referencerow=2,quiet=FALSE)
	{
	numrow <- nrow(mymatrix)
	myrownames <- rownames(mymatrix)
	for (i in 1:numrow)
	{
	rowname <- myrownames[i]
	DiseaseUnexposed <- mymatrix[referencerow,1]
	ControlUnexposed <- mymatrix[referencerow,2]
	if (i != referencerow)

ysaito8015 / file0.txt

Created December 3, 2013 06:23

ubuntu 13.10　で capslock を ctrl　へ変更する ref: http://qiita.com/ysaito8015@github/items/0fd4ec3076293bc372d4

% xmodmap -pke> "~/.xmodmap"