[表示 : 全て 最新50 1-99 101- 201- 301- 401- 501- 601- 701- 801- 2chのread.cgiへ]
Update time : 04/17 13:22 / Filesize : 381 KB / Number-of Response : 837
[このスレッドの書き込みを削除する]
[+板 最近立ったスレ&熱いスレ一覧 : +板 最近立ったスレ/記者別一覧] [類似スレッド一覧]


↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

【R言語】統計解析フリーソフトR 第4章【GNU R】



1 名前:132人目の素数さん mailto:sage [2011/01/09(日) 17:31:49 ]
R は統計計算とグラフィックスのための言語・環境です。
統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。

●関連サイト
The R Project
www.r-project.org/
RjpWiki
www.okada.jp.org/RWiki/
リンク集
www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8

820 名前:132人目の素数さん mailto:sage [2013/04/12(金) 10:13:54.71 ]
ちょっと教えてくれろ。

x <- rnorm(100);hist(x, prob=T, ann=F) ;par(new=T) ;plot(density(x),col="red")
こんな感じの正規分布っぽいグラフでXの範囲指定して、
Y軸0から確率密度の赤線までの面積を表したいんだけど、
ここを色づけする方法が判らない、
例えばXが1~2の間を塗るとき具体的にどんなコードになるか教えて下さい。

821 名前:132人目の素数さん mailto:sage [2013/04/12(金) 15:37:00.07 ]
ペイント使って塗ったらいいじゃん

822 名前:132人目の素数さん mailto:sage [2013/04/13(土) 00:06:14.95 ]
>>820

dat <- with(density(rnorm(100)),data.frame(x,y))

library(ggplot2)

ggplot(data = dat, mapping = aes(x = x, y = y)) +
geom_line() +
layer(data = dat, mapping = aes(x=ifelse(x >1 & x < 2 ,x,0), y=y), geom = "area", geom_params=list(fill="red",alpha=0.5)) +
scale_y_continuous(limits = c(0,max(dat$y)))

823 名前:132人目の素数さん mailto:sage [2013/04/13(土) 10:48:56.27 ]
>>822
ありがとう

824 名前:132人目の素数さん mailto:sage [2013/04/13(土) 23:48:23.04 ]
ここがRの本スレか。
R自体統計プログラムなのかグラフィックソフトウェアなのかよくわからんけど本スレは数学板なのね。
数学嫌いだわ…。Rまで嫌いになりそう

825 名前:132人目の素数さん mailto:sage [2013/04/13(土) 23:52:28.69 ]
むしろ数学界隈以外での R の用途が気になるな。なんでもできるっちゃそうなんだけど。

826 名前:132人目の素数さん mailto:sage [2013/04/14(日) 00:04:02.13 ]
>>825
学生時代数学にはつらい思いでがあって…

数学的能力が全くダメでもいろんな運の巡り合わせでRと向き合わなきゃいけなくなる人もいますんでご容赦を

827 名前:132人目の素数さん mailto:sage [2013/04/14(日) 05:11:16.14 ]
統計学もグラフィックスも数学じゃないよ。数学者はそんなもの相手にせんわ

828 名前:132人目の素数さん [2013/04/14(日) 07:37:32.81 ]
カール・フリードリヒ・ガウス...



829 名前:132人目の素数さん mailto:age [2013/04/14(日) 11:00:26.70 ]
          __ノ)-'´ ̄ ̄`ー- 、_
        , '´  _. -‐'''"二ニニ=-`ヽ、
      /   /:::::; -‐''"        `ーノ
     /   /:::::/           \
     /    /::::::/          | | |  |
     |   |:::::/ /     |  | | | |  |
      |   |::/ / / |  | ||  | | ,ハ .| ,ハ|
      |   |/ / / /| ,ハノ| /|ノレ,ニ|ル' 
     |   |  | / / レ',二、レ′ ,ィイ|゙/   私は只の数ヲタなんかとは付き合わないわ。
