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
381 名前:132人目の素数さん [2012/02/02(木) 16:52:09.32 ] >>380 無理やりの処理だが、こんなんじゃ駄目? mypois <- function(n=50,sum=150,adjtimes=10){ myvec <- rep(0,times=n) mysum <- 0 for(i in 1:adjtimes){ mysum <- sum(myvec) myvec <- myvec + sign(sum-mysum)*rpois(n,abs(sum-mysum)/n) myvec[myvec < 0] <- 0 if( abs(sum-mysum) == 0 ){ break } } return(myvec) } > x <- mypois(50,150) > x [1] 3 2 3 2 5 3 4 5 4 0 1 1 5 4 6 3 5 5 0 4 4 3 4 5 4 7 1 4 1 2 1 2 3 2 1 3 7 2 [39] 0 2 2 3 0 1 3 6 3 4 2 3 > sum(x) [1] 150
382 名前:381 mailto:sage [2012/02/02(木) 17:00:53.33 ] 0から6になってなかったのと、計算の順番がおかしかったので修正した。 mypois <- function(n=50,sum=150,max=6,adjtimes=10){ myvec <- rep(0,times=n) mysum <- 0 for(i in 1:adjtimes){ myvec <- myvec + sign(sum-mysum)*rpois(n,abs(sum-mysum)/n) myvec[myvec < 0] <- 0 myvec[myvec > max] <- max mysum <- sum(myvec) if( abs(sum-mysum) == 0 ){ break } } return(myvec) } > x <- mypois(50,150,6) > x [1] 1 6 1 0 5 3 5 1 2 6 3 3 6 3 3 3 3 3 1 4 2 3 3 0 5 4 3 6 3 5 2 4 4 0 2 1 3 3 [39] 2 4 0 2 1 3 6 4 4 2 1 6 > sum(x) [1] 150
383 名前:381 [2012/02/02(木) 17:33:39.03 ] 青木先生のところの poissondist に嵌めたら、今一つなのは、 0〜6からはみ出たデータを丸めたからかな。上手い方法がないものかねぇ。 > x <- mypois(50,150,6) > x [1] 5 2 1 3 4 4 3 4 3 0 5 5 3 3 1 2 4 4 4 2 5 0 3 4 2 2 3 4 4 1 6 1 1 5 4 4 2 4 [39] 2 5 1 3 2 0 5 2 5 4 1 3 > poissondist(hist(x,plot=F)$counts) ポアソン分布への適合度の検定 data: hist(x, plot = F)$counts X-squared = 9.7821, df = 4, p-value = 0.04426 sample estimates: n lambda 50.00 2.06
384 名前:132人目の素数さん mailto:sage [2012/02/02(木) 18:40:11.25 ] >>381-383 おぉ、十分です。最大は「6ぐらい」なので7とかでもOKですから、 丸めはマイナスになるときだけ必要です。 for文のところはとても勉強になりました。 ありがとうございました。
385 名前:132人目の素数さん mailto:sage [2012/02/06(月) 23:36:20.10 ] こんばんはR初心者です。 一通り調べたのですが、分からない事があるので質問させてください。 ○行ごとのplotを行いたい場合、どのようなコマンドを打てば良いのでしょうか? 例えば、 y1 y2...y10 x1 x2 . . . x100 このような100行×10列のデータフレーム名Xの2列目と4列目をplotしたい場合、 plot(X$y2,X$y4)で1行ごとのplotは出来るのですが、 これをx1,x5,...x100と5行おきにplotしたいのです。 どうかご教授ください。
386 名前:132人目の素数さん mailto:sage [2012/02/07(火) 00:25:49.86 ] >>385 X1 <- X[(0:20)*5+2, ]
387 名前:132人目の素数さん mailto:sage [2012/02/07(火) 08:46:05.39 ] >>385 >このような100行×10列のデータフレーム名Xの2列目と4列目をplotしたい場合、 >plot(X$y2,X$y4)で1行ごとのplotは出来るのですが、 >これをx1,x5,...x100と5行おきにplotしたいのです。 plot(X[seq(1, 100, 5), c(2, 4)]) でどう? [行番号, 列番号] で要素を指定 seq(from, to, by)で1,6,11.16...,100をつくる
388 名前:132人目の素数さん mailto:sage [2012/02/07(火) 22:43:27.83 ] >>385 1,5,...,100 は、最初だけ4行の間隔ですね。 5,10,15,...,100 、あるいは、1,6,11,...,96 ですか? でしたら、私のお薦めも >>387 と同じ方法。
389 名前:132人目の素数さん mailto:sage [2012/02/08(水) 13:38:50.77 ] Emacs ESS でアンダースコアが" <- "になってしまうのが凄く迷惑なんだけど なんであれがデフォルトなんだろう
390 名前:132人目の素数さん mailto:sage [2012/02/08(水) 17:34:58.86 ] タイプ量減るしいいじゃん。
391 名前:132人目の素数さん mailto:sage [2012/02/08(水) 20:31:18.89 ] >>389 アンダースコア2度打ちで、アンダースコアを入力できるというのを 知った上で迷惑と行っているのだろうね。 そもそもRではアンダースコアは使用禁止だったので、 空いているアンダースコアにその機能が割り当てられたという歴史的経緯がある。 デフォルトの理由は、>>390 の言うとおり便利で楽だかから。 いやな人は設定を変更すればいいし。 逆にターミナルでRを使っていて、<-のつもりでアンダースコアを入力して イラッと来ることがある。
392 名前:132人目の素数さん mailto:sage [2012/02/08(水) 22:46:01.26 ] >>391 設定変えるためにググってたらrjpに飛んで、そこで初めて二度打ちやC-q _ でアンダースコアを入力できることを知った アンダースコアを打つためにわざわざテキストモードに変更してた よく調べもしないで愚痴言ってすいません
393 名前:132人目の素数さん mailto:sage [2012/02/13(月) 12:25:46.70 ] emacs essでRコマンダー使うとRコマンダーで外部データを読み込ませた瞬間固まるorz うまくいってるかたいますか?
394 名前:393 mailto:sage [2012/02/13(月) 14:44:06.36 ] すいません自己解決?しました。 emacsのR側で一行でも改行がコマンドを実行すればコマンダー側のフリーズが直りました。 フリーズ自体は何台かのPCで試したところ起こる場合と起こらない場合がありました。 環境に依存するようです。
395 名前:132人目の素数さん [2012/02/13(月) 22:00:37.21 ] >>393 ,394 emacs ess で Rコマンダーを使うメリットって、何? 利点がまったく分からないのだけど。
396 名前:393 mailto:sage [2012/02/14(火) 00:04:31.33 ] >>395 ご指摘のとおりで。 研究仲間がコマンダー使ってて、両方同時に使えると便利だから試しました。 まぁそうじゃなくても使えないよりは使えたほうが少しくらい便利かなと、、、
397 名前:132人目の素数さん mailto:sage [2012/02/15(水) 23:54:50.83 ] Rにはschemeなんかのfoldにあたるような高階関数はないんですか? せっかくベクトル演算言語なので数え上げも高階関数で出来ればと思ったんですが。
398 名前:132人目の素数さん mailto:sage [2012/02/16(木) 09:49:58.30 ] >>397 stat.ethz.ch/R-manual/R-patched/library/base/html/funprog.html プログラミング言語はよくわからないけど、これは参考になるだろうか?
399 名前:132人目の素数さん mailto:sage [2012/02/16(木) 09:57:29.11 ] >>398 たぶん Reduce()がfoldに相当するんじゃね
400 名前:397 mailto:sage [2012/02/16(木) 16:23:10.40 ] >>398 >>399 それです ありがとうございます
401 名前:398 mailto:sage [2012/02/16(木) 16:52:41.50 ] >>400 いまいちイメージがわかないんだけど、面白い使い方があったら、教えてくだされ
402 名前:132人目の素数さん [2012/02/17(金) 05:50:16.64 ] lmによる回帰分析で、切片をあらかじめ固定することはできますか? 切片を取り除くには「-1」を指定するということはわかりますが、ここを好きな値に固定したいのですが・・・
403 名前:132人目の素数さん mailto:sage [2012/02/17(金) 13:00:44.62 ] >>402 ttp://r.789695.n4.nabble.com/R-Forcing-the-intercept-td812422.html
404 名前:397 mailto:sage [2012/02/17(金) 19:52:56.72 ] >>401 おもしろいかどうかわかりませんが、 data <- as.integer(runif(10,1,21)) のようなデータがある時に、 predicate <- function(x) { x < 5 } が真になるような項目を数え上げたいとします。 for文を使うと ret <- 0 for(i in 1:length(data)){ if (predicate(data[i])) { ret <- ret + 1 } } のような形になりますが、高階関数を使うと、 length(Filter(predicate, data)) と書けるようになります。 これは関数プログラミングの技法で、もう少し詳しく学びたいようでしたら programming in haskellなんかが薄くていい本なんじゃないかと思います。
405 名前:132人目の素数さん mailto:sage [2012/02/17(金) 20:15:30.11 ] >>404 >>401 じゃなくて申し訳ないが、 > length(Filter(predicate, data)) [1] 4 > sum(data<5) [1] 4 この場合だと下の方が簡単。しかも、 > system.time(sapply(1:1000, function(x){length(Filter(predicate, data))})) ユーザ システム 経過 0.072 0.000 0.074 > system.time(sapply(1:1000, function(x){sum(data<5)})) ユーザ システム 経過 0.000 0.000 0.002 圧倒的に速い。
406 名前:405 mailto:sage [2012/02/17(金) 20:17:05.57 ] s/下の方が/後者の方が/
407 名前:132人目の素数さん mailto:sage [2012/02/17(金) 22:46:29.44 ] >>404 「おぉ高階関数はこれほどまでに便利なものかと実感できる事例」を期待している。
408 名前:132人目の素数さん mailto:sage [2012/02/18(土) 05:48:26.33 ] rubyで高階関数組み合わせてやってたけど Rに移ってtapplyとかmapplyが便利すぎて禿げた
409 名前:132人目の素数さん mailto:sage [2012/02/18(土) 06:05:41.05 ] >>404 Filter()の例は length(data[data < 5]) みたいなもんか
410 名前:397 mailto:sage [2012/02/18(土) 11:40:22.15 ] >>405 なるほど、そんな書き方が出来るんですか だとしたらそちらの方がよさそうです 速度はおそらく最適化の問題でしょう 下の方がRネイティブっぽいので十分な最適化が行われているんだとおもいます
411 名前:397 mailto:sage [2012/02/18(土) 12:01:56.66 ] >>407 事例はちょっと思いつかないんですが 高階関数はプログラムの中の構造を関数として再利用するためのものです その点では>>404 はあまりいい例ではなくて、 ret<-0 for(i in 1:length(data)){if(predicate(data[i])){ret<-ret+1}else{ret}} について高階関数を利用する例だとすると、 Reduce(function(ret, x){if(predicate(x)){ret+1}else{ret}},data,0) になります。これは外側のループの部分の構造をReduceという関数で再利用していると言う形になります でもどうやらこの手の単純な高階関数はベクトル計算という形でRに組み込まれているようなので あまり出番はなさそうですね
412 名前:132人目の素数さん [2012/02/19(日) 03:55:37.05 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
413 名前:132人目の素数さん mailto:sage [2012/02/19(日) 08:58:51.23 ] 教えてください。 Epiパッケージを使用して、ROC曲線下面積を求めたら0.5未満になりました。 独立変数は二分変数に関係しないと判断しその値(0.425)を記載することは正しいのでしょうか?
414 名前:132人目の素数さん mailto:sage [2012/02/19(日) 17:00:59.94 ] アウトカムの変数は予測因子の数字が大きくなると大きくなると陽性になりますか?もし数が大きい方がアウトカムが陰性ならば、アウトカムの1,0をひっくり返す必要があるかも?
415 名前:132人目の素数さん mailto:sage [2012/02/19(日) 23:14:01.45 ] >>414 ありがとうございます。 明日にでも、いれかえて検討してみます。
416 名前:132人目の素数さん [2012/02/29(水) 23:39:27.42 ] 応答曲面法で最適値を求めるっていうのは出来ないのか? パッケージrsmがそれっぽいんだが
417 名前:132人目の素数さん mailto:sage [2012/03/01(木) 23:08:27.13 ] >>416 なんていうか、もろに、rsm() がそうなんだが? 使い方は www.jstatsoft.org/v32/i07/paper をみてくれ
418 名前:132人目の素数さん [2012/03/03(土) 12:45:30.73 ] Rを2.14.2にしたいんですが、インストールして再起動すると確かに2.14.2にはなるんですが、保存していたオブジェクトが全て消えてしまいます オブジェクトを引き継ぐ方法はないんでしょうか? 今は2.13.1を使ってます
419 名前:132人目の素数さん mailto:sage [2012/03/03(土) 16:40:03.10 ] save(list = ls(all=TRUE), file = "old.RData") で明示的に保存ファイルを作って、 それを新バージョンで load(file = "old.RData") したらどうなる?
420 名前:132人目の素数さん mailto:sage [2012/03/04(日) 18:26:31.31 ] マルチコア対応しましたか〜??
421 名前:132人目の素数さん mailto:sage [2012/03/04(日) 19:53:40.43 ] parallelというパッケージができたね
422 名前:132人目の素数さん mailto:sage [2012/03/04(日) 20:21:12.91 ] >>420 ググると、gotoBLASをマルチコア向けにコンパイルしてRで利用する話が、 詳しく解説されている。 調子に乗ってcuBLASに差し替えられないかと考え、 プログラム板のcudaスレに行ったことがあるが、 あそこの住人はnvccを使わない人間をバカだと考えているみたいだ。
423 名前:420 mailto:sage [2012/03/04(日) 22:56:21.28 ] >>421 >>422 良い情報をありがとうございます。 調べてみます。 シングルスレッド性能でのCPU探しも終わりにしたいですw
424 名前:132人目の素数さん mailto:sage [2012/03/04(日) 23:06:08.66 ] 今なら分散処理じゃ無いかな。IT系のWebサイトじゃあ hadoop とか流行ってるし、 解きたい課題を Map-Reduce に落とし込めればRでの分散処理も可能かと。
425 名前:132人目の素数さん mailto:sage [2012/03/08(木) 19:59:09.03 ] ttp://cran.r-project.org/bin/linux/ubuntu/ preciseはまだまだ来ないね。
426 名前:132人目の素数さん [2012/03/11(日) 23:35:42.80 ] R初心者でpolr関数を使って解析しているのですが、 P値の出し方がわかりません。 お分かりになる方いらっしゃいましたら教えて頂けないでしょうか。 詳細な情報が必要であれば補足します。
427 名前:132人目の素数さん mailto:sage [2012/03/12(月) 16:16:16.41 ] stats.stackexchange.com/questions/7720/how-to-understand-output-from-rs-polr-function-ordered-logistic-regression の二番目の最後のところをみると、正確にはもとまらないから、計算されないようになっている、ということのように思いました。
428 名前:132人目の素数さん mailto:sage [2012/03/13(火) 00:08:39.50 ] www.r-project.org/ R 2.15.0 prerelease versions are currently available. Final release is scheduled for March 30.
429 名前:132人目の素数さん mailto:sage [2012/03/15(木) 07:30:22.52 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
430 名前:132人目の素数さん mailto:sage [2012/03/20(火) 20:20:19.86 ] Rを使ってる最中(emacs+ess) 作ったオブジェクトの数が多くなって、ls()で変数名を見ても、 思い出せない。 こういう状況に対処するにはどうしたら良い? 新しい環境を用意して、そちらに移す?(assignを使って、.GlobalEnvのものを消す?) それとも、特定のファイルに保存する?(複数個のオブジェクトを1ファイルに保存できる?) よい方法などがあれば、教えて。 あと、 emacs+essを使ってるけど、lispとかは分からんので、 emacsの機能での解決方法を示されても反応できないかも。 その時は、ごめんなさい。
431 名前:430 mailto:sage [2012/03/20(火) 20:22:30.66 ] 言葉不足だったので、訂正します。 >>assignを使って、.GlobalEnvのものを消す? assignを使って別の環境にデータを入れて、.GlobalEnvの方は消す?
432 名前:430 mailto:sage [2012/03/20(火) 20:28:40.90 ] >>419 でファイルに保存する方法があったね。見てなかった。 save(list=c("オブジェクト名1","オブジェクト名2"),file="a_path") とかでもいいのか。でも、なんか面倒だな。
433 名前:132人目の素数さん mailto:sage [2012/03/20(火) 21:37:25.25 ] >>430 変数名を見ても思い出せない命名法に問題がありそうだけど あと俺は使っていないけどcomment()を上手に使うとか。
434 名前:430 mailto:sage [2012/03/21(水) 00:50:48.76 ] >>433 こんな感じ? tmp <- sapply(ls(), function(x)comment(eval(parse(text=x)))) tmp[sapply(tmp,is.null)] <- NULL tmp これを関数にしておいて、呼び出せば良いかな? あとは、commentに自動的に参照or作成時間を入れておくとかかな。 どうやったらできるか分からないけど。
435 名前:132人目の素数さん mailto:sage [2012/03/22(木) 22:37:07.40 ] スクリプトをファイルに書き出して基本的にはバッチモードでRを動作させております。 その場合のエラーハンドリングってどうしてますか? stop関数で停止はするのですが、 どの関数のどの行数で停止したかなどがいまいちわかりづらいときがあります。 引数で与えるコメントを変えればいいのでしょうが、いちいち面倒なので。。。 例えばCだと printf("%s %s\n", __FILE__, __LINE__); なんかで対応できるんですが。。 いい方法ありますか?
436 名前:132人目の素数さん mailto:sage [2012/03/26(月) 13:16:29.99 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
437 名前:132人目の素数さん mailto:sage [2012/03/26(月) 17:52:13.93 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
438 名前:132人目の素数さん mailto:sage [2012/03/27(火) 00:36:36.55 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
439 名前:132人目の素数さん mailto:sage [2012/03/27(火) 12:48:45.92 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
440 名前:132人目の素数さん mailto:sage [2012/03/27(火) 20:00:22.88 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
441 名前:132人目の素数さん mailto:sage [2012/03/27(火) 23:36:28.28 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
442 名前:132人目の素数さん mailto:sage [2012/03/29(木) 01:36:02.65 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
443 名前:132人目の素数さん mailto:sage [2012/03/30(金) 22:40:29.56 ] R 2.15.0 is released www.r-bloggers.com/r-2-15-0-is-released/
444 名前:132人目の素数さん [2012/04/03(火) 21:15:56.45 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
445 名前:132人目の素数さん mailto:sage [2012/04/06(金) 18:58:58.10 ] チラ裏 Meadow 2.10(Emacs 21.4.1) + R 2.15.0 で、 ESS 5.14 → ESS 12.04 にしてみたら、 symbol's function definition is void と言われてしまい、 元に戻した今日の夕暮れ。
446 名前:132人目の素数さん mailto:sage [2012/04/10(火) 15:45:27.97 ] ttp://cran.r-project.org/bin/linux/ubuntu/ Preciseはまだ来ないね。4月26日以降になるのかな
447 名前:132人目の素数さん mailto:sage [2012/04/14(土) 11:43:41.70 ] R2.15.0 って、なんか目玉機能あります?
448 名前:132人目の素数さん [2012/04/20(金) 23:36:06.27 ] ソースコードが1つのファイルで2万行ぐらいになってしまったので 複数のファイルに分割したいのですがやり方がわかりません。 助けてください
449 名前:132人目の素数さん mailto:sage [2012/04/20(金) 23:38:20.91 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
450 名前:132人目の素数さん mailto:sage [2012/04/21(土) 07:52:14.18 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
451 名前:132人目の素数さん mailto:sage [2012/04/21(土) 09:49:45.56 ] >>448 2万行って何それ、すごい。 そんなにRの専門家なのに、そんなことを悩むとはw source()を使って階層的にしたら駄目なの?
452 名前:132人目の素数さん [2012/04/21(土) 16:42:07.10 ] >>448 むしろ、何の処理をさせたら、2万行になるのか知りたい。
453 名前:132人目の素数さん [2012/04/21(土) 17:49:35.34 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
454 名前:132人目の素数さん mailto:sage [2012/04/21(土) 18:59:03.29 ] >>448 縦書きか
455 名前:132人目の素数さん [2012/04/24(火) 07:20:26.82 ] SUBSTR<- function(X,S,E){ if(length(unique(nchar(X)))>1)stop("warning") if(unique(nchar(X))<max(S,E)|min(S,E)<=0)stop("warning") temp<-t(sapply(X,function(X){mapply(substr,X,start=S,stop=E)})) dimnames(temp)<-NULL temp }
456 名前:132人目の素数さん [2012/04/24(火) 07:25:35.41 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
457 名前:132人目の素数さん mailto:sage [2012/04/24(火) 10:40:39.43 ] >>455 突然どうしたの? オレならecho X | cut -cS-Eをsystem()でやっちゃうな。
458 名前:132人目の素数さん [2012/04/24(火) 12:16:39.44 ] barplotの質問です cran.md.tsukuba.ac.jp/bin/windows/base/ を使っています。 100個の文字列(labels)と数値(num)を棒グラフに書き出したいのですが barplot(num, names.arg=labels) いざ書き出すと文字列名がすべて表示されません 理想的にはcommons.wikimedia.org/wiki/File:English_letter_frequency_ (alphabetic).svg このように数値と文字列をすべて表示させたいのですが、どうすればよいでしょうか。 また、numでソートした棒グラフとlabelsでソートした棒グラフをそれぞれ作成するには どうすればいいでしょうか。 他のプログラミングでいうとろこの、 ハッシュ構造に文字列と数値を格納してキーと値でそれぞれソートするという手順がよく分かりません。
459 名前:132人目の素数さん mailto:sage [2012/04/24(火) 12:32:11.32 ] >>458 グラフのサイズを多きすくすれば細かい文字も入るよ。 > pdf(file="test.pdf",width=50,height=10) > barplot(round(runif(100),1),names.arg=sapply(1:100,function(x){sample(LETTERS,1)})) > dev.off() ソートは次の例を参考にしてくれ。 > (a <- data.frame(v=round(rnorm(5),1),n=sample(LETTERS,5))) v n 1 -0.1 F 2 -0.4 W 3 1.6 M 4 -0.2 R 5 -0.5 V > a[order(a$v),] v n 5 -0.5 V 2 -0.4 W 4 -0.2 R 1 -0.1 F 3 1.6 M > a[order(a$n),] v n 1 -0.1 F 3 1.6 M 4 -0.2 R 5 -0.5 V 2 -0.4 W
460 名前:459 mailto:sage [2012/04/24(火) 12:37:45.20 ] >>459 の barplot()は > barplot(round(runif(100),1),names.arg=sample(LETTERS,100,replace=TRUE)) に差し替え。最初はアルファベット2文字で書こうとしたのでsapply()が入った。
461 名前:132人目の素数さん mailto:sage [2012/04/24(火) 22:47:48.17 ] __ノ)-'´ ̄ ̄`ー- 、_ , '´ _. -‐'''"二ニニ=-`ヽ、 / /:::::; -‐''" `ーノ / /:::::/ \ / /::::::/ | | | | | |:::::/ / | | | | | | | |::/ / / | | || | | ,ハ .| ,ハ| | |/ / / /| ,ハノ| /|ノレ,ニ|ル' | | | / / レ',二、レ′ ,ィイ|゙/ 私は只の数ヲタなんかとは付き合わないわ。 . | \ ∠イ ,イイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。 | l^,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。 | ` -'\ ー' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて | /(l __/ ヽ、 良い論文の出版を遅らせるお馬鹿な人。 | (:::::`‐-、__ |::::`、 ヒニニヽ、 あなたの論文が Ann of Math に accept される確率は? | / `‐-、::::::::::`‐-、::::\ /,ニニ、\ それとも最近は Inv. Math. の方が上かしら? | |::::::::::::::::::|` -、:::::::,ヘ ̄|'、 ヒニ二、 \ . | /::::::::::::::::::|::::::::\/:::O`、::\ | '、 \ | /:::::::::::::::::::/:::::::::::::::::::::::::::::'、::::\ノ ヽ、 | | |:::::/:::::::::/:::::::::::::::::::::::::::::::::::'、',::::'、 /:\__/‐、 | |/:::::::::::/::::::::::::::::::::::::::::::::::O::| '、::| く::::::::::::: ̄| | /_..-'´ ̄`ー-、:::::::::::::::::::::::::::::::::::|/:/`‐'::\;;;;;;;_| | |/::::::::::::::::::::::\:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/ | /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::O::|::|::::::|:::::::::::::::/
462 名前:132人目の素数さん mailto:sage [2012/04/24(火) 23:01:21.57 ] >>459-460 ありがとうございます。ばっちりできました。
463 名前:132人目の素数さん mailto:sage [2012/04/24(火) 23:52:37.90 ] www.is.titech.ac.jp/~mase/R.html R の公式マニュアル5種類の和訳 このサイトが見れないです
464 名前:459 mailto:sage [2012/04/25(水) 00:40:38.22 ] >>463 10年以上前の情報だぞ。 書店で最新の本を買えよ
465 名前:132人目の素数さん mailto:sage [2012/04/25(水) 00:59:05.13 ] >>464 他にマニュアルの和訳したサイトってないんですか?
466 名前:132人目の素数さん mailto:sage [2012/04/25(水) 14:04:32.25 ] >>465 私の知る限りない。 現状では、マニュアルを読む気概がある人=英語でも大丈夫、という空気になっていて、 初心者や英語が苦手な人は、マニュアルを咀嚼した和書を購入すると思う。
467 名前:132人目の素数さん mailto:sage [2012/04/25(水) 21:29:19.12 ] >>466 ありがとうございます
468 名前:132人目の素数さん mailto:sage [2012/04/27(金) 03:23:30.97 ] sample(c("T", "F"), 1000, replace = TRUE) こうして作成したベルトルからTの個数を数えあげるには どうすればよいでしょうか for(i in 1:1000){}として一つ一つ値を見るしかないのでしょうか。
469 名前:132人目の素数さん mailto:sage [2012/04/27(金) 06:14:45.59 ] >>468 sum(sample(c("T", "F"), 1000, replace = TRUE))
470 名前:132人目の素数さん mailto:sage [2012/04/27(金) 06:17:19.41 ] c("T","F")だから、 sample(c("T", "F"), 1000, replace = TRUE)-> x sum(x=="T") でした。
471 名前:132人目の素数さん mailto:sage [2012/04/27(金) 11:05:43.29 ] >>468 長年のR利用経験から、おっさんが助言する。 "T"や"F"の文字列の使用は、トラブルの原因になる場合があるので、やめておけ。 なかなか気が付きにくいバクの原因になるから、 徹夜でトラブルシューティングするはめになるぞ。 "t"や"f"も使うな。 今回はsample(c("T", "F"), 1000, replace = TRUE)とする代わりに、 sample(c(T, F), 1000, replace = TRUE)とした方がよい。 "T"の代わりにTを、"F"の代わりにFを使うんだ。
472 名前:132人目の素数さん mailto:sage [2012/04/27(金) 14:26:57.77 ] x<-c(2, 3, 5, 4, 7, 9, 10, 2, 3, 4, 5) z<-c(3, 5, 9, 7, 13, 17, 19, 3, 5, 7, 9) y<-c(4, 6, 10, 8, 14, 18, 20, 4, 6, 8, 10) このx,y,zをもとに "要素1" xのベクトル "要素2" yのベクトル "要素3" zのベクトル このような表形式で結果をまとめるにはどうすればいいですか? ベクトルの操作だけでは全体の整形ができないようなのでわかりませんでした
473 名前:132人目の素数さん mailto:sage [2012/04/27(金) 14:31:54.94 ] >>472 「表形式で結果をまとめる」の意味が不明瞭。 ftable()案件なのか、rbind(x,y,z)やcbind(x,y,z)でよいのか、 それとも論文にするの表(TeX)にしたいのか、 分かるように伝えてくれ。
474 名前:473 mailto:sage [2012/04/27(金) 14:36:29.43 ] >>472 >"要素1" xのベクトル >"要素2" yのベクトル >"要素3" zのベクトル これって、例えば、次のような形なんだろうか。 > X <- rbind(x,y,z) > colnames(X) <- month.name[1:ncol(X)] > X January February March April May June July August September October November x 2 3 5 4 7 9 10 2 3 4 5 y 4 6 10 8 14 18 20 4 6 8 10 z 3 5 9 7 13 17 19 3 5 7 9
475 名前:132人目の素数さん mailto:sage [2012/04/27(金) 15:01:57.09 ] x<-c(2, 3, 5, 4, 7, 9, 10, 2, 3, 4, 5) z<-c(3, 5, 9, 7, 13, 17, 19, 3, 5, 7, 9) y<-c(4, 6, 10, 8, 14, 18, 20, 4, 6, 8, 10) f = c("要素1","要素2","要素3") cbind(f,data.frame(rbind(x,y,z))) こういうこと??
476 名前:132人目の素数さん mailto:sage [2012/04/27(金) 16:25:40.38 ] 説明不足ですみません tex形式ではなくプレーンテキストで結果を見ることができれば十分です >>473-475 のrbindやcbindで大丈夫そうです。ありがとうございます。 ただ欲を言えば[,1]のような見出し部分を取り除ければいいのですが、可能でしょうか。
477 名前:132人目の素数さん mailto:sage [2012/04/27(金) 20:26:45.89 ] ?write.table
478 名前:132人目の素数さん mailto:sage [2012/04/27(金) 21:07:02.32 ] >>476 真に必要がないのに、[, 1]などを削るのはいかがかと思うが、 次のようにすればどうだろう。 > x<-c(2, 3, 5, 4, 7, 9, 10, 2, 3, 4, 5) > z<-c(3, 5, 9, 7, 13, 17, 19, 3, 5, 7, 9) > y<-c(4, 6, 10, 8, 14, 18, 20, 4, 6, 8, 10) > hoge <- function(x){ + n <- nrow(x) + l <- rownames(x) + for(i in 1:n){ + cat(paste(l[i], "\t"), fill = FALSE) + cat(x[i, ], fill = TRUE) + } + } > hoge(rbind(x,y,z)) x 2 3 5 4 7 9 10 2 3 4 5 y 4 6 10 8 14 18 20 4 6 8 10 z 3 5 9 7 13 17 19 3 5 7 9
479 名前:132人目の素数さん mailto:sage [2012/04/27(金) 23:02:47.19 ] >>477-478 ありがとうございます。ここまで対応して下さり感謝します。
480 名前:132人目の素数さん mailto:sage [2012/04/28(土) 02:07:58.56 ] rを外部から操作する環境設定について質問します gnupack gnupack.sourceforge.jp/docs/latest/UsersGuide.html のemacsを使って居ます ESS (www.okada.jp.org/RWiki/?ESS )経由でRを利用しているのですが、 割合<-1/3 と書くと 「構文解析中に不正なマルチバイト文字列がありました」 とエラーになります RGuiのR Consoleからはそういうエラーはでません。 割合の"合"が原因のようで、ESSのやりとりではsjisとして文字コードが解釈されてしまうため、R側が適切に認識できないような気がします。 emacsからRをutf-8で操作できれば問題は解消すると思うのですが、どのように設定することになりますか。
481 名前:132人目の素数さん mailto:sage [2012/04/28(土) 05:34:39.27 ] presiceきたー ttp://cran.r-project.org/bin/linux/ubuntu/