Rのggplot2でグラフを作りたい - 300億円欲しい
Rおじさんだよ 描画が綺麗にできるRパッケージ, ggplot2を使いたいです. 紹介のために, ローレンツ方程式の数値解を描画します. ローレンツ方程式とは. http://en.wikipedia.org/wiki/Lorenz_system ローレンツ方程式はカオス的な振る舞いをする有名な問題です. \begin{align*} \frac{\mathrm{d}x}{\mathrm{dt}} &= ax + by \\ \frac{\mathrm{d}y}{\mathrm{dt}} &= b(y-z) \\ \frac{\mathrm{d}z}{\mathrm{dt}} &= -xy + cy -z \end{align*} パラメータを神調整すればカオスアトラクタが見られます. Rで常微分方程式を解く RのdeSolveパッケージは, 適当な微分方程式を指定した解法で解いてくれます
R Internals
R Internals ¶ This is a guide to the internal structures of R and coding standards for the core team working on R itself. This manual is for R, version 4.3.3 (2024-02-29). Copyright © 1999–2023 R Core Team Permission is granted to make and distribute verbatim copies of this manual provided the copyright notice and this permission notice are preserved on all copies. Permission is granted to copy an
俺が思ってたんと違う挙動の遅延評価 - My Life as a Mock Quant
ある程度コードを抽象化しようと言う時に「関数を返す関数」を使うのは良いアイディアだと思うのですが、以下のように書くとハマれるってのでメモ。 ここでは「異なるパラメーター(1〜3)を持つ3つの関数を生成し、それぞれに対して実際の計算を行うものの…期待する結果とはずれる」という話です。 コードはこんな感じ。 #関数返す関数 make.hoge <- function(m){function(x)x*m} #1*x, 2*x, 3*xとなるような関数リストを返してる（つもり） hoges <- sapply(1:3, make.hoge) #1を期待 hoges[[1]](1) #2を期待 hoges[[2]](1) #3を期待 hoges[[3]](1) これを直すにはただ1行、「変数mを強制的に評価する」ようなコードをmake.hoge関数の中に追加すれば良い。これでちゃんと直る！ #関数返
R言語上級ハンドブック - sfchaos's blog
R言語上級ハンドブックを一通り読了しました． R言語上級ハンドブック 作者: 荒引健,石田基広,高橋康介,二階堂愛,林真広出版社/メーカー: シーアンドアール研究所発売日: 2013/09/25メディア: 単行本（ソフトカバー）この商品を含むブログ (9件) を見る 本書は，Rで効率的・効果的にデータ解析や処理を行うための中級〜上級のTips集です．著者はR界隈で著名な方々ばかりで，高度なトピックが平易に解説されています．内容も，実行パフォーマンス(処理速度，メモリ使用量等)，Hadoopや他言語との連携，グラフィクスなど多岐に渡ります．また，knitrやslidifyによるレポート・プレゼン資料作成，shinyによるWebアプリ作成など，比較的新しい話題についても収録されています． Rで分からないことをStack Overflowなどで調べているユーザにとって，本書は座右の書となることで
The R Inferno - Burns Statistics
What is it? Abstract: If you are using R and you think you’re in hell, this is a map for you. A book about trouble spots, oddities, traps, glitches in R. Many of the same problems are in S+. Why would I care? Even if it doesn’t help you with your problem, it might amuse you (and hence distract you from your sorrow). Getting it You can get a free pdf of The R Inferno. Alternatively you can spend a
Welcome · Advanced R.
This is the companion website for “Advanced R”, a book in Chapman & Hall’s R Series. The book is designed primarily for R users who want to improve their programming skills and understanding of the language. It should also be useful for programmers coming to R from other languages, as it explains some of R’s quirks and shows how some parts that seem horrible do have a positive side. Introduction F
pixivのタグ頻度から考えるラブライブのカップリング - 驚異のアニヲタ社会復帰の予備
進撃の巨人を読んだことない人がデータだけでキャラを推測してみる - あんちべ！というのが面白そうだったので、毎度おなじみのラブライブ!でネットワーク解析をしてみる。 pixivの小説からタグを取ってくる。作業は pixiv というフォルダとする。 # Rで tag <- "ラブライブ" page <- 1000 urls <- paste("http://spapi.pixiv.net/iphone/search_novel.php?s_mode=s_tag&p=", seq(page), "&word=", tag, sep="") write.table(urls, "urls.txt", row.names=FALSE, col.names=FALSE, quote=FALSE) こうして出来た urls.txt を wget を使って取得する。 # ターミナルで cd pixiv
zooパッケージを使って行列の欠損値を補間する(R Advent Calendar 2011) - My Life as a Mock Quant
R Advent Calendar 2011 : ATNDに参加中でして、その担当分の記事です。 さて掲題の話ですが、たまに（しょっちゅう？）歯抜けになってる値を含む行列やデータを捌かなきゃいけないことがありまして、その際にはよくRの市販本にあるように x[is.na(x)] <- 0 と欠損値の箇所をある”特定の値”で置いてしまうことが多かったのですが、 もうちょっと行列（あるいはデータフレーム）の特徴を活かした値の補間ができないものかと考えていた所、 普通は時系列処理に使う”zooパッケージ”に入ってるna.locf等の補間系関数を使えばもっと色々とできるじゃないかと思いつきやってみたという内容です。 何はともあれまずはパッケージのインストール＆読み込み install.packages("zoo") library(zoo)ここではサンプルコードとしてちょうど真ん中の要素に欠損値を持
foreach＋doSNOWパッケージを使って、並列処理をやってみた - My Life as a Mock Quant
なんだかweb屋界隈ではHadoopだのMahoutだの象を使うだの使わないだの楽しそうな事をやっていて、私はとてもさみしく、そしてうらやましくもなったわけですが手元にそんな立派な環境なんてないわけで。しかし、そんな私にもマルチコアのパソコンが与えられているのでそれで計算の並列化をやってみた。OSはwindows XP(32bit)。 使うのは統計数理研究所のセミナーでお薦めされていたforeachパッケージ*1。これはRevlution Rの人たちが作っているので今後とも継続的な開発と進化が期待できるとのこと。まずは、doSNOWパッケージも合わせてインストール。 install.packages("foreach") install.packages("doSNOW") foreachパッケージではforループの並列化をしてくれるのですが、処理を並列化させないで書くこともできるのでまず
read.table() でヘッダーをそのまま読み込みたい - Qiita
read.table() でファイルをデータフレームに読み込むとき、ヘッダーに R の変数として使用できない文字があると . に自動的に変換されます。 例えば、下記のようなデータがあったとき、 Name Value(1hours) Value(2hours) Value(3hours) #ofCount
連続する要素の数を数える - 驚異のアニヲタ社会復帰の予備
ベクトル中の連続する要素の連続数は rle で計算できる。 いま、ある塩基配列 s1 <- c("ctgactgagagactttacgtaggagggcccatcgcgccagcggggtgctaaccactagcatgggtaacaccatcacatgctatgtgaaagccctagcggcctgcaaggctgcggggatagttgcgcccacgaaagccaggggactgagaggacgagcggaacctgagggaggccatgagtactctgcccctcctggtgat", "ccccccagaccggaatatgcaaatgtgtccggggccgccgcagatactacctgaccagagacccaaccactccaccccgggctgcctgggaaacagttagacactcccctatcaattcatggctgggatatgggttcgcatggtcct
factorを使ってtableしたら、ないものの数を0としてカウントしてくれる - 驚異のアニヲタ社会復帰の予備
こちらで、ある配列中のAGCTの数を数えている。 要素の数を数えるのは table がよいが、配列にないものはカウントされない。 s0 <- c("A", "G", "G", "C", "C", "C") table(s0) s0 A C G 1 3 2 リンクの行列にapplyで適応すると、長さがまちまちのリストとして結果が返ってくる。 これを回避するには、配列をfactor化する。 base <- c("A", "G", "C", "T") factor(s0, base) [1] A G G C C C Levels: A G C T table(factor(s0, base)) A G C T 1 2 3 0 こうすると、結果が全て長さ4で揃うので、applyした結果は行列になってくれる。
Rにおける遅延評価とスコープ - もうカツ丼はいいよな
http://d.hatena.ne.jp/Rion778/20091210/1260464593 去年書いたんだけど，今ふいに遅延評価ってなんだっけとか思って調べたり上の記事を読み直したりしてみたら割と何書いてあるか分からない理解してなかったっぽいのでもう一度まとめ直し． ただ個人的には納得したものの，内容的には前のと大差ない気がします． まず次の例を見てもらいたい． > ## 一見するとflist[[1]]()が1，flist[[2]]()が2, ... , flist[[10]]()が10を返しそう > flist <- lapply(1:10, function(i){ function() i}) > ## ダメ... > flist[[1]]() [1] 10 > flist[[2]]() [1] 10 別にこれはバグとかではなくて，なぜこのような事が起こるのかというと，遅延評
Overview
The R Journal is the peer-reviewed, open-access scientific journal, published by The R Foundation. It features short to medium length articles (up to 20 pages) covering topics of interest to users and developers of R. For example, description and implementation details of R packages, surveys or reviews of R software related to a particular area, applications of R software or descriptions of teachi
JIN'S PAGE
R、R言語、R環境・・・・・・ Rのダウンロードとインストール リンク集 題名 Chap_01 データ解析・マイニングとR言語 Chap_0２ Rでのデータの入出力 Chap_03 Rでのデータの編集と演算 Chap_04 Rと基本統計量 Chap_05 Rでの関数オブジェクト Chap_06 Rでのデータの視覚化(1) Chap_07 Rでのデータの視覚化(2) Chap_08 Rでのデータの視覚化(3) Chap_09 GGobiとデータの視覚化(Rgobi) Chap_10 Rと確率分布 Chap_11 Rと推定 Chap_12 Rと検定 Chap_13 Rと分散分析 Chap_14 Rと回帰分析 Chap_15 Rと重回帰分析 Chap_16 Rと一般化線形モデル Chap_17 Rと非線形モデル Chap_18 Rと判別分析 Chap_19 Rと樹木モデル Chap_20 WEK
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
R で一部の関数のメッセージの言語が実行条件によって変わる理由とその回避策 - あらびき日記
R Language Definition: 6. Computing on the language
Rcpp
Bio3D in R / FrontPage
Tokyo.Lang.R #0 で「Rのスコープとフレームと環境と」を発表してきました - あらびき日記