.     |   \ ∠イ  ,イイ|    ,`-' |      頭が良くて数学が出来てかっこいい人。それが必要条件よ。
     |     l^,人|  ` `-'     ゝ  |        さらに Ann.of Math に論文書けば十分条件にもなるわよ。
      |      ` -'\       ー'  人          一番嫌いなのは論文数を増やすためにくだらない論文を書いて
    |        /(l     __/  ヽ、           良い論文の出版を遅らせるお馬鹿な人。
     |       (:::::`‐-、__  |::::`、     ヒニニヽ、         あなたの論文が Ann of Math に accept される確率は?
    |      / `‐-、::::::::::`‐-、::::\   /,ニニ、\            それとも最近は Inv. Math. の方が上かしら?
   |      |::::::::::::::::::|` -、:::::::,ヘ ̄|'、  ヒニ二、 \
.   |      /::::::::::::::::::|::::::::\/:::O`、::\   | '、   \
   |      /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ  ヽ、  |
  |      |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、  /:\__/‐、
  |      |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄|
   |     /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_|
   |    |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
    |   /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/

830 名前:132人目の素数さん mailto:sage [2013/04/15(月) 07:52:40.36 ]
あるデータをplot(density(data))でグラフ化してます。
任意のX区間を指定して面積を求めようとしてるのですが、やり方を教えて下さい。

また、正規分布だと、dnorm(c(0,0.2,0.4,0.6,0.8,1))とやればX値に対応したY値がでるみたいですが、
任意のデータをプロットしている場合、X値からY値を取り出すにはどうやればいいのでしょうか?

831 名前:132人目の素数さん mailto:sage [2013/04/15(月) 08:54:48.29 ]
規制されて書き込めなかった。
普通はpolygon()を使う

>>820
> x <- rnorm(100);hist(x, prob=T, ann=F)
> a <- density(x)
> polygon(a$x,a$y, col="red")

832 名前:132人目の素数さん mailto:sage [2013/04/15(月) 09:12:17.89 ]
>>830
こちらも、>>831と同じ要領。
> data <- runif(100)
> plot(density(data))

> a <- density(data)
として、念のために幅を確認して、
> head(diff(a$x))
[1] 0.003151153 0.003151153 0.003151153 0.003151153 0.003151153 0.003151153
横が0.003151153と分かったので、矩形の面積を計算して合計。
> sum(diff(a$x)[1] * a$y )
[1] 1.000883

833 名前:132人目の素数さん mailto:sage [2013/04/15(月) 09:21:38.51 ]
>>825
統計的裏付けのない実感で申し訳ないが、
確かに、10年くらい前は数学・数理統計の関係者がほとんどを占めていたと思う。
しかし、ここ5年ほどで急速に、Rをプログラミング言語としてとらえるIT関係のユーザが増えたし、
心理学研究や医学研究にも地道に普及している。
Rユーザ会の出席者を見渡しても、数学関係者はマイノリティになりつつあると思う。

834 名前:132人目の素数さん mailto:sage [2013/04/15(月) 12:01:36.43 ]
こんな記事もある

「次に来る」プログラミング言語を占ってみる
ttp://www.atmarkit.co.jp/ait/spv/1201/17/news133_2.html

今後の動向を予感させるもう1つの動きは、プログラミング言語「R」が2012年1月のTIOBE PCIの19位に入ったことだ。
(中略)Rが話題になり始めた背景には、最近注目を浴びているある種のコンピュータ技術がある。
大量のデータを効率良く処理する技術だ。「ビッグデータ」という言葉を耳にした方は多いと思う。(中略)
RはGoogleが社内のデータ解析に活用しているなど、ビッグデータを扱う企業において重要なプログラミング言語となっており、今後広い範囲で話題になる可能性がある。

835 名前:132人目の素数さん mailto:sage [2013/04/15(月) 21:12:53.25 ]
後藤大地ってどんなやつかと思ったら小物もいいとこじゃん。聞く耳もつに
値しないね。

836 名前:132人目の素数さん mailto:sage [2013/04/15(月) 23:59:44.25 ]
ちゃんと取材にいってる人だと思うが、どちらかというとPCハードの記者さんだよね。

事実として、企業のコミットが増えてlinuxみたいにRもコミッター増えてんじゃないの?






[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

前100 次100 最新50 [ このスレをブックマーク! 携帯に送る ] 2chのread.cgiへ
[+板 最近立ったスレ&熱いスレ一覧 : +板 最近立ったスレ/記者別一覧](;´∀`)<381KB

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef