楽しい（無駄な）Rアート、その２

R script

楽しい（無駄な）Rアート - こにしき（言葉、日本社会、教育） ←この続編。こころの赴くままに描いてみました。完全オリジナルです。 Script x <- runif(10000) y <- runif(10000) col1 <- ifelse( ((y-.5)^2 + (x-.5)^2 >0.3)&((y>.5)&(x<.5)), "#804000"…

2015-02-10

楽しい（無駄な）Rアート

R script

n <- 2200 X <-runif(n,-1,1); Y <-runif(n,-1,1.5) a <- ifelse( (Y>(X^2)^(1/3)-sqrt(1-X^2))&(Y

2014-05-22

ロト6の当選番号を予想するR関数

R script

いつもロト6が当たらない、とぼやいている人がいたので、いたたまれなくなって、作ってみました。当たるといいですね！ LOTO6 <- function(x,y,z,v,w){ x <- order(runif(43))[1:6] x[order(x)] } 使い方 Rを起動して LOTO6() のカッコの中に、必要事項を記…

2013-07-13

信頼区間をあわせてプロットするR関数

R script

2015年9月7日追記中澤先生のツイートで知りましたが、信頼区間をプロットする関数を含んでいるパッケージがすでに存在するようです。 http://www.inside-r.org/packages/cran/plotrix/docs/plotCI こちらを使うことをお勧めします。 plotrixパッケージにplo…

2012-04-10

Fukuta's U を算出するRスクリプト

R script

UverWorldに似ていると言われた回数から、Fukuta's U を返すRスクリプトを作りました。ご自由にお使い下さい。スクリプト以下をコピー＆ペーストしてください FukutaU <- function(c){ 0.5 / (exp(1/c) - 0.5) } Fukuta's Uの概要 UverWorld のボーカルの…

2012-01-11

SSM職業コード（数字3桁）を職種名に変換するR関数

R script

だれがこんな関数使うんだよのコーナー recode() 関数を利用しているので、 > library(car)で「car パッケージ」を読み込んでおいてください。たぶん、深刻な間違いはしていないと思いますが、たとえば「そんな職業名じゃない！」など、なにかお気づきの点…

2011-12-24

AZC値を算出するRスクリプト（バンドメンバーの独占度を示す指標）

R script

スクリプト以下をコピー＆ペーストしてください azc <- function(m){ sum(max(m/sum(m)) - (m/sum(m)) )/( length(m)-1) } 定義あるバンドにおいて、特定の人物（一般的にボーカル）に人気・知名度・貢献度が一極集中している度合を、「安全地帯」にちなみ…

2011-05-12

雑誌『英語教育』のトレンド分析（および grep()関数について）

R script 雑誌『英語教育』研究ネタ

この記事（雑誌『英語教育』（戦後）の量的傾向をざっと見る - こにしき（言葉、日本社会、教育））の続きです。前述の『英語教育 Fifty』附録の記事タイトルのデータはエクセル形式で提供されていたはずなので、表計算ソフトなどをつかえば、ある用語のトレ…

2010-12-11

有意確率をアステリスクに変換するRスクリプト

R script

hoshi(x) hoshi <- function(x,ns){ X <- NULL for(i in 1:length(x)){ if(x[i] < .001){X[i] <- c("***") }else{ if(x[i] >= .001 & x[i] < .01){X[i] <- c("**") }else{ if(x[i] >= .01 & x[i] < .05){X[i] <- c("*") }else{ if(x[i] >= .05 & x[i] < .1 )…

2010-07-11

日本語文の名詞をランダムに「ヒャッハァー！」に置換するＲスクリプト

R script

RMeCabのインストールが必要 → library(RMeCab) 追記 r-de-rさんに助言いただきました。ありがとうございます。もっと洗練されたコードを提案されたので謹んでコピペ（すみません）します。 hyahhaaa2 <- function(x, y) { b <- unlist(RMeCabC(x)) b[names(…

2010-05-11

名詞の後ろに「（笑）」をつける関数

R script

あらかじめ、形態素解析エンジン和布蕪（MeCab、めかぶ）および、そのR用のパッケージとして、 RMeCab のインストールが必要です。詳細は例えばこちらをご参考にしてください→"RMeCab作成者・石田基広先生による解説（PDF）" （笑）添付関数：RKakkoWarai(…

2010-03-03

リスト形式のオブジェクト重複を削除する

R script

リスト形式をtransactions形式に変換する際、重複があるとエラーを吐き出すらしい。（こちらのページを参考にさせていただきました。多謝。→http://d.hatena.ne.jp/hiru926/20100103/1262508600）というわけで、ヘルプなどを探し回たんだけど、私の「検索…