【R言語】統計解析フリーソフトＲ 第４章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/01/09(日) 17:31:49統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
●関連サイト
The R Project
http://www.r-project.org/
RjpWiki
http://www.okada.jp.org/RWiki/
リンク集
http://www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8
0751１３２人目の素数さん
2013/02/01(金) 11:19:52.70よく分からない。date2numの引数はxだけなのに、なんで、dateとsumのデータフレームを入れるの？
それに、月を7文字目と8文字目ってしているけど、6文字目と7文字目じゃないの？
> date2num("2012/12/05 20:30:00")
以下にエラー if (.mon == 2) c <- 31 else if (.mon == 3) c <- 59 else if (.mon == (from #11) :
TRUE/FALSE が必要なところが欠損値です
[snip]
再現できるように質問しようよ。
もっと言えば、
> as.numeric(as.POSIXlt("2012/12/05 20:30:00"))
[1] 1354707000
とかを使うのは嫌いなの？
0752１３２人目の素数さん
2013/02/01(金) 12:21:01.54自分の力が及ばないせいでデータファイルを作るとき時刻データの前にどうしても消えない半角スペースがあって
それで6番目と7番目ではなく7番目8番目になっていますorz
それと
> as.numeric(as.POSIXlt("2012/12/05 20:30:00"))
[1] 1354707000
これをみて目からうろこでした、自分のやってたことの馬鹿さ加減にあきれました
御助言ありがとうございました
0753１３２人目の素数さん
2013/02/01(金) 13:26:36.91>時刻データの前にどうしても消えない半角スペース
関数を自作するレベルなのになぜw
例えば、
> (a <- c(" 2013-02-01 13:21:17", "2013-01-01 09:21:17", " 2013-03-01 11:21:17"))
[1] " 2013-02-01 13:21:17" "2013-01-01 09:21:17" " 2013-03-01 11:21:17"
こんな感じで、時刻前に空白があったりなかったりするなら、
> library(stringr)
> str_trim(a)
[1] "2013-02-01 13:21:17" "2013-01-01 09:21:17" "2013-03-01 11:21:17"
これでよいのでは？
オプションなしでは、前後の空白を取り除くから、前だけならそのように指定する。
0754１３２人目の素数さん
2013/02/01(金) 13:49:46.73なにからなにまでありがとうございます
0755１３２人目の素数さん
2013/02/01(金) 14:11:14.62>Rの方でtrimするって手もありましたね
別にRの外でも、どこででも空白を削除すれば良いと思う。
例えば下記のようなcsvがあったら、
$ cat tmp.csv
"A","B","C"
"BE DAC","Jan",0.35
"AC BDE","Feb",0.34
" EBACD","Mar",0.66
" DABEC","Apr",0.44
"BEA CD","May",0.84
"CDA EB","Jun",0.3
"EDCAB ","Jul",0.59
"CBEA D","Aug",0.53
"ECDA B","Sep",0.44
"ABDCE ","Oct",0.33
$ sed -e 's/¥" ¥([^ ,]*¥)/¥"¥1/g' tmp.csv
"A","B","C"
"BE DAC","Jan",0.35
"AC BDE","Feb",0.34
"EBACD","Mar",0.66
"DABEC","Apr",0.44
"BEA CD","May",0.84
"CDA EB","Jun",0.3
"EDCAB ","Jul",0.59
"CBEA D","Aug",0.53
"ECDA B","Sep",0.44
"ABDCE ","Oct",0.33
とすると3,4行目にあった文字列前の空白を消せる。
0756１３２人目の素数さん
2013/02/02(土) 10:13:07.27http://cse.naro.affrc.go.jp/takezawa/r-tips/r/43.html
上記サイトの複数の列に対してソートを行う場合は関数を使用して
同じように
a b c
11 1 2 11
12 1 1 12
9 2 4 9
10 2 3 10
7 3 6 7
8 3 5 8
5 5 8 5
6 5 7 6
4 7 9 4
1 9 12 1
2 9 11 2
3 9 10 3
というデータフレームをaを昇順に整列しaで同じものがあった場合はcの小さいものを上にするようにしたいのですが
In Ops.factor(mmm, each) : < 因子に対しては無意味です というエラーが出てうまくいきません
どこが問題でしょうか？
0757１３２人目の素数さん
2013/02/02(土) 10:18:28.45a b c
1 2 11
1 1 12
2 4 9
2 3 10
3 6 7
3 5 8
5 8 5
5 7 6
7 9 4
9 12 1
9 11 2
9 10 3
です
実際にはデータフレーム名test1で
sortlist <- order(test1$a, pmax(test1$a,test1$c))
test1[sortlist,]
以上に実行しました
0758756,757
2013/02/02(土) 10:45:26.080759１３２人目の素数さん
2013/02/02(土) 11:01:54.12複数カラムをキーにしたソートなんて考えたことがなかった。
いったいどんなときに必要になるの？
# 表計算ソフトのように「見せる」目的なら理解できるけど。
0760１３２人目の素数さん
2013/02/02(土) 11:09:52.570761１３２人目の素数さん
2013/02/02(土) 12:03:35.44a,b
3,5
2,3
2,2
2,1
1,3
1,1
a列は必ず降順になっているときに
b列が降順になっていないとき
この例では5行目1,3ですが
これを削除したいのですが 関数化して自動で判別して削除したいです
どのようにしたら良いでしょうか
0762１３２人目の素数さん
2013/02/02(土) 12:16:56.00doByパッケージのorderBy()が使いやすい。
0763１３２人目の素数さん
2013/02/02(土) 12:36:52.26pmax()というのを初めて知った。
bの列を一行下げて頭にInfをいれたものをcとしてデータフレームにいれる、bとcに対してpmax()して、結果をdとしてデータフレームにいれる。dがbと一致しているところを落とす。
でどうだろう？
0764761
2013/02/02(土) 13:18:11.28御助言ありがとうございます
しかし当方最近Rを使いだしたもので何をどうすればいいのか分かりません・・・
あつかましいのですが具体的にどんな感じに書けばいいでしょうか
0765１３２人目の素数さん
2013/02/02(土) 14:23:56.55dat = read.csv(text = "a,b
3,5
2,3
2,2
2,1
1,3
1,1")
dat
dat$c = c(Inf , head(dat$b, -1))
dat
with(dat, pmax(b, c))
dat$d = with(dat, pmax(b, c)) == dat$c
dat
dat[dat$d,]
を想定。関数化されてないけど。もっと大きい例で動くはチェックしていない。
0766１３２人目の素数さん
2013/02/02(土) 15:20:57.47> dat = read.csv(text = "a,b
+ 3,5
+ 2,3
+ 2,2
+ 2,1
+ 1,3
+ 1,1")
> dat[c(TRUE, diff(dat$b) < 0), ]
a b
1 3 5
2 2 3
3 2 2
4 2 1
6 1 1
0767１３２人目の素数さん
2013/02/02(土) 15:35:34.61bが例えば、10, 8, 6, 9, 8, 5とかになっていると、
>>766の方法を1回だけ適用したのでは通用しない。
> d2
a b
1 3 10
2 2 8
3 2 6
4 2 9
5 1 8
6 1 5
> d2[c(TRUE, diff(d2$b) < 0), ]
a b
1 3 10
2 2 8
3 2 6
5 1 8
6 1 5
そのため、十分な回数を繰り返し適用することになる。
> sum(diff(d2$b) > 0)
[1] 1
これが0になるのが、終了条件。
0768１３２人目の素数さん
2013/02/02(土) 19:57:32.38を使えばいいか？
0769１３２人目の素数さん
2013/02/03(日) 12:03:27.26dat = read.csv(text = "a,b
3,5
2,3
2,2
2,1
1,9
1,8
1,1")
dat
dat$cummin <- cummin(dat$b)
dat$correct <- with(dat, b == cummin)
dat
dat[dat$correct == TRUE,]
0770１３２人目の素数さん
2013/02/03(日) 13:34:45.070771１３２人目の素数さん
2013/02/03(日) 22:46:23.58膨大なデータを統計処理して活用するために
RやJavaのhadoopが使える人材が求められてるって書いてあったな
0772１３２人目の素数さん
2013/02/05(火) 08:32:45.91Baseパッケージで何か良い方法ありませんでしょうか？
今は↓のように書いていますが、ablineで線を引くとはみ出してしまう等手間がかかるので、
簡単にかける方法をご存知の方おりましたらお願いします
plot(x=0,y=0,xlim=c(0,6),ylim=c(0,3),lwd=2,ann=F,cex.axis=1,las=1,type="n",axes=F)
axis(side=1,seq(0,6),lab=F,pos=0,tck=0); axis(side=2,seq(0,3),lab=F,tck=0,pos=0)
axis(side=3,at=0:6,lab=F,tck=0,pos=3); axis(side=4,at=0:3,tck=0,lab=F,pos=6)
axis(side=1,at=0,tck=0,pos=0)
>>735-37
レス忘れていました、申し訳ありません
Rconsoleが2つ存在し妙な読み込みがされていたのが原因でした
アドバイスありがとうございました
0773１３２人目の素数さん
2013/02/05(火) 09:04:30.55> plot.new()
> box()
とすると、左下に原点が来るけど。
> points(0, 0)
0774１３２人目の素数さん
2013/02/05(火) 09:17:22.01R的なアプローチとして、複雑な手続きを関数化して簡単にする方法がある。
f772 <- function(){
plot(x=0,y=0,xlim=c(0,6),ylim=c(0,3),lwd=2,ann=F,cex.axis=1,las=1,type="n",axes=F)
axis(side=1,seq(0,6),lab=F,pos=0,tck=0); axis(side=2,seq(0,3),lab=F,tck=0,pos=0)
axis(side=3,at=0:6,lab=F,tck=0,pos=3); axis(side=4,at=0:3,tck=0,lab=F,pos=6)
axis(side=1,at=0,tck=0,pos=0)
}
とスクリプトの最初に書いておけば、
> f772()
以降は、これだけでよい。Baseパッケージのコマンドを探すよりも簡単。
0775１３２人目の素数さん
2013/02/06(水) 03:22:55.78plot(1,1, xaxs = "i", yaxs = "i", xlim = c(0,1), ylim = c(0,1), tck = 0)
?parでbaseの細かいオプションの説明がみられる。
これも役立つかも
www.stat.auckland.ac.nz/~paul/RG2e/chapter3.html
0776１３２人目の素数さん
2013/02/06(水) 06:43:29.05, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/ 私は只の数ヲタなんかとは付き合わないわ。
. | ＼ ∠イ ,ｲイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。
| l＾,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。
| ` -'＼ ｰ' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて
| /(l ＿＿／ ヽ、 良い論文の出版を遅らせるお馬鹿な人。
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、 あなたの論文が Ann of Math に accept される確率は?
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼ それとも最近は Inv. Math. の方が上かしら？
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0777１３２人目の素数さん
2013/02/06(水) 17:31:58.49レスありがとうございます
非常に参考になりました
0778１３２人目の素数さん
2013/02/15(金) 09:19:22.40軸(目盛の位置)がずれてしまうのですが、
揃える方法ってあるのでしょうか？
barplot(1:10, xlim=c(0,10), ylim=c(0,10))
par(new=T)
plot(1:10, col=2, xlim=c(0,10), ylim=c(0,10))
0779１３２人目の素数さん
2013/02/15(金) 10:52:07.52よく分からないが、
> a <- barplot(1:10); points(a, 1:10, col = 2)
こうしたいの？
0780１３２人目の素数さん
2013/02/15(金) 11:18:55.04縦軸の0の位置が、barplot と plot で違うんで揃えたいんです。
0781１３２人目の素数さん
2013/02/15(金) 11:34:08.72いや、そもそも揃っているけど。
> a <- barplot(1:10); points(a, rep(0, 10), col = 2)
0782781
2013/02/15(金) 11:39:38.94あっ、分かった。plot()を無理矢理使っているんだ。
addオプションをサポートしているクラス以外では、plot()を重ねるのは禁止。
想定外の使用形態。
描画を重ねるとき、普通はpoints()やlines()やaxis()を使う。
0783１３２人目の素数さん
2013/02/15(金) 11:44:27.170784781
2013/02/15(金) 11:53:12.37例えば、棒グラフと折れ線グラフを重ねたかったら、次のようにする。
x <- 1:10
n <- length(x)
y <- jitter(rep(10, n), 10)
op <- par(mai = c(1.02, 0.82, 0.82, 0.82))
a <- barplot(x, ylim = c(0, 13))
points(a, y, type = "l", lty = 2)
axis(4, label = (0:6) * 2 / 100, at = (0:6) * 2, las = 2)
par(op)
0785１３２人目の素数さん
2013/03/09(土) 12:17:39.65plyr : 中規模以上のデータで遅すぎ（idataframe使ってもそこまで変わらず）
aggregate : 普通
sqldf : まとも、でもRっぽくない。
他にお勧めがあったら教えて
0786１３２人目の素数さん
2013/03/13(水) 01:03:30.88rawデータをcharacterにしたいです。
chromeの localstorageのデータを整理しようとしています。
ドメインごとのsqlite形式？ファイルがあり、ItemTableというテーブル名で、
key,valueの項目でデータが入っています。
valueのデータ形式がblob型です。
RSQLiteで、
ret <- dbGetGuery(con, "select value ItemTable")
でデータフレームでデータが返って来ました。
class(ret[1,1][[1]])
>> raw
それで、これから先、
ret[1,1][[1]] を characterにしたのですが、
どうすればいいのでしょうか？
0787１３２人目の素数さん
2013/03/13(水) 12:10:10.92内部でSQLiteを使っているソフトが最近多いなと思っていたが、
Chromeもそうなんだ。
再現して検証できない質問なので、適当な回答になるけど、
as.character(ret[1,1][[1]])でどのようなエラーになりますか？
0788786
2013/03/13(水) 12:57:52.15ありがとうございます。
chrome側のconsoleで、 localstorage.hoge = "abc" として、
sqlteデータを入れた時に、
R側でそのblob型のデータ(abcというデータ, 他のsqlite viewerで確認済み) を読むと、
> ret[3,2][[1]]
[1] 32 00 30 00 31 00 33 00 2d 00 30 00 33 00 2d 00 31 00 32 00 54 00 31 00 34
[26] 00 3a 00 33 00 35 00 3a 00 35 00 30 00 2e 00 31 00 34 00 33 00 5a 00
となり、
as.character(ret[3,2][[1]])
[1] "32" "00" "30" "00" "31" "00" "33" "00" "2d" "00" "30" "00" "33" "00" "2d"
[16] "00" "31" "00" "32" "00" "54" "00" "31" "00" "34" "00" "3a" "00" "33" "00"
[31] "35" "00" "3a" "00" "35" "00" "30" "00" "2e" "00" "31" "00" "34" "00" "33"
[46] "00" "5a" "00"
となります。
こういうのを16進数のbyteデータというのですよね？
これをcharacterにencode?してやりたいです。
0789１３２人目の素数さん
2013/03/13(水) 13:54:01.89元データのエンコードが分からないと厳しいと思います。日本語でしょうか。
英語(ASCII)なら、
rawToChar(ret[3,2][[1]])とするとうまく行くかも知れません。
0790789
2013/03/13(水) 15:25:46.98ちゃんと読んでいませんでした。申し訳ない。
これは「abc」という文字列なんですね。
ちなみに"abc"の8進数表記は次の通り
$ echo abc | od -to1c
0000000 141 142 143 012
a b c \n
10進数表記
$ echo abc | od -td1c
0000000 97 98 99 10
a b c \n
16進数表記
$ echo abc | od -tx1c
0000000 61 62 63 0a
a b c \n
>>788 を見る限り、abcじゃなさそうに思える。
16進数と仮定すると、
2 NUL 0 NUL 1 NUL 3 NUL - NUL 3 NUL -
NUL 1 NUL 2 NUL T NUL 1 NUL 4 NUL : NUL 3 NUL
5 NUL : NUL 5 NUL 3 NUL . NUL 1 NUL 4 NUL 3
NUL Z NUL
となるが、NULを飛ばすと、2013-3-12T14:35:53.143Zで、日付時間っぽいな。
これ以上は分からないので、識者の方にパス。
0791１３２人目の素数さん
2013/03/13(水) 15:28:55.56ありがとうございます。
かなり前進しました。
さきほど、データの提示は、間違え(一つずれたものでした)
"abc"というデータを入れたときは、
ret[4,2][[1]]
[1] 61 00 62 00 63 00
rawToChar(ret[4,2][[1]],multiple=T)
[1] "a" "" "b" "" "c" ""
なぜか、空文字が入りますがOKです。
また、日本語の場合ですが、
私のwindows環境の文字コードは、cp932だと思います。
>Sys.getlocale()
[1] "LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.932;......."
それで、データとして "あい" と localstorageに入れて、
RSQLiteでとりだしたところ、
ret[5,2][[1]]
[1] 42 30 44 30
となり、
これは、文字コード表を見ると、
unicode (あい) => UTF16形式 3042 3044 と、
２バイトの間隔で前後入れ替えしたものに見えます。
mapply(function(a,b)paste(b,a,sep=""), x[seq(1,3,2)], x[seq(2,4,2)])
[1] "3042" "3044"
これを、"あい" という文字で、
できれば、UTF8形式でファイルにおとしたいのですが、、
0792１３２人目の素数さん
2013/03/13(水) 15:31:21.55返事が入れ違いになりました。
"abc"ではなく、日付でした。
間違えてました。
0793１３２人目の素数さん
2013/03/13(水) 21:35:55.82下記のようなCSVファイル(,区切り)が存在します。ID、NAME、MONTHは各情報の最初の行にしかデータとして入っておりません。
---------- データ部 ----------
ID, NAME, MONTH, PRICE, NUMBER
AAA1, Apple, 1, 100, 10
, , , 98, 12
, , , 102, 8
AAA2, Orange, 2, 120, 3
, , , 130, 4
, , , 112, 8
------------------------------
これをデータフレームとして読み込みます。
ID NAME MONTH PRICE NUMBER
AAA1 Apple 1 100 10
"" "" "" 98 12
"" "" "" 102 8
AAA2 Orange 2 120 3
"" "" "" 130 4
"" "" "" 112 8
.......
となります。
ここでID NAME MONTHを上の行から補間したいのですが、どのようにすればよいでしょうか?
ID NAME MONTH PRICE NUMBER
AAA1 Apple 1 100 10
AAA1 Apple 1 98 12
AAA1 Apple 1 102 8
AAA2 Orange 2 120 3
AAA2 Orange 2 130 4
........
としたいです。補間さえできれば形式はデータフレームでなくても結構です。よろしくお願いいたします。
0794786
2013/03/13(水) 22:48:40.11http://stackoverflow.com/questions/2776135/last-observation-carried-forward-in-a-data-frame
を参考にするといいよ。
""をNAに変えておかないといけないので、
df[ df == ""] <- NA
としておいて、
library(zoo)
na.locf(df)
まだ、rawデータからcharacterに変えられてない、、、
0795１３２人目の素数さん
2013/03/14(木) 09:59:33.96単純にコピーするだけでは駄目なの？
> sample <- data.frame(matrix(c("AAA1","Apple",1,100,10,
+ "","","",98,12,
+ "","","",102,8,
+ "AAA2","Orange",2,120,3,
+ "","","",130,4,
+ "","","",112,8), ncol = 5, byrow = TRUE))
> names(sample) <- c("ID", "NAME", "MONTH", "PRICE", "NUMBER")
> sample
ID NAME MONTH PRICE NUMBER
1 AAA1 Apple 1 100 10
2 98 12
3 102 8
4 AAA2 Orange 2 120 3
5 130 4
6 112 8
> for(i in 1:(nrow(sample) - 1)){
+ if(sample[i+1, 1] == ""){
+ sample[i+1, 1:3] <- sample[i, 1:3]
+ }
+ }
> sample
ID NAME MONTH PRICE NUMBER
1 AAA1 Apple 1 100 10
2 AAA1 Apple 1 98 12
3 AAA1 Apple 1 102 8
4 AAA2 Orange 2 120 3
5 AAA2 Orange 2 130 4
6 AAA2 Orange 2 112 8
0796793
2013/03/14(木) 22:00:08.24ありがとうございます。
無事うまくいけました。
質問に答えることができればよかったのですが、
残念ながら私にはわかりません。申し訳ありません。
>>795
実際のデータは100万行ほどありますので、
ループを回すと試してはいませんが時間がかかるかと。。。。
0797794
2013/03/15(金) 15:10:27.37いえいえ。
791の疑問は、説明の仕方が悪い気がする。
単純な文字コードの問題なのだろうけど、そこのところがよくわかてっない。
他の言語でやりますのでOKです。
0798１３２人目の素数さん
2013/04/02(火) 00:12:39.59以下の処理を実行するとエラーになります。
library( DMwR )
work00 <- data.frame( kbn=c(rep("A",12),rep("B",3)) , tuikavar=c(rnorm(12,20),rnorm(3,2)) )
SMOTE( kbn ~ . , data=work00 , perc.over=200 , perc.under=100 )
>以下にエラー `colnames<-`(`*tmp*`, value = c("kbn", "tuikavar")) :
>'names' 属性 [2] はベクトル [1] の長さと同じでなければなりません
ドキュメントを読みつくしたつもりなのですが、バグとしか思えません。
詳しい方、このバグは既知のものなのでしょうか？それとも回避方法があるのでしょうか？
0799１３２人目の素数さん
2013/04/02(火) 09:47:57.43その作者のTorgoさんに聞く方がここで聞くよりも早いと思うぞ
0800狢 ◆yEy4lYsULH68
2013/04/02(火) 16:08:36.00> 1 ：西独逸φ ★：2007/08/05(日) 05:47:55 ID:???0
>徳島県警阿南署などは５日未明、東京都足立区千住寿町、
>筑波大学准教授、増田哲也容疑者（５０）を
>県迷惑行為防止条例違反（痴漢行為）容疑で逮捕した。
>
>調べでは、増田容疑者は、４日午後４時２０分ごろから約５０分にわたり、
>ＪＲ牟岐線の列車内で、県内の
>専門学校生の女性（２１）の胸や太ももなどを触った疑い。調べに対し、
>「夏休み期間に、講演活動を兼ね
>て旅行していた。好みの女性だったのでムラムラした」と話しているという。
>
0801１３２人目の素数さん
2013/04/03(水) 23:34:54.430802１３２人目の素数さん
2013/04/04(木) 12:51:51.74R (正確にはR.app GUI)が起動しなくなったw
コンソールの方は使えたので、すぐに対処できたけど、R初心者は要注意だな。
バージョン上げたら、ライブラリのコンパイルをやり直さなくてはいけないので面倒だよな。
configureにオプションを渡さないといけないのもあるし。
> nrow(library()$results)
[1] 267
みんなは、いくつぐらいなの？
0803１３２人目の素数さん
2013/04/04(木) 13:07:53.300804１３２人目の素数さん
2013/04/04(木) 13:11:23.493.1とか3.0.2が出るまでとりあえず様子見だな
0805１３２人目の素数さん
2013/04/04(木) 22:31:14.28[1] 2961
Rcmdrのせいなのか随分と多い。
自分の環境はR version 2.13.0 (2011-04-13)
0806１３２人目の素数さん
2013/04/05(金) 01:02:41.34font = TT MS Gothic
points = 20
background = black
normaltext = white
usertext = LightGreen
highlight = Red
style = normal # Style can be normal, bold, italic
何が悪いんだろう
0807１３２人目の素数さん
2013/04/05(金) 07:49:15.90デフォルトエディタらしいといえばそうなのだが、Rエディタはあまりに低機能過ぎる
モニタが大きい人ならRstudio、ノーパソ派ならNotepad++ & NppToR、あたりが個人的おすすめ
0808１３２人目の素数さん
2013/04/05(金) 09:33:00.80は？ターミナルエミュレータの設定にRは関係ないだろ、と思ったら、
> font = TT MS Gothic
Windowsな方でしたか。
私も回答になっていなくて申し訳ないけど、
背景やフォントカラー、構文色強調を自由に設定できるEmacs+ESSがお薦め。
Emacsは初期学習コストがかかるが、超絶便利。
初めてESSを使ったときは、あまりの便利さに気を失いそうになった。
0809１３２人目の素数さん
2013/04/05(金) 15:31:00.183.0.0で使えているみたいだけど。
RStudio Ver.0.97.336
0810１３２人目の素数さん
2013/04/06(土) 08:59:21.34インスコした後まだ立ち上げる機会がない。
0811１３２人目の素数さん
2013/04/06(土) 13:04:58.94絶対パス: c:\program files\r\r-3.0.0\bin\i386\r.dll
____________________________
____________________________
コンピュータの評価日 2013/04/05 12:12:11
前回の使用 利用不能
起動項目 いいえ
起動済み はい
____________________________
____________________________
少数のユーザー
Norton コミュニティの 50 人未満のユーザーがこのファイルを使いました。
____________________________
ごく新しい
このファイルの更新日はここ 1 週間以内です。
____________________________
高
これは危険度が高いファイルです。
____________________________
脅威の詳細
SONAR 保護はコンピュータ上で疑わしいプログラムの活動を監視します。
____________________________
元ファイル:
r.dll
0812１３２人目の素数さん
2013/04/07(日) 12:09:36.610813１３２人目の素数さん
2013/04/07(日) 17:13:19.15「 (略) 作業効率が上がる（ように錯覚でき）ます． (RjpWikiより)」
0814１３２人目の素数さん
2013/04/08(月) 11:00:43.17逆。試しに使ってみるとか、RCommanderしか使わないとか、
年に数回だけ使うとかそんなレベルなら、
コンソール/RCommanderでもよいが、
ある程度使うなら、試行錯誤するためにも、エディタは必須。
常用しているテキストエディがなければ、微妙だけど。
0815１３２人目の素数さん
2013/04/08(月) 12:42:11.560816１３２人目の素数さん
2013/04/09(火) 19:51:21.47過去にsave.image()しまくった資産があるのでなければ、
基本的に非互換なんて気にしなくてもよいと思うよ。
でも、Ubuntuのリポジトリは2.15.3のままだな。
単に遅れているだけか、何か理由があるのかな。
Debianの方は3.0.0が来た。
0817１３２人目の素数さん
2013/04/10(水) 01:47:48.52因みに2.15だとgplotsが無いみたいでvenn()は使えないみたいでした。
0818１３２人目の素数さん
2013/04/10(水) 08:58:38.30> RSiteSearch("Venn")
でヒットした関数やパッケージは試した上での質問かどうかで回答は変わるけど、
VennDiagramパッケージのような専用パッケージはどう？
0819１３２人目の素数さん
2013/04/10(水) 10:52:02.97>>817-818同じくVennDiagramに一票
0820１３２人目の素数さん
2013/04/12(金) 10:13:54.71x <- rnorm(100);hist(x, prob=T, ann=F) ;par(new=T) ;plot(density(x),col="red")
こんな感じの正規分布っぽいグラフでXの範囲指定して、
Y軸0から確率密度の赤線までの面積を表したいんだけど、
ここを色づけする方法が判らない、
例えばXが1~2の間を塗るとき具体的にどんなコードになるか教えて下さい。
0821１３２人目の素数さん
2013/04/12(金) 15:37:00.070822１３２人目の素数さん
2013/04/13(土) 00:06:14.95dat <- 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)))
0823１３２人目の素数さん
2013/04/13(土) 10:48:56.27ありがとう
0824１３２人目の素数さん
2013/04/13(土) 23:48:23.04Ｒ自体統計プログラムなのかグラフィックソフトウェアなのかよくわからんけど本スレは数学板なのね。
数学嫌いだわ…。Ｒまで嫌いになりそう
0825１３２人目の素数さん
2013/04/13(土) 23:52:28.690826１３２人目の素数さん
2013/04/14(日) 00:04:02.13学生時代数学にはつらい思いでがあって…
数学的能力が全くダメでもいろんな運の巡り合わせでRと向き合わなきゃいけなくなる人もいますんでご容赦を
0827１３２人目の素数さん
2013/04/14(日) 05:11:16.140828１３２人目の素数さん
2013/04/14(日) 07:37:32.810829あぼーん
NGNG0830１３２人目の素数さん
2013/04/15(月) 07:52:40.36任意のX区間を指定して面積を求めようとしてるのですが、やり方を教えて下さい。
また、正規分布だと、dnorm(c(0,0.2,0.4,0.6,0.8,1))とやればX値に対応したY値がでるみたいですが、
任意のデータをプロットしている場合、X値からY値を取り出すにはどうやればいいのでしょうか？
0831１３２人目の素数さん
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")
0832１３２人目の素数さん
2013/04/15(月) 09:12:17.89こちらも、>>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
0833１３２人目の素数さん
2013/04/15(月) 09:21:38.51統計的裏付けのない実感で申し訳ないが、
確かに、10年くらい前は数学・数理統計の関係者がほとんどを占めていたと思う。
しかし、ここ5年ほどで急速に、Rをプログラミング言語としてとらえるIT関係のユーザが増えたし、
心理学研究や医学研究にも地道に普及している。
Rユーザ会の出席者を見渡しても、数学関係者はマイノリティになりつつあると思う。
0834１３２人目の素数さん
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が社内のデータ解析に活用しているなど、ビッグデータを扱う企業において重要なプログラミング言語となっており、今後広い範囲で話題になる可能性がある。
0835１３２人目の素数さん
2013/04/15(月) 21:12:53.25値しないね。
0836１３２人目の素数さん
2013/04/15(月) 23:59:44.25事実として、企業のコミットが増えてlinuxみたいにRもコミッター増えてんじゃないの？
0837１３２人目の素数さん
2013/04/22(月) 15:49:58.990838１３２人目の素数さん
2013/04/22(月) 16:09:28.75ttps://stat.ethz.ch/pipermail/r-announce/2013/000561.html
0839１３２人目の素数さん
2013/04/22(月) 21:01:07.44競馬とオートレースに使ってる。
仕事では使わないなぁ。部下は使ってるみたいでレポート見るだけだね。
0840１３２人目の素数さん
2013/04/22(月) 22:56:40.88もうかりまっか？
0841１３２人目の素数さん
2013/04/22(月) 23:24:49.44オートレースは確実に儲かるよ。ルールやら競技特性やら勉強してパラメータ決めるまで大変だったけど。
競馬はまだキツい。
0842１３２人目の素数さん
2013/04/22(月) 23:26:28.830843１３２人目の素数さん
2013/04/22(月) 23:26:32.86＞競馬で30億円配当に巨額の追徴課税の件
0844１３２人目の素数さん
2013/04/22(月) 23:34:16.04技能を分解して数値化してます。
>>843
派手にやらなきゃ大丈夫。あれがバレたきっかけは株らしいし。
0845１３２人目の素数さん
2013/04/22(月) 23:37:27.43あと技能を分解ってすげーな 何を基準に分けてんのか興味あるｗ
0846１３２人目の素数さん
2013/04/22(月) 23:43:03.50あまり掘り下げるとスレ違いになるのですいません。統計の理屈はギャンブル特性の解析よりある意味簡単ですね。
計算の最後にはsvm使ってますよ。
0847１３２人目の素数さん
2013/04/23(火) 10:17:29.620848１３２人目の素数さん
2013/04/25(木) 17:18:24.32いぢめられるよ？
0849１３２人目の素数さん
2013/04/25(木) 17:33:46.300850１３２人目の素数さん
2013/04/25(木) 19:41:57.32自腹で買う奴いるのか？
会社や学校とか所属組織のカネで買うもんじゃないの？
■ このスレッドは過去ログ倉庫に格納されています