mokjpn

---
title: "Rでe-Statの統計表から都道府県別SMRを描く"
output: html_document
---


このテキストでは、政府統計e-StatのAPIを使って、Rに直接統計表を読み込み、それを処理して特定の疾患（ICD-10に基づいた死因簡単分類で指定）の都道府県別の標準化死亡比(Standard Mortality Ratio, SMR)を算出し、その値によって都道府県を塗り分けた地図をかいてみます。すべてRの中だけで完結し、またプログラムの長さも短くて済んでいることに注目してください。

```{r, echo=FALSE, results='hide'}
appId <- "0000000000000000000000000000000000000000"

mokjpn

from m5stack import lcd
from mqtt import MQTTClient
import time
import ujson
import ussl
import socket
import machine

MACHINE=hex(int.from_bytes(machine.unique_id(), 'big'))[2:]

mokjpn

## 同じ整数の繰り返しである数列に含まれる数を１つずつ使い、目的の数を作るための式を総当たりで求める
## Usage: npuzzle(2,4,24) => (2,2,2,2)を１つずつ使って24を作るにはどうするか．
npuzzle <- function(num, nrep, target, twofuncs = list(
	"*"=function(x,y){x*y},
	"/"=function(x,y){x/y},
	"+"=function(x,y){x+y},
	"-"=function(x,y){x-y},
	"^"=function(x,y){x^y}), onefuncs = list(
	" "=function(x){x},
	"factorial"=function(x){factorial(x)},

mokjpn

## いわゆる目標症例数とリクルートの進捗状況のグラフ

# ダミーの日付データを作る
# リクルート日付は2014-01-01からの２年間からランダムに１００日
dates <- as.Date("2014-01-01") + 365*2*runif(100)

# それぞれの日に１０人までのランダムな人数がリクルートされたとするダミーデータ、100日分
number <- floor(10*runif(100))

# １行１症例のデータセットのイメージ。登録日だけだけど。

mokjpn

## いわゆる目標症例数とリクルートの進捗状況のグラフ

# ダミーの日付データを作る
# リクルート日付は2014-01-01からの２年間からランダムに１００日
dates <- as.Date("2014-01-01") + 365*2*runif(100)

# それぞれの日に１０人までのランダムな人数がリクルートされたとするダミーデータ、100日分
number <- floor(10*runif(100))

# １行１症例のデータセットのイメージ。登録日だけだけど。

mokjpn

## いわゆる目標症例数とリクルートの進捗状況のグラフ

# ダミーの日付データを作る
# リクルート日付は2014-01-01からの２年間からランダムに１００日
dates <- as.Date("2014-01-01") + 365*2*runif(100)

# それぞれの日に１０人までのランダムな人数がリクルートされたとするダミーデータ、100日分
number <- floor(10*runif(100))

# １行１症例のデータセットのイメージ。登録日だけだけど。

mokjpn

## This gist is obsoleted. Refer https://github.com/mokjpn/Define2Validate/ for the new version.

library(R4DSXML)
library(testthat)
library(validate)

define2validate <- function(domain, file="exampleRules.yaml", definexml="Odm_Define.xml") {
  varmd <- subset(getVarMD(definexml), IGD_Name == domain)
  cat(file=file, append=FALSE, "")
  out <- function(...) cat(file=file, append=TRUE, paste(...,"\n",sep=""))

mokjpn

# Read MediaWiki's WikiTable text and convert it into a dataframe

lines <- readLines(file.choose())

columns <- NULL
row <- NULL
df <- data.frame()
ncol <- 0
for(line in lines) {
  if((m <- sub("^\\! *(.*)$", "\\1", line)) != line){

mokjpn

(to name gist)

mokjpn

/* SpeedButton Sketch */

/* NeoPixel */
#include <Adafruit_NeoPixel.h>
// Digital pin to connect NeoPixel
#define PIN            0
// How many NeoPixels are attached to the Arduino?
#define NUMPIXELS      16

Adafruit_NeoPixel pixels = Adafruit_NeoPixel(NUMPIXELS, PIN, NEO_GRB + NEO_KHZ800);