統計解析フリーソフト Ｒ 【第３章】
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2008/10/16(木) 16:33:16統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
0002１３２人目の素数さん
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
0003１３２人目の素数さん
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| ／ ´ ［二＞,,,○|...::::''''
／: :ー ´:ヽ {: :［ﾆ／ヽ:::'''、
／: : : : : : ／￣｀ー ＿＿_ノ／ ／: : :＼
0004１３２人目の素数さん
2008/10/16(木) 23:17:210005１３２人目の素数さん
2008/10/17(金) 22:13:51今、
unko0.dat
unko1.dat
unko2.dat
...という風な名前のファイルが大量にありそれらに似たような処理をＲでしたいとします。
そのため、以下のような関数を作りたいです。
operation <- function(n){
dataframe <- read.table("unko$n.dat")
#データフレームにやりたい処理を書く
}
ここで"unko$n.dat"のように文字列の中で変数を展開させたいのですが
そのようなことは可能でしょうか？よろしくお願いいたします。
0006１３２人目の素数さん
2008/10/18(土) 02:53:02str <- paste("unko",n,".dat",sep="")
でできました。ｽﾚ汚しｽﾏｿ。
0007１３２人目の素数さん
2008/10/19(日) 01:25:33パッケージの中身のソースはどうやってみるのでしょうか。
0008１３２人目の素数さん
2008/10/19(日) 12:52:51そのパッケージのソースをCRANからダウンロードして中をみる以外の方法を聞いているの？
0009１３２人目の素数さん
2008/10/23(木) 23:18:232000×2500のdata.frameの内容をファイルに出すのですが、
時間がかかってしかたありません。
いちいち、sprintfで文字列にして、catで出力するのですが、何か良い方法があればとおもい、
ヒント頂ければと思ってます。
sprintfのフォーマットは"%12.6f".
500万件の出力には、時間必要でしょうかね。
もちろん、ファイルはオープン＆クロースの手順を踏んで、
一つ一つはIO発生させないようにしてるんですが…
テキストファイルへの出力のヒント、頂ければと思います。
0010１３２人目の素数さん
2008/10/24(金) 21:28:35write()でパイプ経由でawkとかperlに渡す場合とどっちが速い？
0011１３２人目の素数さん
2008/10/28(火) 17:33:32思ってるんだが、やっぱRのプログラミングコードとかが
たくさん書いてある本のほうが良いかな。
それとも統計学とRの両方について詳しく説明してる本のほうが良い？
0012１３２人目の素数さん
2008/10/29(水) 14:44:29コーディング経験の有無、統計に関する知識、その辺がすっぽり抜けとる。
0013１３２人目の素数さん
2008/10/29(水) 21:34:29グラフィックボードに日本語を出力しようとしても
全部四角になってしまいます。
解決方法ご存知の方がいらっしゃったら教えてください。
0014１３２人目の素数さん
2008/10/30(木) 11:37:13基礎的なことはネットで資料を探す。
そうすると、Rの基礎みたいな本はほとんどいらないことがわかる。
というかRの使い方という感じで売ってる本だとRのマニュアルのコピー
みたいなので終わってるのが多い。
その上で、Rの使い方＆統計処理のつながりで、お気に入りをあげておくと
「Rによる医療統計学」丸善株式会社
内容は医療統計ではない。タイトルで損をしてるかもしれないが、
このタイトルだから買う人もいるかな。
内容的には、Rの使い方と分析方法がバランスよく書いてあり、使いやすい本。
Rの使い方、確率とか記述統計の話があった後に、１標本、２標本の検定、
分散分析、単・重回帰、ロジスティックモデル、生存分析が載ってる。
「RとS-PLUSによる多変量解析」シュプリンガー
いろいろな分析手法について、Rでどうやればいいかをメインに書いてある。
理論的な話はコラム扱いになっているので、そこは飛ばして読むことも可能。
主成分分析、因子分析、対応分析、クラスター分析、分散分析、重回帰・正準
相関分析、反復測定データの分析（パネルのランダムモデル）、あとはRの
使い方の備忘録がまとめてある。
「統計学：Rを用いた入門書」共立出版
これは現在目を通している所だけど、統計学の基礎をRで経験的に学ぶには
いい本かな。Rの使い方の理解にもいいでしょう。
ブートストラップなどの方法が載ってるのがPCを前提するといいかな。
0015１３２人目の素数さん
2008/10/30(木) 16:58:52R.app だとQuartzGraphicDeviseが日本語に対応していない？ようです。
ですが、ちゃんと調べた訳ではないので、回避方法があるかもしれません。
また、グラフィックデバイスを使わずに直接ファイルに書き込む場合には、
日本語が書き込まれます。
また、Java で動く JGR というGUI環境では、過去にグラフィックデバイス上で
日本語が表示できたとの報告があります。詳しくは、
ttp://www.okada.jp.org/RWiki/?JGR%A4%C7%C6%FC%CB%DC%B8%EC
0016１３２人目の素数さん
2008/11/02(日) 13:53:00Tigar+R 2.7.2だけど別に普通にquartzに日本語が表示されるよ
0017１３２人目の素数さん
2008/11/18(火) 17:45:370018１３２人目の素数さん
2008/11/20(木) 02:03:050019１３２人目の素数さん
2008/11/20(木) 18:31:42にCからRを呼ぶ例がみつからないのだけど
どこかにないのでしょうか
0020１３２人目の素数さん
2008/11/21(金) 01:04:01できる
>>19
R以外の何かはそんなにデータを取り込む手順が面倒でないの？
002118
2008/11/21(金) 05:57:35002218
2008/11/21(金) 05:58:27>>20です
0023１３２人目の素数さん
2008/11/21(金) 22:57:23そんなの自分で調べろ。
> 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'))
0024１３２人目の素数さん
2008/11/22(土) 09:20:38Rcmdrパッケージを導入。
RcommanderからだとGUIでエクセルをファイルを選べる。
0025１３２人目の素数さん
2008/11/22(土) 18:30:190026１３２人目の素数さん
2008/11/26(水) 19:30:170027１３２人目の素数さん
2008/11/29(土) 20:42:49末尾にある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
0028１３２人目の素数さん
2008/11/30(日) 01:56:08過去ログが見られないので質問します。
文字列”123,456,789.123”を数値に変換するにはどうすれば良いのでしょうか？
as.doubleではエラーとなりました。
0029１３２人目の素数さん
2008/11/30(日) 11:35:09専用の関数があった気がするけど、探すよりもカンマを取って変換する方がはやいかも
> 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
0030１３２人目の素数さん
2008/12/04(木) 13:39:40出来ないのでしょうか？
Excel 使わずに、 Calc で済ましているので、
CSVファイルに変換して読み込ませています。
0031１３２人目の素数さん
2008/12/07(日) 22:41:330032１３２人目の素数さん
2008/12/07(日) 23:46:03Rエディタって何？
0033１３２人目の素数さん
2008/12/08(月) 23:20:19で出てくるエディタです。
0034１３２人目の素数さん
2008/12/10(水) 00:41:21edit()って環境変数で指定しているエディタが起動するんじゃなかったっけ？
試しに、今やったら、Emacsが起動した。別のマシンでやったらviだった。
Emacsの色を変えたいなら、手っ取り早いのがcolor-theme。
0035１３２人目の素数さん
2008/12/10(水) 13:26:19http://www.tokyo-tosho.co.jp/books/ISBN978-4-489-02045-2.html
0036１３２人目の素数さん
2008/12/10(水) 20:36:05統計に向いてないのでは？
0037１３２人目の素数さん
2008/12/10(水) 23:40:56逆だと思う。2次元にこだわるのは集計向きだけど、統計向きじゃない。
深い階層構造をもつデータや3次元以上の配列データをどうやって表にするの？
0039１３２人目の素数さん
2008/12/12(金) 01:06:430040１３２人目の素数さん
2008/12/12(金) 19:59:01あおりではないよ。
統計をやる上では、２次元以上のデータを扱わないといけないのが普通。
だから、２次元でないと使いづらいと感じるならば統計処理がわかってない。
それなのに、Ｒの仕様の方に文句つけるってことは、>>36が統計に向いてない
ってこと。
なので、>>36は自分を反省したほうがいいよって内容のコメントです。
あ、>>39が煽りだったか、反応しちゃった。ｗ
0041１３２人目の素数さん
2008/12/12(金) 23:44:40しみじみ時の流れを感じる今日この頃。
0042１３２人目の素数さん
2008/12/13(土) 11:54:59RjpWikiをみると、昔から使ってる人のほうが態度が悪いのが多いと思うが。
昔はよかったってのは嘘じゃないかな？
0043１３２人目の素数さん
2008/12/13(土) 13:13:41とはいつの時代も言われること
0044１３２人目の素数さん
2008/12/13(土) 20:39:22Rjpwikiの始発メンバーであんな回答態度をする人はいなかったけど。
途中参加でとっても辛口な先生1名だけいたけど、全体的にはとっても
雰囲気がよかった。最近の回答態度は酷すぎる。
0045１３２人目の素数さん
2008/12/14(日) 21:43:59「俺達は金もらって仕事で答えてるんじゃない」という意識が強すぎるし、
質問者に対して「学術的態度で臨むこと」を求めすぎ。
ゼミで教授が学部学生をこき下ろすみたいな調子だもんね。
あそこは学者という人種に共通してありがちな人間的欠陥が
よく観察できるコミュニティだと思う。その点では興味深い。
Ｒ自体はとてもいいソフトだと思うけど、あの集団には近寄れない。
結果として自分１人、孤立無援で使うしかない人が大勢生まれてるはず。
004644
2008/12/14(日) 23:16:21人がそれぞれ違うように学者もそれぞれ違うということを理解して欲し
い。それから、今のRjpwikiのノリは、技術系コミュニティに普遍的に見
られるノリをそのまま持ち込んでいるのであって、学者や学術とは関係
がないと思う。むしろ、プロのプログラマがボランティアのコミュニティ
に回答をお願いする態度を求められるような。
「技術系メーリングリスト」で検索すれば、私の言わんとしているとこ
ろが理解してもらえるかもしれない。
0047１３２人目の素数さん
2008/12/15(月) 00:01:57１０年ぐらい前のニュースグループ fj の Java で、学者の罵りあいをはじめて見た。
煽(あお)り、嘲(あざけり)りが、しょっちゅうあったよ。
精神年齢が低いのか、世間と接する機会がないのか、異常な空間だった。
気持ち悪いから、次第に見なくなった。
0048１３２人目の素数さん
2008/12/15(月) 02:48:580049１３２人目の素数さん
2008/12/15(月) 03:35:23demo(persp)
２ばんめのやつ目盛り入ってるけど
005048
2008/12/16(火) 08:21:58ほんとだ。ググってばっかでdemo見てなかった
ありがとうございます
0051１３２人目の素数さん
2008/12/25(木) 16:28:17ランダムに一つの値を返す関数f()を使って、
n個のf()の要素を持ったベクトルを生成したいんですが
どうやってやればいいんでしょうか？
例えばrunif(1)を使って、
n個のrunif(1)を要素に持ったベクトルを作りたいんです。
（この場合はrunif(n)とすれば同じことなんですが。）
runifの定義を見たかったんですがうまく見れなかったんで・・・、
質問がわかりにくくてすいません
0052１３２人目の素数さん
2008/12/25(木) 17:35:41for文で繰り返してf()を呼ぶんじゃあだめなわけ？
0054１３２人目の素数さん
2008/12/25(木) 22:08:06改変した方が早い気がするし、後々便利だと思うのだが…
0055１３２人目の素数さん
2008/12/25(木) 23:37:08こういうこと？
> 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
005651
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でできるんですね。わかりませんでした。
ありがとうございます！
0057１３２人目の素数さん
2008/12/26(金) 21:00:54あなたのやってるのは、ｘを上書きしてるので、ベクトルにならないの。
足すってのは、ベクトルを後ろにのばしていけってことね。
ということでベクトルの追加の仕方をわかってないようなので、
すごく格好悪く書くと
x <- runif(1)
for (i in 1:9) x <- c(x,runif(1))
最初ので、xに値をいれ、次に、その後ろに値を9個追加してる。
大切なのは、x自体に値を加えたものをxに入れることで、ベクトルを
のばしてるってこと。
005851
2008/12/27(土) 01:36:07あ、なるほど、cで繋げばよかったんですね。
ただ、sapplyに比べてかなり効率悪いみたいです・・・(_ _ ；)
自分はかなり大きいデータを扱おうとしてたので
もともとforだと無理だったかもしれません
0059１３２人目の素数さん
2009/01/07(水) 18:50:15GNU R,subset(),&&演算子の使用で結果が変に - gnarl、技術メモ’”%s<marquee>￥
http://d.hatena.ne.jp/gnarl/20090107/1231316900
0060１３２人目の素数さん
2009/01/07(水) 19:51:57勉強半分遊び半分でいじってるんですけどうまく動作しなくて・・・
0061１３２人目の素数さん
2009/01/08(木) 18:39:38乱数で生成したデータでクラスタ分析等をしても実感がわきません。
そこで、実際の実験や調査で得られたデータを使ってみたいのですが、
そういうものはどこかで手に入るものなのでしょうか？
（Rに付属してたりしないんですかね）
0062１３２人目の素数さん
2009/01/08(木) 19:28:14R には有名な Fisher の iris データが付属しています。
あと、省庁が Web で公開している統計データ（都道府県別○○データとか）にも面白いものがあると思います。
006361
2009/01/08(木) 21:31:59data(iris)で呼び出せるとは知りませんでした。
農水省やら環境省を覗いてみましたが統計データがいっぱいありました。
いろいろいじってみようと思います。
0064１３２人目の素数さん
2009/01/21(水) 23:46:2610 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にはこだわってはいません。
0065１３２人目の素数さん
2009/01/22(木) 10:29:45値の方は適当にしたよ。
> 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")
0066１３２人目の素数さん
2009/01/22(木) 16:27:460067１３２人目の素数さん
2009/01/30(金) 00:42:42RとWekaのどちらを学べばよいか迷っています。
どのような違いがあるのでしょうか？
0068１３２人目の素数さん
2009/01/30(金) 15:12:54A
[,1]
1
2
3
4
[,1] [,2]
4 8
2 9
1 0
3 7
というAとBのデータの1列目の一致させて結合したいのですけど、
どうしたらいいでしょうか？
教えてください。
0069１３２人目の素数さん
2009/01/30(金) 20:43:01> 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
0070１３２人目の素数さん
2009/01/30(金) 21:13:39Weka は触ったこともないのだけど、データマイニングに特化しているイメージがある。
R なら、とりあえずなんでもできる。はず。
データマイニングしかしない、というのでなければ、R がいいのでは。
0071１３２人目の素数さん
2009/01/31(土) 08:27:08ヒント merge
0072１３２人目の素数さん
2009/01/31(土) 18:24:05>>71
お答えいただきありがとうございます。
きっとmergeを使う方法がいいのだと思います。
それで一度やってみます。
0073１３２人目の素数さん
2009/02/01(日) 01:15:41merge関数を使う場合に・・・
>merge(data1,data2,all=TRUE,by.x=1,by.y=2)
という風にいれるとおもうんですが、
by.x=やby.y=
には何を意味しているのですか？
0074１３２人目の素数さん
2009/02/01(日) 10:34:510075１３２人目の素数さん
2009/02/05(木) 02:08:45単語の総数が違う二つのコーパスで that の出現頻度に
有意な差があるか調べたいんです
頻度（割合）を出したんですがこのあとどのような検定をすればよいのでしょうか？
値一つでは検定ができないですよね。
0076１３２人目の素数さん
2009/02/05(木) 10:04:58chisq.test()
0077１３２人目の素数さん
2009/02/05(木) 11:08:08回答ありがとうございます。chisq.test()について調べました。
今回の例の場合、引数には２行２列の行列を採用すればよいのでしょうか？
質問ばかりですみません
0078１３２人目の素数さん
2009/02/05(木) 14:19:54現在Rコマンダーというのがありますが半GUIという感じです。
可能でしょうか？
可能なら参考になるものがあれば教えてください
0079１３２人目の素数さん
2009/02/05(木) 14:43:04四分割表が好みなら、それでいいよ。
> 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
0080１３２人目の素数さん
2009/02/05(木) 15:22:46本当にありがとうございます
独学なモノでなかなかわからないんですよね
参考にしてがんばります
0081１３２人目の素数さん
2009/02/06(金) 17:39:17S-plusとか初めから揃っているのではなく、
自分でGUIを設計したいのです。
Rのソースを例えばdll化してVBから読み込むとかは可能でしょうか？
0082１３２人目の素数さん
2009/02/07(土) 00:06:06それよりも、それを責任もってずっとサポートする
ことを考えたら、他のツールも考えていいのでは、
人件費が一番高いことを考えないといけない。
0083１３２人目の素数さん
2009/02/07(土) 03:23:55一生の仕事にするので大丈夫です。
あと、どうしても市販のソフトは現場の方には難しいです。
目的でボタンを選んでいけば結果が出る(統計は裏でやっている)ようなものを作りたいんです。
ExcelからRを使えることもわかりました。試作を作ってみたいと思います。
0084１３２人目の素数さん
2009/02/08(日) 16:34:00SPSSで良いじゃないか。
0085１３２人目の素数さん
2009/02/09(月) 08:49:39現場の人は回帰分析が何に使えるのかさえ知らない。現場を知らなさすぎるよ。
0086１３２人目の素数さん
2009/02/09(月) 11:08:590087１３２人目の素数さん
2009/02/09(月) 20:49:02>目的でボタンを選んでいけば結果が出る(統計は裏でやっている)ようなものを作りたいんです。
ボタンを選んでいくとRのコードを出力するようなものを目指すといいと思う。
そうすれば、教育目的にも使えそう。
あと、ど素人に使えるようにするには、データを解析にかけられるレベルに整形する機能
（データのマージ、欠測値の置き換え、変数変換、欠測値を含む個体の除外など）
を持たせる必要があるのでは？大変だと思う。
0088１３２人目の素数さん
2009/02/09(月) 22:10:36HTMLからRを呼んで結果をWindowに表示するものを作ったことがあるけど、
Rの結果をそのまま貼ったのでは現場ではよろしくないだろう。
出力整形の方が手間がかかると思う。
0089１３２人目の素数さん
2009/02/09(月) 22:43:35＞ボタンを選んでいくとRのコードを出力するようなものを目指すといいと思う。
現場の方はＲとかには興味ないんです。
興味があるのは結局どうすれば改善するの？答えを教えて！
ということだけなので、Ｒのコードすら見せないつもりです。
最終結果のみ出力ですね。
（当然、必要な人のために計算手順、考え方を示したレポートも出力しようとは思っています）
欠損値を含むデータは一式除去します。
（補完は事実ではないので使いたくないしデータは多数あるので）
正規性検定で正規性チェックして
（当然メッセージボックスには正規性チェックＯＫとか処理の意味を記載する）
当然多変量散布図も表示。
（→何をポイントとして見ればよいのかも素人にわかりやすく示す予定）
変数変換はBox-Cox変換でします。（あてはめすぎには注意したいですが。。。）
で（線形）重回帰分析を解く
というところまではＶＢＡで作っていたので。
今後、他の手法もＲで手に入るならば後は全体の分析手順を考えるだけだなと思いまして。
Ｒの本も購入したのでゆっくりやってみます。
0090１３２人目の素数さん
2009/02/09(月) 23:16:15数人ではなく、何十人聞いたことか。
成功させるためには、商売にしてしまうという
考えをしたほうが良いのかもしれない。
0091１３２人目の素数さん
2009/02/09(月) 23:19:51シーケンスなデータがベクトルに入ってるとして、
その極大（小）値を求めることってできる？
0092１３２人目の素数さん
2009/02/09(月) 23:31:28商売？？
もともと社内でそういう仕事やってるので。。。
0093１３２人目の素数さん
2009/02/09(月) 23:40:000094１３２人目の素数さん
2009/02/09(月) 23:42:02それを考える人は結構すくなくて、
社内利用を目的とする人がおおくて
それだと、いつの間にか消えていくもので。
0095１３２人目の素数さん
2009/02/09(月) 23:44:15>Ｒのコードすら見せないつもりです。
R のコードを見せるかどうかは問題ではなく、どうすれば（あなたが）楽をできるか、を考えたつもり。
作成段階でのデバッグとか、保守（不具合が発生した時の対応）とか。
もちろん、ユーザーには解析結果だけ見せればいいだろうけど、R はめちゃくちゃ強力なので、
解析結果の表示の寸前まで R にやらせればいいと思う。
>興味があるのは結局どうすれば改善するの？答えを教えて！
工場の品質管理？
>変数変換はBox-Cox変換でします。（あてはめすぎには注意したいですが。。。）
なんてダーティーな…。
（試行錯誤がある程度許される）品質改善に使うなら、そういうのもありなのかな。
製薬畑で、対数変換しか使ったことがないからよくわからない…。
>>90 >>91
システムとして他の企業に売ったりするつもりが無いなら、オープンソースでもいいかと。
企業内の統計専門家なら実験計画・解析計画と結果の解釈で飯を食えるはずだし。
0096１３２人目の素数さん
2009/02/10(火) 00:53:53自分の楽を考えていては改革は進まないのです。
楽する人には誰も信頼してついてきてくれません。
現場がいかに楽して成果をだせるか、
統計専門家は一段降りて素人と対話することが重要なんです。
わからない専門用語を並べても相手には伝わらないんです。
統計はあくまで手段にすぎないのです。
統計がなくても目的が達成できれば、統計すらいらないのです。
統計を使うことが目的ではなく、現場をよくすることが目的だからです。
本来は現場の方が自主的に統計を学べばよいのですが、現場は現在の仕事で手一杯状態。
そこを補佐するのが私の任務です。
あと、わかっているとは思いますが、
Box-Cox変換の中には対数変換も含まれていますよ。
変換は、平方根変換、対数変換が主でしょうが、
プログラム内ではBox-Cox使う方が値を変更するだけなのでコード記載が楽なのです。
モジュール一つで済みますから。
＞飯を食えるはずだし。
飯を食えるとか発想が保守的すぎる。
現場が望むものなら不況であれ何であれやろうってなりますよ。
改革派の俺にはわからん考え方ですな。。。
0097１３２人目の素数さん
2009/02/10(火) 01:15:52それを「業務」としてやっているなら>>95のいうように飯をくっていることに変わりはない。
まあどんなに簡単に使えるツールを作っても業務が改善されなければ現場からは
役に立たないと言われるだろう。難しいぞ。
0098１３２人目の素数さん
2009/02/10(火) 01:21:41自分しかできないから存在価値があるのです。
誰でもできることばかりやっていては
いきなり首を切られても文句は言えません。
0099１３２人目の素数さん
2009/02/10(火) 01:28:46で考え方は違うし、行動に差がでるよ。
前者は、飯をくうために、できる限り楽して金をかせぐことを考えてしまう。
後者は、目的を達成することに全勢力を注ぎ込む。結果的に周りからも認められ金を得られる。
0100１３２人目の素数さん
2009/02/10(火) 01:53:47シーケンスなデータがベクトルに入ってるとして、
その極大（小）値を求めることってできる？
0101１３２人目の素数さん
2009/02/10(火) 10:21:21>>83
応援だけする（笑。ただ、そういうものなら裏でR走らせるよりも、
自分でコードした方が楽かもしれない、とも思う。
0103１３２人目の素数さん
2009/02/10(火) 13:16:18Rpadがすでに用意されているんだから、車輪の再発明はいらないと思う。
0104１３２人目の素数さん
2009/02/10(火) 13:53:28あ、ごめん、勝手に最大最小に置き換わってたw
0105１３２人目の素数さん
2009/02/10(火) 19:20:18確実にR使う方が楽だよ。
0106１３２人目の素数さん
2009/02/10(火) 21:57:330107１３２人目の素数さん
2009/02/10(火) 23:28:47手段が目的化するのは普通のことだし、そもそも業務ならこんなところで
質問してないで仕事しろって感じ。ｗ
0108１３２人目の素数さん
2009/02/11(水) 00:53:12>楽する人には誰も信頼してついてきてくれません。
楽をしろ、というのは、性能・品質を妥協しろ、とか、素人にわかりにくい出力をそのまま出せ、という意味ではない。
同じ動作をするシステムを作るなら、なるべく楽をして作るべき（無駄に苦労するのは社会人として失格）だ、ということ。
そのための方策として、結果出力の寸前までを（他の言語に投げずに） R に任せることを勧めている。
あと、デバッグ・保守を楽にできないのだとしたら、それは稼働率も信頼性も低いダメなシステムを作ったということ。
プログラミングに限っていうと、一般に楽につくったシステムはデバッグ・保守も楽な「良い」システムになるから、一石二鳥。
>>98
>誰でもできることばかりやっていては
>いきなり首を切られても文句は言えません。
あなたは別の理由で職を失う心配をした方がいいかも。主張がちょっぴりうつ病患者っぽいよ。
0109１３２人目の素数さん
2009/02/11(水) 01:33:44通常ならなにもやらないやつが切られるのにな(笑)
ぼく、仕事はお遊びじゃないんだよ。
0110１３２人目の素数さん
2009/02/11(水) 01:36:15完全に慢性病。
ご愁傷さま。
0111１３２人目の素数さん
2009/02/11(水) 02:00:01詭弁。口がうまいというか世渡り上手というか。
自分が不利にならないように話の解釈を変えている。
口がうまい分実力がなくても大丈夫。
実力者についてわたり歩けばいいんだから。
口がうまい人なりの仕事の仕方はある。
否定はしない。
0112１３２人目の素数さん
2009/02/11(水) 11:28:320113１３２人目の素数さん
2009/02/11(水) 13:15:02誰も動こうとしないから楽なことは楽だが会社が潰れないか心配だ。
0114１３２人目の素数さん
2009/02/11(水) 14:50:50衰退に気づかない『ゆでかえる』かもしれない。
俺には恥ずかしながら現状を具体策をもって行動できるパワーはない。
だからその課題解決力が羨ましく思う。
反対するやつは必ずいる。負けずに行動してほしい。
0115１３２人目の素数さん
2009/02/11(水) 15:23:12ただ違うことはわかるという評論家のようなやつはいる。
違うことすらわからず、
人の性格、年齢等まったく異なる話にもっていき攻める幼稚なやつさえいる。
子供か大人かは話の論点がそもそも違う。
会社に利益をもたらすかどうかが重要。
優秀さとは、自分のアタマで考えることができ、行動できる人だと思う。
がんばってほしいと思う。
0116１３２人目の素数さん
2009/02/11(水) 16:13:54統計も統計ツールも現場も知らないお子様。
恥ずかしい限り。
> Rを他の人が使いやすいように完全GUIから入力したいのです。
> 現在Rコマンダーというのがありますが半GUIという感じです
> 目的でボタンを選んでいけば結果が出る(統計は裏でやっている)ようなものを作りたいんです。。
>SPSSは統計をわかっている人向け。
>現場の人は回帰分析が何に使えるのかさえ知らない。現場を知らなさすぎるよ。
0117１３２人目の素数さん
2009/02/11(水) 16:33:30あと>>116
具体策何か提示しないと
お子様にも負けていることを自分で言っているようなものだと思う。
0118１３２人目の素数さん
2009/02/11(水) 16:38:01いや、>>116も夢と表現しているところからすると
内心はお子様の目指している方向の正しさを認めているようだ。
しかし、実現させたいが自分にはできないといったところじゃないか。
素直にはなれない頑固さというところか。
0119１３２人目の素数さん
2009/02/11(水) 16:59:320120１３２人目の素数さん
2009/02/11(水) 17:05:52統計をわかってないと、そもそも「目的でボタンを選んで」ができない。ｗ
「現場の人は回帰分析が何に使えるのかさえ知らない」のなら、そんな人が
回帰分析をしてはいけない。
ある程度の不便さは、使う側を選別するために必要だと思う。
0121１３２人目の素数さん
2009/02/11(水) 17:12:010122１３２人目の素数さん
2009/02/11(水) 17:16:14無理じゃなくて。
できるとすれば、具体的にどうすればよいと思うか？
否定だけでは今のところみな評論家のままだよ。
俺同様、この点でお子様に負けていると思うんだよ。
0123１３２人目の素数さん
2009/02/11(水) 17:18:37＞ある程度の不便さは、使う側を選別するために必要だと思う。
これ彼が言ってた。
統計専門家は現場を意識していない。
素人と向き合っていないってことだね。
変われない専門家じゃないかよ。
0124１３２人目の素数さん
2009/02/11(水) 17:28:04いや、統計学については、中卒、高卒でもわかるようにノウハウは作られてる。
そのために先人がいろいろとテキストを書いてきてるし、きちんとしてる所は
社員に対してそのためのトレーニングもしてる。
問題は、ソフトをある程度以上に簡単化することは、統計の誤用を引き起こす
原因になるってこと。
私はそっちの方がリスクが高いと思ってる。
ということで、統計をできる人間を入れて、やらせればいいんであって、
みんなにやらせる必要はないよ。その方が無駄ってのが私の考え。
0125１３２人目の素数さん
2009/02/11(水) 17:29:32＞統計をわかってないと、そもそも「目的でボタンを選んで」ができない。ｗ
これは違うと思う。
例えば株なら○％で上がる株を探したいとか目的がある。
それは株の知識があればできる。いわゆる現場の人ならできるんだよね。
＞「現場の人は回帰分析が何に使えるのかさえ知らない」のなら、そんな人が
＞回帰分析をしてはいけない。
＞ある程度の不便さは、使う側を選別するために必要だと思う。
当然、現状のＳＰＳＳ等のソフトのような場合、わからずに使っては怪我をする。
ならこんなのはどうだろう。指導者が横にいれば使ってもよい。だろ？
彼は指導者の部分をプログラムに組み込もうとしていると思うんだよ。
素人が使っても正しく回帰分析を使うことができている。
しかも意識せずにというのが使用者にストレスがなくていいね。
0126１３２人目の素数さん
2009/02/11(水) 17:29:37分析にたえないデータってあるよね。
データの取り方から計画しないといけないことを考えると、ソフトだけ自動化
してどうするの？ってのがある。
0127１３２人目の素数さん
2009/02/11(水) 17:35:48>例えば株なら○％で上がる株を探したいとか目的がある。
>それは株の知識があればできる。いわゆる現場の人ならできるんだよね。
知識があれば上みたいなことは考えない。
こういう喩えを出す段階で、>>125の知識が危ないでしょ。
>>125は、○％で上がる株を予測できるシステムを組めるの？
（まあ、いまだとありませんとかでそうだがｗ）
また、素人に正しく回帰分析をさせるような指導込みのプログラムを組んだとしよう。
すると、「それはできません」「その処方は間違っています」と指導される。
これだけじゃ素人には改善できないので、なぜダメなのか理由も説明することになる
だろう。
でも、そんなこと言われたら、「意識せずにストレスにならずに使える」には
ならんでしょ。
その理由自体が読んでもわからないかもしれないのに。
0128１３２人目の素数さん
2009/02/11(水) 17:37:55それアンケートの話じゃないかな。
多分彼がやりたいのは製造品質っぽい。
計画部分の重要性が薄れてくる分野かも。
0129１３２人目の素数さん
2009/02/11(水) 18:02:09＞すると、「それはできません」「その処方は間違っています」と指導される。
それは自動でやってくれたらいいよね。
間違わないように指導してくれるんだから、駄目だというだけでは指導にならない。
＞知識があれば上みたいなことは考えない。
＞また、素人に正しく回帰分析をさせるような指導込みのプログラムを組んだとしよう。
君の考える指導込みのプログラムって何かな？
知識があるなら言えるってことだよね。
0130１３２人目の素数さん
2009/02/11(水) 18:08:55具体的に聞き出していったらSPSSそのもので笑わしてくれる気がしてならない。
既存のツールだって十分に統計素人にも使えるように設計されているぞ？
EnterpriseMinerもClementineもSPSSも、そしてRもな。
細かいパラメータの設定なんて一切なしの、自然言語で全て入力できるような
UIを考えているのか？
まあ、書き込み内容を見ていると統計知識も情報技術も無いのがすぐ解るけど。
0131１３２人目の素数さん
2009/02/11(水) 18:11:57ある程度のリテラシーは意思の疎通を大いに助けるから。
その時間すらとれないほどのブラック企業なら、この不景気で明日にでも倒産しておかしくない。
>>125
>素人が使っても正しく回帰分析を使うことができている。
よくご存じですねw
でも、意識が高い連中なら、統計の入門書を買って読むぐらいのことは既にしていて、本当の素人ではないと思う。
「マンガでわかる統計学」、「マンガでわかる統計学 回帰分析編」とか、売れまくってるし。
大卒のエンジニアなら、数学の才能は下手をすると彼を超えるだろうし。
0132１３２人目の素数さん
2009/02/11(水) 18:14:560133１３２人目の素数さん
2009/02/11(水) 18:15:47突っ込みたい一部分だけを取り出すのは、ネットではよくあることだが、
>>126の君が引用した部分の後も引用してもらえないかな。
>これだけじゃ素人には改善できないので、なぜダメなのか理由も説明することになる
>だろう。
>でも、そんなこと言われたら、「意識せずにストレスにならずに使える」には
>ならんでしょ。
>その理由自体が読んでもわからないかもしれないのに。
要は指導付きプログラムを作ったとして、どこまで指導するのかってこと。
結局、教科書と同じになったら意味がないんだけど。
でだ、>>129は逆質問で返してきてるけど、君の持ち出した指導込みプログラム
ってのは具体性がないので、何も言ってないのと一緒だよ。
まずは、どういう風に指導するプログラムなら、ド素人でもわかるのか言いだしっぺ
が例をだしてもらえる。
私の意見は、まったくのど素人にはムリ、専門家を育てて配置するか、それなりに
勉強しろってことなので、君の方から示してくれないと答えようがないんだけど。
私の目からうろこが落ちるような、すばらしい発想なんだと思うけどね。
0134１３２人目の素数さん
2009/02/11(水) 18:17:40＞でも、意識が高い連中なら、統計の入門書を買って読むぐらいのことは既にしていて、本当の素人ではないと思う。
＞「マンガでわかる統計学」、「マンガでわかる統計学 回帰分析編」とか、売れまくってるし。
マクロとミクロの話がごっちゃ。
素人は話にならん。本当に統計の誤用とはまさにこのこと。
0135１３２人目の素数さん
2009/02/11(水) 18:18:02Ｒバンザイ！
0136１３２人目の素数さん
2009/02/11(水) 18:20:37思考停止。慢性病患者だね。
自分で考えようとしない。
0137１３２人目の素数さん
2009/02/11(水) 18:22:27とふと思った。
0138１３２人目の素数さん
2009/02/11(水) 18:27:36統計＋プログラミング＋現場知識が必要になるね。
0139１３２人目の素数さん
2009/02/11(水) 18:28:25実際、統計専門家は現場知らないから難しいんだよ！
0140１３２人目の素数さん
2009/02/11(水) 18:29:59優秀な現場のおっちゃんが来た感じ
0141１３２人目の素数さん
2009/02/11(水) 18:30:400142１３２人目の素数さん
2009/02/11(水) 18:31:06解説よろしく！
「マンガでわかる統計学」の英語版を口語表現の参考に読んだことしかないのだけど、
普通に記述統計・推測統計の基礎を書いたものだったよ？
0143１３２人目の素数さん
2009/02/11(水) 18:32:16また統計素人がでてきた。。。
0144１３２人目の素数さん
2009/02/11(水) 18:32:54相手に「自分で考えよとしない」と書いてるけど、君は考えてるの？
できるって言うなら案を出してよ。
できるといってる方が案を出すものなのに、できないといってる人（私）に案を
出せっておかしくない？ｗ
できるって言うなら、いい案を出して、私を納得させてよ。
何も案を出さずに、
指導込みのプログラムにすればいい！
その内容は誰でもわかるものでー！
ど素人でも間違いなく使えるー！
なんていってても、単なるバカだよ。
私にすれば、>>136は現実を見ることができず、大きなことを言ってるだけで
何も作り上げることのできない、理想主義者に思える。
0145１３２人目の素数さん
2009/02/11(水) 18:36:17ソフトはある程度汎用にしておいて、あとはユーザーが自分の置かれた状況に
合わせて勉強するのが一番早い。
0146１３２人目の素数さん
2009/02/11(水) 18:36:29う ぜ ぇ
IDがないから何人で議論してるか分からん
0147１３２人目の素数さん
2009/02/11(水) 18:38:39できないという理由は、現実、誰にも証明できないから
できないとは誰も言えないはずなんだよ。
しかし君はできないと主張してる。
その時点で話す価値のない対象になってるってことに早く気づけ。
0148１３２人目の素数さん
2009/02/11(水) 18:40:590149１３２人目の素数さん
2009/02/11(水) 18:42:570150１３２人目の素数さん
2009/02/11(水) 18:49:19彼は既にＶＢＡで一度プログラムしてるようだな。
>>138の統計＋プログラムの力はあるようだ。
それで使われなかったこともあったんだろう。
現場知識の現状も把握しているようだ。
0151１３２人目の素数さん
2009/02/11(水) 18:51:38できないことの証明はできないにしても、それが「だからできる」って
証明にはならんのよ。
まったく、案を示さないままに、上記のようなことを言ってる段階で痛い。
数学スレにいるのに、頭悪いと思う。
それとだ、私のできないってのは、他の案のほうが効率的なのに、安くできるのに
無意味ってのも入ってる。
元々の人はSPSSですら、ド素人にはムリと書いているので、あれ以上のUIおよび
簡略化の上で、>>147は現場に合わせた指導まで入れるつもりだよね。
素直に、勉強会開いて統計を教えた方が早いよ。
0152１３２人目の素数さん
2009/02/11(水) 18:58:52本人ではないが少なからずできるよ。
わざわざできる理由を聞いてる方が無知に見える。
インテリジェントCAEが一例。
0153１３２人目の素数さん
2009/02/11(水) 19:00:56>>147だけど、元の人の発想については、業務を限定して、その作業が定型化されて
いるならばUIを制限することで、作業を効率化することは可能だと思う。
まあ、マクロ組めばいいって気もするが。
でも、>>125の場合は、そもそも適用する対象からしておかしいし、そのための
手法も夢を見すぎ、しかも具体案がないので、突っ込んでる。
0154１３２人目の素数さん
2009/02/11(水) 19:03:32単に、世の中で統計の本が売れまくってるから、彼の職場でも統計の本を読む奴がいておかしくない、ってことだよ。
ひょっとして、彼の職場の社員集団は日本全体の無作為サンプリングではない、とか考えた？
なんでもかんでも統計学的な言説だと解釈してはいかんよ。
>>150
>現場知識の現状も把握しているようだ。
なぜ”彼”の言うことを鵜呑みにできるのか、わからないね。
>>152
横レスだが、インテリジェントCAE を使うのは素人じゃないと思う。
0155１３２人目の素数さん
2009/02/11(水) 19:04:20正直、そんな低レベルな話をされても困る。
現場のことも実は何も理解していないから
現場の人間に受け入れてもらえなかったんだよ。
頭でっかちのルーキーには良くある事。
自分もそうだった。
0156１３２人目の素数さん
2009/02/11(水) 19:18:37う ぜ ぇ
0157１３２人目の素数さん
2009/02/11(水) 19:27:50>>151、>>153だが、>>152が見えてなかった。
>>151でも書いたけど、業務が定型化していて、やることが限られているなら
いけると思うよ。
でも、インテリジェント機能にしても、元の使い道がわからんとどうしようもない。
ある程度わかってる（本当は相当に精通してる）から、そのサポートで効率化できる
んであって、ド素人が使うものじゃないでしょ。
ド素人はサポートされると混乱する．．．
結局は、少しは勉強してド素人を脱出した方が早い。
0158１３２人目の素数さん
2009/02/11(水) 20:05:35現場を理解してないのは俺らも同じだろ？
現場に統計使え、おまいら教えるから勉強しろと。
けど数学得意な俺らにすれば理解はたやすいが得意でないのが多い。
0159１３２人目の素数さん
2009/02/11(水) 20:09:46君は一生成功しないな。
現場は統計はど素人だが現場の専門家だよ。
俺らは現場の素人。
現場からしたら少しは現場を勉強しろとなる。
どちらかが歩み寄らねばならない。
いがみ合ってもお互いにとって損だ。
0160１３２人目の素数さん
2009/02/11(水) 20:13:14どんな現場なの？
ほんとに統計勉強したことない人ばっかりなの？
そういう会社で統計使うんかな？
0161１３２人目の素数さん
2009/02/11(水) 20:19:06現場知らない人1人発見
0162１３２人目の素数さん
2009/02/11(水) 20:22:32統計知識がある君はどこに就職したい？
友達はどこに就職したいと言っている？
統計人材が流れる先は大方、アクチュアリ、リサーチャーじゃないかな。
0163１３２人目の素数さん
2009/02/11(水) 20:28:420164１３２人目の素数さん
2009/02/11(水) 20:36:43逆に統計を専門に扱わない、あくまでも手段の一つ、ただ使った方がより効率的な仕事。
こんな現場ではこのスレのような問題が起きる。
0165１３２人目の素数さん
2009/02/11(水) 21:16:430166１３２人目の素数さん
2009/02/11(水) 21:25:27具体案はそれぞれで考えるしかない。
まして現場知識がないのにこのスレで具体案がでるわけがないと思う。
各個人が現場と話し合って作り上げることが必要になる。
0167１３２人目の素数さん
2009/02/11(水) 21:28:540168１３２人目の素数さん
2009/02/11(水) 21:32:530169１３２人目の素数さん
2009/02/11(水) 21:36:300170１３２人目の素数さん
2009/02/11(水) 21:42:010171１３２人目の素数さん
2009/02/11(水) 21:43:060172１３２人目の素数さん
2009/02/11(水) 21:46:340173１３２人目の素数さん
2009/02/11(水) 21:52:220174１３２人目の素数さん
2009/02/11(水) 21:54:51ストレス社会嫌・・・
0175１３２人目の素数さん
2009/02/11(水) 21:57:08製薬会社の統計担当者だけど、薬事法、生理学、病理学、メジャーな薬や開発中の
新薬の作用メカニズムを新人研修で教わったし、その後も自主的に勉強してるよ。
生物学は高校でも選択してなかったから、たしかに大変だったけど、
営業マン（MR）相手なら作用メカニズムを教えることぐらいはできるよ。
逆に、医学・薬学・獣医学の出身者に対する新人研修には統計学が含まれるし、
工場でも研究所でも自主的な統計の勉強会をやっているよ。
どちらかが歩み寄るのではなく、どちらも歩み寄るのが当然だと思う。
職場のパートナーがやっている仕事に関する知識が素人レベル、
なんてあまりに非効率的で、冗談にもならない。どの業界・業種でもきっとそう。
>>167
彼自身がどういう業種・現場で働いているか、ということを明かさないのだから、仕方ない。
よほどマイナーな業種でなければ、統計・現場知識の両方で彼を凌駕する者はいると思う。
明かせなくなってしまったのかもしれないけど。
>>171
同感。でも、統計に限らず、学問っていうのは手段。
0176１３２人目の素数さん
2009/02/11(水) 22:06:53生物統計、医学統計ってのがあるぐらいだ。
0177１３２人目の素数さん
2009/02/11(水) 22:09:030178１３２人目の素数さん
2009/02/11(水) 22:23:24彼は本当に初歩の初歩的な質問しかしていないぞ。
周りが話を大きくし過ぎ。
> S-plusとか初めから揃っているのではなく、
> 自分でGUIを設計したいのです。
> Rのソースを例えばdll化してVBから読み込むとかは可能でしょうか？
0179１３２人目の素数さん
2009/02/11(水) 22:30:460180１３２人目の素数さん
2009/02/11(水) 22:35:420181１３２人目の素数さん
2009/02/11(水) 22:37:20彼はいずこへ！？
0182１３２人目の素数さん
2009/02/11(水) 22:38:110183１３２人目の素数さん
2009/02/11(水) 22:38:410184１３２人目の素数さん
2009/02/11(水) 22:39:33ExcelのVBAで多変量解析を駆使し、投資用ソフト開発出来る方いましたら500万円で買います。
ロジックは私が用意しますからプログラミングだけお願いします。
ご連絡お待ちしております。
[email protected]
0185１３２人目の素数さん
2009/02/11(水) 22:40:24きいてらんね〜や・・・・と
仕組みを作り始めてる
0186１３２人目の素数さん
2009/02/11(水) 22:42:53あとは実行するだけだろ、彼ならできそうだ。
0187１３２人目の素数さん
2009/02/11(水) 22:43:530188１３２人目の素数さん
2009/02/11(水) 22:44:23たぶん俺ら煽り隊かな。
0189１３２人目の素数さん
2009/02/11(水) 22:45:18その彼と、今回話を大きくして、現場とか言ってるのは別人じゃないのか？
今回の彼が持ち出した現場の例が「株」だしなー．．．
0190１３２人目の素数さん
2009/02/11(水) 22:45:410191１３２人目の素数さん
2009/02/11(水) 22:53:54Excel の VBA、っていうきつい縛りは何のため？エンドユーザーにソフトを売ることを考えてるのかな？
ロジックの特許を押さえるか、ガチガチの契約を結ばないと丸ごと持ち逃げされるよ。特許がとれたら、このスレで報告してね。
0192１３２人目の素数さん
2009/02/11(水) 23:06:360193１３２人目の素数さん
2009/02/11(水) 23:10:33ＶＢＡでプログラムした彼ならできるのでは。。
しかし彼はいない。。
0194１３２人目の素数さん
2009/02/11(水) 23:22:000195１３２人目の素数さん
2009/02/11(水) 23:28:19いや、崇高な理念に燃えていたので、こんなところで時間をつぶしてはいないはずさ！！
0196１３２人目の素数さん
2009/02/11(水) 23:31:33> S-plusとか初めから揃っているのではなく、
> 自分でGUIを設計したいのです。
> Rのソースを例えばdll化してVBから読み込むとかは可能でしょうか？
0197１３２人目の素数さん
2009/02/11(水) 23:44:59それにしても、彼とか彼女とかわけわからんw
コテ付けてくれ。
0198１３２人目の素数さん
2009/02/11(水) 23:48:34エキスパートシステムが瓦解したことを知らない世代が増えたってことだな。
0199１３２人目の素数さん
2009/02/11(水) 23:52:23あんまりLispをいじめないで下さい(><)
0200１３２人目の素数さん
2009/02/11(水) 23:55:080201１３２人目の素数さん
2009/02/12(木) 00:06:42エキスパートシステムはいろんな業界で生きていると思う。
アパレルショップに置いてあるファッションコーディネートシステムとか。
0202１３２人目の素数さん
2009/02/12(木) 00:23:07いや、崇高な理念というのは、何もわかってない現場のために、使いやすいソフト
開発しようと言う部分なんだが。
0203１３２人目の素数さん
2009/02/12(木) 00:34:52自らデータを漁って自ら最適解を見つけ出すソフトを開発するよ。
具体案はナッシングだけど。
0204１３２人目の素数さん
2009/02/12(木) 00:37:500205１３２人目の素数さん
2009/02/12(木) 00:38:26「何かわかんないけど、これの言うとおりにして、こことここに表示される数値を
写しておけばいいらしいよ」となり
「何か、出てきた数字がまずいみたいなんすけど」「あ、それは写すときに適当に
なおしておけばいいよ。p＝0.02と書いとけば問題ないみたいだ」となったり、
「何か、計算してくれないんすけど」「あ、それはデータに問題があるってこと
だから、それらしい数字を足しておけばいいよ」とかになるんだよな。
現場はいそがしいから、手取り足取りしてくれるソフトつくっても、手順が面倒だと
思われて、逆にすっとばされることがある．．．
0206１３２人目の素数さん
2009/02/12(木) 00:38:530207１３２人目の素数さん
2009/02/12(木) 00:41:310208１３２人目の素数さん
2009/02/12(木) 00:42:34アクションにつながってないから当然無意味。
まさにSPSS使ってる現場だな。
0209１３２人目の素数さん
2009/02/12(木) 00:43:53それいいアイデアだよね！
で 具体的にどうしようか？
0210１３２人目の素数さん
2009/02/12(木) 00:54:48社内プレゼンでつかいたいなら>>205でいんじゃね？
予測して儲けたい、とかいうなら、手取り足取りソフトに頼るな、っつ〜話しだし。
使う手法が決まってるなら、薄いGUIラッパ書いて終わりじゃん。
そもそも、彼の質問ってなんだったっけ？
0211１３２人目の素数さん
2009/02/12(木) 00:55:490212１３２人目の素数さん
2009/02/12(木) 00:58:400213１３２人目の素数さん
2009/02/12(木) 01:07:01男が続いたから、次こそは女の子だと思って頑張るんじゃないぞ！
0214１３２人目の素数さん
2009/02/12(木) 01:07:45「今日は安全日なの」を信じちゃダメだ。
0215１３２人目の素数さん
2009/02/12(木) 01:10:34品質管理や生産プラントの収量向上が目的なら、もうちょっとマシになると思う。
システムが吐く出力が「行程 xx の温度を 10℃上げてみてください。」とか、
「データ不足により、意味のある答えはえられませんでした。」とか、
「多重共線性が生じています。彼を呼んで変数選択をさせてください。」とか、
「内部エラーです。彼を呼びつけて叱り飛ばしてください。」とか。
それぐらいなら作れそうな気もする…
>>208
アクションにつながらない作業なんて、まっとうな企業には無い。
マーケティングデータや学術研究データの統計解析は、
その後の営業戦略策定や研究結果の考察に直結する。
それらに SPSS や SAS は思いっきり使われてる。
そして、>>205 みたいなデタラメをさせないために統計専門家がいたりいなかったりする。
0216１３２人目の素数さん
2009/02/12(木) 01:38:28>>89とか読むと
>現場の方はＲとかには興味ないんです。
>興味があるのは結局どうすれば改善するの？答えを教えて！
>ということだけなので、Ｒのコードすら見せないつもりです。
>最終結果のみ出力ですね。
>（当然、必要な人のために計算手順、考え方を示したレポートも出力しようとは思っています）
>
>欠損値を含むデータは一式除去します。
>（補完は事実ではないので使いたくないしデータは多数あるので）
データが多数あって、欠損があるのは除去しても問題ない環境らしいので、
データ不足は最初からおきなさそう。
データが多数ある、結果だけとか書いてあるし、>>125が彼だとすると、
製品の品質管理とかじゃなく、株とか景気の予測屋さんで、経済とか財務
データの分析かな。
0217１３２人目の素数さん
2009/02/12(木) 01:44:390218１３２人目の素数さん
2009/02/12(木) 01:55:110219sage
2009/02/12(木) 02:10:160220１３２人目の素数さん
2009/02/12(木) 15:25:40株の話をしてる時点で違う。
0221１３２人目の素数さん
2009/02/12(木) 22:41:02株屋さんなどの金融関係なら、「現場」とか「改善」なんて言わないと思いますよ？
どう考えても、品質管理とかプロセス関連。
>>220
いや、>>125 は間違いなく彼だね。本人しか知り得ない事実↓を知っているから。
>素人が使っても正しく回帰分析を使うことができている。
>しかも意識せずにというのが使用者にストレスがなくていいね。
客観的には、彼の認識している世界における事実、っぽいけど。
0222１３２人目の素数さん
2009/02/12(木) 23:25:48それは俺でも想像できる。
俺でもそう答える。
0223１３２人目の素数さん
2009/02/12(木) 23:29:39推測と事実を混同してはぶれる。
事実はなんだ？
0224１３２人目の素数さん
2009/02/12(木) 23:47:30う ぜ ぇ
0225１３２人目の素数さん
2009/02/13(金) 02:37:31まあ、例を出す時に株を出しちゃったけど、やってるのは品質管理とかかもね。
株の現場を知ってるにしては、稚拙な喩えだったしね。
しかし、たまにあらわれて「うぜぇ」って言ってるのは誰なんだろうね？
理想は人間を狭量にするのかな．．．
0226１３２人目の素数さん
2009/02/13(金) 07:21:27仕事できねぇ典型だな。
0227１３２人目の素数さん
2009/02/13(金) 07:52:25このスレを建設的に進めたい。
個人の悔しさをぶつけるのは他でやってくれないか。
0228１３２人目の素数さん
2009/02/14(土) 08:30:21ものすごくレスが増えてて驚いたが、イタいコが複数人湧いてただけか
>>225と、あとは分かんね
この板にもID欲しいな
0229１３２人目の素数さん
2009/02/14(土) 10:27:29いまさら書いてる君も痛い子。
0230１３２人目の素数さん
2009/02/15(日) 14:28:13志と実力がつりあってないから「イタい」と批判されるわけで、
志に見合った実力を身につけて、見返してやればいいさ。
その過程で、実は批判者の側が正しいことが分かったりもするけれど、
それもまた人生。
0231１３２人目の素数さん
2009/02/15(日) 21:39:310232１３２人目の素数さん
2009/02/15(日) 22:52:02>>225は俺だが、>>231とは別人の件
（これは本当ね）
0233１３２人目の素数さん
2009/02/15(日) 22:53:20>>225は俺だが、>>230とは別人の件
（これは本当ね）
0234１３２人目の素数さん
2009/02/15(日) 23:42:410235１３２人目の素数さん
2009/02/16(月) 08:53:57このような人に、逆に格上の人を挙げてみてくださいと言えば、有名人しか挙げられないんじゃないかな。
部分的に見れば大人ですら小学生に負けてる知識もある。
勝ち負けの発想は和を生まないし幸せをもたらさない。
どんな相手でも敬っていきたいものです。
みんなで協力して一つのものを作り上げていきたいものです。
0236１３２人目の素数さん
2009/02/16(月) 11:07:50Ｒ好きのどこかの教授みたいだな。
0237KingGold ◆3waIkAJWrg
2009/02/16(月) 21:54:140238１３２人目の素数さん
2009/02/16(月) 22:09:070239１３２人目の素数さん
2009/02/17(火) 03:00:170240１３２人目の素数さん
2009/02/18(水) 01:17:44King→王様
Gold→一位
どんだけ自分大好きなんだよ。
ネタ仕込みに座布団一枚。
0241１３２人目の素数さん
2009/02/18(水) 02:46:10あまりのウザサに、俺も含め昔はいた優良回答者もだんだんいなくなってる気がします。
何とかしてください。
0242１３２人目の素数さん
2009/02/18(水) 08:45:08肩書きもあるんだから
恥ずかしいので出てこないで欲しい
そのうち誰かばれちゃうぞ
0243１３２人目の素数さん
2009/02/18(水) 10:31:102チャンネルのノリで書き込んでいるcapetapetaという盆栽好きのニートに限っては、質問者の方に
問題があると思うけど。回答者を馬鹿にするにもほどがある。もちろん、スルーできない大人も悪いが。
0244１３２人目の素数さん
2009/02/18(水) 11:08:09離れていく人も（の方が？）多い事実に
気が付いていないようす
Sもそうだった
0245１３２人目の素数さん
2009/02/22(日) 19:30:31>macでRを使ってるんですが、
>グラフィックボードに日本語を出力しようとしても
>全部四角になってしまいます。
昨日からRを使い始めたのですが、par(family="Osaka")とかで
対応できますか？
ttp://oku.edu.mie-u.ac.jp/~okumura/stat/first.html
0246１３２人目の素数さん
2009/02/23(月) 12:46:46ttp://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
0247１３２人目の素数さん
2009/02/27(金) 20:27:46ここで三つの軸の範囲を自分で設定するにはどうすればよいのでしょうか？
よろしくお願いします。
> 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)
0248１３２人目の素数さん
2009/02/27(金) 21:28:31persp()にxlim,ylim,zlimっていう引数があるけど？
0249１３２人目の素数さん
2009/02/27(金) 23:02:25それでうまくいきました。ありがとうございました。
0250１３２人目の素数さん
2009/02/28(土) 07:57:35hist2dでcountsをヒストグラム表示しましたが、これを頻度にして表示する方法はないでしょうか？
(1, 5, 4, 2, 1, 7)というcountを(0.05, 0.25, 0.2, 0.1, 0.05, 0.35)と頻度に換算して表示したいのです。
0251１３２人目の素数さん
2009/02/28(土) 15:34:45h2d$counts を h2d$counts/sum(h2d$counts) にするだけだと思いますが。
0252１３２人目の素数さん
2009/02/28(土) 15:42:40ちょっと自分で考えれば分かる質問をするクズや、それに丁寧に回答する優しいクズもいることも知って感心した
0253１３２人目の素数さん
2009/02/28(土) 16:16:220254１３２人目の素数さん
2009/02/28(土) 19:06:00よく分かっていれば質問もしてこないだろうし
場所も場所だし。これが現実だろう。
どうしていきたい？
いい解決案でもあるかな？
自分の考えているレベルのスレッドか
ブログでも建てるのも手じゃない？
0255１３２人目の素数さん
2009/02/28(土) 19:46:49ありがとうございました。その通りで行けましたmOm
0256１３２人目の素数さん
2009/02/28(土) 20:29:260257１３２人目の素数さん
2009/03/02(月) 13:45:00http://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
どっちのが古いとかってありましたっけ？
0258１３２人目の素数さん
2009/03/02(月) 14:52:34復活したな
0259１３２人目の素数さん
2009/03/02(月) 17:56:23どうも中央農業総合研究センターのほうが新しいみたいね。
0260１３２人目の素数さん
2009/03/05(木) 23:01:16gtools.dllがR2.8.1のフォルダに入っているという・・・
何で？
0261１３２人目の素数さん
2009/03/06(金) 19:51:04アンチウイルスソフトは何ですか？
アンチウイルスソフトのサポートページを見ると情報が得られるかもしれません。
0262１３２人目の素数さん
2009/03/07(土) 02:33:36単刀直入に難しいですか？
参考までに
高校時代は数学は２Bまでやりました（あまり覚えてないですが）
経済商系統の学部の学生です
統計学は多少わかります
0263１３２人目の素数さん
2009/03/07(土) 11:10:54基本はコマンド入力だからね。
0264１３２人目の素数さん
2009/03/07(土) 19:33:190265１３２人目の素数さん
2009/03/07(土) 23:06:21アセンブリ言語？
0266１３２人目の素数さん
2009/03/08(日) 01:07:58文系だと「コマンド入力」という言葉自体を何それ？と思っても不思議じゃないが。
0267１３２人目の素数さん
2009/03/08(日) 12:31:03r-tips.pdf を見ているのですが、よく分かりません。
何か変数に ~ を付けてみても、
> x <- 1:10
> ~x
~x
となるだけで、何をするものなのか分かりません。
0268１３２人目の素数さん
2009/03/08(日) 14:45:590269１３２人目の素数さん
2009/03/08(日) 19:21:58~ って演算子じゃないんじゃないか。
モデル式のとき = の代わりに使うしかできないと思うけど。
0270１３２人目の素数さん
2009/03/08(日) 20:55:24例えば、
> グループXの標準偏差<-
と入力した時点でXじゃなくてYの間違いだと気づいたとします。
そうして左矢印キーで赤いキャレットを戻してXのところに合わせてBackSpaceキーで消そうとすると、
「X」ではなくて「プ」が消えます。
ちなみに環境はWindows VistaでRのバージョンは2.8.1で日本語入力システムはATOK2007です。
0271１３２人目の素数さん
2009/03/09(月) 01:24:15なるほど。= の代わりか。ありがとうございます。
ただ、どうも演算子ではあるようなのです。
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)
0272１３２人目の素数さん
2009/03/09(月) 19:45:52> plot(group ~ extra, data=sleep)
> plot(extra ~ group, data=sleep)
も面白い。
モデル式同様関数の引数の中でしか使わないのかな。
0273１３２人目の素数さん
2009/03/10(火) 11:30:55help('~')
0274１３２人目の素数さん
2009/03/10(火) 21:16:37~xが何なのか見たければstr(~x)とすればよい
~xがクラスformulaを返していることが分かるでしょう。
~はクラスformulaを返す演算子ということではないでしょうか。
0275１３２人目の素数さん
2009/03/17(火) 16:00:25Warning message:
In gsub(perl = TRUE, "(?<!\\\\)((\\\\\\\\)*)\"", "\\1\\\\\"", "f") :
perl=TRUE は UTF-8 ロケールに対してのみ完全実装されています
という警告文が出てきて、反映されるのは、編集直後の表示だけなのですが、どうすれば、その後もデータを反映させれますか？
エラー文からして文字コードがいけないと思うのですが、文字コードはどこで変えればいいのかわかりません。
0276１３２人目の素数さん
2009/03/17(火) 18:06:20環境とかは秘密ですか？
> 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' 中に因子水準を加えました
こちらでは特に問題ありません。あなたの問題が再現する方法を提示した方がよいよ。
0277１３２人目の素数さん
2009/03/17(火) 19:53:080278１３２人目の素数さん
2009/04/04(土) 20:09:02あとはここだけだな。
データマネジメントも苦手か。
0279１３２人目の素数さん
2009/04/05(日) 09:49:37それってものすごいメモリ消費になってるんじゃ・・・？
0280１３２人目の素数さん
2009/04/05(日) 11:37:290281１３２人目の素数さん
2009/04/05(日) 12:59:31それってなんか違うような気がする。
SASからRを呼び出すイメージ？
0282１３２人目の素数さん
2009/04/05(日) 13:06:09ｃとかjavaとかね
0283１３２人目の素数さん
2009/04/05(日) 13:45:44使い勝手悪すぎ！
ならcで全て処理を書くよ！
java･･･全部Rで行くよ。
0284１３２人目の素数さん
2009/04/05(日) 15:58:08臨機応変に出来なきゃ
ｃでかけるならR要らないだろ
出来ないから使ってるんじゃないかな
そういうものだよ
0285１３２人目の素数さん
2009/04/05(日) 23:33:02awkか？perlか？それともSASか？
備え付けのデータセットしか分析していないとかいうオチはないよな？
0286１３２人目の素数さん
2009/04/06(月) 11:27:24データ加工の前処理は、日常的に使っているsed tr paste cut sort
uniq colrm nkfあたりで十分だけど。そのあとは、すべてRの中で処
理しちゃうな。たまにperlも使うけど。
0287１３２人目の素数さん
2009/04/06(月) 22:06:16なんとストイックな。でももう少し楽したいとは思わないの？
あと扱うレコード数は１万件ぐらいがMAXでしょ？
0288１３２人目の素数さん
2009/04/07(火) 08:38:51メモリ次第だろJK
0289１３２人目の素数さん
2009/04/07(火) 11:05:04私が日頃使うデータは、そんなタイプ。何重にも階層化されたlist構造
をもつので、sapplyを入れ子にしまくり。これをperlで書いたら大変な
ことになるよ。
0290１３２人目の素数さん
2009/04/07(火) 11:55:32わける必要があり、289の言うのは取り込めた後で、
285の言わんとするところは（おそらく）取り込む前の
ことを指していると見ることができそう。
たしかにperlで書いても出力後の果は、またテキストだし
値の管理もできない。
取り込めた後なら289が言うようにRは強い。
元のデータが煩雑な場合にRで取り込みやすいように加工するには
286の形になるだろう。
データの管理を簡易に任せたい、あるいは、汎用コンピュータからの
日本語のEBCDICや他のコードが入り混じっている、あるいは、極端に
複雑なレイアウトで大容量なら、SASを使用るメリットは大きい。
どの点を言っているのか、どの程度まで必要かによって向き不向きがあるので、
一言では、な〜んとも答えきれないんだろうな〜
誰かデータ準備してくれる人いれば、あるいはDBで管理されているデータだけ
を使用するのでよければ、Rだけで済むかもしれないし。ふ〜〜
0291１３２人目の素数さん
2009/04/07(火) 14:21:10加工前と後を分けるのではなくて、どこまで加工した段階で、Rに渡すのかという
タイミングの話をしていたのでないの？例えば、めちゃくちゃ複雑な階層構造をもつ
XMLがあったとして、perlなどで十分に加工してからRに取り込むのか、またはXMLを
とりあえず取り込んでからRの中で編集加工するのか、とか。
0292１３２人目の素数さん
2009/04/07(火) 14:32:08ちなみに、上記で言った意味ではXMLは複雑な構造までいかない。
きれいな構造だし。
めっちゃめっちゃなのがあるんだよ。
0293１３２人目の素数さん
2009/04/07(火) 18:22:39気をつけてる点とかってどんなのありますか？
私はt<-12345とか死ぬほどやって頭にきた挙句、自分が使う変数は必ず
.t <- 12345
とか、関数なら
my.func <- function(...) ...
とかドット始まりか自分プレフィックス始まりにするようにしたんですが、
なんでもドット始まりとか正直ウザー、です。このあたりのR的慣習とか
ベストプラクティスってどこかで紹介されてたりするんでしょうか？
0294１３２人目の素数さん
2009/04/07(火) 21:18:25命名するとき、"_"（アンダースコア）を最初につけるって
いう風習をどっかで学んでから、シェルスクリプトを書く
ときにはそうしている。
Rで_は全くお薦めでないが。.（ドット）は隠しファイルみ
たいでやだな。実際、隠し関数（内部関数とも言う）には、
ドットから始まるように命名しているのも多いし。
あとは、何がRの変数名で使えるんだったけ？Namespaceをう
まく使えばよいのかもしれないけど、
Internal name space support functions. Not intended to be called
directly.
って書いてあるからユーザは普通使わないだろうな。
0295１３２人目の素数さん
2009/04/07(火) 21:38:49(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
0296１３２人目の素数さん
2009/04/07(火) 22:04:59> 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 はちゃんと呼ばれるようだ。
0297１３２人目の素数さん
2009/04/08(水) 00:38:50あれ？とりあえず判り易い例としてtを出す前にt<-12345で試して
これなら問題ないなと確認したつもりだったのにorz
まあ、"t"限定でなく、要はうっかりぶつかったり混乱しそうな名前付けとかを
回避するためのネーミングルールとか、何か慣習的にやってるものありますか？ってことで。変数はぶつからないから気にしなくてもダイジョブ？
0298１３２人目の素数さん
2009/04/08(水) 01:15:54ドットは慣例があるからたとえば、
is.xxx なら、xxxであるかどうか判定する関数だし
read.xxx なら、xxx形式のデータの読みこみだし
じゃあwrite.xxxは？ plot.xxx? わかるでしょ。直感的に。
0299１３２人目の素数さん
2009/04/08(水) 01:25:21一文字変数は極力避けたほうがいいよ。
そして、先頭ドットは隠れ関数/変数の意味合いがあるよ。
0300１３２人目の素数さん
2009/04/08(水) 01:33:42Rではネーム空間と関数空間が別ってことでしょう。
これはlisp系の影響を受けた言語ならあながち不思議ではないし、普通に同じ変数と
同じ関数の名前があるなんてざらなのです。Cの系統から見れば不思議なのですがね。
仕様をみると、Rって記述はCやjavascriptに似ておりますが、中身は実質lispの方言という
印象なんでね。
0301１３２人目の素数さん
2009/04/08(水) 01:51:53;(defun foo (vector)
; (let ((length (length vector))
; …))
みたいなことを差して 300を書いてます。一つめのlengthは変数だし二つめは
関数だしね。Rで書けば
> foo <- function(vector)
+ {
+ length<-length(vector)
+ …
+ }
こんなんかな。
0302１３２人目の素数さん
2009/04/08(水) 02:07:30お、言語に詳しい人が来たみたいね。
だけどそれは違う。R の名前空間は、変数と関数とで分かれていない。
影響を受けているのはLISPの中でもScheme。
もしCommon Lispみたいに、名前空間が分かれているなら、
> x <- function() print("x")
> x()
[1] "x"
> x <- 1
> x()
エラー： 関数 "x" を見つけることができませんでした
このエラーは起こらないはずでしょ？
変数のスロットが一つしかなくて、function() print("x") が 1で上書きされ
てしまうから、こういうことになる。
0303１３２人目の素数さん
2009/04/08(水) 02:16:26なるほど、てっきりわかれてるのかとおもったけど違うのか。
0304１３２人目の素数さん
2009/04/08(水) 02:35:18>>302 の x はかち合って、c や t がかち合わないのは、c や t がユーザ
定義変数のスコープの外側にある、組み込み変数スコープで定義されているから。
関数と変数で名前空間の分かれた言語(Common Lisp, Rubyなど)は、名前の衝突
は起こりにくいが、変数に関数を代入する手続きが煩雑になり、高階関数を書
きづらい。
関数と変数の名前空間が分かれていない言語(Scheme, JavaScript, Pythonな
ど)は、関数を変数の一種のように扱えて、高階関数が書きやすいが、関数名と
変数名がぶつかってしまう。だから例えばSchemeだと、リストを収める変数に
list という名前は使えず、lis とか書くのが慣例。
つうことで R は、>>296の通り、関数呼び出しの文脈で、ある変数に束縛され
ているオブジェクトが実行可能でなければ、その外のスコープに同名の変数を
探しに行くことにした。
このおかげで、変数に関数を自在にセットする便利さと、組み込みの関数名と
変数名がかち合わない便利さの両方を得たわけ。
0305１３２人目の素数さん
2009/04/08(水) 11:11:050306１３２人目の素数さん
2009/04/12(日) 19:01:44具体的には以下のようなファイルを読み込む場合です。
1,2,3,4,5,
,,8,9,10,
11,12,,,
scan("cliptable",fill=TRUE,na.strings="")
ではうまく読み込めません。
欠損値が詰められます。
0307１３２人目の素数さん
2009/04/12(日) 20:10:160308１３２人目の素数さん
2009/04/12(日) 20:32:44そんなわけで、
scan("cliptable",fill=TRUE,na.strings="",sep="\t")
でいけました。
お騒がせしました。
read.table系は処理コストが高いのであまり好きではありません。
0309１３２人目の素数さん
2009/04/12(日) 20:34:27できるのはわかってたのね。 ふーん
0310１３２人目の素数さん
2009/04/12(日) 22:05:41固有値分解しないと駄目？
0311１３２人目の素数さん
2009/04/13(月) 00:02:310312１３２人目の素数さん
2009/04/13(月) 00:15:12最小多項式を求めるとか
0313１３２人目の素数さん
2009/04/15(水) 12:59:35＝ 統計解析フリーソフト Ｒ 【第２章】 ＝
http://science6.2ch.net/test/read.cgi/math/1152449095/
＝ 統計解析フリーソフト Ｒ ＝
http://science4.2ch.net/test/read.cgi/math/1062650510/
0314１３２人目の素数さん
2009/04/19(日) 23:04:19重要な問題だと思うんだけどなぁ。
0315１３２人目の素数さん
2009/04/20(月) 00:11:49なんであんな糞わかりにくい本が5つ星なんだよ
0316１３２人目の素数さん
2009/04/27(月) 01:51:53他にまともな本があるならともかく
0317１３２人目の素数さん
2009/05/03(日) 15:04:27適当なSQL鯖にRODBCでつなげばいいんじゃね？
0318１３２人目の素数さん
2009/05/04(月) 16:08:05model{
for ( i in 1:n ){
x[i] ~ dnorm ( b0 , tau )
}
b0 ~ dnorm (0 , 0.001 )
tau ~ dnorm(0 ,0.001)
}
単純にデータｘがあって、平均と標準偏差を求めたいだけです。
0319１３２人目の素数さん
2009/05/04(月) 20:30:25そう言うもの(別ツール）に頼るかしかないのかってこと
を言っていると思うよ。
0320１３２人目の素数さん
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
0321１３２人目の素数さん
2009/05/05(火) 12:10:110322１３２人目の素数さん
2009/05/07(木) 22:33:39知っている方がいましたら教えていただけると助かります。
0323１３２人目の素数さん
2009/05/08(金) 01:08:47heatmap(z)
で、heatmapを描いてくれるのですが、
縦軸横軸に個々の要素の名前（名前が付いていない場合は数字）
が付きます。この名前、もしくは数字を消すにはどうした良い
ですか？教えていただけると大変助かります。
0324１３２人目の素数さん
2009/05/08(金) 13:40:32糖鎖とタンパク質との反応をheatmapで描いている人？
heatmapエンター
でheatmapの定義内容が表示されるので、当該部分を編集して
my.heatmap <- function(....
と新コマンドを作成したらどうですか？
あと、?axis を読んでおく。
0325１３２人目の素数さん
2009/05/09(土) 10:32:40データ数が少ないときデータの平均値はt分布に
従う、と書いてある。笑わせてもらいました。
0326１３２人目の素数さん
2009/05/09(土) 11:53:222つの変数の両方が合致したときのみデータフレームを連結してくれるような。
0327１３２人目の素数さん
2009/05/17(日) 12:00:220328１３２人目の素数さん
2009/05/17(日) 12:11:27平均とかと一緒に
0329１３２人目の素数さん
2009/05/17(日) 12:19:080330１３２人目の素数さん
2009/05/17(日) 13:35:15あまりに素人質問過ぎて同一人物なのがミエミエ。
最大値を求めるのはmax、最小値を求めるのがmin、
summaryやfivenumでその他の統計量と同時に出力。
cutは区分値を求めるものであり最大値・最小値を求めるものではない。
0331328
2009/05/17(日) 14:02:120332１３２人目の素数さん
2009/05/17(日) 14:33:130333１３２人目の素数さん
2009/05/17(日) 17:50:34してくれば？
0334１３２人目の素数さん
2009/05/17(日) 18:18:220335１３２人目の素数さん
2009/05/17(日) 22:30:280336１３２人目の素数さん
2009/05/18(月) 01:50:05???
どーいうこと?
…体重でもなんでもいいけど、例えば四分位数求めたい時にcut使う(?)(間違ってたらスマソ)
0337１３２人目の素数さん
2009/05/18(月) 09:29:37>>328
>>329
素人かどうか、同一人物かどうかはどうでもいいけど、取りあえずRをインストールして
みたらどうだろう。少なくともRを使ったことのある人の書き込みとは到底思えない。
>>327の質問の内容・仕方では、>>330の答え以上は期待できないよ。
fivenum()なんて知らんかった。
0338１３２人目の素数さん
2009/05/18(月) 16:43:39Rの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上で使用しています。
0339１３２人目の素数さん
2009/05/18(月) 16:51:50この問題の答えがabθ/2になるのですが、途中計算がわからないので教えてくれませんか？
∫の後ろの括弧は積分範囲です
0340１３２人目の素数さん
2009/05/18(月) 16:52:430341338
2009/05/18(月) 16:56:49Excel 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)
うーん、何なんでしょ。
0342１３２人目の素数さん
2009/05/18(月) 17:48:28> (4.8+2.5)/2 - (4.9+2.4)/2
[1] -4.440892099e-16
0344388
2009/05/19(火) 10:26:48ttp://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)
}
としました。お騒がせしました。
0345１３２人目の素数さん
2009/05/19(火) 18:28:12ちょっとお聞きしたいのですが、クリップボードにあるcsvもしくはタブ区切りのファイルを
データフレームとして取り込みたいと思ったらどうするのが一番手っ取り早いでしょうか？
0346１３２人目の素数さん
2009/05/20(水) 11:14:14Windowsユーザではないが、放置されてかわいそうなので、
> RSiteSearch("clipboard")
であなたの代わりに調べてあげたよ。
ttp://finzi.psych.upenn.edu/R/library/psych/html/read.clipboard.html
0347１３２人目の素数さん
2009/05/21(木) 09:27:21まずヘルプ回りの関数類の使い方を憶える事と、パッケージの概念を理解するところからですね、ややこしいです……
0348１３２人目の素数さん
2009/05/26(火) 00:17:17なんつーか見る者を絶望させる言語ですね……
ヘルプの日本語訳(古いバージョンのでもいい)があるといいんですけど……
0349１３２人目の素数さん
2009/05/26(火) 01:03:240350１３２人目の素数さん
2009/05/26(火) 09:56:02なにぃぃぃぃ、そいえばステップ実行はしたけれど、そこにはそれらしい情報なかったな……
0351１３２人目の素数さん
2009/05/26(火) 09:58:390352１３２人目の素数さん
2009/05/26(火) 12:37:13まぁなんだ。頑張れ！
身近に教えてくれる人はいないの？
0353１３２人目の素数さん
2009/05/27(水) 12:30:260354１３２人目の素数さん
2009/05/27(水) 12:57:03というだけの話。コンピュータ開発言語と同列に考えたら駄目だよ。
>>349
エラーが出たときにエラーが出た行で止まってくれるから、発生行番号が
不明でも一目瞭然なんだけど。わざわざoptions()でechoを止めているの？
>>353
意味が分からん。Rの問題ではなくて、開発環境（テキストエディタなど）
の問題じゃないのか？C-c C-cで文脈を判断して、パラグラフなり関数定義なりをまとめて実行してくれるから、まとめてデバッグできるけど。、
0355１３２人目の素数さん
2009/05/28(木) 00:04:030356１３２人目の素数さん
2009/05/28(木) 11:10:30>エラーが出たときにエラーが出た行で止まってくれる
そうか、エラーが出た時はコマンドプロンプトにスクリプトをコピペして
動かしてみれば良かったんだな。
0357328
2009/05/28(木) 12:34:10C-cはCtrl+cって言う意味。
MはAltかEsc
0358１３２人目の素数さん
2009/05/28(木) 12:35:21エディタから、いちいちコマンドプロンプトにコピペしていたら、疲れない？
文明の利器を使いなよ。エディタから、カーソル行、または選択複数行をRに
送って実行する機能があるだろ。
0359１３２人目の素数さん
2009/05/28(木) 12:36:01ごめんほかのスレの質問と勘違いした。
xyzzyスレだと思った。
0361１３２人目の素数さん
2009/05/28(木) 13:14:11Rを一度使ったことあったからそのときに2chブラウザのこのスレをお気にスレに入ってただけ。
てっきりxyzzyの初心者の質問かと思って。
0362１３２人目の素数さん
2009/05/28(木) 13:33:08例は良くないけど、
基本は、昔の汎用コンピュータでの
ラインバイラインでのインタプリタと
同様と考えるといい。
0363１３２人目の素数さん
2009/06/02(火) 00:39:32データの保存について質問です
ワークスペースを消す際に作業スペースを保存しますか？
という質問が出ますがここではいを選択すればこれまでの作業内容が保存されるのでしょうか？
前回この操作をしたにも関わらず新しくRを立ち上げても
以前にセーブされたワークスペースを復帰します
という文が出てくるだけで以前の作業内容が反映されません
作業スペースの読み込みコマンドからR Workspace (.RData)を選択して開こうとしても
load("C:\\Users\\kobayashi\\Documents\\.RData")
の一行が返されるだけです
他にどのような操作をすればよいのでしょうか？
0364１３２人目の素数さん
2009/06/02(火) 10:10:03初心者ならではの初々しい質問だな。
.RDataは現在Rが作業しているフォルダに保存されるので、違う場所（フォルダ）で
Rを起動すればそのときに自動的に読み込まれる.RDataは別物になる。現在地（カレン
トディレクトリ）とは別の場所の.RDataを読み込みたいときには、適宜自分でload()を使う。
loadやsage.imageのヘルプを読もう。
> 以前の作業内容が反映されません
library()の結果は保存されない。保存されるのはRオブジェクト（関数、リスト、
データフレーム、行列など）のみ。
ls()とすると、前回保存したRオブジェクトが今もあるのかどうか確認できる。
最初はいろいろととまどうこともあるだろうけど、頑張れ小林君w
0365363
2009/06/03(水) 11:36:14ls()で確認したところ以前保存した内容がオブジェクトとして保存されていることが
分かりました
ところでこのオブジェクトをワークスペース上に呼び出すにはどのように
したらよいのでしょうか？
ヘルプを読んだ限りでは
load("all.Rdata", .GlobalEnv)
のように（ファイル名は自分のものに変えています）すれば出てくるのかと思ったのですが
以下にエラー load("普及率.Rdata", .○○○) : オブジェクト '.○○○' がありません
のようなメッセージが出てしましました
他にload(○○○）load("○○○")などもためしてみましたがうまく行きません
どこが間違っているのでしょうか？
何度もすいません
0366363
2009/06/03(水) 11:39:13そのまま名前を入れれば良かったのですね
364さん本当にありがとうございました
0367１３２人目の素数さん
2009/06/03(水) 12:28:38どういたしまして。
数時間とか作成に時間がかかるRオブジェクトは保存した方がよいけど、
外部ファイルを読み込んでから数秒、数十秒でできあがるRオブジェク
トはいちいち保存しないで、Rスクリプトの方を保存した方がメンテナ
ンス性が高いと思う。
Rのインターフェイスと直接的に対話的操作をするのではなく、テキスト
エディタを通じて、Rスクリプトの現在行、選択行、パラグラフをRに送
って実行する形態で作業した方が効率が上がると思うよ。Rと通信できる
テキストエディタはいくつかある。Windows版を使っていないので、いい
加減なことを言っているかもしれないけど、最近のWindows版Rは、Rとや
りとりする内蔵エディタが組み込まれていなかったけ？
0368１３２人目の素数さん
2009/06/04(木) 01:43:10ある動物に取り付け、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以降のデータをくっつける
といった作業が必要なのでしょうか？それはあまりに手間がかかる気がしたので質問させて頂きました。
0369１３２人目の素数さん
2009/06/04(木) 03:04:350370１３２人目の素数さん
2009/06/04(木) 05:49:33http://academy6.2ch.net/test/read.cgi/psycho/1012801769/
0371１３２人目の素数さん
2009/06/04(木) 10:44:49読み込み時／後に、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)
0372１３２人目の素数さん
2009/06/04(木) 11:17:06> 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における日付型の扱いを勉強しよう
0373１３２人目の素数さん
2009/06/20(土) 00:13:21いう情報がネットにないのですが、
なぜなんでしょう？
特にMTS法はmahalanobis()関数
使って何とかならないものでしょうか？
出来る方やり方教えてください。(´・ω・｀)
0374１３２人目の素数さん
2009/06/23(火) 23:18:50データﾌﾚｰﾑ上の変数の値を、別途定義したベクトルの数値に
置き換えたいのですがどうにもうまくいきません。
for文でまわしてコピーするしかないのでしょうか？
0375１３２人目の素数さん
2009/06/24(水) 01:36:50x[ , 3] <- y
みたいではできない？
0376１３２人目の素数さん
2009/06/24(水) 04:07:19質問が抽象すぎて答えにくい。
> (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
こういうこと？
0377374
2009/06/24(水) 20:58:37自作関数で計算したベクトルをデータフレームに代入しようと
してたんだけど、返値のベクトルをnumeric(0)で初期化してたのが
まずかったみたい。
numeric(0)での初期化をやめたら375さん、376さんの方式で
できました。
ありがとー。
0378１３２人目の素数さん
2009/06/25(木) 17:35:08lkappaでカッパ係数簡単に計算したいんだが、library(psy)が読み込めないから、てで計算するしかないんですか？
0379１３２人目の素数さん
2009/06/25(木) 18:32:40psyパッケージをインストールしていないだけでは？
> RSiteSearch("kappa")
でカッパ係数を計算する関数がざくざくでてくる。あちこちのパッケージ作者が
自分で独自に実装しているみたいで、車輪の再発明の典型だな。
0380１３２人目の素数さん
2009/06/25(木) 18:49:04サンキュです！
0381１３２人目の素数さん
2009/06/26(金) 17:51:55それで379さんの書いたものを読み込ませサイトに接続したのですが、どうやってパッケージをインストールするのかわかりません
0382378
2009/06/26(金) 18:17:50試しにやってみたら、普通にインストールできるけど。
好みの問題かもしれないけど、管理者権限で実行した方が無難かも。
っていうか、ちょっと調べたり本を読めばすぐに分かるレベルの質問をしても相手にする人は少ないと思うよ。
スルーできていないのは私だけみたいだし。
> 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)
0383379
2009/06/26(金) 18:21:280384１３２人目の素数さん
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
・・・・
ニーズありますよね？自作でも良いのでもしありましたらご紹介下さい。
0385１３２人目の素数さん
2009/07/20(月) 14:55:38できずに悩んでいます。
> d <- data.frame(rbind(c(1,2,3),c(4,5,6)))
> d
X1 X2 X3
1 1 2 3
2 4 5 6
これに対して、各行ごとにカラムの最小値を選んで
X1 X2 X3 X4
1 1 2 3 1
2 4 5 6 4
のような結果を得たいのですが
> d$X4 <- min(d$X1, d$X2, d$X3)
だと
> d
X1 X2 X3 X4
1 1 2 3 1
2 4 5 6 1
のようになってしまいます。そうなる理由は分かったのですが、ではどうすれば
望みの結果が得られるのか、というところで詰まっております。
参照すべきドキュメントのポインタでもかまわないのでどなたかお教えいただけ
ますでしょうか。
0386１３２人目の素数さん
2009/07/21(火) 00:07:50> d <- data.frame(rbind(c(1,2,3),c(4,5,6)))
> cbind(d,X4=apply(d,1,min))
X1 X2 X3 X4
1 1 2 3 1
2 4 5 6 4
apply関数群を勉強してくれ
0388１３２人目の素数さん
2009/07/22(水) 07:01:370389１３２人目の素数さん
2009/07/22(水) 10:15:48ほぉ、では、pminを使って、>>385の課題を解いてみて。
0390１３２人目の素数さん
2009/07/23(木) 01:25:18memory.sizeで回避できるの？
教えてエロイ人。
0391１３２人目の素数さん
2009/07/23(木) 10:48:47>>390
9MBという判断はどうやって判断したの？
> x <- runif(10000000)
> save(x,file="tmp.gz")
> system("ls -lh tmp.gz")
-rw-r--r-- 1 hoge hoge 51M 2009-07-23 10:46 tmp.gz
runif(10000000)で51MBになったよ。
0392１３２人目の素数さん
2009/07/23(木) 20:49:56時間がかかるのですが。。。
0393１３２人目の素数さん
2009/07/23(木) 20:51:52必要ない
0394１３２人目の素数さん
2009/07/23(木) 22:02:21dat <- scan( dat.csv , what=list(a="",b=0,c=0,d=0,e=0) , sep="," )
こんな読込みをした際に、エラーが発生しました。
エラーメッセージに、「9MB以上のベクトルは生成できません。」とありました。
0395１３２人目の素数さん
2009/07/24(金) 10:14:49scan()は日頃使わないから、外すかもしれんが、そもそもwhatの中身はそれでいいの？
If 'what' is a
list, it is assumed that the lines of the data file are
records each containing 'length(what)' items ('fields') and
the list components should have elements which are one of the
first six types listed or 'NULL'
ってヘルプには書いてあるけど。
１レコードが複数行とか特別の事情がない限りscan()を直接使わないで、read.csv()で十分では。
0396１３２人目の素数さん
2009/07/25(土) 00:21:19のSTLデータを作成するために、建物形状のポリゴン(凹の場合もあり)の
三角形分割をしたいのだけど、gpclibのtriangulateで分割すると、
各節点からのx方向へのスキャン・ラインを使った、無駄な三角分割をする。
adsのtriangulateならそんなことないのだけど、ポリゴンによってはmemory
not mappedで死んでしまう。
他に良いlibraryや手はない?
0397１３２人目の素数さん
2009/07/26(日) 03:51:140398１３２人目の素数さん
2009/08/15(土) 11:01:32libmpi not found. exiting...
ということなので、どのパッケージか調べてみると
$ apt-cache search libmpi |grep dev
libmpich-mpd1.0-dev - mpich static libraries and development files
libmpich-shmem1.0-dev - mpich static libraries and development files
libmpich1.0-dev - mpich static libraries and development files
3種類もあった。どれを入れたらよいのか判断がつかないので、助言をお願いします。
0399１３２人目の素数さん
2009/08/16(日) 11:39:37どうか手助けをお願いします。
glm()の前にtableで度数分布表を作ると度数0のマスがあってオッズ比がぶっ飛んでしまう変数が出てしまいます。
wikiとかですとサンプルを増やせ、デザインが悪いってことになるでしょうが、サンプルは現実的には増やせません。
そこでイェーツの補正みたいなことはglm()では出来ないでしょうか？
独立変数は2つ(仮にI,II)、それぞれ3つのカテゴリデータ(仮にA,B,C)です。IとIIの交互作用?といいますかjoint odds ratioを出したいのですが、
つまり最終的には、
IA:IIA OR=1.0(Reference)
IA:IIB OR=
IA:IIC OR=
IB:IIA OR=
.
.
IC:IIC OR=
といった結果を得たいと考えております。
試しに
x <- glm(Y ~ I+II+I:II,,,,)
y <- glm(Y ~ I*II,,,,)
をおこなったところ、出力される項目が異なります。
xだとIA:IIB, IA:IIC, IB:IIA, IC:IIAが出てきませんし、
yだとIA:IIAがreferenceにならないみたいです。
xの場合、IB,IC,IIB,IICはそれぞれIIAかIAで調整済みとなっているのでしょうか。
それともとても大きな勘違いをしています？。
長くなってすいませんがよろしくお願いします。
0400１３２人目の素数さん
2009/08/16(日) 18:02:130401399
2009/08/19(水) 11:01:11現在、前述のオッズ比の算出で手こずっています。Epitoolsとかいうパッケージを使えばIA:IIAをオッズ比1にして計算してくれるのですが、glmのオッズ比とは違う様子。
0402１３２人目の素数さん
2009/08/19(水) 17:56:490403１３２人目の素数さん
2009/08/19(水) 23:52:27大きく勘違いしてないか？
ロジスティック回帰なのに度数分布？
Rだと生データのままかけるんだけど。
それともやりたいのは対数線型モデル？
ちなみに対数線型モデルでも最尤法で解くんだから
０が出てきてもたくさんでなければ解ける。
0404399
2009/08/20(木) 11:28:56そもそもロジスティック回帰自体を良く分かっていないので自分でも常識的なことを間違えているんだろうなと思ってます
データは、
>df
ID 結果 要因A 要因B
001 1 0 1
002 0 0 0
003 0 0 0
004 1 1 2
005 1 2 2
.
と全部で300件くらい。
GLM<-glm(結果~要因A+要因B,family=binomial(logit),data=df)としています。
ORはexp(coef(GLM))で出ると思うのですが、このORについて理解できていません。
(結果~要因A+要因B) 要因Aの"0"を基準にした要因Aの"1"と"2"のORと要因B"0"を基準にした要因Bの"1"と"2"のORだと思うのですが、
(結果~要因A*要因B) とすると要因A1:要因B1、要因A2:要因B1、要因A1:要因B2、要因A2:要因B2が何に対するORなのか、また、
要因A0:要因B1とか要因A0:要因B2とか要因A1:要因B0とかは何で出てこないのか。
(結果~要因A:要因B) とすると全ての組み合わせが出てくるが、一体なにを基準したORなのか。
あたりが分かっていません。
で、とりあえずepitoolsとかいうパッケージを使って分割表(3次のarray?)を作ってORを計算してみているところです。
最終的には要因Aと要因Bの単点での結果に対する影響はこれくらいだけど要因Aと要因Bが組み合わされると影響はこれくらいだよというのが言いたいです。
分かりにくくてすいません、こんがらがってきました。
0405１３２人目の素数さん
2009/08/21(金) 02:49:15交互作用の自由度はその例だと(3-1)x(3-1)=4なわけでだから4つしか出ません。
他の組み合わせ（たとえば要因A0:要因B1など）はすべて0と基準化されているということです。
0406１３２人目の素数さん
2009/08/22(土) 14:45:18seqでbyを0にして作ろうとするとエラーになってしまいます。
for文を使って作るしかないのでしょうか？
0407１３２人目の素数さん
2009/08/22(土) 17:53:21rep(2,5)
0408１３２人目の素数さん
2009/08/22(土) 20:25:01ありがとうございます。たすかりました。
0409１３２人目の素数さん
2009/08/23(日) 03:43:29どうしたらよいのでしょうか？
0410１３２人目の素数さん
2009/08/23(日) 20:44:24ソフトの使い方程度で質問しないでください
0411１３２人目の素数さん
2009/08/24(月) 02:12:53申し訳ありません。
初めてここ見たばかりですし、Ｒも最近始めたばかりの初心者だったもので。
もう少し自分でしっかり勉強してから出直します。
0412１３２人目の素数さん
2009/08/24(月) 12:02:17>>410はスルーの方向で。ここは便所の落書きなので何を書き込んでもよいよ。
初心者の質問やいまさら人に聞けない質問も書き込んでもいい。ただし、その
質問に回答する人が現れるかかどうかは別問題という認識は必要。
回答が欲しい質問は、ここじゃなくてRjpwikiへ。質問をするときの注意事項をよ
く吟味してから書き込むこと。
>>403
> Rだと生データのままかけるんだけど。
ダウト。生データも使えるが正しい。
0413１３２人目の素数さん
2009/08/25(火) 02:09:30だから生データがあるのにわざわざ度数分布にする必要はないだろ。
0414１３２人目の素数さん
2009/08/25(火) 03:03:17おめえ結果が>>410と全く一緒なのわかってるよなw
0415１３２人目の素数さん
2009/09/02(水) 21:30:57変数の型が途中で変えられるのに非常に違和感を感じます。
C言語なんかでは不可能です。
関数が値渡しだからかのう？？？
頭が整理されていないので、チンプンカンプンな質問になっていますが、
どなたかご説明いただけますか？
x <- c(2,2,3,1,5)
x <- factor(x)
↑numericがfactorに！
0416１３２人目の素数さん
2009/09/02(水) 22:19:50Rは型安全じゃないし動的型だから。
違和感を感じるって言われても困るから。
0417１３２人目の素数さん
2009/09/03(木) 01:08:18ググル為のキーワードを下さい。
0418１３２人目の素数さん
2009/09/03(木) 10:37:22別の変数を同じ変数名で上書きができると考えたらいいんじゃね？
x <- c(2,2,3,1,5)
y <- factor(x)
としたとき、xとyは別々の変数だけど、
x <- y
として、yをxで上書きできる。
y <- factor(x); x <- y
を１文にすると、
x <- factor(x)
となる。
>>415が言っている論法だと、
> x <- c(2,2,3,1,5)
> x <- function(i){return(i+1)}
> class(x)
[1] "function"
↑numericがfunctionに！
と言っているのと大差がない。
0419418
2009/09/03(木) 10:45:18結果と結論は同じだけど。
0420１３２人目の素数さん
2009/09/03(木) 20:23:26ttp://www.okada.jp.org/RWebRef/R-lang.jp.v110.html
これによると、Cソース上ではRのオブジェクトはSEXPRECへのポインタ*SEXP。
定義は、src/include/Rinternals.h にある、
SEXPREC {
SEXPREC_HEADER;
union {
struct primsxp_struct primsxp; （略）
} u;
} SEXPREC, *SEXP;
マクロSEXPREC_HEADERを追っていくと、全てのR上の変数をenumしている、
SEXPTYPE に行き着く。ココに変数のR上での型を入れてあるようだ。
Rの変数は全てオブジェクトらしいから、genericな型としてSEXPを使っている
のでは？generic型のやり取りになるから、全ての変数を入れ替え可能となる。
ちゃんと読んでないから知らんけど。
0421１３２人目の素数さん
2009/09/04(金) 08:42:44もしそうなら、Schemeの変数ってCのvoid *見たいなもんだから、型が変わるとかそもそも意味が無い。
(Schemeでは、値には型があるが変数には型は無い)
0422１３２人目の素数さん
2009/09/23(水) 22:55:36WindowsのGUIメニューで「パッケージのインストール」を実行すると
警告： リポジトリー http://cran.md.tsukuba.ac.jp/bin/windows/contrib/2.9 に対する索引にアクセスできません
警告： リポジトリー http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.9 に対する索引にアクセスできません
以下にエラー install.packages(NULL, .libPaths()[1L], dependencies = NA, type = type) :
如何なるパッケージも指定されませんでした
と出てきます。どういうことでしょうか・・・。
0423422
2009/09/23(水) 22:58:380424１３２人目の素数さん
2009/09/26(土) 21:39:03最適化というのは，ホワイトノイズがよりホワイトノイズらしく（ガウシアン）なるようなNを求めたいということです．
0425１３２人目の素数さん
2009/10/11(日) 11:58:53SASからの完全移行は無理だな･･･
0426１３２人目の素数さん
2009/10/12(月) 16:29:50決まってるやんか。
0427１３２人目の素数さん
2009/10/13(火) 02:07:490428１３２人目の素数さん
2009/10/13(火) 02:35:4050Gでもサクサクです。
0429１３２人目の素数さん
2009/10/14(水) 23:40:16引き続いてのシェッフェの多重比較で有意差がでなくなるのがわかりません
0430１３２人目の素数さん
2009/10/15(木) 00:20:35これも関数として定義されているからなのでしょうか？
具体的には、クラス"term"に対する"["に動作が予想外で戸惑っています。
0431１３２人目の素数さん
2009/10/15(木) 01:27:59下請けの関数に渡されます。plot()やsummary()も同様。
0432１３２人目の素数さん
2009/10/15(木) 05:51:36なるほど総称関数になっているんですね。
では、総称関数"["の内容を見るにはどうすれば良いのでしょうか？
例えば、総称関数"plot"だったら、引数を与えずに plot と入力すれば中身が見えますよね？
"["は [ と入力してもエラーになります。
0433KY
2009/10/23(金) 22:05:28行列も計算も出来るじゃん！
俺のTI-89でも出来ないでもないけど
入力が('A`)超ﾏﾝﾄﾞｸｾ
P.S., 変数削除するときの「本気ですか？」
ワロタｗｗｗ
0434学生
2009/10/24(土) 21:48:5815128 12 0 0 1
15128 6 1 0 0
10584 2 0 0 1
24352 5 0 0 4
15128 10 0 1 1
24352 1 1 0 1
10584 9 0 0 1
19742 3 0 1 3
Rを学んでいる学生です。例えばこのようなデータの場合、死者＋重傷者の事故と軽傷者の事故を分けたいのですがどのようにRで解析すれば良いでしょうか??
先生から出された宿題なんですが全くわからなくて…。見づらいと思いますがよろしくお願いします。
0435１３２人目の素数さん
2009/10/24(土) 22:38:530436１３２人目の素数さん
2009/10/24(土) 23:38:57>死者＋重傷者の事故と軽傷者の事故を分けたい
問題が理解できません。「分ける」とは、どういう意味でしょうか？
どういう結果を得たら、「分かれた」と言えるのでしょうか？
R ではなく、手計算（手集計）を行うとしたら、どうすればいいのかわかっていますか？
0437学生
2009/10/24(土) 23:43:32死者＋重傷者を１軽傷者をー１という感じに分けたいんです。
どうにかなるものなのですかね??
0438学生
2009/10/24(土) 23:46:150439KY
2009/10/25(日) 04:35:410440KY
2009/10/25(日) 04:57:47データ・フレーム名がtrafficAccidentsなら
trafficAccidents[c(3,4)]
と
trafficAccidents[c(3,5)]
でいいんじゃね？
0441１３２人目の素数さん
2009/10/25(日) 23:40:030442１３２人目の素数さん
2009/10/26(月) 03:02:19gplots
gdata
gmodels
gtools
…の四つのパッケージをインストールする必要がある、と書かれています。
しかし、その四つを確実にインストールしているのに
> plotmeans(w ~ g, data=mydata3)
エラー： 関数 "plotmeans" を見つけることができませんでした
というメッセージが出ます。
他に何かインストールしないといけないパッケージがあるんですか？
検索しても出てきません。
0443442
2009/10/26(月) 03:13:39インストールしてたんですけど読み込んでなかったです
どうか忘れてください
0444１３２人目の素数さん
2009/10/28(水) 07:30:06一般化加法モデルとか
0445１３２人目の素数さん
2009/11/06(金) 21:55:27A <- 1
を
A=1
って書いたらスクリプトを公開するときに、恥ずかしいですか？
0446１３２人目の素数さん
2009/11/06(金) 21:56:20出力のところに出てくる
[1] って、なんか意味あるんですか？ 邪魔なんですけど。
0447１３２人目の素数さん
2009/11/07(土) 17:01:12「各言語の流儀にすら合わせる事が出来ない、柔軟性や学習意欲に欠けた人間です」と自己紹介することになるだけだから、大した事無いでしょ。
0448１３２人目の素数さん
2009/11/07(土) 22:13:45ひどいｗｗ
0449１３２人目の素数さん
2009/11/08(日) 09:44:59plot(x,y, col='red')
を
plot(x,y, col <- 'red')
と書けないのはなぜ？
0450１３２人目の素数さん
2009/11/08(日) 16:40:49<-はインスタンスを生成すると言う、明確なメッセージが有る。
お前さんの例は、インスタンスの属性の指定だから。
そこら辺を意識して書き分けられないのは、オブジェクト指向性が分かってないと言う事。
0451１３２人目の素数さん
2009/11/08(日) 20:36:57その <- っていうのは S の用法なんですかね？
Rを作るときにやっぱ＝も使えるようにしちゃえーとか。
なんでもいいけど、Rっていうネーミングは、検索かけると違うものがいっぱいひっかかるので
ちょっと考えてほしかったぞｗ
0452１３２人目の素数さん
2009/11/15(日) 10:48:28> なんでもいいけど、Rっていうネーミングは、検索かけると違うものがいっぱいひっかかるので
> ちょっと考えてほしかったぞｗ
本当に。
Rって、名前付けのセンスが悪いから、ずっと無視してきた。
けど、最近Scilabから移行するか思案中。
最近は、信号処理分野でも、統計っぽいことをすることが増えてきたので。
0453１３２人目の素数さん
2009/11/15(日) 17:36:53Rの残念なところは、Scilabと違って、
>A=1+2
と打っても
3
と答え出してくれないので、メモリ電卓代わりには使えないって事ですよねｗｗ
まぁそのまま 1+2 って打てばいいんですが。
scilabなら
A=1+2;
はエコーバックしない。
A=1+2
はエコーバックする、と区別がついていたので。
信号処理するなら、微分方程式とか状態変数とかラプラス変換とか出来ないと辛いんじゃないですか？
Rで出来ましたっけ？
0454１３２人目の素数さん
2009/11/16(月) 21:55:14optimを使って，目的変数を最大化（最小化）する計算を速くしたいのです．
4コアより8コアにしたほうが計算は速いですか？
それとも，クロック数に依存するのでしょうか？
教えてください．
0455１３２人目の素数さん
2009/11/22(日) 04:30:58ttp://www.amazon.co.jp/dp/4274067831/
0456１３２人目の素数さん
2009/11/22(日) 06:46:12目的が言語なのか
処理なのか
中途半端ともいえる
0458１３２人目の素数さん
2009/11/26(木) 12:49:40ttp://brainarray.mbni.med.umich.edu/brainarray/rgpgpu/
話が出ないですけど、どんなもんでしょ？
0459１３２人目の素数さん
2009/11/27(金) 19:23:28どうやったら確かめられる？
0460１３２人目の素数さん
2009/11/30(月) 01:22:07単精度と倍精度で誤差の出るような演算をしてみるとか？
0461１３２人目の素数さん
2009/11/30(月) 03:02:51> (A <- 1 + 2)
[1] 3
0462１３２人目の素数さん
2009/11/30(月) 23:26:36どうにかならないものでしょうか…
0463１３２人目の素数さん
2009/12/01(火) 00:08:08ペースト後、
WordのフォントをMS ゴシックに変えてみたらどうですか？
0464１３２人目の素数さん
2009/12/01(火) 12:54:53ありがとうございます！
かなり悩んでいたのにあっさり解決してびっくりです。
0465１３２人目の素数さん
2009/12/01(火) 23:27:07信頼区間の幅が0.06以下になるには何個抽出すればいいですか？
0466１３２人目の素数さん
2009/12/02(水) 22:32:20初めての板だからさっきは別のスレに書いてしまいました。
ＳＡＳユーザーです。最近は東大もＳＡＳ辞めてしまいましたが、
Ｒって使い勝手はどうですか？
0467１３２人目の素数さん
2009/12/03(木) 01:00:510468１３２人目の素数さん
2009/12/03(木) 02:21:280469１３２人目の素数さん
2009/12/03(木) 11:26:47検索がまるで役に立たない検索結果ばかりで泣きそうになる･ﾟ･(ﾉД`)･ﾟ･
0470１３２人目の素数さん
2009/12/03(木) 19:18:59かな
0471１３２人目の素数さん
2009/12/03(木) 21:46:31R 統計
でOK
0472１３２人目の素数さん
2009/12/05(土) 15:55:58CRAN-RとかCRANRなら、被らないと思うんでそっちに統一して欲しいです
0473１３２人目の素数さん
2009/12/07(月) 19:38:19って何？
某インストーラーでlib=に指定されてたけど、そのまま試したら「そんなもの無い」って怒られた。
lib=指定をなくして再度インストーラー試したら
rgument 'lib' is missing: using '$HOME/R/x86_64-pc-linux-gnu-library/2.8'
ってなった。無事にインストールされたか心配。
0474１３２人目の素数さん
2009/12/08(火) 00:13:170475１３２人目の素数さん
2009/12/10(木) 00:41:34> /Library/Frameworks/R.framework/Versions/2.6/Resources/RLib64
MacOSXだね
> rgument 'lib' is missing: using '$HOME/R/x86_64-pc-linux-gnu-library/2.8'
今度はLinux？
これでWindowsユーザなら笑ってしまう
0476１３２人目の素数さん
2009/12/20(日) 03:49:050477１３２人目の素数さん
2009/12/21(月) 12:48:27もちろんいいよ
0478１３２人目の素数さん
2009/12/24(木) 18:17:15> expression(C^o)
expression(C^o)
> expression(^oC)
エラー： 予想外の '^' です ( "expression(^" の)
superscriptの直前に文字が必要みたい
0479１３２人目の素数さん
2009/12/24(木) 21:07:460480１３２人目の素数さん
2009/12/24(木) 21:12:32どこのどなたか存じ上げませんが、誠にありがとうございました。
解決しました。
0481１３２人目の素数さん
2010/01/01(金) 23:20:22コンソールのフォントを変えても、立ち上げ直すたびに貧弱なフォントに
戻ってしまいます。
どうすればよいのかご教授ください。
0482１３２人目の素数さん
2010/01/09(土) 23:09:51> 信号処理するなら、微分方程式とか状態変数とかラプラス変換とか出来ないと辛いんじゃないですか？
> Rで出来ましたっけ？
信号処理に使っているけど、Scilabに微分方程式とかラプラス変換をやらすことはないなあ。
最終的にデジタルシステムに落とすのが目的なので、
微分方程式を解くとかはあらかじめMapleとかでやれば十分。
って、Rでも odesolveライブラリ使えば、微分方程式解けるんじゃないの？
0483１３２人目の素数さん
2010/01/12(火) 19:21:093年ぐらい前に初めてRをインストールした時、同じ症状だったなぁ。
あきらめて Terminal から R を呼んで使ってるよ。
0484１３２人目の素数さん
2010/01/12(火) 19:29:21何か初期設定ファイルとかあるんじゃないの。でも、コンソールのフォントにこだわりがあるなんて
さすがMacゆーざって感じですな。グラフィックディバイスのフォントは高品位である必要があるので、
とことんこだわるべきだが。
ちなみに当方は、MacOSX Leopardでemacs+essからRを使っているから、コンソールのフォントは、
emacsの方で設定できる。
0485１３２人目の素数さん
2010/01/14(木) 00:43:52質問させてください。今130万×130万くらいの行列があります。
ただ要素は、ほとんどが0で一行に数字が入っているのは多分32までで、
負の数もありません。
これでオライリーの「集合知プログラミング」10章の非負値行列因子分解をしたいのですが、
あまりにも大きいのでとりあえずデータベース板に行きましたが、統計ソフトR
という名前が出てきたので、ここに来ましたが、SparseMというのを使えばいいんでしょうか。
何も分からない素人なのでどこから手を付ければいいか教えてください
0486１３２人目の素数さん
2010/01/14(木) 10:24:21NMFNパッケージ
> library(NMFN) #NMFNパッケージの読み込み
> X <- matrix(1:12,3,4) #サンプルの行列を作成
> X
[,1] [,2] [,3] [,4]
[1,] 1 4 7 10
[2,] 2 5 8 11
[3,] 3 6 9 12
> ss <- nnmf(X,3) # 3 factors
Multiplicative Update Algorithm
Iter = 50 relative error = 0.001003956 diff = 52.33016 eucl dist = 0.0009181338
Iter = 100 relative error = 0.0008185111 diff = 0.01446467 eucl dist = 0.0006126808
Iter = 150 relative error = 0.0006704387 diff = 0.01154965 eucl dist = 0.0004113811
[snip]
Iter = 1000 relative error = 2.457965e-05 diff = 0.0004078737 eucl dist = 5.550067e-07
> sx <- nnmf(X,3,'nnmf_als') # 3 factors
Alternating Least Squares Algorithm
Iter = 50 relative error = 8.113168e-17 diff = 35.19087 eucl dist = 7.999543e-30
Execution finishes at iteration = 50
130万x130万行ができるかどうかは、あなたの環境次第かもしれない。
Windowsだとメモリの限界がすぐに来るかもしれないので、LinuxやMacOSXの方がよいかも。
0487１３２人目の素数さん
2010/01/14(木) 18:52:09レスありがとうございます。そのパッケージが標準で付いてくるものなのか、
別途インストールが必要なものかググっても分からなかったのですが、
とりあえずRをubuntuにインストールしようとしましたが以下のエラーが出ました。
$ sudo apt-get install r-base
パッケージリストを読み込んでいます... 完了
依存関係ツリーを作成しています
状態情報を読み取っています... 完了
インストールすることができないパッケージがありました。おそらく、あり得
ない状況を要求したか、(不安定版ディストリビューションを使用しているの
であれば) 必要なパッケージがまだ作成されていなかったり Incoming から移
動されていないことが考えられます。
単純な操作を行っただけなので、このパッケージは単にインストールできない
可能性が高いです。そのため、このパッケージへのバグレポートを送ってくだ
さい。
以下の情報がこの問題を解決するために役立つかもしれません:
以下のパッケージには満たせない依存関係があります:
r-base: 依存: r-base-core (>= 2.10.1-2karmic0) しかし、インストールされようとしていません
依存: r-recommended (= 2.10.1-2karmic0) しかし、インストールされようとしていません
E: 壊れたパッケージ
Synapticも試しましたが、やはり依存関係でエラーが出ます。ubuntuは8.04なのですが、
これが古いということでしょうか。
0488１３２人目の素数さん
2010/01/14(木) 20:37:36Ubuntuは初めてですか、お客さん？
karmicはKarmic Koala (Ubuntu 9.10)のこと。で、あなたのUbuntuはHardy Heron。
apt-lineをkarmicに向けちゃったら、そらエラーになる。karmicって書いちゃったところを
hardyに書き直すとどうだい？
0489１３２人目の素数さん
2010/01/14(木) 20:50:36Linuxでは、追加パッケージを入れるとき、aptでインストールできるもの以外は、Rの中で、
> install.packages("NMFN")
という感じでやるのだけど、インストール先が~/でよいならユーザ権限でRを起動しても良いけど、
私なら、
$ sudo R
とroot権限でRを起動してパッケージを入れる。ネットワークインストールがいやなら、
自分でダウンロードしてきて、
$ sudo R CMD INSTALL hoge.tar.gz
としてもよい（たぶん）。
あと、重要なことだけど、追加パッケージは、ソースの形でダウンロードされ、そのままコンパイルされる。
コンパイルできる環境や必要な外部ライブラリの開発ファイルをあらかじめ入れておく必要がある。
他のディストリビューションならほぼ問題ないが、Ubuntuは開発環境を自分で入れる必要があったと思う（たぶん）。
0490１３２人目の素数さん
2010/01/15(金) 00:34:18よくわからんのだが、130万×130万の行列って、１要素に8バイト使うと、
1.3e6*1.3e6*8/1e12=13テラになるんだが、それでいいのか？
0491１３２人目の素数さん
2010/01/16(土) 06:51:20>>488
すみませんレスに気が付かずVer.UPしてしまいました。その途中で画面に出てる文字でそのことに気が付きました。
>>489
とりあえずinstall.packages("NMFN") とやって見てエラーメッセージを読んで考えようと思ったら通ってしまったようです。
>>490
そうですね。もとの行列は0のところは何も書かないで、x行,y列の値をx,y,値とファイルに書いているだけなので
そんなにかさばらないのですが、非負値行列因子分解をする過程で、0以外の値の入る要素が増える可能性はあるだろうし、
考慮する必要がありますね… 13テラのHDは今準備できませんし。
0492１３２人目の素数さん
2010/01/16(土) 12:47:19その名もずばりSparseMなるパッケージもあるみたい。
使ったことが無いから、どんなものかは知らないけど。
0493１３２人目の素数さん
2010/01/16(土) 20:43:430494485
2010/01/17(日) 06:38:11でSparseMは疎行列をあつかうパッケージなのですよね。SparseMを使うと自分で
非負値行列因子分解のアルゴリズムを実装する必要があるというわけですよね。
0495１３２人目の素数さん
2010/01/17(日) 11:23:24wilcox.testで、「'x'は数値でなければなりません」になります。
データは、data<-read.csv("data.csv",header=TRUE) で読込み、
BとCの項目を、wilcox.test（ data["B"],data["C"])で渡しています。
データ数はB,Cとも400程度で、5桁以下の整数値のみです。
どのように渡せばいいにでしょうか？
0496495
2010/01/18(月) 09:52:21mode()でみると、listになっています。
as.numeric()を掛けると、(list)オブジェクトは'double'に変換できません。で
as.integer（）だと、(list)オブジェクトは'integer'に変換できません。になります。
数値と指定して読み込むには、どうしたらいいのでしょう？
0497１３２人目の素数さん
2010/01/18(月) 10:10:56自助努力をしているみたいなので、ヒントを差し上げよう。
> dat <- data.frame(B=1:5,C=6:10)
> dat
B C
1 1 6
2 2 7
3 3 8
4 4 9
5 5 10
> class(dat["B"])
[1] "data.frame"
> class(dat[,1])
[1] "integer"
> wilcox.test(dat[,1],dat[,2])
Wilcoxon rank sum test
data: dat[, 1] and dat[, 2]
W = 0, p-value = 0.007937
alternative hypothesis: true location shift is not equal to 0
0498１３２人目の素数さん
2010/01/18(月) 11:06:48有難う御座いました。
Vector?と思い込んでいましたが、
array?だったのですね。
助かりました。
0499１３２人目の素数さん
2010/01/18(月) 19:22:44では上に表示されます…どうにか対処できるでしょうか？
あとグラフのlabel内で、平方メートルの単位はm^2と表記するしか
ないですよね…最近Ｒの不便さを感じます。
0500１３２人目の素数さん
2010/01/18(月) 19:44:44表のタイトルが下？？図じゃなくて？Rで表を作ることもあるけど、出力してから
普通自分で編集しない？
Rのグラフィックを論文投稿用に使うとき、基本的にmainは使わない。
講義用とかならmainを使った方がよいけど。
それから、m^2は文字通り「m^2」と書くという意味？なぜ、superscriptを使わないの？
0501１３２人目の素数さん
2010/01/18(月) 19:55:50> d2 <- 0.45
> plot(rnorm(100), main=bquote(Family ~ Gaussian ~ ";" ~ D^2 == .(d2)), ylab=expression(m^2))
0502１３２人目の素数さん
2010/01/18(月) 23:13:08wilcox.test(data$B, data$C) で桶
どうしても ["B"] とかを使いたいなら、
wilcox.test(unlist(data["B"]), unlist(data["C"]))
理由は、
> class(dat["B"])
[1] "data.frame"
> class(dat$B])
[1] "integer"
だから
>>499
main 使わずに sub 使えばいいんじゃね
0504１３２人目の素数さん
2010/01/26(火) 22:37:14どこに何をいれたか忘れてしまったためです。
0505１３２人目の素数さん
2010/01/27(水) 00:02:18ls() じゃダメ? 変数の中身は見れないけど
0506１３２人目の素数さん
2010/01/27(水) 00:13:45essを使っているならess-rdired.elとか
他にも似たような機能を持つRのGUIアプリはいくつかあったような。
0507１３２人目の素数さん
2010/01/27(水) 19:56:05ありがとうございます。できました。
0508１３２人目の素数さん
2010/01/29(金) 20:27:41その曲線を任意の範囲での積分（つまり、確率）
ってできますか？
0509１３２人目の素数さん
2010/01/31(日) 14:20:110510１３２人目の素数さん
2010/02/01(月) 09:20:19コンピュータ言語じゃないという話もあるから 「GNU R」 で検索させてくれ
0511１３２人目の素数さん
2010/02/01(月) 10:30:11できる/できないの話ならできるんだろうけど、面倒くさいな。
任意の範囲をマウスで指定するならもっと面倒。
0512１３２人目の素数さん
2010/02/01(月) 17:44:58どうやってやるんですか？近似曲線を求めるとか？
任意の範囲はマウスではなく数字指定で十分です。
0513１３２人目の素数さん
2010/02/01(月) 18:07:08> plot(density(c(-20,rep(0,98),20)), xlim = c(-4,4))
> x <- density(c(-20,rep(0,98),20))
> str(x)
List of 7
$ x : num [1:512] -23.1 -23 -22.9 -22.8 -22.7 ...
$ y : num [1:512] 4.46e-05 5.77e-05 7.40e-05 9.45e-05 1.20e-04 ...
$ bw : num 1.02
$ n : int 100
$ call : language density.default(x = c(-20, rep(0, 98), 20))
$ data.name: chr "c(-20, rep(0, 98), 20)"
$ has.na : logi FALSE
- attr(*, "class")= chr "density"
上の例で言えば、矩形の面積を511個計算して合計すると全体、部分が欲しければ、
その部分の矩形の面積を合計。0.1刻みじゃ嫌だ！というなら、もっと細かくすればよい。
0514１３２人目の素数さん
2010/02/01(月) 23:46:29多くの場合、半角スペースを検索できない
0515１３２人目の素数さん
2010/02/02(火) 07:41:25なるほど！density自体にそんな結果が含まれていたのですね。
plotで使ってただけだから気付かなかった。
ありがとうございます。大変助かりました。
0516１３２人目の素数さん
2010/02/02(火) 23:39:17横軸を年と月がざっくり判るような目盛でplotするにはどうすれば
0517１３２人目の素数さん
2010/02/03(水) 10:05:48example(format.Date)
example(axis.Date)
0518１３２人目の素数さん
2010/02/09(火) 22:03:33計算することができますか？
0519１３２人目の素数さん
2010/02/09(火) 23:08:030520１３２人目の素数さん
2010/02/10(水) 00:09:27金融商品だけに魔法じゃなくて錬金術、とか。
0521１３２人目の素数さん
2010/02/10(水) 13:23:19「できる/できない」を問う質問は不毛だよ。
>>519みたいな答えや「できます」という答えが返ってくるだけ。
0522１３２人目の素数さん
2010/02/15(月) 17:17:32http://science6.2ch.net/test/read.cgi/math/1265261354/
STATAのスレ立てた人なにスレ沈めてんだ
盛り上げましょうって宣言したじゃないか
0523１３２人目の素数さん
2010/02/16(火) 18:58:130524１３２人目の素数さん
2010/02/17(水) 17:31:340525１３２人目の素数さん
2010/02/18(木) 17:14:25まじめに答えてみました。
標準誤差(standard error)seを求める．
標準誤差は，標本平均の標準偏差(standard deviation)のことを意味します．
中心極限定理: 平均μ，標準偏差σが決まる母集団からであれば，どんな分布をした母集団からでも，
大きさnの無作為標本に基づく標本平均は，nが十分に大きいとき，平均μ，標準偏差(標準誤差)σ/√nの正規分布に近づく.
中心極限定理より標本平均は，平均μ，標準偏差(標準誤差)sea=σ/√nの正規分布に従うことが分かります．
しかし，母平均, 母分散は不明なため次のように推定します．
母平均μの推定=標本平均m
母分散σの推定=不偏分散(unbiased variance)uv
不偏分散の公式より
> uv <- (m ? DATA$a)^2 ) / (length(DATA$a) - 1) #length(DATA$a)は標本数．
se=√(uv/n)より
> (se <- sqrt(uv/length(DATA$a)) #sqrt() は()内の正の平方根を計算する．sqrtは，square rootの略．
Rに用意された標準偏差を求めるsd()コマンド(関数)を用いて
> (se <- sd(DATA$a)/sqrt(length(DATA$a))
0526１３２人目の素数さん
2010/02/18(木) 17:21:09スレを立てた本人です。
済みません。
簡単に沈んでしまうものなんですね。
２ｃｈ初心者なもので。
今のところSTATAをいじっている余裕がないので
またいじり始めたら建て直すと思います。
ちなみに、
525=526で、
？の文字化けは-(マイナス)です。当たり前ですが。
0527１３２人目の素数さん
2010/02/18(木) 17:43:21・スレッドのレス数が10以下
・レスが一週間書き込まれていない
二つの条件が満たされるとスレッドが沈むので
そうならないようにしておけば沈まないよ
0528１３２人目の素数さん
2010/02/18(木) 19:23:38ありがとう. 解決しました.
0529１３２人目の素数さん
2010/02/18(木) 20:58:180530１３２人目の素数さん
2010/02/18(木) 21:00:51また、どのような処理を行っているのかを調べる方法は、
CやFortranで書かれたプログラムを解析する以外にあるのでしょうか？
0531１３２人目の素数さん
2010/02/20(土) 07:40:520532１３２人目の素数さん
2010/02/23(火) 16:22:27具体的なデータは、
「数字,10000文字の日本語文字列,数字,数字」
の４列で百万行程度あります。
解析したいのは、１列目と３列目です。
２列目の文字列の容量が大きすぎるため、
Ｒに単純には取りこめません。
そのため、２列目を除いて、１列目と３列目だけＲに取り入れるいいアイデアが
あったら、ご教授いただけないでしょうか？
（私が考え付くのは、１００行づつ取り入れる２行目を除去。これを一万回繰り返し、
最後に併合くらいです。）
普通だったら、
（１）awkやcutで一列だけ切り取る。
（２）切り出したものをＲで解析する。
という手順でおこないたいのですが、
解析環境が「windows＋Ｒ」のみで、
新しいソフトのインストール禁止。
データの取り出しも禁止。
という環境化なので、
Ｒですべてを処理するしかない状況です...。
0533１３２人目の素数さん
2010/02/23(火) 20:26:22それだけ厳しい条件だと、>>532の方法くらいしか思いつかない。
0534１３２人目の素数さん
2010/02/23(火) 20:43:13R sed:lineでエディタの現在行を送ろうしてしても、Rguiの方が一行改行されるだけ。
Tinn-R 2.3.4.4
R 2.10.1
WinXP
何が足りないのだろう。Rprofile.siteもTinn-Rのメニューから作成・保存したのに。
0535１３２人目の素数さん
2010/02/23(火) 20:48:15でどう？
１００万行までできるかわからんけど
0536１３２人目の素数さん
2010/02/23(火) 20:52:50awkならインストールしなくてもフラッシュメモリなどから読み込めばいいのでは。
あと、jscriptやvbscriptとかcmdとか。
Rってデータを分割して読み込む方法あるのかな。
0537１３２人目の素数さん
2010/02/23(火) 23:20:57データベースから読めばええやん
0538532
2010/02/24(水) 01:37:26ありがとうございます。
>>536 javascriptを使う発想はなかったです。頑張ればブラウザでデータ整理できますね。
今後、もう、そのアイデアいざって時につかわさせていただきます（ＲなくてもＯＫだし）。
>>535 ベスト回答。この答えを待っていました！ 美しいし短いし最高です。
scanもいじってたんですが、NULLまではたどりつけなくてあきらめてました。
先ほどまで、私は苦肉の策でC言語かよって感じに
file<-file("test.txt","r")
s<-0
d1<-0
d2<-0
i<-1
repeat{
s<-strsplit(readLines(f,1),",")
d1[i]<-as.numeric(s[[1]][1])
d2[i]<-as.numeric(s[[1]][3])
i<-i+1
}
とかしてましたが、今後は>>535方式で行こうと思います。
2chやっぱりすごいですね。ありがとうございました！
0539532
2010/02/24(水) 02:58:53535の方法で132万行×126列のデータを２列のみ３分程度でＲに読み込めました。
ありがとうございました。
0540１３２人目の素数さん
2010/02/25(木) 21:31:50例えば、タネを撒いた時刻は出芽数に影響を与えるか、っていう感じです
これって、どうにかRで扱えないでしょうか？
0541１３２人目の素数さん
2010/02/26(金) 00:12:300542１３２人目の素数さん
2010/02/26(金) 12:40:32Rを使う/使わない以前に、分析の原理・手法を検討するのが先。
bud~Poisson(mu)
log(bud)=alpha+beta*s(time,k)+e
とか、まずモデルを検討して、検討した結果、採用候補の分析手法を
実行するRライブラリ/コマンドが存在するかどうか尋ねられたらいかがでしょうか。
0543１３２人目の素数さん
2010/02/26(金) 16:25:21方向の関数としてのカウントデータって
なんかいい解析法ないですかね？
0544１３２人目の素数さん
2010/02/26(金) 16:27:220545１３２人目の素数さん
2010/02/27(土) 07:48:34適当にロジスティック関数を使えばいいのではないかと
0546１３２人目の素数さん
2010/02/27(土) 21:41:36Tinn-Rの1.19.4.7ならOK。
Tinn-R：2.ｘ以降はWindowsでは、動かなかった。
原因は不明・・。
誰か知ってたら教えて・・。
0547１３２人目の素数さん
2010/03/04(木) 15:41:220548１３２人目の素数さん
2010/03/04(木) 19:01:28汚くしようと思えば汚くできるし、綺麗にしようと思えば綺麗になる。そんなのあなた次第。
画像のサイズ（解像度）をコントロールする方法と、ベクタ画像とラスタ画像の違いを勉強しなさい
0549１３２人目の素数さん
2010/03/06(土) 09:18:060550１３２人目の素数さん
2010/03/06(土) 18:54:510551１３２人目の素数さん
2010/03/09(火) 01:27:59ありがとうございます！
ラスタ画像って、ちょっと難しい言葉だけど、ググッテ見たらなんとなく意味が分かりました！
ラスタ画像はコピペッした後で、画像サイズをコントロールすると、ピクセルとピクセルの間が補完されてしまうから、汚くなっちゃうんですね！
今度から、ベクタ画像とラスタ画像の違いを勉強して使います！！
それにしても、ベクタ画像とラスタ画像って、まさに僕が知りたい情報だったんですが、548さんは、どうして僕の気持ちが分かったのでしょうか？
まさにエスパーですね。今度から548さんをエスパーイシヤ人と呼ばせてください！
２チャンネルはすごいですね！
感動した！！
0552１３２人目の素数さん
2010/03/09(火) 08:42:490553１３２人目の素数さん
2010/03/10(水) 11:56:360554548
2010/03/10(水) 13:40:02>>551
そこまで調べたなら、次はベクタ画像を出力するRのグラフィックディバイス
例えば、ps(), pdf(), svg(), cairo_pdf(), cairo_ps()あたりを使って一度ファイルに
出力して、emfに変換すればよいと思うよ。
svg(), cairo_pdf(), cairo_ps()は私の環境では日本語不可でした。
0555１３２人目の素数さん
2010/03/10(水) 15:46:420556１３２人目の素数さん
2010/03/12(金) 07:28:02どのぐらいすごいのでしょうか
0557１３２人目の素数さん
2010/03/12(金) 10:52:58普通数万から１０万円超えるような統計ソフトと同等かそれ以上の機能を持ちながら無料
これ最強
0558１３２人目の素数さん
2010/03/12(金) 12:34:01"a_5"
から変数
a[5]
の値を返す関数みたいなものはないのでしょうか
0559１３２人目の素数さん
2010/03/12(金) 15:37:59あなたはだまされています。もしくは数式処理と統計処理を混同しています。
>>558
そんな便利な関数はない。自分で書く。
下記にヒントをしるす
> a <- month.abb
> b <- "a_5"
> s <- strsplit(b,'_')
> s
[[1]]
[1] "a" "5"
> eval(parse(text=paste(s[[1]][1],'[',s[[1]][2],']')))
[1] "May"
0560１３２人目の素数さん
2010/03/12(金) 15:47:10->
=
で何か違いはあるのでしょうか？
0561１３２人目の素数さん
2010/03/12(金) 18:05:050562１３２人目の素数さん
2010/03/12(金) 20:48:16オレも
M = read.csv( .... )
って書きたいんだが、周りの無言の圧力により
M <- read.csv( .... )
って、仕方なく書いている。
0563１３２人目の素数さん
2010/03/12(金) 21:41:22Rのパッケージャによっては、
M = read.csv( .... )
を使う人もいるし、また、R関連本（洋書）を見れば、わりと"="派な人がいる。
私自身は、関数のオプションと区別するためや
"=="と混同しないように、代入には <- を使う。
0564１３２人目の素数さん
2010/03/12(金) 23:25:18internはないけどevalはあるのですね
ありがとうございました
0565１３２人目の素数さん
2010/03/13(土) 05:14:450566１３２人目の素数さん
2010/03/13(土) 10:28:030567１３２人目の素数さん
2010/03/13(土) 10:51:02>>560
0568１３２人目の素数さん
2010/03/13(土) 23:55:11The R Tips第2版（オーム社）より
> for (i in 1:10) {
+ assign(paste("x", i, sep=""), i) # 変数xiにiを代入
+ }
> x1
[1] 1
> x2
[1] 2
0569１３２人目の素数さん
2010/03/14(日) 16:11:000570１３２人目の素数さん
2010/03/15(月) 17:03:31[,20100101] [,20100101]
[,user1] 0 1
[,user2] 2 3
> x <- read.table(file)
> x
X..20100101. X..20100101..1
[,user1] 0 1
[,user2] 2 3
列のヘッダに付いてる X.. は何なのでしょう。。。
0571１３２人目の素数さん
2010/03/15(月) 17:16:250572１３２人目の素数さん
2010/03/15(月) 17:45:19[user1,] 0 1 3 2
[user2,] 2 3 4 0
↓
X..ab. X..cd. X..ef. X..gh.
[user1,] 0 1 3 2
[user2,] 2 3 4 0
うーむ。。。
[R.app GUI 1.31 (5537) i386-apple-darwin9.8.0] なんですが、
列のヘッダはこれが標準なんでしょうか。。。
0573１３２人目の素数さん
2010/03/15(月) 18:27:54Rコンソールウィンドーの一番上に出てくるのが、本当の R のバージョン。
でな、ピリオドも同じ。特殊？文字や数値からは始められない。
プログラム言語なので、そういうモノ。
0574１３２人目の素数さん
2010/03/15(月) 20:15:500575１３２人目の素数さん
2010/03/20(土) 12:28:050576１３２人目の素数さん
2010/03/21(日) 15:50:520577１３２人目の素数さん
2010/03/21(日) 21:22:13正にjordan標準形を求めたいのです。
0578１３２人目の素数さん
2010/03/23(火) 10:37:07library(rms)
> (a <- diag(1:3))
[,1] [,2] [,3]
[1,] 1 0 0
[2,] 0 2 0
[3,] 0 0 3
> -matinv(a.inv1, 2:3, negate=FALSE)
[,1] [,2] [,3]
[1,] 1 0.0 0.0000000
[2,] 0 0.5 0.0000000
[3,] 0 0.0 0.3333333
attr(,"rank")
[1] 2
attr(,"swept")
[1] TRUE TRUE TRUE
0579１３２人目の素数さん
2010/03/23(火) 12:51:05http://atnd.org/events/3454
0580１３２人目の素数さん
2010/03/23(火) 12:53:13R 掲示板
http://www.okada.jp.org/RWiki/?R%B7%C7%BC%A8%C8%C4
0581578
2010/03/23(火) 13:11:35> library(rms)
> example(matinv)
>>580
こんなところにRjpwikiのurlを貼らないで。
0582１３２人目の素数さん
2010/03/23(火) 13:21:450583１３２人目の素数さん
2010/03/23(火) 14:42:25IBM?
individual-based modelぐらいしか思いつかない。
Rは、商用パッケージと違ってユーザ登録がないので、
ユーザの専門分野構成比について、誰も客観的な答えが出せない。
0584580
2010/03/24(水) 11:49:59なんで？
Spamが来るから？
このサイト、書きこみにはパスワードがかかってるよね？
0585１３２人目の素数さん
2010/03/24(水) 23:49:21平均600、標準偏差100の、対数正規分布に従う乱数を発生するにはどうすればよいのでしょうか？
> log( rlnorm( 100 , meanlog=600 , sdlog=100 ) )では、Infが多発してしまいます。
0586１３２人目の素数さん
2010/03/25(木) 00:15:21> log( rlnorm( 100 , meanlog=600 , sdlog=100 ) )では意味合いがまったく違いますね。
アドバイスお待ちしています。
0587１３２人目の素数さん
2010/03/25(木) 00:50:43> mean(log(rlnorm(100,600,100)),trim=1000)
[1] 618.7184
> mad(log(rlnorm(100,600,100)))
[1] 94.40666
0588１３２人目の素数さん
2010/03/25(木) 01:09:28ご回答ありがとうございます。
しかし、log(rlnorm(100,600,100))だと、
対数をとった結果が、平均600、標準偏差が100の正規分布に従うという意味ですよね？
私が求めたいのは、平均600、標準偏差が100に従う対数正規分布そのものです。
（やりたいことは、所得の乱数を発生させて各種シミュレーションを行うことです。
所得の分布に対数正規分布を仮定しています。）
引き続きアドバイスをお願いします。
0589１３２人目の素数さん
2010/03/25(木) 08:27:09とりあえず「データが対数正規分布で、もとの値での平均値を問題」をキーワードに検索してみれば？
0590１３２人目の素数さん
2010/03/27(土) 19:05:10普通の正規分布に従う乱数の発生はできますか？
0591１３２人目の素数さん
2010/03/29(月) 00:20:50レスありがとうございます。
検索したページを読んでみましたが･･･
γ(600)=exp(μ+σ^2/2)、
δ(100)=μ^2・{exp(σ^2)-1}
として、uとσを求めればよいのでしょうか？…良くわかりません。
>>590
rnorm(100,600,100) です。
これぐらいならば問題なく大丈夫です。
0592１３２人目の素数さん
2010/03/29(月) 12:38:31Ｒに乗り換えるといいことありますか？
pythonで出来なくてRにできることって何ですか？
0593１３２人目の素数さん
2010/03/29(月) 18:59:39ない。「Cにできなくてpythonにできること」がないのと同様。
0594１３２人目の素数さん
2010/03/29(月) 23:32:35Pythonで Rモジュールを利用することが出来るとぐぐったことがある。
Pythonを使っているのなら、調べてみてはいかが？
（Pythonは使ったことがないので、未確認ですが・・）
0595１３２人目の素数さん
2010/04/08(木) 20:06:11コンサルやMBA(笑)みたいな胡散臭い奴らが学んでるし。
理工系のエンジニアとしては彼らと一緒にされたくないので
統計学を学ぶのには嫌悪感がある。
0596１３２人目の素数さん
2010/04/09(金) 00:10:46ん？>>595は単に文系レベルの統計学しか知らないという事では？
0597１３２人目の素数さん
2010/04/09(金) 08:22:140598１３２人目の素数さん
2010/04/09(金) 14:16:15は他スレでも同じことを書いていたね。
きっと、コンピューターはインテリが使うものだという感じがするので嫌い。
とか言うんだろう。
0599１３２人目の素数さん
2010/04/09(金) 21:29:04理系の研究者で統計を知らない奴はいないよ
0600１３２人目の素数さん
2010/04/10(土) 00:00:57温度の定義を知ってるかい?
0601１３２人目の素数さん
2010/04/10(土) 08:22:26>>595は理系ぶってる文系が同族嫌悪してるんだから！
0602１３２人目の素数さん
2010/04/10(土) 09:46:00「物理屋」を自称していても、意外と統計力学の基礎がすっぽり
抜け落ちてる人がいるんだよな。熱力でストップしてるの。
0603１３２人目の素数さん
2010/04/10(土) 10:22:02だから、物理「屋」なんだろうね。
0604１３２人目の素数さん
2010/04/10(土) 11:02:12物理系なら統計力学は必須だろ
0605１３２人目の素数さん
2010/04/12(月) 00:43:56統計を知らない物理屋はいないだろう。
0606１３２人目の素数さん
2010/04/12(月) 15:03:350607１３２人目の素数さん
2010/04/17(土) 10:37:57【toto】サッカーくじBIGスレ part10
http://yutori7.2ch.net/test/read.cgi/soccer/1266583447/
0608１３２人目の素数さん
2010/04/17(土) 13:18:100609１３２人目の素数さん
2010/04/18(日) 08:12:280610１３２人目の素数さん
2010/04/18(日) 13:11:10最近は文系どころか理系もひどい
ttp://www.itmedia.co.jp/enterprise/articles/1004/09/news008.html
ttp://www.itmedia.co.jp/enterprise/articles/1004/09/news008_2.html
0611１３２人目の素数さん
2010/04/18(日) 15:20:02名前は知らなくていいけど、概念を説明しても納得しない
AかつBでなければならない、という法律に違反している状態は
AとBの少なくとも一方が否定されれば満たされるのに、
AでなくかつBでない状態だと信じてそのまま判決を出す裁判官は実在する
0612１３２人目の素数さん
2010/04/18(日) 16:16:29a < -3 AND 5 > a
を否定するときは
-3 <= a OR a <= 5
っていう話
0613１３２人目の素数さん
2010/04/18(日) 16:17:36a < -3 AND 5 < a
です
orz
0614１３２人目の素数さん
2010/04/18(日) 16:18:49例を挙げるなら
a < -3 OR 5 < a
を否定するときは
-3 <= a AND a <= 5
っていう話
です
吊ってくる
0615１３２人目の素数さん
2010/04/18(日) 16:22:310616１３２人目の素数さん
2010/04/18(日) 16:48:280617１３２人目の素数さん
2010/04/18(日) 20:34:10マーフィーの法則
バカと喧嘩はするな、はたからはどっちがバカだか分らない。
0618１３２人目の素数さん
2010/04/18(日) 22:50:200619１３２人目の素数さん
2010/04/22(木) 21:32:54delta2 <- function(x,y){
return(log(y/exp(2*log(x)) + 1))
}
myu <- function(x,y){
return(-0.5 * log(y/exp(2*log(x)) +1) + log(x) )
}
xに600、ｙに100^2を放り込んでみたら？
0620１３２人目の素数さん
2010/04/22(木) 21:37:52一応答え書いておくと
> mean(rlnorm(100000,6.38323,0.1655264))
[1] 599.9213
> sd(rlnorm(1000000,6.38323,0.1655264))
[1] 99.98724
0621１３２人目の素数さん
2010/04/25(日) 14:06:440622１３２人目の素数さん
2010/04/25(日) 16:34:390623１３２人目の素数さん
2010/04/27(火) 08:58:29頑張れ！R
http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1139981306
0624１３２人目の素数さん
2010/04/27(火) 15:09:23ん。エクセルの対抗馬としてRはないかな。
エクセルの対抗馬としては、グーグルドキュメントや、
オープンオフィス、Macの表計算ソフトかな。
Rではなく、アクセスなら、データベースソフトとしての対抗馬になったかしらん。
0625１３２人目の素数さん
2010/04/27(火) 15:35:54想像するに、その質問者は、
先生から、（統計解析を勉強するなら）今からエクセルをやっても確実に負ける。Rでなら勝てると言われました。
この（）の部分を聞き流したのではないだろうか。
0626１３２人目の素数さん
2010/04/27(火) 20:30:34アドインは
http://aoki2.si.gunma-u.ac.jp/Hanasi/excel/
だし、
http://www.forest.impress.co.jp/article/2007/10/11/excel2007patch.html
なんてのがあるから、統計処理にExcelは向かないね。
0627１３２人目の素数さん
2010/05/03(月) 22:25:39という本を読めば？
0628１３２人目の素数さん
2010/05/07(金) 00:11:170629１３２人目の素数さん
2010/05/09(日) 16:36:37コンソールに出力された情報をそのままテキストファイルに出力するのがかなり面倒なんですね。
sinkだと計算結果のみで、命令文は出力されない。また、テキストファイルと、コンソール画面同時には表示されない。
capture.outputだと表とかは崩れて出力されてしまう。
http://www.okada.jp.org/RWiki/?R%BD%D0%CE%CF%A4%CE%B5%AD%CF%BF
ここには色々方法は書かれているんだが、プラットフォームに左右されない良い方法というのは無さそうな印象。
皆さんどうされていますか？
0630１３２人目の素数さん
2010/05/09(日) 17:19:53RjWikiの最初に書かれている方法だけあって、Unixでteeを使う方法は素晴らしいですね。
しかし、レクチャーに使いたいと思っていて、プラットフォームに関係のない解決方法、
なるべくならRアプリケーションそのものから解決したいんです。
outputtext <- capture.output(
コマンド１
コマンド２
・・・
コマンドｘ
)
sink("outputtext.txt")
outputtext
sink()
こんな感じで実行ファイルを作ってみたんですが、エラーになります。
良い方法ありますか？
sinkだと、表が崩れるという問題は残りますが。
0631１３２人目の素数さん
2010/05/10(月) 10:27:50Rを適当なインターフェイスを通じて使わずに直接使うからでは？
例えばテキストエディタから使っていれば、Rの実行画面をそのまま
テキストファイルとして保存すればコンソールの出力を丸ごと保存
できるよ。
0632１３２人目の素数さん
2010/05/10(月) 13:48:40Unix系、私はMacなんでターミナルからですが、上記のR|tee R.logで上手くいきました。
テキストエディタというと思いつくのは、Emacs + ESSですが、
それだと講義でEmacsの使い方まで教えないといけないかなと思い躊躇しています。
でも他には方法がなさそうですね。
まぁ、もともと個人的には、結果を見たいところでわざとエラーで止めて、
logファイルを見る使い方はしていなかったので、ちょっと考えてみます。
0633１３２人目の素数さん
2010/05/10(月) 14:37:14テキストエディタは例えばの話で、RのGUIインターフェイスはいろいろとある。
マルチプラットフォームならJGRとか。今、JGRを試してみたけど、保存のアイコンをクリックしたら、
コンソールの出力を丸ごとファイルに保存できたよ（保存するときに、全部、コマンドだけ、出力だけ
を選ぶみたい）。
他のインターフェイスでも同様に出来るのではないだろうか。
Windows版のR GUIも実は保存というメニューがあるかも(手元にないので確認できないが)。
0634１３２人目の素数さん
2010/05/10(月) 15:43:25いろいろありがとうございました。
かなりバカな質問をしていたようです。
MacのGUIインターフェイスのCarbon Emacsでもちゃんと保存のメニューがありました。
コンソールの容量も大きいようで（限界は不明ですが）、長い出力にも十分耐えられそうです。
別のソフトのレクチャーをR用に直すのがRの独学に一番良いかと思いやっています。
またお世話になると思いますが宜しくお願い致します。
0635１３２人目の素数さん
2010/05/11(火) 02:28:09コンソールの保存を手動（マウス操作）ではなく、GUIインターフェイスのR上から、
コマンドで行う方法はあるのでしょうか？
今考えているのは、Rのコマンドのみを記録したcommand.Rファイルにlogの保存の指定も
出来るように書ければ、便利かと思ったのですが。
command.Rファイルを実行すれば、自動的にlog.txtが作成される。
0636１３２人目の素数さん
2010/05/11(火) 10:18:15コンソールの出力を保存する方法は、ターミナルソフトウェアの仕事であって、
Rの仕事ではないから、R自体にそのようなことをさせることは出来ないし、
させるべきでもないと思う。そもそも標準入力、標準出力の概念や「小さなツールを
組み合わせて大きな仕事する」哲学に反する。標準出力に出力された後は、
OS依存のアクセサリーツールが、ファイルに保存するなりなんなり対応すればよい。
また、GUIにこだわらず、単に「Rのコマンドのみを記録したcommand.Rファイルにlogの保存
の指定も出来る」のみの要件であれば、R CMD BATCH command.R log.txtで足りる。
0637１３２人目の素数さん
2010/05/11(火) 15:14:16なるほど。とても参考になります。
ちなみにMac OSXのターミナルからは、
command.Rを用意すれば、
R --save <command.R >commandlog.txt
でcommandlog.txtにログを全てとることができました。
--saveと--no-saveの違いを理解していませんが。
もうちょっと考えてみます。
0638１３２人目の素数さん
2010/05/12(水) 02:44:43combn(5,2)を実行すると
組み合わせが全て表示されます。
この組み合わせを表示するのではなく、
5C2の計算結果の10を得るためにはどうしたらいいですか？
0639１３２人目の素数さん
2010/05/12(水) 03:24:39でどうよ？
0640１３２人目の素数さん
2010/05/12(水) 07:38:29なるほどです。Number of columnなんですね。
勉強になりました。
しかし、今回の場合は
combn(52,13)とかだともう数字が大きすぎて、
エラーになってしまいます。
結局下記リンクのように自作関数を作るのがいいみたいでした。
http://www.geocities.co.jp/HiTeens/5433/R/comb.html
0641１３２人目の素数さん
2010/05/13(木) 12:15:31プログラムを書くために参考となる文献や書籍はありませんかね？
0642１３２人目の素数さん
2010/05/13(木) 17:13:160643１３２人目の素数さん
2010/05/13(木) 17:49:46ML
0644１３２人目の素数さん
2010/05/14(金) 23:11:22要素の大きさの順番ではなく要素自体をソートして返す
orderのコマンドってありますか？
行列・配列・リストでいうところの(x[order(x[,1]),])のようなかんじで
0645636
2010/05/15(土) 05:21:42皆さんはどのようにRを使っていますか？
基本として、
コマンドのみを記したcommand.Rファイルと
その実行結果を表示するテキストファイルを出力するように使いたいのです。
出来れば一度設定してしまえば、MacやPCのようなプラットフォームに関係なくできる環境を
構築したいと思っています。
皆さんの作業環境を教えて頂ければと思います。
上記は講義で、一般的な作業環境構築について話したいので質問しています。
個人的には、矩形処理が得意なEmacs上で、ESSを使っていますが、Rの講義で
Emacsの説明までするのは、やりすぎかなと思っています。
0646１３２人目の素数さん
2010/05/15(土) 10:31:37>>636を書いたのは私だけど、なりすましされたのは初めてだ。
電卓として使うならlittlerを使う。本格的に使うならESS+Sweaveを使う。
講義に関しては、Emacsの講義をすでに履修している状態でもなければ、
Emacsを使うのはやり過ぎだろう。Emacsに限らず、Rに用意された
インターフェイス以外を使うのは混乱の元だと思うが。Windows/MacOSXなら
そのままGUIを、Linuxなら適当なターミナルソフトを使うべきだと思う。
0647１３２人目の素数さん
2010/05/15(土) 12:09:05このスレを埋めた後、R言語で検索できるようになったら
0648635
2010/05/16(日) 10:20:20どうも済みません。なりすましするつもりは無かったのですが、間違えました。申し訳ない。
Sweaveはかなり本格的な感じですね。
http://stat.epfl.ch/webdav/site/stat/shared/Regression/EPFL-Sweave-powerdot.pdf
0649１３２人目の素数さん
2010/05/17(月) 20:56:020650１３２人目の素数さん
2010/05/20(木) 02:35:236行のうち1、3、5行目だけを使いたいのですが、指定の仕方が分かりません。
他の関数も使わないと指定はできないのでしょうか？
0651１３２人目の素数さん
2010/05/20(木) 10:29:07困っている度が低そうで答える気に余りならないが、他の関数が"["だったら、
許容範囲なのかな？
> ff <- tempfile()
> cat(file=ff,paste(1:6,"234567890123",sep=""),sep="\n")
> system(paste("cat",ff))
1234567890123
2234567890123
3234567890123
4234567890123
5234567890123
6234567890123
> read.fortran(ff, c("F5","F5","I3"))[c(1,3,5),]
V1 V2 V3
1 12345 67890 123
3 32345 67890 123
5 52345 67890 123
0652635
2010/06/03(木) 03:40:13例えば
a1, a2, a3, a4 ・・・a10という名前の変数がるときに、
この変数に入った値を一度に表示する方法。
Stataなら
>display a*
でOK
forを使って通し番号のついた変数を自動生成する。
Rではこれでいいですか？
for (i in 1:10){assign(paste("a",i,sep=""),代入する値)}
代入する値に通し番号がついた変数名（例b1,b2..）を使う場合にはどうなりますか？
for (i in 1:10){assign(paste("day",i,sep=""),b`i'}
Stataなら上記のように`i'で添字を操作できます。
解決方法としては、grep()など、正規表現が使える関数を利用することも
考えましたが、上手く扱えません。
例
>apropos("a*") #これでは変数名の羅列で変数の中身がでてきません。
>for (i in 1:10) cat(get(paste("a", i, sep="")),"\n")
この方法ではなんとか目的は達成できますが、もっと簡単な方法はないでしょうか？
宜しくお願い致します。
0653１３２人目の素数さん
2010/06/03(木) 09:34:55リスト使えばいいよ。
0654１３２人目の素数さん
2010/06/03(木) 10:12:21>>653の言うとおりlistを使えばよいと思う。
> a <- lapply(1:10,function(x){runif(2)})
> a[[1]]
[1] 0.02104710 0.54894781
> a[[2]]
[1] 0.3112974 0.5590326
> a[[10]]
[1] 0.1201111 0.7264913
という乱数が2個ずつ入った変数があったとしたら、
> a
[[1]]
[1] 0.02104710 0.54894781
[[2]]
[1] 0.3112974 0.5590326
[[3]]
[1] 0.09237121 0.45438306
[中略]
[[10]]
[1] 0.1201111 0.7264913
で中身が表示される。通し番号 1, 2, 3, ... 10の代わりに、[[1]], [[2]], [[3]],...,[[10]]
と考えればOKじゃないかな。
0655635
2010/06/03(木) 13:56:58>>654
オーなーるほど。
早速の回答ありがとうございました。
変数の扱い方自体の考え方が違うんですね。
しかし、Rは懐が深いですね。
0656635
2010/06/03(木) 17:17:10少し疑問があるのですが、
リストとデータフレームの違いがイマイチ良く理解できていません。
データフレームはいわゆるエクセルなどのスプレッドシート上のデータ群と
理解しています。
リストはデータフレームよりもフレキシブルなんでしょうか？
リストからデータフレームに変換するのには次の方法で良いのでしょうか？
>x <- list(4:6,7:9,10:12)
>x <-data.frame(x)
宜しくお願い致します。
0657１３２人目の素数さん
2010/06/03(木) 17:48:52そのあたりの超入門的なことは、入門書を一読してから尋ねた方がいいと思う。
購入がいやならオンラインでもいろいろとある。
> リストからデータフレームに変換するのには次の方法で良いのでしょうか？
こんなことを聞かれても、欲しいデータフレームがどういうものか見せてくれないと、
答えようがないし、list構造の概念を理解していないからこそ発せられる質問だと思うよ。
> as.data.frame(x)
X4.6 X7.9 X10.12
1 4 7 10
2 5 8 11
3 6 9 12
> as.data.frame(unlist(x))
unlist(x)
1 4
2 5
3 6
4 7
5 8
6 9
7 10
8 11
9 12
どっちもデータフレームだよ。
0658635
2010/06/03(木) 18:58:57オンラインのものは、- 統計解析ソフトR の備忘録PDF -とか、Rjpwikiとかは一通り読みました。
まだちゃんとは読んでいませんが、英語のR本も何冊か持っています。
今のところの体感認識では、
データフレームは、Rデータエディタで表示すると、きれいにスプレッドシートとして表示される。
リストは、そうでない。
また、リストはリストそのものも要素に扱えるということから、
雑多なデータを寄せ集めて表現しているものという認識で間違っていないでしょうか？
つまり、
データフレーム形式では、行・列の２次元でデータに意味付けされている。
リストでは、そういう縛りがない。
いろいろ丁寧に教えていただき勉強になっています。
ありがとうございます。
0659635
2010/06/03(木) 19:18:16リストはデータフレームそのものも要素になりうるんですね。
リストにデータフレームを1から5まで読み込んで、ひとつのデータフレームにまとめたいときに
>newdataframe<-as.data.frame(original_list)
>newdataframe<-merge(original_list[[1:5]])
上記二つの方法が出来ることがわかりました。
二番目の方法では、列データがまるまる重複するデータは自動的に削除されるので
二番目の方法の方が扱いやすいのかと理解しています。
他の方法などあれば、また教えて下さい。
0660635
2010/06/04(金) 12:34:01merge()関数は二つのデータフレームしか結合できないようなので、
>newdataframe<-original_list[[1]]
>for(i in 2:5){newdataframe <- merge(newdataframe, original_list[[i]])}
が正しいようでした。
0661１３２人目の素数さん
2010/06/06(日) 11:29:560662１３２人目の素数さん
2010/06/11(金) 03:23:470663１３２人目の素数さん
2010/06/12(土) 12:16:16下記コマンドでエクセルからデータを取り込みました。
library(foreign)
db <- read.xls("sample.xls", sheet=1)
すると
str(db)
'data.frame': 35 obs. of 38 variables:
$ NO : int 1 2 3 4 5 6 7 8 9 10 ...
$ ID : Factor w/ 35 levels "AA001","AB002",..: 29..
で確かめると要素が文字列（つまりchr）であることが期待された変数IDは、
上記のように因子レベルになっています。
数値を因子に直す方法は、factor(変数名,levels=,labels=c("",""))を使うことは分かったのですが、
因子を文字列に直す方法が分かりません。
宜しくお願い致します。
0664１３２人目の素数さん
2010/06/13(日) 04:29:43as.characterでよくね？
0665１３２人目の素数さん
2010/06/15(火) 00:07:362時間掛かっているのに終わりません。
こんなに時間が掛かるものなのでしょうか？
0666１３２人目の素数さん
2010/06/15(火) 09:58:16そんなのマシンによるだろ。
依存ライブラリのコンパイルもしているだろうから、
依存ライブラリが予め入っていないなら+αの時間がかかるし
0667１３２人目の素数さん
2010/06/16(水) 02:40:13> c1<-rnorm(1024);
> x<-data.frame(c1=c1);
> write.csv(x,file="/tmp/data.csv");
> system("cat /tmp/data.csv | ssh [email protected] \"command line\"");
[email protected]'s password:
>
これで一応転送は出来てるのですが、最初のwrite.csvをやらないようにしたい
のですが、やり方がわかりません。出力先をパイプにすればいいのでしょうか？
どうかよろしくお願いします。
0668１３２人目の素数さん
2010/06/16(水) 02:55:12とかやると出来そうな雰囲気もあるのですが、パイプに出力をうまくかけません。
0669１３２人目の素数さん
2010/06/17(木) 11:42:06rownameが付いたデータがあります。
as.vectorでもas.listでもfalseになるので型がよく判りませんが、
数値の入った一次元配列？です。
データが1以下の行は使わないので、0とかNAにしましたが、
そのあとの計算で邪魔になってします。
データが1以下の行を削除するには、どうしたらいいのでしょうか？
調べても判らなかったので、ここを見ろ、でも有り難いです。
0670１３２人目の素数さん
2010/06/17(木) 13:04:09> a <- (1:4)/2
> a
[1] 0.5 1.0 1.5 2.0
> a[a>1]
[1] 1.5 2.0
こういうこと？
0672１３２人目の素数さん
2010/06/23(水) 14:34:4010行2列のマトリックスAと20行2列のマトリックスBを使って
回帰分析を行いたいんですが、エラーが出ます。
具体的には、A,Bの一列目を従属変数(濃度)、二列目を独立変数(スペクトル)にしてます。
関数lmを使ってAで作った検量線に、関数predictを使ってBの二列目を代入したい
(スペクトルの値から濃度を予測したい)と考えているのですが、
「データ数が違う」という旨の警告が出てできません。
分かりにくいかもしれませんが、宜しくお願いします。
0673１３２人目の素数さん
2010/06/23(水) 15:33:28Rjpwikiの方が親切に回答してもらえると思うぞ。
> A <- data.frame(matrix(1:20,nrow=10))
> head(A)
X1 X2
1 1 11
2 2 12
3 3 13
4 4 14
5 5 15
6 6 16
> B <- data.frame(matrix(runif(40),nrow=20))
> head(B)
X1 X2
1 0.6353921 0.2679418
2 0.2917369 0.3124313
3 0.2857031 0.2795522
4 0.1951349 0.4510532
5 0.4426567 0.4894139
6 0.5981174 0.3142484
> predict(lm(X1~X2,A),B)
1 2 3 4 5 6 7 8
-9.732058 -9.687569 -9.720448 -9.548947 -9.510586 -9.685752 -9.718648 -9.481105
9 10 11 12 13 14 15 16
-9.477939 -9.628372 -9.697301 -9.711742 -9.717192 -9.791366 -9.910725 -9.650475
17 18 19 20
-9.903119 -9.638157 -9.547609 -9.051490
ってこと？具体的なコード（再現性のあるサンプルコード）を示してくれないと、
なんで出来ないのかこちらには分からない。
0674１３２人目の素数さん
2010/07/04(日) 01:35:420675名無しさん＠そうだ選挙に行こう
2010/07/11(日) 17:05:40anova(fit)
で出てくるF値は、t値の二乗になっていないですよね。
このF値は統計学的にどういう意味のあるF値なんですか。
0676名無しさん＠そうだ選挙に行こう
2010/07/11(日) 20:31:09Rtoolsを利用して通常のtar→makeまではできたのですが、UNIXでの./configure --enable-R-shlib のようなコマンドがわからず困っています…
0677名無しさん＠そうだ選挙に行こう
2010/07/11(日) 20:39:550678名無しさん＠そうだ選挙に行こう
2010/07/11(日) 21:41:450679１３２人目の素数さん
2010/07/13(火) 10:02:28少なくともあなたがコンピュータに無縁なのはよく分かった
0680１３２人目の素数さん
2010/07/14(水) 21:58:23factanalを使って1因子の因子分析をしているのですが、初期解を求めようとすると
「以下にエラー solve.default(cv) :
システムは数値的に特異です：条件数の逆数 = 2.13584e-17 」
というエラーが出てしまいます。
prcompを使うと正しく計算されます。
しかし「前回課題で使用したデータ」を「factanalを使って1因子分析」するのが条件なので、
関数やデータをいじるわけにもいきません。
どうしたら解決できるのでしょうか。
0681１３２人目の素数さん
2010/07/14(水) 22:06:130682１３２人目の素数さん
2010/07/14(水) 22:18:12Rに比べると怪しげに見える。
よく知らんけど、Rはgccと同じくらいのレベルで標準を確立したりしていないか？
0683１３２人目の素数さん
2010/07/14(水) 23:13:21それってなぜエラーが出るか考えさせる問題ってことはないよね？
0684１３２人目の素数さん
2010/07/15(木) 11:33:10まだプロジェクトが続いていたのか。
昔試したときには、GUI化される前のSPSS Baseのコピーって
かんじだったけど。公式ページを見たら、psppもGUI化されたみたいだな。
0685１３２人目の素数さん
2010/07/17(土) 20:20:27例えば、
x<-1:10
which(x<5&&x>2)
ってしたら
integer(0)
ってでてくる
0686１３２人目の素数さん
2010/07/17(土) 20:25:440687１３２人目の素数さん
2010/07/20(火) 10:16:24&と&&は意味が違うので注意って入門書に書いていなかった？
0688１３２人目の素数さん
2010/08/01(日) 12:35:35>GUI化
結構いい線いってると思うけど、本物と比べると出来ないことが多すぎて
代わりにはなりませんな。
0689１３２人目の素数さん
2010/08/01(日) 13:09:12次にRscriptを呼び出して結果を出力して、
後始末して、ということをしてるけど、
batファイルというのがどうも不格好なので、
rファイルだけで全部やらせられない?
■ このスレッドは過去ログ倉庫に格納されています