統計解析フリーソフト　Ｒ　【第３章】

1 ：１３２人目の素数さん：2008/10/16(木) 16:33:16

R は統計解析とグラフィックスのためのフリーソフトです。
統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。

2 ：１３２人目の素数さん：2008/10/16(木) 16:34:58

★ＲプロジェクトのＨＰ

「The R Project」
http://www.r-project.org/
「CRAN」（日本ミラーサイト・Ｒ本体はここから手に入ります）
http://cran.md.tsukuba.ac.jp/

◆日本語ドキュメント

「R の日本語文章 (pdf 版)」
http://buran.u-gakugei.ac.jp/~mori/LEARN/R/
「Ｒ−Ｔｉｐｓ　ＰＤＦ版」
http://cse.naro.affrc.go.jp/takezawa/r-tips.pdf

3 ：１３２人目の素数さん：2008/10/16(木) 18:56:36

　　　　　　 |: : : : : : : :,イ　／'7 /: : : : : : :／　＼　＼
　　　　　　 |: : : : : :／: {　ゝ::/　7: : : : ／　　　　＼　ヽ
　　　　　 　|:__; ;- ´: : : ｀ ´: :｀:´○::／　 ',　　　　　ヽ　,
　　　　　　 'ヽ: : : : : : : : : : : : : :.iイ　　 i　|ヽ __|　　i i, ',`､　　　３げっと
　　　　　　　　｀ - : :_:_:_:_: : イ」_|_|,　　 ﾄ,イ´,　|| |　| .|　,ヽ
　　　　　　　　　　　 /　|'／´|V,,ゝﾄ, 　ｲ ﾚ 　}ﾉ |.|　i|i　　, ＼
　　　　　　 　 　 　 / 　 |　　{ O::::`, ヽﾉ　 ／⌒'7 ﾉｲ) 　 , 　＼
　　　　　　　（⌒ヽ/　/　＼ |ヽヾﾆｿ　,　 　　 　/／ ﾚ　　 i＼　＼ _,,,--,,,_
　　　　　　 　 ｀{ニ;(´＼　| ヽﾄ.　　　 ＿ノ　　 ノ|　 / 　 |　|　｀ ー|　　　　＼
　　　　　　　／{＿ﾌ7　 }　|　 | ＞ ,,＿＿,, ＜　/　/　　 |　|　　　 ｀´⌒ヽ　 ヽ
　　　　 　 ／´ ヽ フ　 ノ　 |　|'＼_:_::ゝ／´: : :V　/　|　 ||　|　　　　　　　} 　　}
　　　　　　　　　 7i二二it ﾅ＼|　［:ニ［|＼-'-/／　 /　/'| ,|　　 く´⌒｀ｰ　　ノ
　　　　　　　　　 |　|　　| ロ　 ヽ　/:/ |;ト′ |￣''ヽ/__./　|/　 r ヽヽ_, ‐-‐‐´
　　　　　 　 　　 |　|　　 ー-ｯ--く∠__|:》___　ー´　〈〈｀i　　...:::ヽノ
　　　 　 　 　 　 |　i|　　　／　　　　　 ´　　［二＞,,,○|...::::''''
　　　　　　　　　 ／: :ー ´:ヽ　　　　　　　　{: :［ﾆ／ヽ:::'''、
　　　　　　　　／: : : : : : ／￣｀ー　＿＿_ノ／　　 ／: : :＼

4 ：１３２人目の素数さん：2008/10/16(木) 23:17:21

ちんこ乙

5 ：１３２人目の素数さん：2008/10/17(金) 22:13:51

すみません質問です。最近Ｒを使い始めました。
今、
unko0.dat
unko1.dat
unko2.dat
...という風な名前のファイルが大量にありそれらに似たような処理をＲでしたいとします。
そのため、以下のような関数を作りたいです。

operation <- function(n){

dataframe <- read.table("unko$n.dat")

#データフレームにやりたい処理を書く

}

ここで"unko$n.dat"のように文字列の中で変数を展開させたいのですが
そのようなことは可能でしょうか？よろしくお願いいたします。

6 ：１３２人目の素数さん：2008/10/18(土) 02:53:02

事故解決しました。

str <- paste("unko",n,".dat",sep="")

でできました。ｽﾚ汚しｽﾏｿ。

7 ：１３２人目の素数さん：2008/10/19(日) 01:25:33

最近Ｒを使い始めたのですが質問があります。
パッケージの中身のソースはどうやってみるのでしょうか。

8 ：１３２人目の素数さん：2008/10/19(日) 12:52:51

>>7
そのパッケージのソースをCRANからダウンロードして中をみる以外の方法を聞いているの？

9 ：１３２人目の素数さん：2008/10/23(木) 23:18:23

初めてRを使うものです。
2000×2500のdata.frameの内容をファイルに出すのですが、
時間がかかってしかたありません。
いちいち、sprintfで文字列にして、catで出力するのですが、何か良い方法があればとおもい、
ヒント頂ければと思ってます。
sprintfのフォーマットは"%12.6f".
500万件の出力には、時間必要でしょうかね。
もちろん、ファイルはオープン＆クロースの手順を踏んで、
一つ一つはIO発生させないようにしてるんですが…
テキストファイルへの出力のヒント、頂ければと思います。

10 ：１３２人目の素数さん：2008/10/24(金) 21:28:35

>>9
write()でパイプ経由でawkとかperlに渡す場合とどっちが速い？

11 ：１３２人目の素数さん：2008/10/28(火) 17:33:32

Rを勉強しようと思って本でも買ってみようかと
思ってるんだが、やっぱRのプログラミングコードとかが
たくさん書いてある本のほうが良いかな。
それとも統計学とRの両方について詳しく説明してる本のほうが良い？

12 ：１３２人目の素数さん：2008/10/29(水) 14:44:29

何がしたいかによるんだから、具体的に書かないとレス貰えないよ。
コーディング経験の有無、統計に関する知識、その辺がすっぽり抜けとる。

13 ：１３２人目の素数さん：2008/10/29(水) 21:34:29

macでRを使ってるんですが、
グラフィックボードに日本語を出力しようとしても
全部四角になってしまいます。

解決方法ご存知の方がいらっしゃったら教えてください。

14 ：１３２人目の素数さん：2008/10/30(木) 11:37:13

>>11
基礎的なことはネットで資料を探す。
そうすると、Rの基礎みたいな本はほとんどいらないことがわかる。
というかRの使い方という感じで売ってる本だとRのマニュアルのコピー
みたいなので終わってるのが多い。

その上で、Rの使い方＆統計処理のつながりで、お気に入りをあげておくと
「Rによる医療統計学」丸善株式会社
内容は医療統計ではない。タイトルで損をしてるかもしれないが、
このタイトルだから買う人もいるかな。
内容的には、Rの使い方と分析方法がバランスよく書いてあり、使いやすい本。
Rの使い方、確率とか記述統計の話があった後に、１標本、２標本の検定、
分散分析、単・重回帰、ロジスティックモデル、生存分析が載ってる。

「RとS-PLUSによる多変量解析」シュプリンガー
いろいろな分析手法について、Rでどうやればいいかをメインに書いてある。
理論的な話はコラム扱いになっているので、そこは飛ばして読むことも可能。
主成分分析、因子分析、対応分析、クラスター分析、分散分析、重回帰・正準
相関分析、反復測定データの分析（パネルのランダムモデル）、あとはRの
使い方の備忘録がまとめてある。

「統計学：Rを用いた入門書」共立出版
これは現在目を通している所だけど、統計学の基礎をRで経験的に学ぶには
いい本かな。Rの使い方の理解にもいいでしょう。
ブートストラップなどの方法が載ってるのがPCを前提するといいかな。

15 ：１３２人目の素数さん：2008/10/30(木) 16:58:52

>>13
R.app だとQuartzGraphicDeviseが日本語に対応していない？ようです。
ですが、ちゃんと調べた訳ではないので、回避方法があるかもしれません。
また、グラフィックデバイスを使わずに直接ファイルに書き込む場合には、
日本語が書き込まれます。

また、Java で動く JGR というGUI環境では、過去にグラフィックデバイス上で
日本語が表示できたとの報告があります。詳しくは、
ttp://www.okada.jp.org/RWiki/?JGR%A4%C7%C6%FC%CB%DC%B8%EC

16 ：１３２人目の素数さん：2008/11/02(日) 13:53:00

>>13
Tigar+R 2.7.2だけど別に普通にquartzに日本語が表示されるよ

17 ：１３２人目の素数さん：2008/11/18(火) 17:45:37

RってExcelファイルとかのデータ取り込む手順がメンドイから嫌だ。

18 ：１３２人目の素数さん：2008/11/20(木) 02:03:05

plot()とかでx軸の目盛りを0,pi/2,pi,3pi/2,2piみたいにできますか？

19 ：１３２人目の素数さん：2008/11/20(木) 18:31:42

RcppTemplate
にCからRを呼ぶ例がみつからないのだけど
どこかにないのでしょうか

20 ：１３２人目の素数さん：2008/11/21(金) 01:04:01

>>18
できる
>>19
R以外の何かはそんなにデータを取り込む手順が面倒でないの？

21 ：18：2008/11/21(金) 05:57:35

どうやって？

22 ：18：2008/11/21(金) 05:58:27

安価忘れた
>>20です

23 ：１３２人目の素数さん：2008/11/21(金) 22:57:23

>>22
そんなの自分で調べろ。
> x <- (1:10)/5*pi; plot(x,sin(x),type='o',axes=FALSE); axis(1,c(0,pi/2,pi,3*pi/2,2*pi),c('0','pi/2','pi','3pi/2','2pi'))

24 ：１３２人目の素数さん：2008/11/22(土) 09:20:38

>>17
Rcmdrパッケージを導入。
RcommanderからだとGUIでエクセルをファイルを選べる。

25 ：１３２人目の素数さん：2008/11/22(土) 18:30:19

RODBCとかread.xls(file.choose())とかが出てきていない件

26 ：１３２人目の素数さん：2008/11/26(水) 19:30:17

うるさい。

27 ：１３２人目の素数さん：2008/11/29(土) 20:42:49

scale関数の結果はただのmatrixクラスですか？
末尾にあるattr情報はなんなんでしょうか？？


> scale(c(1,2,3,4,5))
[,1]
[1,] -1.2649111
[2,] -0.6324555
[3,] 0.0000000
[4,] 0.6324555
[5,] 1.2649111
attr(,"scaled:center")
[1] 3
attr(,"scaled:scale")
[1] 1.581139

28 ：１３２人目の素数さん：2008/11/30(日) 01:56:08

以前にもこんな質問があったような気がしますが、
過去ログが見られないので質問します。

文字列”123,456,789.123”を数値に変換するにはどうすれば良いのでしょうか？

as.doubleではエラーとなりました。

29 ：１３２人目の素数さん：2008/11/30(日) 11:35:09

>>28
専用の関数があった気がするけど、探すよりもカンマを取って変換する方がはやいかも
> a <- "123,456,789.123"
> a
[1] "123,456,789.123"
> b <- as.numeric(gsub(',','',a))
> b
[1] 123456789
> options(digits=12)
> b
[1] 123456789.123

30 ：１３２人目の素数さん：2008/12/04(木) 13:39:40

OpenOffice Calc のファイルを直接読み込んむことは
出来ないのでしょうか？
Excel 使わずに、 Calc で済ましているので、
CSVファイルに変換して読み込ませています。

31 ：１３２人目の素数さん：2008/12/07(日) 22:41:33

Rエディタのフォント（フォアグランド・バックグランド）を設定するにはどうすればいいのでしょうか？

32 ：１３２人目の素数さん：2008/12/07(日) 23:46:03

>>31
Rエディタって何？

33 ：１３２人目の素数さん：2008/12/08(月) 23:20:19

edit(x)

で出てくるエディタです。

34 ：１３２人目の素数さん：2008/12/10(水) 00:41:21

>>33
edit()って環境変数で指定しているエディタが起動するんじゃなかったっけ？
試しに、今やったら、Emacsが起動した。別のマシンでやったらviだった。
Emacsの色を変えたいなら、手っ取り早いのがcolor-theme。

35 ：１３２人目の素数さん：2008/12/10(水) 13:26:19

豊田秀樹の新刊、データマイニング入門−Rで学ぶ最新データ解析
http://www.tokyo-tosho.co.jp/books/ISBN978-4-489-02045-2.html

36 ：１３２人目の素数さん：2008/12/10(水) 20:36:05

Rって表形式じゃないから凄く見にくいんですけど
統計に向いてないのでは？

37 ：１３２人目の素数さん：2008/12/10(水) 23:40:56

>>36
逆だと思う。2次元にこだわるのは集計向きだけど、統計向きじゃない。
深い階層構造をもつデータや3次元以上の配列データをどうやって表にするの？

38 ：１３２人目の素数さん：2008/12/11(木) 21:47:43

>>36
それはね、>>36が統計に向いてないの。

39 ：１３２人目の素数さん：2008/12/12(金) 01:06:43

こんな過疎スレで煽りレス入れて楽しい？

40 ：１３２人目の素数さん：2008/12/12(金) 19:59:01

>>39
あおりではないよ。
統計をやる上では、２次元以上のデータを扱わないといけないのが普通。
だから、２次元でないと使いづらいと感じるならば統計処理がわかってない。
それなのに、Ｒの仕様の方に文句つけるってことは、>>36が統計に向いてない
ってこと。

なので、>>36は自分を反省したほうがいいよって内容のコメントです。

あ、>>39が煽りだったか、反応しちゃった。ｗ

41 ：１３２人目の素数さん：2008/12/12(金) 23:44:40

Rユーザの層が拡大するにつれ、>>40のようなユーザも出現するようになったのだなと
しみじみ時の流れを感じる今日この頃。

42 ：１３２人目の素数さん：2008/12/13(土) 11:54:59

>>41
RjpWikiをみると、昔から使ってる人のほうが態度が悪いのが多いと思うが。
昔はよかったってのは嘘じゃないかな？

43 ：１３２人目の素数さん：2008/12/13(土) 13:13:41

「昔はよかった」
とはいつの時代も言われること

44 ：１３２人目の素数さん：2008/12/13(土) 20:39:22

>>42
Rjpwikiの始発メンバーであんな回答態度をする人はいなかったけど。
途中参加でとっても辛口な先生1名だけいたけど、全体的にはとっても
雰囲気がよかった。最近の回答態度は酷すぎる。

45 ：１３２人目の素数さん：2008/12/14(日) 21:43:59

同感。
「俺達は金もらって仕事で答えてるんじゃない」という意識が強すぎるし、
質問者に対して「学術的態度で臨むこと」を求めすぎ。
ゼミで教授が学部学生をこき下ろすみたいな調子だもんね。
あそこは学者という人種に共通してありがちな人間的欠陥が
よく観察できるコミュニティだと思う。その点では興味深い。
Ｒ自体はとてもいいソフトだと思うけど、あの集団には近寄れない。
結果として自分１人、孤立無援で使うしかない人が大勢生まれてるはず。

46 ：44：2008/12/14(日) 23:16:21

私も大学教員の端くれだけど、「あの集団」の中の人かも知れないけど、
人がそれぞれ違うように学者もそれぞれ違うということを理解して欲し
い。それから、今のRjpwikiのノリは、技術系コミュニティに普遍的に見
られるノリをそのまま持ち込んでいるのであって、学者や学術とは関係
がないと思う。むしろ、プロのプログラマがボランティアのコミュニティ
に回答をお願いする態度を求められるような。
「技術系メーリングリスト」で検索すれば、私の言わんとしているとこ
ろが理解してもらえるかもしれない。

47 ：１３２人目の素数さん：2008/12/15(月) 00:01:57

今は見ていないけど。
１０年ぐらい前のニュースグループ fj の Java で、学者の罵りあいをはじめて見た。
煽(あお)り、嘲(あざけり)りが、しょっちゅうあったよ。
精神年齢が低いのか、世間と接する機会がないのか、異常な空間だった。
気持ち悪いから、次第に見なくなった。

48 ：１３２人目の素数さん：2008/12/15(月) 02:48:58

流れをぶったぎってすみませんがpersp()で作った図には目盛をいれられないのでしょうか？いろいろ検索して見たのですが目盛が入った図形を見かけません。

49 ：１３２人目の素数さん：2008/12/15(月) 03:35:23

>>48
demo(persp)
２ばんめのやつ目盛り入ってるけど

50 ：48：2008/12/16(火) 08:21:58

>>49
ほんとだ。ググってばっかでdemo見てなかった
ありがとうございます

51 ：１３２人目の素数さん：2008/12/25(木) 16:28:17

初心者です。
ランダムに一つの値を返す関数f()を使って、
n個のf()の要素を持ったベクトルを生成したいんですが
どうやってやればいいんでしょうか？
例えばrunif(1)を使って、
n個のrunif(1)を要素に持ったベクトルを作りたいんです。
（この場合はrunif(n)とすれば同じことなんですが。）
runifの定義を見たかったんですがうまく見れなかったんで・・・、
質問がわかりにくくてすいません

52 ：１３２人目の素数さん：2008/12/25(木) 17:35:41

>>51
for文で繰り返してf()を呼ぶんじゃあだめなわけ？

53 ：１３２人目の素数さん：2008/12/25(木) 21:11:17

>>51

>>52さんの説明に追加しておくと、後ろにｎ回足していけば要素ｎのベクトルになるからね。

54 ：１３２人目の素数さん：2008/12/25(木) 22:08:06

質問の意味がよく分かってないんだけど、素直に f(n) が長さnの配列を返すように
改変した方が早い気がするし、後々便利だと思うのだが…

55 ：１３２人目の素数さん：2008/12/25(木) 23:37:08

>>51
こういうこと？
> n <- 100
> sapply(1:n,function(x){runif(1)})
  [1] 0.818340010 0.651594423 0.127955778 0.171274880 0.563699459 0.384024996 0.029357331 0.016647415
  [9] 0.101343935 0.133003978 0.747144830 0.744196239 0.003289461 0.959092745 0.629223455 0.383857910
 [17] 0.037938089 0.006460069 0.773271457 0.150311993 0.196387192 0.836505357 0.412903159 0.095657646
 [25] 0.553321234 0.236131052 0.497607133 0.593558011 0.467087544 0.579318360 0.809348525 0.286566273
 [33] 0.023595126 0.898196327 0.829228607 0.051619953 0.442055013 0.116213767 0.310974789 0.256060483
 [41] 0.046852290 0.139351573 0.945665213 0.977257256 0.320223742 0.892197528 0.051279639 0.224508460
 [49] 0.634450256 0.776448613 0.048133771 0.008337753 0.330989671 0.897959799 0.160857094 0.760585248
 [57] 0.570918484 0.606187105 0.284859592 0.779051957 0.005772743 0.434992523 0.847252370 0.102974025
 [65] 0.977090385 0.820228452 0.456311142 0.405863907 0.626831524 0.559654842 0.301375986 0.118507505
 [73] 0.420883286 0.217497744 0.951507267 0.539087156 0.572545737 0.401707581 0.037067011 0.820561413
 [81] 0.796349923 0.014069674 0.416268822 0.466225379 0.514395235 0.907431938 0.528704567 0.993907087
 [89] 0.518415753 0.076568432 0.716981598 0.427751862 0.571849909 0.692945911 0.214970149 0.687965385
 [97] 0.634484604 0.064504115 0.484531934 0.719008875

56 ：51：2008/12/26(金) 16:54:57

返事遅れてすいません。
レスどうもですm(_ _ m)

>>52
最初for文でやろうとしたんですがうまくできなかったんです。
n<-100
x<-for(i in n) runif(1)
x
としてもベクトルになってもらえなかったので。

>>53
足す、とはどうすればいいんでしょうか・・・？

>>54
本当はそういう風にrunif(1)から自分でrunif(n)みたいなのを作りたかったんですが
ベクトルを返り値にするには関数をどう定義すればいいのかわからなくて。


>>55
おおーこんな感じです！
んーsapplyでできるんですね。わかりませんでした。
ありがとうございます！

57 ：１３２人目の素数さん：2008/12/26(金) 21:00:54

>>56

あなたのやってるのは、ｘを上書きしてるので、ベクトルにならないの。
足すってのは、ベクトルを後ろにのばしていけってことね。

ということでベクトルの追加の仕方をわかってないようなので、
すごく格好悪く書くと

x <- runif(1)
for (i in 1:9) x <- c(x,runif(1))

最初ので、xに値をいれ、次に、その後ろに値を9個追加してる。
大切なのは、x自体に値を加えたものをxに入れることで、ベクトルを
のばしてるってこと。

58 ：51：2008/12/27(土) 01:36:07

>>57
あ、なるほど、cで繋げばよかったんですね。
ただ、sapplyに比べてかなり効率悪いみたいです・・・(_ _ ；)
自分はかなり大きいデータを扱おうとしてたので
もともとforだと無理だったかもしれません

59 ：１３２人目の素数さん：2009/01/07(水) 18:50:15

質問です。これってなぜですかねえ？

GNU R,subset(),&&演算子の使用で結果が変に - gnarl、技術メモ’”%s<marquee>￥
http://d.hatena.ne.jp/gnarl/20090107/1231316900

60 ：１３２人目の素数さん：2009/01/07(水) 19:51:57

SEMのパッケージで構造方程式モデリングやってる面白い論文ありますか？
勉強半分遊び半分でいじってるんですけどうまく動作しなくて・・・

61 ：１３２人目の素数さん：2009/01/08(木) 18:39:38

Rの使い方に慣れる為にテキスト片手に勉強しているのですが、
乱数で生成したデータでクラスタ分析等をしても実感がわきません。

そこで、実際の実験や調査で得られたデータを使ってみたいのですが、
そういうものはどこかで手に入るものなのでしょうか？
（Rに付属してたりしないんですかね）

62 ：１３２人目の素数さん：2009/01/08(木) 19:28:14

>>61
R には有名な Fisher の iris データが付属しています。
あと、省庁が Web で公開している統計データ（都道府県別○○データとか）にも面白いものがあると思います。

63 ：61：2009/01/08(木) 21:31:59

さっそくありがとうございました。
data(iris)で呼び出せるとは知りませんでした。
農水省やら環境省を覗いてみましたが統計データがいっぱいありました。
いろいろいじってみようと思います。

64 ：１３２人目の素数さん：2009/01/21(水) 23:46:26

はじめまして。例えば2007/10から2008/9のデータが以下のようにあるとして、
10 0.5
11 0.4
12 0.3
1 0.4
2 0.5
3 0.6
4 0.7
5 0.8
6 0.9
7 0.8
8 0.7
9 0.6
（最初のカラムは年月日の月）
っていうのから折れ線グラフや散布図を作りたいんだけど、
簡単にRcmdrでやってみたら月がソートされてしまう。
データ順通り表示するにはどうしたらいいの？
なお、Rcmdrにはこだわってはいません。

65 ：１３２人目の素数さん：2009/01/22(木) 10:29:45

>>64
値の方は適当にしたよ。
> Oct07 <- as.Date("2007-10-1")
> time <- seq(Oct07, len=12, by="1 month")
> df.data <- data.frame(time = time, value = rpois(12, 3))
> plot(df.data$time, df.data$value, type="l")

66 ：１３２人目の素数さん：2009/01/22(木) 16:27:46

age

67 ：１３２人目の素数さん：2009/01/30(金) 00:42:42

統計フリーソフトとして、
RとWekaのどちらを学べばよいか迷っています。
どのような違いがあるのでしょうか？

68 ：１３２人目の素数さん：2009/01/30(金) 15:12:54

Rはほぼ初心者なんですけど・・・

A
[,1]
1
2
3
4


[,1] [,2]
4 8
2 9
1 0
3 7

というAとBのデータの1列目の一致させて結合したいのですけど、
どうしたらいいでしょうか？
教えてください。

69 ：１３２人目の素数さん：2009/01/30(金) 20:43:01

Aの並びは任意？それも昇順？
> B <- cbind(c(4,2,1,3),c(8,9,0,7))
> B[c(1,2,3,4),]
[,1] [,2]
[1,] 4 8
[2,] 2 9
[3,] 1 0
[4,] 3 7
> B[order(B[,1]),]
[,1] [,2]
[1,] 1 0
[2,] 2 9
[3,] 3 7
[4,] 4 8

70 ：１３２人目の素数さん：2009/01/30(金) 21:13:39

>>67
Weka は触ったこともないのだけど、データマイニングに特化しているイメージがある。
R なら、とりあえずなんでもできる。はず。
データマイニングしかしない、というのでなければ、R がいいのでは。

71 ：１３２人目の素数さん：2009/01/31(土) 08:27:08

>68
ヒント merge

72 ：１３２人目の素数さん：2009/01/31(土) 18:24:05

>>69
>>71
お答えいただきありがとうございます。
きっとmergeを使う方法がいいのだと思います。
それで一度やってみます。

73 ：１３２人目の素数さん：2009/02/01(日) 01:15:41

なんどもすいません。
merge関数を使う場合に・・・
>merge(data1,data2,all=TRUE,by.x=1,by.y=2)
という風にいれるとおもうんですが、
by.x=やby.y=
には何を意味しているのですか？

74 ：１３２人目の素数さん：2009/02/01(日) 10:34:51

help(merge)

75 ：１３２人目の素数さん：2009/02/05(木) 02:08:45

初心者です

単語の総数が違う二つのコーパスで　that　の出現頻度に
有意な差があるか調べたいんです

頻度（割合）を出したんですがこのあとどのような検定をすればよいのでしょうか？
値一つでは検定ができないですよね。

76 ：１３２人目の素数さん：2009/02/05(木) 10:04:58

>>75
chisq.test()

77 ：１３２人目の素数さん：2009/02/05(木) 11:08:08

>>76
回答ありがとうございます。chisq.test()について調べました。

今回の例の場合、引数には２行２列の行列を採用すればよいのでしょうか？
質問ばかりですみません

78 ：１３２人目の素数さん：2009/02/05(木) 14:19:54

Rを他の人が使いやすいように完全GUIから入力したいのです。
現在Rコマンダーというのがありますが半GUIという感じです。

可能でしょうか？
可能なら参考になるものがあれば教えてください

79 ：１３２人目の素数さん：2009/02/05(木) 14:43:04

>>77
四分割表が好みなら、それでいいよ。
> a <- matrix(c(20,55,11,43),ncol=2)
> colnames(a) <- LETTERS[1:2]
> row.names(a) <- c("that","others")
> a
A B
that 20 11
others 55 43

>>78
むしろ、S-plus

80 ：１３２人目の素数さん：2009/02/05(木) 15:22:46

>>79
本当にありがとうございます

独学なモノでなかなかわからないんですよね
参考にしてがんばります

81 ：１３２人目の素数さん：2009/02/06(金) 17:39:17

78です。
S-plusとか初めから揃っているのではなく、
自分でGUIを設計したいのです。
Rのソースを例えばdll化してVBから読み込むとかは可能でしょうか？

82 ：１３２人目の素数さん：2009/02/07(土) 00:06:06

S-PLUSは両方できるよ

それよりも、それを責任もってずっとサポートする
ことを考えたら、他のツールも考えていいのでは、
人件費が一番高いことを考えないといけない。

83 ：１３２人目の素数さん：2009/02/07(土) 03:23:55

>>82
一生の仕事にするので大丈夫です。
あと、どうしても市販のソフトは現場の方には難しいです。
目的でボタンを選んでいけば結果が出る(統計は裏でやっている)ようなものを作りたいんです。

ExcelからRを使えることもわかりました。試作を作ってみたいと思います。

84 ：１３２人目の素数さん：2009/02/08(日) 16:34:00

なんだろうこの子供は。

SPSSで良いじゃないか。

85 ：１３２人目の素数さん：2009/02/09(月) 08:49:39

SPSSは統計をわかっている人向け。
現場の人は回帰分析が何に使えるのかさえ知らない。現場を知らなさすぎるよ。

86 ：１３２人目の素数さん：2009/02/09(月) 11:08:59

強者が妬まれるのはこの世の常

87 ：１３２人目の素数さん：2009/02/09(月) 20:49:02

>>83
>目的でボタンを選んでいけば結果が出る(統計は裏でやっている)ようなものを作りたいんです。
ボタンを選んでいくとRのコードを出力するようなものを目指すといいと思う。
そうすれば、教育目的にも使えそう。

あと、ど素人に使えるようにするには、データを解析にかけられるレベルに整形する機能
（データのマージ、欠測値の置き換え、変数変換、欠測値を含む個体の除外など）
を持たせる必要があるのでは？大変だと思う。

88 ：１３２人目の素数さん：2009/02/09(月) 22:10:36

教育は目的ではないだろ。ｗ
HTMLからRを呼んで結果をWindowに表示するものを作ったことがあるけど、
Rの結果をそのまま貼ったのでは現場ではよろしくないだろう。
出力整形の方が手間がかかると思う。

89 ：１３２人目の素数さん：2009/02/09(月) 22:43:35

>>87
＞ボタンを選んでいくとRのコードを出力するようなものを目指すといいと思う。
現場の方はＲとかには興味ないんです。
興味があるのは結局どうすれば改善するの？答えを教えて！
ということだけなので、Ｒのコードすら見せないつもりです。
最終結果のみ出力ですね。
（当然、必要な人のために計算手順、考え方を示したレポートも出力しようとは思っています）

欠損値を含むデータは一式除去します。
（補完は事実ではないので使いたくないしデータは多数あるので）
正規性検定で正規性チェックして
（当然メッセージボックスには正規性チェックＯＫとか処理の意味を記載する）
当然多変量散布図も表示。
（→何をポイントとして見ればよいのかも素人にわかりやすく示す予定）

変数変換はBox-Cox変換でします。（あてはめすぎには注意したいですが。。。）
で（線形）重回帰分析を解く

というところまではＶＢＡで作っていたので。

今後、他の手法もＲで手に入るならば後は全体の分析手順を考えるだけだなと思いまして。
Ｒの本も購入したのでゆっくりやってみます。

90 ：１３２人目の素数さん：2009/02/09(月) 23:16:15

そのような（ほとんど同じ）ことを事を考える人多いね。
数人ではなく、何十人聞いたことか。
成功させるためには、商売にしてしまうという
考えをしたほうが良いのかもしれない。

91 ：１３２人目の素数さん：2009/02/09(月) 23:19:51

流れ切って悪いが質問

シーケンスなデータがベクトルに入ってるとして、
その極大（小）値を求めることってできる？

92 ：１３２人目の素数さん：2009/02/09(月) 23:31:28

>>90
商売？？
もともと社内でそういう仕事やってるので。。。

93 ：１３２人目の素数さん：2009/02/09(月) 23:40:00

外販が商売なの？

94 ：１３２人目の素数さん：2009/02/09(月) 23:42:02

だとすれば、可能性有るかもしれない。
それを考える人は結構すくなくて、
社内利用を目的とする人がおおくて
それだと、いつの間にか消えていくもので。

95 ：１３２人目の素数さん：2009/02/09(月) 23:44:15

>>89
>Ｒのコードすら見せないつもりです。
R のコードを見せるかどうかは問題ではなく、どうすれば（あなたが）楽をできるか、を考えたつもり。
作成段階でのデバッグとか、保守（不具合が発生した時の対応）とか。
もちろん、ユーザーには解析結果だけ見せればいいだろうけど、R はめちゃくちゃ強力なので、
解析結果の表示の寸前まで R にやらせればいいと思う。

>興味があるのは結局どうすれば改善するの？答えを教えて！
工場の品質管理？

>変数変換はBox-Cox変換でします。（あてはめすぎには注意したいですが。。。）
なんてダーティーな…。
（試行錯誤がある程度許される）品質改善に使うなら、そういうのもありなのかな。
製薬畑で、対数変換しか使ったことがないからよくわからない…。

>>90 >>91
システムとして他の企業に売ったりするつもりが無いなら、オープンソースでもいいかと。
企業内の統計専門家なら実験計画・解析計画と結果の解釈で飯を食えるはずだし。

96 ：１３２人目の素数さん：2009/02/10(火) 00:53:53

>>95
自分の楽を考えていては改革は進まないのです。
楽する人には誰も信頼してついてきてくれません。
現場がいかに楽して成果をだせるか、
統計専門家は一段降りて素人と対話することが重要なんです。
わからない専門用語を並べても相手には伝わらないんです。

統計はあくまで手段にすぎないのです。
統計がなくても目的が達成できれば、統計すらいらないのです。
統計を使うことが目的ではなく、現場をよくすることが目的だからです。
本来は現場の方が自主的に統計を学べばよいのですが、現場は現在の仕事で手一杯状態。
そこを補佐するのが私の任務です。

あと、わかっているとは思いますが、
Box-Cox変換の中には対数変換も含まれていますよ。
変換は、平方根変換、対数変換が主でしょうが、
プログラム内ではBox-Cox使う方が値を変更するだけなのでコード記載が楽なのです。
モジュール一つで済みますから。

＞飯を食えるはずだし。
飯を食えるとか発想が保守的すぎる。
現場が望むものなら不況であれ何であれやろうってなりますよ。
改革派の俺にはわからん考え方ですな。。。

97 ：１３２人目の素数さん：2009/02/10(火) 01:15:52

若いな…。
それを「業務」としてやっているなら>>95のいうように飯をくっていることに変わりはない。
まあどんなに簡単に使えるツールを作っても業務が改善されなければ現場からは
役に立たないと言われるだろう。難しいぞ。

98 ：１３２人目の素数さん：2009/02/10(火) 01:21:41

難しいから仕事なのです。
自分しかできないから存在価値があるのです。
誰でもできることばかりやっていては
いきなり首を切られても文句は言えません。

99 ：１３２人目の素数さん：2009/02/10(火) 01:28:46

目的が『飯をくうために』。結果『飯をくっている』
で考え方は違うし、行動に差がでるよ。

前者は、飯をくうために、できる限り楽して金をかせぐことを考えてしまう。
後者は、目的を達成することに全勢力を注ぎ込む。結果的に周りからも認められ金を得られる。

100 ：１３２人目の素数さん：2009/02/10(火) 01:53:47

流れ切って悪いが質問

シーケンスなデータがベクトルに入ってるとして、
その極大（小）値を求めることってできる？

101 ：１３２人目の素数さん：2009/02/10(火) 10:21:21

min(), max() って話ではないのか？？？

>>83
応援だけする（笑。ただ、そういうものなら裏でR走らせるよりも、
自分でコードした方が楽かもしれない、とも思う。

102 ：100：2009/02/10(火) 12:20:16

>>101
いや、極値を全部取り出したかったんだが
まぁ自分で実装するか

103 ：１３２人目の素数さん：2009/02/10(火) 13:16:18

RgtkやR-Tcl/Tkが用意されているんだから、GUIインターフェイスはいかようにも作成できるし、HTMLやWebアプリにしたければ、
Rpadがすでに用意されているんだから、車輪の再発明はいらないと思う。

104 ：１３２人目の素数さん：2009/02/10(火) 13:53:28

>>102
あ、ごめん、勝手に最大最小に置き換わってたw

105 ：１３２人目の素数さん：2009/02/10(火) 19:20:18

>>110
確実にR使う方が楽だよ。

106 ：１３２人目の素数さん：2009/02/10(火) 21:57:33

>>110に期待

107 ：１３２人目の素数さん：2009/02/10(火) 23:28:47

やっぱり若いな。
手段が目的化するのは普通のことだし、そもそも業務ならこんなところで
質問してないで仕事しろって感じ。ｗ

108 ：１３２人目の素数さん：2009/02/11(水) 00:53:12

>>96
>楽する人には誰も信頼してついてきてくれません。
楽をしろ、というのは、性能・品質を妥協しろ、とか、素人にわかりにくい出力をそのまま出せ、という意味ではない。
同じ動作をするシステムを作るなら、なるべく楽をして作るべき（無駄に苦労するのは社会人として失格）だ、ということ。
そのための方策として、結果出力の寸前までを（他の言語に投げずに） R に任せることを勧めている。

あと、デバッグ・保守を楽にできないのだとしたら、それは稼働率も信頼性も低いダメなシステムを作ったということ。
プログラミングに限っていうと、一般に楽につくったシステムはデバッグ・保守も楽な「良い」システムになるから、一石二鳥。

>>98
>誰でもできることばかりやっていては
>いきなり首を切られても文句は言えません。
あなたは別の理由で職を失う心配をした方がいいかも。主張がちょっぴりうつ病患者っぽいよ。

109 ：１３２人目の素数さん：2009/02/11(水) 01:33:44

鬱の患者が多い駄目会社で常人扱いされるあながが一番よくわかる症状らしい。
通常ならなにもやらないやつが切られるのにな(笑)
ぼく、仕事はお遊びじゃないんだよ。

110 ：１３２人目の素数さん：2009/02/11(水) 01:36:15

>>107
完全に慢性病。
ご愁傷さま。

111 ：１３２人目の素数さん：2009/02/11(水) 02:00:01

>>108
詭弁。口がうまいというか世渡り上手というか。
自分が不利にならないように話の解釈を変えている。
口がうまい分実力がなくても大丈夫。
実力者についてわたり歩けばいいんだから。
口がうまい人なりの仕事の仕方はある。
否定はしない。

112 ：１３２人目の素数さん：2009/02/11(水) 11:28:32

完全に子供だろ。社会人１年目か２年目の。

113 ：１３２人目の素数さん：2009/02/11(水) 13:15:02

私の会社にはこれだけ意思のある人はいないな。
誰も動こうとしないから楽なことは楽だが会社が潰れないか心配だ。

114 ：１３２人目の素数さん：2009/02/11(水) 14:50:50

俺らは何かを忘れているのかもしれない。
衰退に気づかない『ゆでかえる』かもしれない。
俺には恥ずかしながら現状を具体策をもって行動できるパワーはない。
だからその課題解決力が羨ましく思う。
反対するやつは必ずいる。負けずに行動してほしい。

115 ：１３２人目の素数さん：2009/02/11(水) 15:23:12

どこにでも行動できない何したらよいかわからない
ただ違うことはわかるという評論家のようなやつはいる。
違うことすらわからず、
人の性格、年齢等まったく異なる話にもっていき攻める幼稚なやつさえいる。
子供か大人かは話の論点がそもそも違う。

会社に利益をもたらすかどうかが重要。
優秀さとは、自分のアタマで考えることができ、行動できる人だと思う。
がんばってほしいと思う。

116 ：１３２人目の素数さん：2009/02/11(水) 16:13:54

夢だけ見てるあほだろ。

統計も統計ツールも現場も知らないお子様。
恥ずかしい限り。

> Rを他の人が使いやすいように完全GUIから入力したいのです。
> 現在Rコマンダーというのがありますが半GUIという感じです

> 目的でボタンを選んでいけば結果が出る(統計は裏でやっている)ようなものを作りたいんです。。

>SPSSは統計をわかっている人向け。
>現場の人は回帰分析が何に使えるのかさえ知らない。現場を知らなさすぎるよ。

117 ：１３２人目の素数さん：2009/02/11(水) 16:33:30

知識は後からつけられるが、意識の高さ、考え方は後からつけるのは難しいからな。

あと>>116
具体策何か提示しないと
お子様にも負けていることを自分で言っているようなものだと思う。

118 ：１３２人目の素数さん：2009/02/11(水) 16:38:01

>>117
いや、>>116も夢と表現しているところからすると
内心はお子様の目指している方向の正しさを認めているようだ。
しかし、実現させたいが自分にはできないといったところじゃないか。
素直にはなれない頑固さというところか。

119 ：１３２人目の素数さん：2009/02/11(水) 16:59:32

自分のアタマで考えることができ、行動できる人

120 ：１３２人目の素数さん：2009/02/11(水) 17:05:52

俺は>>116じゃないけど、無理な方向だとは思うよ。

統計をわかってないと、そもそも「目的でボタンを選んで」ができない。ｗ

「現場の人は回帰分析が何に使えるのかさえ知らない」のなら、そんな人が
回帰分析をしてはいけない。

ある程度の不便さは、使う側を選別するために必要だと思う。

121 ：１３２人目の素数さん：2009/02/11(水) 17:12:01

統計解析あ〜る

122 ：１３２人目の素数さん：2009/02/11(水) 17:16:14

>>120
無理じゃなくて。
できるとすれば、具体的にどうすればよいと思うか？
否定だけでは今のところみな評論家のままだよ。
俺同様、この点でお子様に負けていると思うんだよ。

123 ：１３２人目の素数さん：2009/02/11(水) 17:18:37

>>120
＞ある程度の不便さは、使う側を選別するために必要だと思う。
これ彼が言ってた。
統計専門家は現場を意識していない。
素人と向き合っていないってことだね。
変われない専門家じゃないかよ。

124 ：１３２人目の素数さん：2009/02/11(水) 17:28:04

>>122
いや、統計学については、中卒、高卒でもわかるようにノウハウは作られてる。
そのために先人がいろいろとテキストを書いてきてるし、きちんとしてる所は
社員に対してそのためのトレーニングもしてる。

問題は、ソフトをある程度以上に簡単化することは、統計の誤用を引き起こす
原因になるってこと。
私はそっちの方がリスクが高いと思ってる。

ということで、統計をできる人間を入れて、やらせればいいんであって、
みんなにやらせる必要はないよ。その方が無駄ってのが私の考え。

125 ：１３２人目の素数さん：2009/02/11(水) 17:29:32

>>120
＞統計をわかってないと、そもそも「目的でボタンを選んで」ができない。ｗ
これは違うと思う。
例えば株なら○％で上がる株を探したいとか目的がある。
それは株の知識があればできる。いわゆる現場の人ならできるんだよね。


＞「現場の人は回帰分析が何に使えるのかさえ知らない」のなら、そんな人が
＞回帰分析をしてはいけない。
＞ある程度の不便さは、使う側を選別するために必要だと思う。

当然、現状のＳＰＳＳ等のソフトのような場合、わからずに使っては怪我をする。
ならこんなのはどうだろう。指導者が横にいれば使ってもよい。だろ？

彼は指導者の部分をプログラムに組み込もうとしていると思うんだよ。
素人が使っても正しく回帰分析を使うことができている。
しかも意識せずにというのが使用者にストレスがなくていいね。

126 ：１３２人目の素数さん：2009/02/11(水) 17:29:37

ついでに書いておくと、このデータを分析してくれって持ってこられても
分析にたえないデータってあるよね。

データの取り方から計画しないといけないことを考えると、ソフトだけ自動化
してどうするの？ってのがある。

127 ：１３２人目の素数さん：2009/02/11(水) 17:35:48

>>125
>例えば株なら○％で上がる株を探したいとか目的がある。
>それは株の知識があればできる。いわゆる現場の人ならできるんだよね。

知識があれば上みたいなことは考えない。
こういう喩えを出す段階で、>>125の知識が危ないでしょ。
>>125は、○％で上がる株を予測できるシステムを組めるの？
（まあ、いまだとありませんとかでそうだがｗ）

また、素人に正しく回帰分析をさせるような指導込みのプログラムを組んだとしよう。
すると、「それはできません」「その処方は間違っています」と指導される。
これだけじゃ素人には改善できないので、なぜダメなのか理由も説明することになる
だろう。

でも、そんなこと言われたら、「意識せずにストレスにならずに使える」には
ならんでしょ。
その理由自体が読んでもわからないかもしれないのに。

128 ：１３２人目の素数さん：2009/02/11(水) 17:37:55

>>126
それアンケートの話じゃないかな。
多分彼がやりたいのは製造品質っぽい。
計画部分の重要性が薄れてくる分野かも。

129 ：１３２人目の素数さん：2009/02/11(水) 18:02:09

>>127
＞すると、「それはできません」「その処方は間違っています」と指導される。
それは自動でやってくれたらいいよね。
間違わないように指導してくれるんだから、駄目だというだけでは指導にならない。

＞知識があれば上みたいなことは考えない。
＞また、素人に正しく回帰分析をさせるような指導込みのプログラムを組んだとしよう。
君の考える指導込みのプログラムって何かな？
知識があるなら言えるってことだよね。

130 ：１３２人目の素数さん：2009/02/11(水) 18:08:55

どんなUIをイメージしているんだ？

具体的に聞き出していったらSPSSそのもので笑わしてくれる気がしてならない。
既存のツールだって十分に統計素人にも使えるように設計されているぞ？
EnterpriseMinerもClementineもSPSSも、そしてRもな。

細かいパラメータの設定なんて一切なしの、自然言語で全て入力できるような
UIを考えているのか？

まあ、書き込み内容を見ていると統計知識も情報技術も無いのがすぐ解るけど。

131 ：１３２人目の素数さん：2009/02/11(水) 18:11:57

>>124 も指摘してるけど、年間に何時間か割いてもらって社内で統計の勉強会を開いた方がいいと思う。
ある程度のリテラシーは意思の疎通を大いに助けるから。
その時間すらとれないほどのブラック企業なら、この不景気で明日にでも倒産しておかしくない。

>>125
>素人が使っても正しく回帰分析を使うことができている。
よくご存じですねw
でも、意識が高い連中なら、統計の入門書を買って読むぐらいのことは既にしていて、本当の素人ではないと思う。
「マンガでわかる統計学」、「マンガでわかる統計学 回帰分析編」とか、売れまくってるし。
大卒のエンジニアなら、数学の才能は下手をすると彼を超えるだろうし。

132 ：１３２人目の素数さん：2009/02/11(水) 18:14:56

Rのお陰で統計に対する敷居が低くなったのが良く解る。 ・д・

133 ：１３２人目の素数さん：2009/02/11(水) 18:15:47

>>129
突っ込みたい一部分だけを取り出すのは、ネットではよくあることだが、
>>126の君が引用した部分の後も引用してもらえないかな。

>これだけじゃ素人には改善できないので、なぜダメなのか理由も説明することになる
>だろう。
>でも、そんなこと言われたら、「意識せずにストレスにならずに使える」には
>ならんでしょ。
>その理由自体が読んでもわからないかもしれないのに。

要は指導付きプログラムを作ったとして、どこまで指導するのかってこと。
結局、教科書と同じになったら意味がないんだけど。

でだ、>>129は逆質問で返してきてるけど、君の持ち出した指導込みプログラム
ってのは具体性がないので、何も言ってないのと一緒だよ。
まずは、どういう風に指導するプログラムなら、ド素人でもわかるのか言いだしっぺ
が例をだしてもらえる。

私の意見は、まったくのど素人にはムリ、専門家を育てて配置するか、それなりに
勉強しろってことなので、君の方から示してくれないと答えようがないんだけど。

私の目からうろこが落ちるような、すばらしい発想なんだと思うけどね。

134 ：１３２人目の素数さん：2009/02/11(水) 18:17:40

>>131
＞でも、意識が高い連中なら、統計の入門書を買って読むぐらいのことは既にしていて、本当の素人ではないと思う。
＞「マンガでわかる統計学」、「マンガでわかる統計学 回帰分析編」とか、売れまくってるし。

マクロとミクロの話がごっちゃ。
素人は話にならん。本当に統計の誤用とはまさにこのこと。

135 ：１３２人目の素数さん：2009/02/11(水) 18:18:02

>>132
Ｒバンザイ！

136 ：１３２人目の素数さん：2009/02/11(水) 18:20:37

>>133
思考停止。慢性病患者だね。
自分で考えようとしない。

137 ：１３２人目の素数さん：2009/02/11(水) 18:22:27

全世界に現場がある会社だったらかなり厳しいな。
とふと思った。

138 ：１３２人目の素数さん：2009/02/11(水) 18:27:36

これを実現するには
統計＋プログラミング＋現場知識が必要になるね。

139 ：１３２人目の素数さん：2009/02/11(水) 18:28:25

>>138
実際、統計専門家は現場知らないから難しいんだよ！

140 ：１３２人目の素数さん：2009/02/11(水) 18:29:59

>>132
優秀な現場のおっちゃんが来た感じ

141 ：１３２人目の素数さん：2009/02/11(水) 18:30:40

結局インテリジェントＣＡＥってことだな。

142 ：１３２人目の素数さん：2009/02/11(水) 18:31:06

>>134
解説よろしく！
「マンガでわかる統計学」の英語版を口語表現の参考に読んだことしかないのだけど、
普通に記述統計・推測統計の基礎を書いたものだったよ？

143 ：１３２人目の素数さん：2009/02/11(水) 18:32:16

>>142
また統計素人がでてきた。。。

144 ：１３２人目の素数さん：2009/02/11(水) 18:32:54

>>136
相手に「自分で考えよとしない」と書いてるけど、君は考えてるの？
できるって言うなら案を出してよ。
できるといってる方が案を出すものなのに、できないといってる人（私）に案を
出せっておかしくない？ｗ

できるって言うなら、いい案を出して、私を納得させてよ。

何も案を出さずに、
指導込みのプログラムにすればいい！
その内容は誰でもわかるものでー！
ど素人でも間違いなく使えるー！
なんていってても、単なるバカだよ。

私にすれば、>>136は現実を見ることができず、大きなことを言ってるだけで
何も作り上げることのできない、理想主義者に思える。

145 ：１３２人目の素数さん：2009/02/11(水) 18:36:17

>>138
ソフトはある程度汎用にしておいて、あとはユーザーが自分の置かれた状況に
合わせて勉強するのが一番早い。

146 ：１３２人目の素数さん：2009/02/11(水) 18:36:29

う　ぜ　ぇ



IDがないから何人で議論してるか分からん

147 ：１３２人目の素数さん：2009/02/11(水) 18:38:39

>>144
できないという理由は、現実、誰にも証明できないから
できないとは誰も言えないはずなんだよ。

しかし君はできないと主張してる。
その時点で話す価値のない対象になってるってことに早く気づけ。

148 ：１３２人目の素数さん：2009/02/11(水) 18:40:59

統計知識も情報技術も現場も知らないガキであることがわかりました。

149 ：１３２人目の素数さん：2009/02/11(水) 18:42:57

まじプロファイラーだな。

150 ：１３２人目の素数さん：2009/02/11(水) 18:49:19

上に書いてある。
彼は既にＶＢＡで一度プログラムしてるようだな。
>>138の統計＋プログラムの力はあるようだ。
それで使われなかったこともあったんだろう。

現場知識の現状も把握しているようだ。

151 ：１３２人目の素数さん：2009/02/11(水) 18:51:38

>>147
できないことの証明はできないにしても、それが「だからできる」って
証明にはならんのよ。
まったく、案を示さないままに、上記のようなことを言ってる段階で痛い。

数学スレにいるのに、頭悪いと思う。

それとだ、私のできないってのは、他の案のほうが効率的なのに、安くできるのに
無意味ってのも入ってる。
元々の人はSPSSですら、ド素人にはムリと書いているので、あれ以上のUIおよび
簡略化の上で、>>147は現場に合わせた指導まで入れるつもりだよね。

素直に、勉強会開いて統計を教えた方が早いよ。

152 ：１３２人目の素数さん：2009/02/11(水) 18:58:52

>>151
本人ではないが少なからずできるよ。
わざわざできる理由を聞いてる方が無知に見える。

インテリジェントCAEが一例。

153 ：１３２人目の素数さん：2009/02/11(水) 19:00:56

>>150

>>147だけど、元の人の発想については、業務を限定して、その作業が定型化されて
いるならばUIを制限することで、作業を効率化することは可能だと思う。
まあ、マクロ組めばいいって気もするが。

でも、>>125の場合は、そもそも適用する対象からしておかしいし、そのための
手法も夢を見すぎ、しかも具体案がないので、突っ込んでる。

154 ：１３２人目の素数さん：2009/02/11(水) 19:03:32

>>134
単に、世の中で統計の本が売れまくってるから、彼の職場でも統計の本を読む奴がいておかしくない、ってことだよ。
ひょっとして、彼の職場の社員集団は日本全体の無作為サンプリングではない、とか考えた？
なんでもかんでも統計学的な言説だと解釈してはいかんよ。

>>150
>現場知識の現状も把握しているようだ。
なぜ”彼”の言うことを鵜呑みにできるのか、わからないね。

>>152
横レスだが、インテリジェントCAE を使うのは素人じゃないと思う。

155 ：１３２人目の素数さん：2009/02/11(水) 19:04:20

>>150
正直、そんな低レベルな話をされても困る。

現場のことも実は何も理解していないから
現場の人間に受け入れてもらえなかったんだよ。

頭でっかちのルーキーには良くある事。
自分もそうだった。

156 ：１３２人目の素数さん：2009/02/11(水) 19:18:37

う　ぜ　ぇ

157 ：１３２人目の素数さん：2009/02/11(水) 19:27:50

>>152

>>151、>>153だが、>>152が見えてなかった。

>>151でも書いたけど、業務が定型化していて、やることが限られているなら
いけると思うよ。

でも、インテリジェント機能にしても、元の使い道がわからんとどうしようもない。
ある程度わかってる（本当は相当に精通してる）から、そのサポートで効率化できる
んであって、ド素人が使うものじゃないでしょ。
ド素人はサポートされると混乱する．．．

結局は、少しは勉強してド素人を脱出した方が早い。

158 ：１３２人目の素数さん：2009/02/11(水) 20:05:35

>>155
現場を理解してないのは俺らも同じだろ？
現場に統計使え、おまいら教えるから勉強しろと。
けど数学得意な俺らにすれば理解はたやすいが得意でないのが多い。

159 ：１３２人目の素数さん：2009/02/11(水) 20:09:46

>>157
君は一生成功しないな。
現場は統計はど素人だが現場の専門家だよ。
俺らは現場の素人。
現場からしたら少しは現場を勉強しろとなる。

どちらかが歩み寄らねばならない。
いがみ合ってもお互いにとって損だ。

160 ：１３２人目の素数さん：2009/02/11(水) 20:13:14

ちなみに、現場で活躍中の方に質問なんだけど
どんな現場なの？
ほんとに統計勉強したことない人ばっかりなの？
そういう会社で統計使うんかな？

161 ：１３２人目の素数さん：2009/02/11(水) 20:19:06

>>160
現場知らない人1人発見

162 ：１３２人目の素数さん：2009/02/11(水) 20:22:32

>>160
統計知識がある君はどこに就職したい？
友達はどこに就職したいと言っている？

統計人材が流れる先は大方、アクチュアリ、リサーチャーじゃないかな。

163 ：１３２人目の素数さん：2009/02/11(水) 20:28:42

コンサルティング会社ってことですか

164 ：１３２人目の素数さん：2009/02/11(水) 20:36:43

就職したい＝統計を専門に活用する仕事

逆に統計を専門に扱わない、あくまでも手段の一つ、ただ使った方がより効率的な仕事。

こんな現場ではこのスレのような問題が起きる。

165 ：１３２人目の素数さん：2009/02/11(水) 21:16:43

結局なんの具体案も書けないスイーツの夢物語だったな。

166 ：１３２人目の素数さん：2009/02/11(水) 21:25:27

統計の現場は多数ある。
具体案はそれぞれで考えるしかない。
まして現場知識がないのにこのスレで具体案がでるわけがないと思う。

各個人が現場と話し合って作り上げることが必要になる。

167 ：１３２人目の素数さん：2009/02/11(水) 21:28:54

誰も現場を知らないことがよくわかった。

168 ：１３２人目の素数さん：2009/02/11(水) 21:32:53

一体誰だよ現場、現場って言い出したのは・・・

169 ：１３２人目の素数さん：2009/02/11(水) 21:36:30

たぐってみると、質問者の方のようですね

170 ：１３２人目の素数さん：2009/02/11(水) 21:42:01

いや傍観者もいるよ。

171 ：１３２人目の素数さん：2009/02/11(水) 21:43:06

よくよく考えると統計って単なる手段だよなと気付く。

172 ：１３２人目の素数さん：2009/02/11(水) 21:46:34

ふむふむ

173 ：１３２人目の素数さん：2009/02/11(水) 21:52:22

統計をもっといろんな分野に使ってもらえば俺らの仕事の領域が広がるんじゃないかな

174 ：１３２人目の素数さん：2009/02/11(水) 21:54:51

あー、もう一回学生に戻ってダラダラ勉強したいな
ストレス社会嫌・・・

175 ：１３２人目の素数さん：2009/02/11(水) 21:57:08

>>158 >>159
製薬会社の統計担当者だけど、薬事法、生理学、病理学、メジャーな薬や開発中の
新薬の作用メカニズムを新人研修で教わったし、その後も自主的に勉強してるよ。
生物学は高校でも選択してなかったから、たしかに大変だったけど、
営業マン（MR）相手なら作用メカニズムを教えることぐらいはできるよ。
逆に、医学・薬学・獣医学の出身者に対する新人研修には統計学が含まれるし、
工場でも研究所でも自主的な統計の勉強会をやっているよ。

どちらかが歩み寄るのではなく、どちらも歩み寄るのが当然だと思う。
職場のパートナーがやっている仕事に関する知識が素人レベル、
なんてあまりに非効率的で、冗談にもならない。どの業界・業種でもきっとそう。

>>167
彼自身がどういう業種・現場で働いているか、ということを明かさないのだから、仕方ない。
よほどマイナーな業種でなければ、統計・現場知識の両方で彼を凌駕する者はいると思う。
明かせなくなってしまったのかもしれないけど。

>>171
同感。でも、統計に限らず、学問っていうのは手段。

176 ：１３２人目の素数さん：2009/02/11(水) 22:06:53

製薬は統計活用が進んでるからな。
生物統計、医学統計ってのがあるぐらいだ。

177 ：１３２人目の素数さん：2009/02/11(水) 22:09:03

統計担当ってのがある時点でうらやましいよな。

178 ：１３２人目の素数さん：2009/02/11(水) 22:23:24

oioi、最初から読み直せ。

彼は本当に初歩の初歩的な質問しかしていないぞ。
周りが話を大きくし過ぎ。

> S-plusとか初めから揃っているのではなく、
> 自分でGUIを設計したいのです。
> Rのソースを例えばdll化してVBから読み込むとかは可能でしょうか？

179 ：１３２人目の素数さん：2009/02/11(水) 22:30:46

現場、とか言い出して話を大きくしたのは彼だと思う。

180 ：１３２人目の素数さん：2009/02/11(水) 22:35:42

いや、彼女の方こそ。

181 ：１３２人目の素数さん：2009/02/11(水) 22:37:20

たぶんおもしろがって煽っていただけの周りだと思うが。
彼はいずこへ！？

182 ：１３２人目の素数さん：2009/02/11(水) 22:38:11

いや、彼女の方こそ。

183 ：１３２人目の素数さん：2009/02/11(水) 22:38:41

それほどでも、てへっ！

184 ：１３２人目の素数さん：2009/02/11(水) 22:39:33

すみませんｽﾚﾁですが、
ExcelのVBAで多変量解析を駆使し、投資用ソフト開発出来る方いましたら500万円で買います。
ロジックは私が用意しますからプログラミングだけお願いします。
ご連絡お待ちしております。
motorhead@55mail.cc

185 ：１３２人目の素数さん：2009/02/11(水) 22:40:24

>>181
きいてらんね〜や・・・・と
仕組みを作り始めてる

186 ：１３２人目の素数さん：2009/02/11(水) 22:42:53

まあやり方はわかったんだから
あとは実行するだけだろ、彼ならできそうだ。

187 ：１３２人目の素数さん：2009/02/11(水) 22:43:53

俺はいったい誰と戦っていたんだ。。。

188 ：１３２人目の素数さん：2009/02/11(水) 22:44:23

>>187
たぶん俺ら煽り隊かな。

189 ：１３２人目の素数さん：2009/02/11(水) 22:45:18

>>178-179
その彼と、今回話を大きくして、現場とか言ってるのは別人じゃないのか？

今回の彼が持ち出した現場の例が「株」だしなー．．．

190 ：１３２人目の素数さん：2009/02/11(水) 22:45:41

煽りたいから煽り隊！

191 ：１３２人目の素数さん：2009/02/11(水) 22:53:54

>>184
Excel の VBA、っていうきつい縛りは何のため？エンドユーザーにソフトを売ることを考えてるのかな？
ロジックの特許を押さえるか、ガチガチの契約を結ばないと丸ごと持ち逃げされるよ。特許がとれたら、このスレで報告してね。

192 ：１３２人目の素数さん：2009/02/11(水) 23:06:36

こういうとき自分でプログラムできたらなーって思う。

193 ：１３２人目の素数さん：2009/02/11(水) 23:10:33

>>191
ＶＢＡでプログラムした彼ならできるのでは。。
しかし彼はいない。。

194 ：１３２人目の素数さん：2009/02/11(水) 23:22:00

彼はずっと張り付いてるジャンｗ

195 ：１３２人目の素数さん：2009/02/11(水) 23:28:19

>>194
いや、崇高な理念に燃えていたので、こんなところで時間をつぶしてはいないはずさ！！

196 ：１３２人目の素数さん：2009/02/11(水) 23:31:33

この程度で崇高な理念ｗ

> S-plusとか初めから揃っているのではなく、
> 自分でGUIを設計したいのです。
> Rのソースを例えばdll化してVBから読み込むとかは可能でしょうか？

197 ：１３２人目の素数さん：2009/02/11(水) 23:44:59

すごい数のレス付いてるから、何か大きなアップデートでもあったのかと思って覗いてみれば…

それにしても、彼とか彼女とかわけわからんw
コテ付けてくれ。

198 ：１３２人目の素数さん：2009/02/11(水) 23:48:34

いつの間にかこんなに盛り上がったのか。ｗ
エキスパートシステムが瓦解したことを知らない世代が増えたってことだな。

199 ：１３２人目の素数さん：2009/02/11(水) 23:52:23

>>198
あんまりLispをいじめないで下さい(><)

200 ：１３２人目の素数さん：2009/02/11(水) 23:55:08

まあ別にLispが悪いわけじゃーない

201 ：１３２人目の素数さん：2009/02/12(木) 00:06:42

人工知能でエキスパートシステムを作る目論見は瓦解してるけど、
エキスパートシステムはいろんな業界で生きていると思う。
アパレルショップに置いてあるファッションコーディネートシステムとか。

202 ：１３２人目の素数さん：2009/02/12(木) 00:23:07

>>196
いや、崇高な理念というのは、何もわかってない現場のために、使いやすいソフト
開発しようと言う部分なんだが。

203 ：１３２人目の素数さん：2009/02/12(木) 00:34:52

じゃあおれは現場の人間の代わりに自ら問題意識を持って、
自らデータを漁って自ら最適解を見つけ出すソフトを開発するよ。

具体案はナッシングだけど。

204 ：１３２人目の素数さん：2009/02/12(木) 00:37:50

彼のは統計手法の組み合わせだけでなんとかなりそうだが、君のは統計外の知識が必要になりそうね。

205 ：１３２人目の素数さん：2009/02/12(木) 00:38:26

まあ、現場が何も知らなくても使えるソフトを作ると

「何かわかんないけど、これの言うとおりにして、こことここに表示される数値を
写しておけばいいらしいよ」となり

「何か、出てきた数字がまずいみたいなんすけど」「あ、それは写すときに適当に
なおしておけばいいよ。p＝0.02と書いとけば問題ないみたいだ」となったり、

「何か、計算してくれないんすけど」「あ、それはデータに問題があるってこと
だから、それらしい数字を足しておけばいいよ」とかになるんだよな。

現場はいそがしいから、手取り足取りしてくれるソフトつくっても、手順が面倒だと
思われて、逆にすっとばされることがある．．．

206 ：１３２人目の素数さん：2009/02/12(木) 00:38:53

子供を作ってまともな教育を与えれば生体コンピューターになるんでね？

207 ：１３２人目の素数さん：2009/02/12(木) 00:41:31

それってコストパフォーマンス悪そうじゃない？

208 ：１３２人目の素数さん：2009/02/12(木) 00:42:34

>>205
アクションにつながってないから当然無意味。
まさにSPSS使ってる現場だな。

209 ：１３２人目の素数さん：2009/02/12(木) 00:43:53

>>206
それいいアイデアだよね！
で 具体的にどうしようか？

210 ：１３２人目の素数さん：2009/02/12(木) 00:54:48

で、現場は統計を使って何をしたいの？
社内プレゼンでつかいたいなら>>205でいんじゃね？
予測して儲けたい、とかいうなら、手取り足取りソフトに頼るな、っつ〜話しだし。

使う手法が決まってるなら、薄いGUIラッパ書いて終わりじゃん。
そもそも、彼の質問ってなんだったっけ？

211 ：１３２人目の素数さん：2009/02/12(木) 00:55:49

まず子供をつくるところからだな。

212 ：１３２人目の素数さん：2009/02/12(木) 00:58:40

子作りの現場で注意することはありますか？

213 ：１３２人目の素数さん：2009/02/12(木) 01:07:01

男が生まれるか、女が生まれるかの確率を自分の子供に適用しちゃダメってことだ。
男が続いたから、次こそは女の子だと思って頑張るんじゃないぞ！

214 ：１３２人目の素数さん：2009/02/12(木) 01:07:45

>>212
「今日は安全日なの」を信じちゃダメだ。

215 ：１３２人目の素数さん：2009/02/12(木) 01:10:34

>>205
品質管理や生産プラントの収量向上が目的なら、もうちょっとマシになると思う。
システムが吐く出力が「行程 xx の温度を 10℃上げてみてください。」とか、
「データ不足により、意味のある答えはえられませんでした。」とか、
「多重共線性が生じています。彼を呼んで変数選択をさせてください。」とか、
「内部エラーです。彼を呼びつけて叱り飛ばしてください。」とか。
それぐらいなら作れそうな気もする…

>>208
アクションにつながらない作業なんて、まっとうな企業には無い。
マーケティングデータや学術研究データの統計解析は、
その後の営業戦略策定や研究結果の考察に直結する。
それらに SPSS や SAS は思いっきり使われてる。
そして、>>205 みたいなデタラメをさせないために統計専門家がいたりいなかったりする。

216 ：１３２人目の素数さん：2009/02/12(木) 01:38:28

>>205、>>215

>>89とか読むと

>現場の方はＲとかには興味ないんです。
>興味があるのは結局どうすれば改善するの？答えを教えて！
>ということだけなので、Ｒのコードすら見せないつもりです。
>最終結果のみ出力ですね。
>（当然、必要な人のために計算手順、考え方を示したレポートも出力しようとは思っています）
>
>欠損値を含むデータは一式除去します。
>（補完は事実ではないので使いたくないしデータは多数あるので）

データが多数あって、欠損があるのは除去しても問題ない環境らしいので、
データ不足は最初からおきなさそう。

データが多数ある、結果だけとか書いてあるし、>>125が彼だとすると、
製品の品質管理とかじゃなく、株とか景気の予測屋さんで、経済とか財務
データの分析かな。

217 ：１３２人目の素数さん：2009/02/12(木) 01:44:39

彼を主成分分析にかけてみてください＞＜

218 ：１３２人目の素数さん：2009/02/12(木) 01:55:11

彼女もお願いします＞＜

219 ：sage：2009/02/12(木) 02:10:16

彼女を彼で内挿してください＞＜

220 ：１３２人目の素数さん：2009/02/12(木) 15:25:40

>>125は彼じゃないだろ。
株の話をしてる時点で違う。

221 ：１３２人目の素数さん：2009/02/12(木) 22:41:02

>>216
株屋さんなどの金融関係なら、「現場」とか「改善」なんて言わないと思いますよ？
どう考えても、品質管理とかプロセス関連。

>>220
いや、>>125 は間違いなく彼だね。本人しか知り得ない事実↓を知っているから。
>素人が使っても正しく回帰分析を使うことができている。
>しかも意識せずにというのが使用者にストレスがなくていいね。

客観的には、彼の認識している世界における事実、っぽいけど。

222 ：１３２人目の素数さん：2009/02/12(木) 23:25:48

>>221
それは俺でも想像できる。
俺でもそう答える。

223 ：１３２人目の素数さん：2009/02/12(木) 23:29:39

>>221
推測と事実を混同してはぶれる。
事実はなんだ？

224 ：１３２人目の素数さん：2009/02/12(木) 23:47:30

う　ぜ　ぇ

225 ：１３２人目の素数さん：2009/02/13(金) 02:37:31

>>221
まあ、例を出す時に株を出しちゃったけど、やってるのは品質管理とかかもね。
株の現場を知ってるにしては、稚拙な喩えだったしね。

しかし、たまにあらわれて「うぜぇ」って言ってるのは誰なんだろうね？
理想は人間を狭量にするのかな．．．

226 ：１３２人目の素数さん：2009/02/13(金) 07:21:27

だめだこいつ。
仕事できねぇ典型だな。

227 ：１３２人目の素数さん：2009/02/13(金) 07:52:25

>>225
このスレを建設的に進めたい。
個人の悔しさをぶつけるのは他でやってくれないか。

228 ：１３２人目の素数さん：2009/02/14(土) 08:30:21

今北

ものすごくレスが増えてて驚いたが、イタいコが複数人湧いてただけか
>>225と、あとは分かんね

この板にもID欲しいな

229 ：１３２人目の素数さん：2009/02/14(土) 10:27:29

>>228
いまさら書いてる君も痛い子。

230 ：１３２人目の素数さん：2009/02/15(日) 14:28:13

自分は「イタい子」って嫌いじゃないけどな。
志と実力がつりあってないから「イタい」と批判されるわけで、
志に見合った実力を身につけて、見返してやればいいさ。

その過程で、実は批判者の側が正しいことが分かったりもするけれど、
それもまた人生。

231 ：１３２人目の素数さん：2009/02/15(日) 21:39:31

>>225と>>230の文体が似ている件

232 ：１３２人目の素数さん：2009/02/15(日) 22:52:02

>>231

>>225は俺だが、>>231とは別人の件
（これは本当ね）

233 ：１３２人目の素数さん：2009/02/15(日) 22:53:20

間違ってるし．．．

>>225は俺だが、>>230とは別人の件
（これは本当ね）

234 ：１３２人目の素数さん：2009/02/15(日) 23:42:41

必死ｗ

235 ：１３２人目の素数さん：2009/02/16(月) 08:53:57

これらの言動の根本原因は自分より格下だと考える発想からきてるようですね。
このような人に、逆に格上の人を挙げてみてくださいと言えば、有名人しか挙げられないんじゃないかな。

部分的に見れば大人ですら小学生に負けてる知識もある。
勝ち負けの発想は和を生まないし幸せをもたらさない。
どんな相手でも敬っていきたいものです。
みんなで協力して一つのものを作り上げていきたいものです。

236 ：１３２人目の素数さん：2009/02/16(月) 11:07:50

>これらの言動の根本原因は自分より格下だと考える発想からきてるようですね。
Ｒ好きのどこかの教授みたいだな。

237 ：KingGold ◆3waIkAJWrg ：2009/02/16(月) 21:54:14

敬うことも蔑むこともない、しかし悪を除外する。

238 ：１３２人目の素数さん：2009/02/16(月) 22:09:07

なんだそりゃ？

239 ：１３２人目の素数さん：2009/02/17(火) 03:00:17

悪なんかいないけど。

240 ：１３２人目の素数さん：2009/02/18(水) 01:17:44

>>237
King→王様
Gold→一位
どんだけ自分大好きなんだよ。
ネタ仕込みに座布団一枚。

241 ：１３２人目の素数さん：2009/02/18(水) 02:46:10

久しぶりに見たら、RjpwikiのQ&Aで回答しないで文句だけ言ってる回答者が相変わらずウザイです。
あまりのウザサに、俺も含め昔はいた優良回答者もだんだんいなくなってる気がします。
何とかしてください。

242 ：１３２人目の素数さん：2009/02/18(水) 08:45:08

kingも大学で教えているんだし
肩書きもあるんだから
恥ずかしいので出てこないで欲しい
そのうち誰かばれちゃうぞ

243 ：１３２人目の素数さん：2009/02/18(水) 10:31:10

>>241
2チャンネルのノリで書き込んでいるcapetapetaという盆栽好きのニートに限っては、質問者の方に
問題があると思うけど。回答者を馬鹿にするにもほどがある。もちろん、スルーできない大人も悪いが。

244 ：１３２人目の素数さん：2009/02/18(水) 11:08:09

自分で広めているつもりだが
離れていく人も（の方が？）多い事実に
気が付いていないようす
Sもそうだった

245 ：１３２人目の素数さん：2009/02/22(日) 19:30:31

>>13

>macでRを使ってるんですが、
>グラフィックボードに日本語を出力しようとしても
>全部四角になってしまいます。

昨日からRを使い始めたのですが、par(family="Osaka")とかで
対応できますか？
ttp://oku.edu.mie-u.ac.jp/~okumura/stat/first.html

246 ：１３２人目の素数さん：2009/02/23(月) 12:46:46

ここ見れば
ttp://www.okada.jp.org/RWiki/?%BD%E9%B5%E9%A3%D1%A1%F5%A3%C1%20%A5%A2%A1%BC%A5%AB%A5%A4%A5%D6%289%29#n4ae5f25

247 ：１３２人目の素数さん：2009/02/27(金) 20:27:46

hist2d()とpersp()で二変量データの密度ヒストグラムを描かせています。
ここで三つの軸の範囲を自分で設定するにはどうすればよいのでしょうか？
よろしくお願いします。

> h2d <- hist2d(x,y,show=FALSE, same.scale=TRUE, nbins=c(20,30))
> persp( h2d$x, h2d$y, h2d$counts,
+ ticktype="detailed", theta=60, phi=30,
+ expand=0.5, shade=0.5, col="cyan", ltheta=-30)

248 ：１３２人目の素数さん：2009/02/27(金) 21:28:31

>>247
persp()にxlim,ylim,zlimっていう引数があるけど？

249 ：１３２人目の素数さん：2009/02/27(金) 23:02:25

>>248
それでうまくいきました。ありがとうございました。

250 ：１３２人目の素数さん：2009/02/28(土) 07:57:35

昨夜のものです。
hist2dでcountsをヒストグラム表示しましたが、これを頻度にして表示する方法はないでしょうか？
(1, 5, 4, 2, 1, 7)というcountを(0.05, 0.25, 0.2, 0.1, 0.05, 0.35)と頻度に換算して表示したいのです。

251 ：１３２人目の素数さん：2009/02/28(土) 15:34:45

>>250
h2d$counts を h2d$counts/sum(h2d$counts) にするだけだと思いますが。

252 ：１３２人目の素数さん：2009/02/28(土) 15:42:40

一週間ほど前の荒れ具合を見て、このスレにいるのはどうしようもないクズばかりだと感じていたが、
ちょっと自分で考えれば分かる質問をするクズや、それに丁寧に回答する優しいクズもいることも知って感心した

253 ：１３２人目の素数さん：2009/02/28(土) 16:16:22

>>252屑乙

254 ：１３２人目の素数さん：2009/02/28(土) 19:06:00

>>252
よく分かっていれば質問もしてこないだろうし
場所も場所だし。これが現実だろう。

どうしていきたい？
いい解決案でもあるかな？

自分の考えているレベルのスレッドか
ブログでも建てるのも手じゃない？

255 ：１３２人目の素数さん：2009/02/28(土) 19:46:49

>>251
ありがとうございました。その通りで行けましたmOm

256 ：１３２人目の素数さん：2009/02/28(土) 20:29:26

√(256) = 16 才は思春期

257 ：１３２人目の素数さん：2009/03/02(月) 13:45:00

おおお？ R-tips.pdf のページがなくなってる？！
http://cse.naro.affrc.go.jp/takezawa/r-tips.pdf
http://cse.naro.affrc.go.jp/takezawa/r-tips/r.html

東工大のはまだあるみたい。
http://www.is.titech.ac.jp/~shimo/class/doc/r-tips.pdf

どっちのが古いとかってありましたっけ？

258 ：１３２人目の素数さん：2009/03/02(月) 14:52:34

http://cse.naro.affrc.go.jp/takezawa/r-tips/r.html
復活したな

259 ：１３２人目の素数さん：2009/03/02(月) 17:56:23

よかったよかった。
どうも中央農業総合研究センターのほうが新しいみたいね。

260 ：１３２人目の素数さん：2009/03/05(木) 23:01:16

R2.8.1入れて使いはじめたらアンチウイルスソフトが反応。
gtools.dllがR2.8.1のフォルダに入っているという・・・
何で？

261 ：１３２人目の素数さん：2009/03/06(金) 19:51:04

>>260
アンチウイルスソフトは何ですか？
アンチウイルスソフトのサポートページを見ると情報が得られるかもしれません。

262 ：１３２人目の素数さん：2009/03/07(土) 02:33:36

今度これ使い始めるんですが、
単刀直入に難しいですか？

参考までに
高校時代は数学は２Bまでやりました（あまり覚えてないですが）
経済商系統の学部の学生です
統計学は多少わかります

263 ：１３２人目の素数さん：2009/03/07(土) 11:10:54

統計学のこともだけどコンピュータが得意かも問題。
基本はコマンド入力だからね。

264 ：１３２人目の素数さん：2009/03/07(土) 19:33:19

たいしたコマンド入力でもネーヨ

265 ：１３２人目の素数さん：2009/03/07(土) 23:06:21

コンピュータに得意とかいえるほどのコマンド入力って
アセンブリ言語？

266 ：１３２人目の素数さん：2009/03/08(日) 01:07:58

ともに得意な人の発言だね。
文系だと「コマンド入力」という言葉自体を何それ？と思っても不思議じゃないが。

267 ：１３２人目の素数さん：2009/03/08(日) 12:31:03

すみません、R の演算子で ~ ってどういう機能なんでしょう？
r-tips.pdf を見ているのですが、よく分かりません。
何か変数に ~ を付けてみても、
> x <- 1:10
> ~x
~x
となるだけで、何をするものなのか分かりません。

268 ：１３２人目の素数さん：2009/03/08(日) 14:45:59

RでLU分解を行う関数ってありますか？

269 ：１３２人目の素数さん：2009/03/08(日) 19:21:58

>>267
~ って演算子じゃないんじゃないか。
モデル式のとき = の代わりに使うしかできないと思うけど。

270 ：１３２人目の素数さん：2009/03/08(日) 20:55:24

既出かもしれませんが、RGuiでコンソールで日本語入力した時にキャレットの位置がずれるんですが仕様ですか？
例えば、
> グループXの標準偏差<-
と入力した時点でXじゃなくてYの間違いだと気づいたとします。
そうして左矢印キーで赤いキャレットを戻してXのところに合わせてBackSpaceキーで消そうとすると、
「X」ではなくて「プ」が消えます。
ちなみに環境はWindows VistaでRのバージョンは2.8.1で日本語入力システムはATOK2007です。

271 ：１３２人目の素数さん：2009/03/09(月) 01:24:15

>>269
なるほど。= の代わりか。ありがとうございます。

ただ、どうも演算子ではあるようなのです。

http://cran.r-project.org/doc/contrib/manuals-jp/R-lang.jp.v110.pdf
(P9)
~ チルダ、公式で使われ、単項演算子としても二項演算子としても使え
る

また、こういう使い方は、= とはちょっと違うように見えます。
http://cse.naro.affrc.go.jp/takezawa/r-tips.pdf
P122
> plot(~ group, data=sleep)
> plot(~ extra, data=sleep)

272 ：１３２人目の素数さん：2009/03/09(月) 19:45:52

plotにそんな使い方があるとは知らなかった。
> plot(group ~ extra, data=sleep)
> plot(extra ~ group, data=sleep)
も面白い。
モデル式同様関数の引数の中でしか使わないのかな。

273 ：１３２人目の素数さん：2009/03/10(火) 11:30:55

>>267
help('~')

274 ：１３２人目の素数さん：2009/03/10(火) 21:16:37

>>267
~xが何なのか見たければstr(~x)とすればよい
~xがクラスformulaを返していることが分かるでしょう。
~はクラスformulaを返す演算子ということではないでしょうか。

275 ：１３２人目の素数さん：2009/03/17(火) 16:00:25

データをeditで編集したら、
Warning message:
In gsub(perl = TRUE, "(?<!\\\\)((\\\\\\\\)*)\"", "\\1\\\\\"", "f") :
perl=TRUE は UTF-8 ロケールに対してのみ完全実装されています
という警告文が出てきて、反映されるのは、編集直後の表示だけなのですが、どうすれば、その後もデータを反映させれますか？

エラー文からして文字コードがいけないと思うのですが、文字コードはどこで変えればいいのかわかりません。

276 ：１３２人目の素数さん：2009/03/17(火) 18:06:20

>>275
環境とかは秘密ですか？
> x <- data.frame(x=runif(5),y=LETTERS[1:5])
> x
x y
1 0.3677259 A
2 0.9520577 B
3 0.5019692 C
4 0.2780160 D
5 0.5301348 E
> edit(x)
x y
1 0.3677259 日本語
2 0.9520577 問題なし
3 0.5019692 C
4 0.2780160 D
Warning message:
In edit.data.frame(x) : 'y' 中に因子水準を加えました

こちらでは特に問題ありません。あなたの問題が再現する方法を提示した方がよいよ。

277 ：１３２人目の素数さん：2009/03/17(火) 19:53:08

環境はXP ProでRのver. 2.8.1です。

278 ：１３２人目の素数さん：2009/04/04(土) 20:09:02

Rってデータ量５M位から、もう動きが怪しくならないか？

あとはここだけだな。

データマネジメントも苦手か。

279 ：１３２人目の素数さん：2009/04/05(日) 09:49:37

>>278
それってものすごいメモリ消費になってるんじゃ・・・？

280 ：１３２人目の素数さん：2009/04/05(日) 11:37:29

大容量を扱うときって外部プログラムからRのライブラリ使うような感じの方が効率よくないかな

281 ：１３２人目の素数さん：2009/04/05(日) 12:59:31

>>280
それってなんか違うような気がする。
SASからRを呼び出すイメージ？

282 ：１３２人目の素数さん：2009/04/05(日) 13:06:09

まあそれでも良いけど
ｃとかjavaとかね

283 ：１３２人目の素数さん：2009/04/05(日) 13:45:44

>>282
使い勝手悪すぎ！
ならcで全て処理を書くよ！

java･･･全部Rで行くよ。

284 ：１３２人目の素数さん：2009/04/05(日) 15:58:08

でもメモリ足りないんだろ
臨機応変に出来なきゃ
ｃでかけるならR要らないだろ
出来ないから使ってるんじゃないかな
そういうものだよ

285 ：１３２人目の素数さん：2009/04/05(日) 23:33:02

おい、お前ら。結局データ加工は何でやってるんよ？
awkか？perlか？それともSASか？

備え付けのデータセットしか分析していないとかいうオチはないよな？

286 ：１３２人目の素数さん：2009/04/06(月) 11:27:24

>>285
データ加工の前処理は、日常的に使っているsed tr paste cut sort
uniq colrm nkfあたりで十分だけど。そのあとは、すべてRの中で処
理しちゃうな。たまにperlも使うけど。

287 ：１３２人目の素数さん：2009/04/06(月) 22:06:16

>>286
なんとストイックな。でももう少し楽したいとは思わないの？
あと扱うレコード数は１万件ぐらいがMAXでしょ？

288 ：１３２人目の素数さん：2009/04/07(火) 08:38:51

>>287
メモリ次第だろJK

289 ：１３２人目の素数さん：2009/04/07(火) 11:05:04

データの形態によってはRの中でデータ加工をした方が楽な場合がある。
私が日頃使うデータは、そんなタイプ。何重にも階層化されたlist構造
をもつので、sapplyを入れ子にしまくり。これをperlで書いたら大変な
ことになるよ。

290 ：１３２人目の素数さん：2009/04/07(火) 11:55:32

データを取り込む前の加工と取り込んだ後の加工とに
わける必要があり、289の言うのは取り込めた後で、
285の言わんとするところは（おそらく）取り込む前の
ことを指していると見ることができそう。
たしかにperlで書いても出力後の果は、またテキストだし
値の管理もできない。
取り込めた後なら289が言うようにRは強い。

元のデータが煩雑な場合にRで取り込みやすいように加工するには
286の形になるだろう。
データの管理を簡易に任せたい、あるいは、汎用コンピュータからの
日本語のEBCDICや他のコードが入り混じっている、あるいは、極端に
複雑なレイアウトで大容量なら、SASを使用るメリットは大きい。

どの点を言っているのか、どの程度まで必要かによって向き不向きがあるので、
一言では、な〜んとも答えきれないんだろうな〜
誰かデータ準備してくれる人いれば、あるいはDBで管理されているデータだけ
を使用するのでよければ、Rだけで済むかもしれないし。ふ〜〜

291 ：１３２人目の素数さん：2009/04/07(火) 14:21:10

>>290
加工前と後を分けるのではなくて、どこまで加工した段階で、Rに渡すのかという
タイミングの話をしていたのでないの？例えば、めちゃくちゃ複雑な階層構造をもつ
XMLがあったとして、perlなどで十分に加工してからRに取り込むのか、またはXMLを
とりあえず取り込んでからRの中で編集加工するのか、とか。

292 ：１３２人目の素数さん：2009/04/07(火) 14:32:08

文章読んでるとそういう意味でもなさそうな気もするんだけどな。

ちなみに、上記で言った意味ではXMLは複雑な構造までいかない。
きれいな構造だし。
めっちゃめっちゃなのがあるんだよ。

293 ：１３２人目の素数さん：2009/04/07(火) 18:22:39

みなさんRでのコーディングルールというか、そんな真面目なものじゃないけど
気をつけてる点とかってどんなのありますか？

私はt<-12345とか死ぬほどやって頭にきた挙句、自分が使う変数は必ず

　.t <- 12345

とか、関数なら

　my.func <- function(...) ...

とかドット始まりか自分プレフィックス始まりにするようにしたんですが、
なんでもドット始まりとか正直ウザー、です。このあたりのR的慣習とか
ベストプラクティスってどこかで紹介されてたりするんでしょうか？

294 ：１３２人目の素数さん：2009/04/07(火) 21:18:25

そういえばさぁ、シェルスクリプトの中で自分用の変数を
命名するとき、"_"（アンダースコア）を最初につけるって
いう風習をどっかで学んでから、シェルスクリプトを書く
ときにはそうしている。
Rで_は全くお薦めでないが。.（ドット）は隠しファイルみ
たいでやだな。実際、隠し関数（内部関数とも言う）には、
ドットから始まるように命名しているのも多いし。

あとは、何がRの変数名で使えるんだったけ？Namespaceをう
まく使えばよいのかもしれないけど、
Internal name space support functions. Not intended to be called
directly.
って書いてあるからユーザは普通使わないだろうな。

295 ：１３２人目の素数さん：2009/04/07(火) 21:38:49

t にデータを代入しても、関数tの挙動は変わらない件について。
(R version 2.8.1 (2008-12-22) on Fedora 10)

> t(matrix(c(0, 0)))
     [,1] [,2]
[1,]    0    0

> t <- 12345

> t(matrix(c(0, 0)))
     [,1] [,2]
[1,]    0    0

296 ：１３２人目の素数さん：2009/04/07(火) 22:04:59

cでも同様になる。
> c <- 10000

> c(1,2,3,4)
[1] 1 2 3 4

Rの処理系は、文脈が関数呼び出しの場合、ある変数に束縛されたオブジェクト
が実行可能かどうか判断して、実行可能でないなら外のスコープを見に行くよ
うだ。

この環境で、文脈が変数参照なら、こうなる。
> c
[1] 10000
> t
[1] 12345

> c(c, t)
[1] 10000 12345
すげえ、もう訳わかんねえ。

?c や ?t はちゃんと呼ばれるようだ。

297 ：１３２人目の素数さん：2009/04/08(水) 00:38:50

>>295
あれ？とりあえず判り易い例としてtを出す前にt<-12345で試して
これなら問題ないなと確認したつもりだったのにorz

まあ、"t"限定でなく、要はうっかりぶつかったり混乱しそうな名前付けとかを
回避するためのネーミングルールとか、何か慣習的にやってるものありますか？ってことで。変数はぶつからないから気にしなくてもダイジョブ？

298 ：１３２人目の素数さん：2009/04/08(水) 01:15:54

>>294
ドットは慣例があるからたとえば、
is.xxx なら、xxxであるかどうか判定する関数だし
read.xxx なら、xxx形式のデータの読みこみだし
じゃあwrite.xxxは？ plot.xxx? わかるでしょ。直感的に。

299 ：１３２人目の素数さん：2009/04/08(水) 01:25:21

上のは293あてね。
一文字変数は極力避けたほうがいいよ。
そして、先頭ドットは隠れ関数/変数の意味合いがあるよ。

300 ：１３２人目の素数さん：2009/04/08(水) 01:33:42

>>296,>>295
Rではネーム空間と関数空間が別ってことでしょう。
これはlisp系の影響を受けた言語ならあながち不思議ではないし、普通に同じ変数と
同じ関数の名前があるなんてざらなのです。Cの系統から見れば不思議なのですがね。
仕様をみると、Rって記述はCやjavascriptに似ておりますが、中身は実質lispの方言という
印象なんでね。

301 ：１３２人目の素数さん：2009/04/08(水) 01:51:53

ごめん ネーム空間じゃなくて変数空間だわ。
;(defun foo (vector)
; (let ((length (length vector))
; …))

みたいなことを差して 300を書いてます。一つめのlengthは変数だし二つめは
関数だしね。Rで書けば

> foo <- function(vector)
+ {
+ length<-length(vector)
+ …
+ }
こんなんかな。

302 ：１３２人目の素数さん：2009/04/08(水) 02:07:30

>>300-301
お、言語に詳しい人が来たみたいね。
だけどそれは違う。R の名前空間は、変数と関数とで分かれていない。
影響を受けているのはLISPの中でもScheme。

もしCommon Lispみたいに、名前空間が分かれているなら、

> x <- function() print("x")
> x()
[1] "x"
> x <- 1
> x()
エラー： 関数 "x" を見つけることができませんでした

このエラーは起こらないはずでしょ？

変数のスロットが一つしかなくて、function() print("x") が 1で上書きされ
てしまうから、こういうことになる。

303 ：１３２人目の素数さん：2009/04/08(水) 02:16:26

>>302
なるほど、てっきりわかれてるのかとおもったけど違うのか。

304 ：１３２人目の素数さん：2009/04/08(水) 02:35:18

ここから先は、言語オタクのための説明。

>>302 の x はかち合って、c や t がかち合わないのは、c や t がユーザ
定義変数のスコープの外側にある、組み込み変数スコープで定義されているから。

関数と変数で名前空間の分かれた言語(Common Lisp, Rubyなど)は、名前の衝突
は起こりにくいが、変数に関数を代入する手続きが煩雑になり、高階関数を書
きづらい。

関数と変数の名前空間が分かれていない言語(Scheme, JavaScript, Pythonな
ど)は、関数を変数の一種のように扱えて、高階関数が書きやすいが、関数名と
変数名がぶつかってしまう。だから例えばSchemeだと、リストを収める変数に
list という名前は使えず、lis とか書くのが慣例。

つうことで R は、>>296の通り、関数呼び出しの文脈で、ある変数に束縛され
ているオブジェクトが実行可能でなければ、その外のスコープに同名の変数を
探しに行くことにした。

このおかげで、変数に関数を自在にセットする便利さと、組み込みの関数名と
変数名がかち合わない便利さの両方を得たわけ。

305 ：１３２人目の素数さん：2009/04/08(水) 11:11:05

rpy2使ってる人いますか？

306 ：１３２人目の素数さん：2009/04/12(日) 19:01:44

欠損値""を欠損値として読み込むにはどうすればいいでしょうか？
具体的には以下のようなファイルを読み込む場合です。

1,2,3,4,5,
,,8,9,10,
11,12,,,

scan("cliptable",fill=TRUE,na.strings="")
ではうまく読み込めません。
欠損値が詰められます。

307 ：１３２人目の素数さん：2009/04/12(日) 20:10:16

read.csv("cliptable",header=FALSE) やったらあかんの？

308 ：１３２人目の素数さん：2009/04/12(日) 20:32:44

Excelのクリップボードからの読込みでした。

そんなわけで、
scan("cliptable",fill=TRUE,na.strings="",sep="\t")

でいけました。
お騒がせしました。

read.table系は処理コストが高いのであまり好きではありません。

309 ：１３２人目の素数さん：2009/04/12(日) 20:34:27

>>308
できるのはわかってたのね。 ふーん

310 ：１３２人目の素数さん：2009/04/12(日) 22:05:41

行列のべき乗って、適当な関数ないの？

固有値分解しないと駄目？

311 ：１３２人目の素数さん：2009/04/13(月) 00:02:31

ha?

312 ：１３２人目の素数さん：2009/04/13(月) 00:15:12

>>310
最小多項式を求めるとか

313 ：１３２人目の素数さん：2009/04/15(水) 12:59:35

今更だが前スレを貼る。

＝　統計解析フリーソフト　Ｒ　【第２章】　＝
http://science6.2ch.net/test/read.cgi/math/1152449095/

＝　統計解析フリーソフト　Ｒ　＝
http://science4.2ch.net/test/read.cgi/math/1062650510/

314 ：１３２人目の素数さん：2009/04/19(日) 23:04:19

結局、データマネジメントにはどの言語（ツール）を併用するのが良いのかね？

重要な問題だと思うんだけどなぁ。

315 ：１３２人目の素数さん：2009/04/20(月) 00:11:49

R-tips、アマゾンでべた褒めすぎだろ
なんであんな糞わかりにくい本が5つ星なんだよ

316 ：１３２人目の素数さん：2009/04/27(月) 01:51:53

>>315
他にまともな本があるならともかく

317 ：１３２人目の素数さん：2009/05/03(日) 15:04:27

>>314
適当なSQL鯖にRODBCでつなげばいいんじゃね？

318 ：１３２人目の素数さん：2009/05/04(月) 16:08:05

このbugs式はどこが間違っていますか？

model{

for ( i in 1:n ){

x[i] ~ dnorm ( b0 , tau )

}

b0 ~ dnorm (0 , 0.001 )
tau ~ dnorm(0 ,0.001)

}

単純にデータｘがあって、平均と標準偏差を求めたいだけです。

319 ：１３２人目の素数さん：2009/05/04(月) 20:30:25

そこに入れるまでの過程とか、
そう言うもの(別ツール）に頼るかしかないのかってこと
を言っていると思うよ。

320 ：１３２人目の素数さん：2009/05/05(火) 02:20:31

Ｒに、最も当てはまりの良い分布を自動的に選択してくれるパッケージってありますか？

イメージとしてはこんなのです。結構便利で需要もありそうですよね？

x <- rnorm( 100 , 10 , 4 )
result <- FitDistribution( x , type="continuous" )

summary( result )
[[1]]
[1] normal distribution
[2] parameter：mean= 9.92 sd=4.1
[3] AIC：421.3

[[2]]
[1]logistic distribution
[2]parameter：location=9.92 , scale=3
[3] AIC：501.2

321 ：１３２人目の素数さん：2009/05/05(火) 12:10:11

322 ：１３２人目の素数さん：2009/05/07(木) 22:33:39

どなたかブートストラップ法とクロスバリデーションの講習会をご存じないでしょうか？
知っている方がいましたら教えていただけると助かります。

323 ：１３２人目の素数さん：2009/05/08(金) 01:08:47

z <- matrix()
heatmap(z)

で、heatmapを描いてくれるのですが、
縦軸横軸に個々の要素の名前（名前が付いていない場合は数字）
が付きます。この名前、もしくは数字を消すにはどうした良い
ですか？教えていただけると大変助かります。

324 ：１３２人目の素数さん：2009/05/08(金) 13:40:32

>>323
糖鎖とタンパク質との反応をheatmapで描いている人？
heatmapエンター
でheatmapの定義内容が表示されるので、当該部分を編集して
my.heatmap <- function(....
と新コマンドを作成したらどうですか？
あと、?axis を読んでおく。

325 ：１３２人目の素数さん：2009/05/09(土) 10:32:40

テキストデータの統計科学入門　という本に
データ数が少ないときデータの平均値はt分布に
従う、と書いてある。笑わせてもらいました。

326 ：１３２人目の素数さん：2009/05/09(土) 11:53:22

mergeの2変数をkeyにできる関数ってないのかな？
2つの変数の両方が合致したときのみデータフレームを連結してくれるような。

327 ：１３２人目の素数さん：2009/05/17(日) 12:00:22

統計ソフトRを用いて、CUTを使い例えば体重の最大値と最小値を求める場合、どうすれば…？

328 ：１３２人目の素数さん：2009/05/17(日) 12:11:27

最大値最小値って簡単にまとめて出てこなかったっけ？
平均とかと一緒に

329 ：１３２人目の素数さん：2009/05/17(日) 12:19:08

cut使い、最大値と最小値求める場合は簡単に出来る？

330 ：１３２人目の素数さん：2009/05/17(日) 13:35:15

>>327=>>328=>>329

あまりに素人質問過ぎて同一人物なのがミエミエ。

最大値を求めるのはmax、最小値を求めるのがmin、
summaryやfivenumでその他の統計量と同時に出力。

cutは区分値を求めるものであり最大値・最小値を求めるものではない。

331 ：328：2009/05/17(日) 14:02:12

おれも素人だけど別人です。

332 ：１３２人目の素数さん：2009/05/17(日) 14:33:13

糞でる

333 ：１３２人目の素数さん：2009/05/17(日) 17:50:34

>>332
してくれば？

334 ：１３２人目の素数さん：2009/05/17(日) 18:18:22

糞したら負けかなと思ってる

335 ：１３２人目の素数さん：2009/05/17(日) 22:30:28

じゃ、栓をしときなさい。

336 ：１３２人目の素数さん：2009/05/18(月) 01:50:05

>>327、328、329
???
どーいうこと?

…体重でもなんでもいいけど、例えば四分位数求めたい時にcut使う(?)(間違ってたらスマソ)

337 ：１３２人目の素数さん：2009/05/18(月) 09:29:37

>>327
>>328
>>329
素人かどうか、同一人物かどうかはどうでもいいけど、取りあえずRをインストールして
みたらどうだろう。少なくともRを使ったことのある人の書き込みとは到底思えない。
>>327の質問の内容・仕方では、>>330の答え以上は期待できないよ。
fivenum()なんて知らんかった。

338 ：１３２人目の素数さん：2009/05/18(月) 16:43:39

ちょっとお尋ねします。

Rのround関数はJIS方式のため、例えば3.65を丸めれば3.7では無く3.6になる、という所は理解できたのですが、
直前に3.65を求める計算によって結果が変わるケースに出くわしました。

> round(3.65,1)
[1] 3.6　　#想定通り

> b<-(4.8+2.5)/2
> b
[1] 3.65
> round(b,1)
[1] 3.6 　　#ここは想定通り

> a<-(4.9+2.4)/2
> a
[1] 3.65
> round(a,1)
[1] 3.7　　#なぜ？

Rjp wikiのJIS式四捨五入ページなど調べてみましたが、どうにも分かりません。
ご教授いただけますと幸いです。

Rのバージョンは2.81で、Windows MCE上で使用しています。

339 ：１３２人目の素数さん：2009/05/18(月) 16:51:50

S=b/a∫(acosθ→a)√(a^2-x^2)dx+acosθbsinθ/2
この問題の答えがabθ/2になるのですが、途中計算がわからないので教えてくれませんか？
∫の後ろの括弧は積分範囲です

340 ：１３２人目の素数さん：2009/05/18(月) 16:52:43

すいません誤爆しました

341 ：338：2009/05/18(月) 16:56:49

追記です

Excel VBAのRound関数もJIS式だった事を思い出したので、念のため試してみましたが、
いずれも3.6になりました。

Cells(1, 1) = Round(3.65, 1)

b = (4.8 + 2.5) / 2
Cells(2, 1) = Round(b, 1)

a = (4.9 + 2.4) / 2
Cells(3, 1) = Round(a, 1)

うーん、何なんでしょ。

342 ：１３２人目の素数さん：2009/05/18(月) 17:48:28

ヒント
> (4.8+2.5)/2 - (4.9+2.4)/2
[1] -4.440892099e-16

343 ：338：2009/05/18(月) 17:52:45

>>342
ありがとうございます。浮動小数点演算誤差なんですかね？
ズルっぽい雰囲気も漂いますが、2段階で丸めるようにします。

344 ：388：2009/05/19(火) 10:26:48

結局、2段階丸めでは近似値の問題
ttp://www.okada.jp.org/RWiki/?%A3%D1%A1%F5%A3%C1%20%28%BE%E5%B5%E9%BC%D4%A5%B3%A1%BC%A5%B9%29#v8aea385
を解決できなかったため、

ttp://hp.vector.co.jp/authors/VA008683/QA_JISROUND.htm
を移植して

round10 <- function(x, n) {
　　x<-x*10^n
　　if((x %% 1) < 0.5){
　　x<-trunc(x)
　}else if((x %% 1) > 0.5){
　　x<-ceiling(x)
　}else if(trunc(x) %% 2 == 0){
　　x<-trunc(x)
　}else{
　　x<-ceiling(x)
　}
　　y<-x/10^n
　return(y)
}

としました。お騒がせしました。

345 ：１３２人目の素数さん：2009/05/19(火) 18:28:12

今日初めてWin版のインストールして、ヘルプ見ながら操作しています。
ちょっとお聞きしたいのですが、クリップボードにあるcsvもしくはタブ区切りのファイルを
データフレームとして取り込みたいと思ったらどうするのが一番手っ取り早いでしょうか？

346 ：１３２人目の素数さん：2009/05/20(水) 11:14:14

>>345
Windowsユーザではないが、放置されてかわいそうなので、
> RSiteSearch("clipboard")
であなたの代わりに調べてあげたよ。

ttp://finzi.psych.upenn.edu/R/library/psych/html/read.clipboard.html

347 ：１３２人目の素数さん：2009/05/21(木) 09:27:21

ありがとうございます、できました。
まずヘルプ回りの関数類の使い方を憶える事と、パッケージの概念を理解するところからですね、ややこしいです……

348 ：１３２人目の素数さん：2009/05/26(火) 00:17:17

統計用語の英語は全然わからない、ヘルプまともに読めない、言語もトリッキーな仕様
なんつーか見る者を絶望させる言語ですね……
ヘルプの日本語訳(古いバージョンのでもいい)があるといいんですけど……

349 ：１３２人目の素数さん：2009/05/26(火) 01:03:24

真の絶望は、エラーが出た際にエラー発生行数すら教えてくれない事じゃないかとw

350 ：１３２人目の素数さん：2009/05/26(火) 09:56:02

>>349
なにぃぃぃぃ、そいえばステップ実行はしたけれど、そこにはそれらしい情報なかったな……

351 ：１３２人目の素数さん：2009/05/26(火) 09:58:39

一週間にして心が折れそうだ(ToT)

352 ：１３２人目の素数さん：2009/05/26(火) 12:37:13

>>351
まぁなんだ。頑張れ！
身近に教えてくれる人はいないの？

353 ：１３２人目の素数さん：2009/05/27(水) 12:30:26

それって一行一行デバッグしろってことか

354 ：１３２人目の素数さん：2009/05/27(水) 12:57:03

Rは対話的な統計解析環境で、スクリプトを書いて一気に処理することもできる
というだけの話。コンピュータ開発言語と同列に考えたら駄目だよ。

>>349
エラーが出たときにエラーが出た行で止まってくれるから、発生行番号が
不明でも一目瞭然なんだけど。わざわざoptions()でechoを止めているの？

>>353
意味が分からん。Rの問題ではなくて、開発環境（テキストエディタなど）
の問題じゃないのか？C-c C-cで文脈を判断して、パラグラフなり関数定義なりをまとめて実行してくれるから、まとめてデバッグできるけど。、

355 ：１３２人目の素数さん：2009/05/28(木) 00:04:03

C-cってなに？

356 ：１３２人目の素数さん：2009/05/28(木) 11:10:30

>>354
>エラーが出たときにエラーが出た行で止まってくれる

そうか、エラーが出た時はコマンドプロンプトにスクリプトをコピペして
動かしてみれば良かったんだな。

357 ：328：2009/05/28(木) 12:34:10

大文字CはCtrl
C-cはCtrl+cって言う意味。
MはAltかEsc

358 ：１３２人目の素数さん：2009/05/28(木) 12:35:21

>>356
エディタから、いちいちコマンドプロンプトにコピペしていたら、疲れない？
文明の利器を使いなよ。エディタから、カーソル行、または選択複数行をRに
送って実行する機能があるだろ。

359 ：１３２人目の素数さん：2009/05/28(木) 12:36:01

>355
ごめんほかのスレの質問と勘違いした。
xyzzyスレだと思った。

360 ：354：2009/05/28(木) 12:38:30

>>359
C-cの説明はそれであってるよ。

xyzzyからRを使っているの？

361 ：１３２人目の素数さん：2009/05/28(木) 13:14:11

xyzzyはいつもエディタに使ってる。
Rを一度使ったことあったからそのときに2chブラウザのこのスレをお気にスレに入ってただけ。
てっきりxyzzyの初心者の質問かと思って。

362 ：１３２人目の素数さん：2009/05/28(木) 13:33:08

>>356
例は良くないけど、
基本は、昔の汎用コンピュータでの
ラインバイラインでのインタプリタと
同様と考えるといい。

363 ：１３２人目の素数さん：2009/06/02(火) 00:39:32

最近Rを使い始めたものです
データの保存について質問です
ワークスペースを消す際に作業スペースを保存しますか？
という質問が出ますがここではいを選択すればこれまでの作業内容が保存されるのでしょうか？
前回この操作をしたにも関わらず新しくRを立ち上げても
以前にセーブされたワークスペースを復帰します
という文が出てくるだけで以前の作業内容が反映されません
作業スペースの読み込みコマンドからR Workspace (.RData)を選択して開こうとしても
load("C:\\Users\\kobayashi\\Documents\\.RData")
の一行が返されるだけです
他にどのような操作をすればよいのでしょうか？

364 ：１３２人目の素数さん：2009/06/02(火) 10:10:03

>>363
初心者ならではの初々しい質問だな。

.RDataは現在Rが作業しているフォルダに保存されるので、違う場所（フォルダ）で
Rを起動すればそのときに自動的に読み込まれる.RDataは別物になる。現在地（カレン
トディレクトリ）とは別の場所の.RDataを読み込みたいときには、適宜自分でload()を使う。

loadやsage.imageのヘルプを読もう。

> 以前の作業内容が反映されません
library()の結果は保存されない。保存されるのはRオブジェクト（関数、リスト、
データフレーム、行列など）のみ。
ls()とすると、前回保存したRオブジェクトが今もあるのかどうか確認できる。

最初はいろいろととまどうこともあるだろうけど、頑張れ小林君w

365 ：363：2009/06/03(水) 11:36:14

364さんありがとうございます
ls()で確認したところ以前保存した内容がオブジェクトとして保存されていることが
分かりました
ところでこのオブジェクトをワークスペース上に呼び出すにはどのように
したらよいのでしょうか？
ヘルプを読んだ限りでは
load("all.Rdata", .GlobalEnv)
のように（ファイル名は自分のものに変えています）すれば出てくるのかと思ったのですが
以下にエラー load("普及率.Rdata", .○○○) : オブジェクト '.○○○' がありません
のようなメッセージが出てしましました
他にload(○○○）load("○○○")などもためしてみましたがうまく行きません
どこが間違っているのでしょうか？
何度もすいません

366 ：363：2009/06/03(水) 11:39:13

自己解決しました
そのまま名前を入れれば良かったのですね
364さん本当にありがとうございました

367 ：１３２人目の素数さん：2009/06/03(水) 12:28:38

>>366
どういたしまして。

数時間とか作成に時間がかかるRオブジェクトは保存した方がよいけど、
外部ファイルを読み込んでから数秒、数十秒でできあがるRオブジェク
トはいちいち保存しないで、Rスクリプトの方を保存した方がメンテナ
ンス性が高いと思う。

Rのインターフェイスと直接的に対話的操作をするのではなく、テキスト
エディタを通じて、Rスクリプトの現在行、選択行、パラグラフをRに送
って実行する形態で作業した方が効率が上がると思うよ。Rと通信できる
テキストエディタはいくつかある。Windows版を使っていないので、いい
加減なことを言っているかもしれないけど、最近のWindows版Rは、Rとや
りとりする内蔵エディタが組み込まれていなかったけ？

368 ：１３２人目の素数さん：2009/06/04(木) 01:43:10

こんばんは。最近Rを使い始めた新参者です。自己解決しないので質問させていただきます。

ある動物に取り付け、gpsから得られたn×８列のデータ
Time Date Longitude Latitude Altitude GES speed Ang
13:39:06.25 21/09/2007 141.980460 39.308037 -7.4 35.6 1.44 333.8
13:39:06.75 21/09/2007 141.980458 39.308032 -6.8 35.6 0.21 336.5
13:39:07.25 21/09/2007 141.980460 39.308031 -6.1 35.6 0.35 339.0
・
・
をx=read.table("C:x.txt",header=T)
で取り込むことはできるのですが、各列の型が
factor,factor,numeric,numeric,numeric,numeric,numeric,numeric
となり、Timeの部分をnumericで取り込む事ができません。

colClassesを用いて、
x=read.table("C:x.txt",header=T,colClasses=c("numeric",
+ "factor","numeric","numeric","numeric","numeric","numeric","numeric"))
にて問題解決を試みるも、
「以下にエラー scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
scan 関数は 'a real' を期待したのに、得られたのは '13:39:06.25' でした」
というエラー文が出てしまいます。

Timeの"："が問題だと思っています。Timeの列をnumericに取り込む事はできますでしょうか？

もしできないのであれば、
ファイルxの第一列を取り出す→":"毎に分ける→13時以降の端数を小数点に直す→TimeとDate以降のデータをくっつける
といった作業が必要なのでしょうか？それはあまりに手間がかかる気がしたので質問させて頂きました。

369 ：１３２人目の素数さん：2009/06/04(木) 03:04:35

S-PLUSの質問はどこですればいいんだ？

370 ：１３２人目の素数さん：2009/06/04(木) 05:49:33

>>369
http://academy6.2ch.net/test/read.cgi/psycho/1012801769/

371 ：１３２人目の素数さん：2009/06/04(木) 10:44:49

>>368
読み込み時／後に、POSIXlt をつかえばいんでね？ help(ISOdate)

> old.digits <- getOption("digits.secs")
> options(digits.secs = 3)
> dateObj <- strptime("13:39:06.25 21/09/2007", "%H:%M:%OS %d/%m/%Y")
> dateObj
[1] "2007-09-21 13:39:06.25"

これで時間オブジェクトになるから後は煮るなり焼くなり。最後にオプション戻す
>options(digits.secs = old.digits)

372 ：１３２人目の素数さん：2009/06/04(木) 11:17:06

>>368
> x <- read.table("test.dat",header=T,colClasses=c("character","character","numeric","numeric","numeric","numeric","numeric","numeric"))
> x
         Time       Date Longitude Latitude Altitude  GES speed   Ang
1 13:39:06.25 21/09/2007  141.9805 39.30804     -7.4 35.6  1.44 333.8
2 13:39:06.75 21/09/2007  141.9805 39.30803     -6.8 35.6  0.21 336.5
3 13:39:07.25 21/09/2007  141.9805 39.30803     -6.1 35.6  0.35 339.0
> options("digits.secs"=2)
> strptime(paste(x$Time,x$Date),"%H:%M:%OS %d/%m/%Y")
[1] "2007-09-21 13:39:06.25" "2007-09-21 13:39:06.75"
[3] "2007-09-21 13:39:07.25"
> cbind(DateTime=strptime(paste(x$Time,x$Date),"%H:%M:%OS %d/%m/%Y"),x[,-(1:2)])
                DateTime Longitude Latitude Altitude  GES speed   Ang
1 2007-09-21 13:39:06.25  141.9805 39.30804     -7.4 35.6  1.44 333.8
2 2007-09-21 13:39:06.75  141.9805 39.30803     -6.8 35.6  0.21 336.5
3 2007-09-21 13:39:07.25  141.9805 39.30803     -6.1 35.6  0.35 339.0
Rにおける日付型の扱いを勉強しよう

373 ：１３２人目の素数さん：2009/06/20(土) 00:13:21

Ｒでパラメータ設計とかMTS法(Mahalanobis Taguchi System)とか
いう情報がネットにないのですが、
なぜなんでしょう？

特にMTS法はmahalanobis()関数
使って何とかならないものでしょうか？

出来る方やり方教えてください。(´・ω・｀)

374 ：１３２人目の素数さん：2009/06/23(火) 23:18:50

こんばんは、もしよければ教えて下さい。

データﾌﾚｰﾑ上の変数の値を、別途定義したベクトルの数値に
置き換えたいのですがどうにもうまくいきません。

for文でまわしてコピーするしかないのでしょうか？

375 ：１３２人目の素数さん：2009/06/24(水) 01:36:50

>>374
x[ , 3] <- y
みたいではできない？

376 ：１３２人目の素数さん：2009/06/24(水) 04:07:19

>>374
質問が抽象すぎて答えにくい。
> (dat <- data.frame(x=runif(5),y=LETTERS[1:5],stringsAsFactors=FALSE))
          x y
1 0.9966249 A
2 0.3132333 B
3 0.6491857 C
4 0.7796462 D
5 0.3020311 E
> (z <- month.abb[1:5])
[1] "Jan" "Feb" "Mar" "Apr" "May"
> dat[c(1,3,4),2] <- z[c(1,3,4)]
> dat
          x   y
1 0.9966249 Jan
2 0.3132333   B
3 0.6491857 Mar
4 0.7796462 Apr
5 0.3020311   E
こういうこと？

377 ：374：2009/06/24(水) 20:58:37

回答thx

自作関数で計算したベクトルをデータフレームに代入しようと
してたんだけど、返値のベクトルをnumeric(0)で初期化してたのが
まずかったみたい。

numeric(0)での初期化をやめたら375さん、376さんの方式で
できました。

ありがとー。

378 ：１３２人目の素数さん：2009/06/25(木) 17:35:08

Rで、カッパ係数出すために色々数値入力して、library(psy)を読み込ませようとしたんだが、“「psy」というパッケージは存在しません”と表示された。
lkappaでカッパ係数簡単に計算したいんだが、library(psy)が読み込めないから、てで計算するしかないんですか？

379 ：１３２人目の素数さん：2009/06/25(木) 18:32:40

>>378
psyパッケージをインストールしていないだけでは？

> RSiteSearch("kappa")
でカッパ係数を計算する関数がざくざくでてくる。あちこちのパッケージ作者が
自分で独自に実装しているみたいで、車輪の再発明の典型だな。

380 ：１３２人目の素数さん：2009/06/25(木) 18:49:04

>>379
サンキュです！

381 ：１３２人目の素数さん：2009/06/26(金) 17:51:55

378さんに関連してですが、僕もpsy(？)を読み込ませたらパッケージが存在しないと表示されました。

それで379さんの書いたものを読み込ませサイトに接続したのですが、どうやってパッケージをインストールするのかわかりません�ｫ

382 ：378：2009/06/26(金) 18:17:50

>>381
試しにやってみたら、普通にインストールできるけど。
好みの問題かもしれないけど、管理者権限で実行した方が無難かも。
っていうか、ちょっと調べたり本を読めばすぐに分かるレベルの質問をしても相手にする人は少ないと思うよ。
スルーできていないのは私だけみたいだし。

> install.packages("psy",repos='http://cran.md.tsukuba.ac.jp/')
install.packages("psy", repos = "http://cran.md.tsukuba.ac.jp/") 中で警告がありました：
引数 'lib' が欠けています：/usr/local/lib/R/site-library を使います
URL 'http://cran.md.tsukuba.ac.jp/src/contrib/psy_0.7.tar.gz' を試しています
Content type 'application/x-gzip' length 27221 bytes (26 Kb)
開かれた URL
==================================================
downloaded 26 Kb

* Installing *source* package ‘psy’ ...

[中略]

** building package indices ...
* DONE (psy)

383 ：379：2009/06/26(金) 18:21:28

すまん。私は378じゃなくて379だった。

384 ：１３２人目の素数さん：2009/06/30(火) 22:44:25

こんなイメージのパッケージ（関数）ありませんか？

使用例：　（連続変数である、一変量確率変数xの、AIC基準で最も当てはまりの良い分布を検索。
Fitting.Distribution( x , "continuouts" , "AIC")

アウトプット例：

RANK：１
NAME：Normal Distribution
parameter：(mean=3.1 , sd=1.2 )
AIC：324

RANK：２
NAME：Log Normal Distribution
parameter：(shape1=4.2 , shape2=2.3)
AIC：411
・・・・


ニーズありますよね？自作でも良いのでもしありましたらご紹介下さい。