読者です 読者をやめる 読者になる 読者になる

ルビスコくん空を飛ぶ

博士学生(D1)の日記です.日記に加えてプログラミング(R)のメモ書きもします.

R

パッケージ'plotflow'を使って2軸表示したときのバグ

R

以前「'plotflow'パッケージを使えばggplot2で2軸表示ができる」という記事を書いたのですが,よくわからないバグ(?)が発生するようになってしまいました。以前にはなかったのですが。 #データフレームの準備 dd <- 1:10 * 10ee <- 1:10 / 20ff <- 1:10dd…

回帰直線を引いて同時に回帰式も図に表示させる

R

わかりやすくするために図の中に式を埋め込みたいけど,式をRで調べていちいちパワポやイラストレータなどで文字を打ち込むのは時間の無駄なのでRで(ggplotで)やってしまいましょう.今回もデータセットはairqualityを使います. この中のOzoneとTempの回…

Rで素数判定プログラム

R

気になる自然数を入力するとその数が素数かどうかを判別してくれるプログラムを,思いつきでつくってみました.その理屈はいたって簡単.自然数Xを2から順番に割っていって,初めて割り切れた数がXなら素数ですから"X is a prime number!",それ以外なら"X i…

read.csv()の段階で列名を指定する

R

col.names = で指定できる.read.csv(x, header = F, skip = 1, col.names = c("T-Yamada", "T-Okada", "T-Kanemoto"))

2種類のカテゴリデータにそれぞれ色とプロットのマーカーを割り当てて散布図を描きたい

R

「変数XとYの関係に,別の変数ZとWが与えている影響をみたいとき」はXとYをプロットしZとWで色とマーカーを変えた図を出すという方法がある(ただし図が煩雑にはなるのは想像に難くない).こんな図は,描画パッケージggplot2を使えば簡単にできる. 例に使…

パッケージggplot2による2軸表示

R

先日,Rのパッケージggplot2のグラフ2軸機能 - ルビスコくん空を飛ぶでggplot2ではグラフの2軸は表示できないという記事を書いたが,どうやらできるようになったらしい.「plotflow」パッケージをインストールするために,以下のようなコマンドを打ち込む. …

文字列から特定の文字列を削除する

R

"hogehoge.txt"という文字列から".txt"を削除する方法.list.files()などでファイル名を取り出したのはいいが,拡張子は削除したい,というときに便利. gsub()を使う. 一つ目の引数に検索文字列,次の引数に置換後の文字列,最後に置換したい文字列を入力…

データセットにおける数値のカテゴリ化

R

データセットのなかのある変数でデータを区切り,区間化してカテゴリに変換する方法.Rにデフォルトで組み込まれているデータセットairquality(ニューヨークの大気状態観測値.1973 年5~9月の1日ごとに4つの変数がある.変数は下図の通り)を用いる.今回は…

Rの関数作成

R

Rを使ってデータ解析をしていると,「標準誤差を出す」とか「何らかの曲線で回帰して,r2値・係数・切片を求める」とか同じ作業が多い.そこでいちいち res_a <- summary(lm(y1 ~ x1))$coef[1] #直線による回帰式の傾きを得る res_b <- summary(lm(y1 ~ x1))…

Rのggplot2でのエラーバーの表示

R

Rのパッケージggplot2で平均±標準偏差をプロットする方法のメモ. データセットはirisを使う.まず,doByパッケージのsummaryBy()でデータをまとめる.この関数は,自分が指定した変数(下のコードではSpecies)についてクラス分けをして,クラスごと&指定…

Rのパッケージggplot2のグラフ2軸機能

R

ggplot2は,エレガントで効率的に図を生成してくれるRのパッケージである.最近はデフォルトのplot()関数よりもggplot2を使って図を描くことが多くなってきた. だが,そんなggplot2にも弱点があるようで,y軸を二つに分けることができないらしい.どういう…

source()関数

R

RのコマンドラインでRスクリプトを実行する,あるいはRスクリプトの中で別のスクリプトを実行したい場合,source()関数を使う. そこで,いちいち別のスクリプトがある場所にカレントディレクトリをsetwd()で移さなくても,source()の中にパスごと書いてしま…

Rのggplot2

R

Rメモ forループでまわしながらggplot2でグラフを大量生産する作業をしていたのだが,何度やってもグラフが指定したフォルダに保存されないという問題が生じた.結局,以前作った同じ名前のファイルを削除したら保存できた(このことに気づくのに1時間ほど…