Rで描画回りははまるので、なるべく使わないようにしているのですが、 やっぱりRが便利なのでつい使ってしまってはまってしまいます。 良く忘れるので備忘録にメモしておきます。 ベクターの場合・・・cutsを使う。 x <- data.frame( x=c(141.1,141.1,141.1,…

次年度予算を積算しておりますが、例えば以下のようなテキスト入力用のフォームから数値のみを抽出して合計する関数を検討。 ハードディスク 20千円 研究打合せ旅費（東京−仙台、1回、2泊3日、1人）50千円 研究打合せ旅費(東京−京都、1回、2泊3日、1人)50千…

Rをあたらしくインストールした際に、以前のバージョンで使用していたパッケージを一括でインストールしたい。 #一つ前のRバージョンでパッケージがインストールされているフォルダを確認しメモする Sys.getenv("R_LIBS_USER") #最新バージョンのRを起動して…

久々に描画周りを使うとわすれていることが多くて悩ましい。 library(sp) #ポイントデータの場合 data(meuse) coordinates(meuse)=~x+y spplot(meuse, "cadmium", do.log = TRUE, key.space=list(x=0.2,y=0.9,corner=c(0,1)), cuts=c(1, 2, 3, 4, 10, 20),sc…

https://cran.r-project.org/web/packages/oce/vignettes/oce.htmlあまり日本語の情報がないので、どなたかまとめていただけると助かります。 CTDデータの読み込み library(oce) #デモデータ x <- read.oce(system.file("extdata", "ctd.cnv", package="oce"…

久々にSOC関連の文献を見ていたら、そこそこ進展があるようです。 というかこのグループの論文のプロダクティビティとIFの高さはすごい、こんなに読み切れない。 一番の進展はCASCADEというプロジェクトが2012年ごろからはじまって2017年までつづくようです…

よく困って忘れたりするので、ここに気づいたことをメモしていきたいと思います。関連するR-jpwiki など http://cse.naro.affrc.go.jp/takezawa/r-tips/r/18.html http://r.livedocs.net/vector/nan.html is.na()でNAを操作できないこと、 logical(0)は対応…

緯度 子午線弧長の平均約40008km。 平均的な長さは 1度 約111km 1分 約1.85km 1秒 約30.9m GRS 80楕円体表面上の代表地点 緯度00度上 1秒 約30.7m 1度 約111km 緯度35度上 1秒 約30.8m 1度 約111km 緯度90度上 1秒 約31.0m 1度 約112km http://ja.wikipedia.…

査読者にGAMの%Deviance explainedが50%に近いということはランダムな要因を入れた場合と変わらないのではないかという指摘をうける。それは勘違いじゃないかとおもって、しらべてみる。 計算方法 これ。 (1- Residual Deviance/Null Deviance)*100 GLMでもd…

Example 13. read.csv.sql and read.csv2.sql read.csv.sqlはread.csv＋sqlの引数が使えるsqldfのインターフェース 1)SQLiteのimport facility via RSQLite 2)provided SQL statement 直接RSQLiteで最初に読むのでRが扱えるサイズにまでフィルタリングすれば…

library(sqldf) #省略名と値のデータを作成 Abbr <- data.frame( Species = rep(levels(iris$Species),3) ,Abbr = rep(c("S", "Ve", "Vi"),3) ,Petal.Width = c(rep(0.1,3),rep(1.5,3),rep(2,3)) ) str(iris) #（なぜか.が_になるのを無視して以下のsql式を…