タグ

ブックマーク / okumuralab.org/~okumura (9)

  • 東京都のCOVID-19

    東京都の情報は2021-05-24以降は毎日16:45に東京都福祉保健局からPDFで発表されているが,手作業での更新のようで,うまくスクレープできないこともあった。そうこうしているうちに厚労省が16:00ごろ全国の都道府県についてオープンデータを更新することになったので、そちらを取得して COVID-tokyo.csv を作成し,自動ツイートしている(描画部分のコード:COVID-tokyo.py)。そのグラフ: 実効再生産数を,(直近7日間の陽性数/n日前7日間の陽性数)**(5/n)で近似したもの。東洋経済オンラインでは n=7 としているが,ここでは 1 から 7 までの7のグラフを示した。 全区間にわたってプロットしたもの: 縦軸を対数目盛にしたもの: 奥村 晴彦 Last modified: 2022-12-31 16:26:39 JST

  • ライフゲーム

    統計とは関係ないが,JavaScriptで書いたライフゲームをRに移植してみた。9行で書けた。 a = b = matrix(0, nrow=50, ncol=50) a[24:26,25] = a[25,24] = a[26,26] = 1 for (g in 1:300) { b[2:49,2:49] = a[1:48,1:48] + a[1:48,2:49] + a[1:48,3:50] + a[2:49,1:48] + a[2:49,3:50] + a[3:50,1:48] + a[3:50,2:49] + a[3:50,3:50] a = ifelse(b != 2, b == 3, a) heatmap(a, Rowv=NA, Colv=NA, scale="none", col=0:1, main=g) }

    Nyoho
    Nyoho 2016/06/27
  • アイスクリームの売れ方

    問題 2016年1月17日の入試センター試験「数学①」の「数学I」(全問必答)の第4問および「数学I・数学A」の第2問(必答問題)で,2003〜2012年の120か月の各世帯の1日あたりアイスクリーム平均購入額と,東京の月別気象データの相関を調べる問題が出題された。 データを調べる アイスクリーム平均購入額については,「総務省統計局(2013)『家計調査年報』」が出典と書かれているが,Googleで検索して出てきた家計調査年報(家計収支編)平成25年(2013年)を見てもどこにあるのかわからない。 そこで,e-Stat のキーワード検索で「アイスクリーム」を入れる。一覧から「家計調査」→「家計収支編」をクリック,左で「データベース」を選び「表題一覧表示」をクリック,「月次」があるのは上の三つだけなので,最新のもの(一番上)の「DB」を選ぶ。 「1/5 表章項目」で「選択」をクリック,「金額

  • 猛暑日の増加

    [2018-07-20] as.POSIXct() を as.Date() に変えました。 東京管区気象台の東京における雷日数や真夏日等の日数の変化というページを見ると,猛暑日(日最高気温35℃以上)の日数が1876年以来著しく増加していることがわかる。ところが,平均気温を見ると,地球温暖化・都市化の影響で徐々に増加しているが,激しく変わっているようには見えない。平均ではなくて猛暑日の日数を見るべきであるという議論もある。しかし,「猛暑日」(日最高気温35℃以上)や「真夏日」(日最高気温30℃以上)のような分布の裾の日数は,何℃で切るかによって印象が大きく異なる。気温変化を誇張することにならないか。 「東京」の観測地点は2014年12月2日に約900m離れた地点に移転している。詳しくは「東京」の観測地点の移転について(PDF,2014年11月14日,気象庁観測部)を参照されたい。 東京の日

    Nyoho
    Nyoho 2015/08/09
    気象データ から 奥村先生のページ
  • 統計グラフの色

    [TODO] R 4.0.0 以降の色についてはそのうち書きます。とりあえず palette() のヘルプをご覧ください。palette("Okabe-Ito") とするとsafe colorsになります。→ A New palette() for R,Color Universal Design,Essentials of color in R。あと Colorspace 2.0 とその JSS paper。 はじめに 統計グラフに色を付けることは広く行われています。しかし,色は万人に共通のものではありません。日人男性の5%,白人男性の8%は,RGB(赤緑青)のうち赤と緑の区別がうまくできません。その内訳は1:3で赤の感受性がないP型(1型,protanopia)と緑の感受性がないD型(2型,deuteranopia)に分かれます。青を感じない人や,RGBのうち2色以上を感じない人もい

  • Rの初歩

    Rのインストール R(アール)は統計・データ解析・統計グラフ作成のためのオープンソースソフトです。CRAN(「シーラン」または「クラン」、Comprehensive R Archive Network)のミラーサイト(日では統数研 https://cran.ism.ac.jp/ など)からダウンロードし、インストールしてください。Windows用、Mac用、Linux用があります。 最近では下図のような RStudio(アール・スタジオ、英語の発音により忠実に書けばアール・ステューディオまたはアール・ストゥーディオ)という統合開発環境(IDE)を介してRを使う人が増えています。Download RStudio Desktop のページからダウンロードしてインストールします。 RStudio。左上ペインがソースエディタ(この図では隠れている)、左下がRコンソール、右上がグローバル環境のオブ

    Nyoho
    Nyoho 2012/05/20
  • 中学生のコンピュータの問題 | Okumura's Blog

    Twitterで教えていただいた中学の技術・家庭科副読を買ってみた。以下はこれに載っている問題の一部である。 文字,数字,図形,静止画像,動画像,音声などを(情報)という。 漢字変換は,(日語変換)システムで行っている。 インターネットに接続するには,プロバイダと契約し,メールアドレス以外に何を取得する必要があるか。(パスワード) まだまだある。Twitterを検索されたい。

    Nyoho
    Nyoho 2010/12/14
    FEP って出てくるのはおっさんだなー (わしも浮かんだ
  • PDFのフォント制限は無意味 | Okumura's Blog

    情報処理学会全国大会の原稿送信要項には,原稿(PDF)の日フォントWindowsではMS明朝またはMSゴシック,Macでは細明朝または中ゴシック・平成明朝または平成角ゴシック,英字フォントはArial,Century,Times,Times New Roman,Helvetica,Symbolに限定するように書かれている。TeXではこれだけに限定することは至難の業である。ちょっとした数式を書けば数式用フォントが使われてしまう。そもそもなぜ「これ以外のフォントを使用されると,DVD-ROMを利用する環境によっては文字化けや表示エラーを起こすことがあります」なのか理解ができない。 情報処理学会の研究報告原稿(PDFファイル)作成についてには,LaTeXのスタイルとWordのテンプレートが用意されており,PDFファイルには必ずフォントを埋め込むよう指示がある。この「必ずフォントを埋め込む」

    Nyoho
    Nyoho 2010/03/08
  • Webで数式を簡単に使う方法 | Okumura's Blog

    以前Webで数式を書く方法について書いたが,今なら Google Chart Tools のAPIを使うほうが簡単。例: <img src="http://chart.apis.google.com/chart?cht=tx&amp;chl=x=\frac{-b\pm\sqrt{b^2-4ac}}{2a}" alt="" /> URL中に使えない文字は%16進2桁で表記する。例えばスペースは試した限りではそのまま使えたが用心するなら%20とする。詳しくはGoogleの解説 Mathematical (TeX) Formulas 参照。 問題点:ピクセルサイズの指定はできるが,標準の2倍の大きさにしたいという指定はできないみたい。

    Nyoho
    Nyoho 2010/03/04
  • 1