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

1 名前：１３２人目の素数さん mailto:sage [2011/01/09(日) 17:31:49 ] R は統計計算とグラフィックスのための言語・環境です。 統計計算で重宝するデータ型や、複数要素を処理する演算や関数、 解析結果を表示するグラフィックなど、多彩な機能を提供します。 ●関連サイト The R Project www.r-project.org/ RjpWiki www.okada.jp.org/RWiki/ リンク集 www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8 53 名前：１３２人目の素数さん mailto:sage [2011/04/28(木) 13:50:34.33 ] >>52 nonrandomパッケージのVignettesを読んでからもう一度質問してください。 ttp://cran.r-project.org/web/packages/nonrandom/vignettes/nonrandom.pdf 54 名前：１３２人目の素数さん mailto:sage [2011/04/28(木) 17:38:07.68 ] >>53 ありがとうございます。 55 名前：１３２人目の素数さん mailto:sage [2011/04/30(土) 09:11:31.24 ] nattyきたー ttp://cran.r-project.org/bin/linux/ubuntu/ 56 名前：１３２人目の素数さん mailto:sage [2011/04/30(土) 13:11:01.00 ] >>55 Releaseのファイルサイズが0でまだ使えないね $ w3m -dump cran.r-project.org/bin/linux/ubuntu/natty/ |grep -oE '[a-zA-Z0-9_\.-]*amd64\.deb' littler_0.1.3-1natty0_amd64.deb python-rpy_1.0.3-17natty0_amd64.deb python-rpy_1.0.3-18natty0_amd64.deb r-base-core-dbg_2.13.0-1natty0_amd64.deb r-base-core_2.13.0-1natty0_amd64.deb r-cran-class_7.3-3-1natty0_amd64.deb r-cran-cluster_1.13.3-1natty0_amd64.deb r-cran-foreign_0.8.44-1natty0_amd64.deb r-cran-kernsmooth_2.23-4-1natty0_amd64.deb r-cran-lattice_0.19-26-1natty0_amd64.deb r-cran-mass_7.3-11-1natty0_amd64.deb r-cran-matrix_0.999375-50-1natty0_amd64.deb r-cran-mgcv_1.7-6-1natty0_amd64.deb r-cran-nlme_3.1.100-1natty0_amd64.deb r-cran-nnet_7.3-1-1natty0_amd64.deb r-cran-rodbc_1.3-2-1natty0_amd64.deb r-cran-rpart_3.1.50-1natty0_amd64.deb r-cran-spatial_7.3-2-1natty0_amd64.deb r-cran-survival_2.36-8-1natty0_amd64.deb r-mathlib_2.13.0-1natty0_amd64.deb 今あるパッケージはこれだけ。 57 名前：１３２人目の素数さん mailto:sage [2011/05/03(火) 20:04:04.55 ] ニューラルネットのパッケージnnetについて質問させてください。 すでにデータAを学習したネットワークにさらに別のデータBを学習させるにはどうすればいいのでしょうか？先にAとBを結合して読ませるのは時間がかかるために困っています。 他のパッケージなら出来るという情報もお待ちしています。 58 名前：１３２人目の素数さん mailto:sage [2011/05/05(木) 21:13:39.01 ] MySQLからダウンロードしたデータをRで計算して MySQLに返しつつ、gnuplotで描画しようと思っています RubyとかからMySQL、R、gnuplotを操作すべきか、 RからMySQLとgnuplotを操作すべきかで迷ってます みなさんならどんな風にしますか？ 59 名前：１３２人目の素数さん mailto:sage [2011/05/06(金) 13:30:23.98 ] >>58 このスレ的には、Rの中からMySQLにアクセスして、 Rで描画する方向が第一選択だよな。 Rのグラフィック能力とgnuplotってほぼ同等と思っていたけど、 Rを使っているのにgnuplotを加えて使うのは理由は？ >>57 反応がないから、Rjpwikiの方に書かれてはいかが。 私はnnetパッケージを使わないので分からない。 60 名前：57 mailto:sage [2011/05/06(金) 22:06:00.16 ] >>59 そうしてみます。ありがとうございました。 61 名前：１３２人目の素数さん mailto:sage [2011/05/06(金) 23:00:08.59 ] gnuplotの方が慣れているというだけです。 Rで描けるようになればR+MySQLなのでシンプルでいいですね。 ちょっとやってみます。ありがとうございました。 62 名前：１３２人目の素数さん mailto:sage [2011/05/09(月) 23:39:55.72 ] pi<-seq(0,1,0.0001) y_n<-60 n<-100 posterior<-((pi^y_n)*(1-pi)^(n-y_n))/beta(y_n+1,n-y_n+1) の95％信頼区間を求めたいんですがその関数がよくわからないです誰かわかる人いらっしゃいますか？ 63 名前：１３２人目の素数さん mailto:sage [2011/05/10(火) 09:08:37.12 ] >>62 Bayesian credible intervalは信頼区間じゃないよ 64 名前：１３２人目の素数さん mailto:sage [2011/05/10(火) 15:28:06.31 ] >>63　だね．　bayequentistかな 65 名前：１３２人目の素数さん mailto:sage [2011/05/14(土) 21:58:27.95 ] for (i in 1:5) {sprintf("%s", ticker[i,1])} だと上手く出力されないのはどうしてでしょう？ 素人の質問ですみません 66 名前：１３２人目の素数さん mailto:sage [2011/05/14(土) 22:42:04.29 ] >>65 > sprintf("%s", ticker[,1]) > sprintf("%s", ticker[1:5,1]) > as.character(ticker[,1]) > as.character(ticker[1:5,1]) などと1行で書くのがR流じゃないかと。 ループ使った場合でも > for (i in 1:5) {a[i]<-sprintf("%s", ticker[i,1])} > a とやれば配列aが表示される。 67 名前：１３２人目の素数さん mailto:sage [2011/05/15(日) 17:30:01.07 ] >>65 for文の中ではprint()が必要です。 for (i in 1:5) {print(sprintf("%s", ticker[i,1]))} 普通は>>66の言うとおり、 as.character(ticker[1:5,1])だけどね。 68 名前：１３２人目の素数さん [2011/05/25(水) 00:02:27.53 ] データ配列の最初から20個ずつをまとめ、 それぞれの20個のデータの平均値と標準偏差を求めるプログラムstatsample(m,data)を作成する。 プログラムソースが分かる人は教えてください。 69 名前：１３２人目の素数さん mailto:sage [2011/05/25(水) 09:29:38.65 ] >>68 とりあえず、たたき台を示すので、データ数が20で割り切れないをどうするのかを 考えて、その処理を入れて、その他のエラー処理を入れたら完成だと思う。 statsample <- function(x){ n <- length(x) m <- floor(n/20) for (i in 0:(m-1)){ print(paste("平均値:",mean(x[(1+i*20):((i+1)*20)]),"標準偏差",sd(x[(1+i*20):((i+1)*20)]))) } } > statsample(runif(100)) [1] "平均値: 0.409618140733801 標準偏差 0.284202820077949" [1] "平均値: 0.449437574110925 標準偏差 0.286492195081962" [1] "平均値: 0.555866089393385 標準偏差 0.26665116632303" [1] "平均値: 0.28123870670097 標準偏差 0.226341301484355" [1] "平均値: 0.503801094950177 標準偏差 0.282596081360767" 70 名前：１３２人目の素数さん mailto:sage [2011/05/25(水) 19:42:32.88 ] >>69 datamatrix<-matrix(data,20,length(data)/20) apply(datamatrix,2,mean) apply(datamatrix,2,sd) とやるのが普通じゃないか。 データ数が20で割り切れない場合のwarningも表示されるし。 71 名前：１３２人目の素数さん [2011/05/25(水) 19:45:58.57 ] 農家の果物の良品、不良品の検査や工場での製品の良品、不良品のチェックなどをシュミレーションにより実験する場合を考える。 このシュミレーションを行うために、不良品率ｘ％でn個のサンプルの中に不良品を発生させるプログラムsampling(x,y)を以下の手順で作成する。 (1)不良品の出現確率ｘ/100とする。発生させた一様乱数がｘ/100より小さければ、不良品(1)、そうでなければ、良品(0)とする。 (2)不良品出現確率ｘ/100とした試行をn回繰り返し、不良品が出現した場合の数を求める。 　その数がsampling(ｘ,n)の出力とする。 以下の手順をプログラムし、関数sampling(ｘ,n)を作る。 誰か教えてください。 72 名前：１３２人目の素数さん mailto:sage [2011/05/25(水) 21:51:26.02 ] #include "/dev/console" 73 名前：１３２人目の素数さん [2011/05/26(木) 00:30:46.84 ] ＃ 74 名前：１３２人目の素数さん [2011/05/26(木) 01:25:55.19 ] >>69 75 名前：１３２人目の素数さん mailto:sage [2011/05/26(木) 09:18:07.75 ] >>71 宿題？ ヒントだけ > runif(10) < 5/100 [1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE 76 名前：１３２人目の素数さん mailto:sage [2011/05/26(木) 22:53:09.94 ] rbinom(1,n,x/100) 77 名前：１３２人目の素数さん [2011/05/31(火) 02:28:47.96 ] sample関数について、中でどうやってドローしているんだろうと、 見ようとしたら、 .Internal( sample.... となってて、中身が見えません。 他の言語で作成しているのだろうと思いましたけど、 sample関数のソースって、どこにあるんでしょうか？ 78 名前：７７ [2011/05/31(火) 02:32:14.67 ] 自己解決しました。 CRANからソースを取得しました。 79 名前：あんでぃは非存在 ◆AdkZFxa49I [2011/06/06(月) 14:26:34.02 ] あんでぃ 80 名前：１３２人目の素数さん [2011/06/07(火) 04:41:09.16 ] 乱数のところ、ソースみてもわからん。 81 名前：１３２人目の素数さん mailto:sage [2011/06/07(火) 22:07:44.63 ] 乱数だからな 82 名前：１３２人目の素数さん mailto:sage [2011/06/07(火) 22:18:53.87 ] メルセンヌツイスターとかいうやつか 83 名前：１３２人目の素数さん mailto:sage [2011/06/17(金) 13:02:41.63 ] すいません、初心者です 株価の日次データが数年分あるんですが グラフのX軸の目盛りを年次にしたい場合はどうすればいいでしょう？ 84 名前：１３２人目の素数さん mailto:sage [2011/06/17(金) 13:23:39.86 ] >>83 a <- ts(1:2000,start=c(2011,6),frequency=365) plot(a) 85 名前：１３２人目の素数さん mailto:sage [2011/06/17(金) 13:35:41.47 ] >>84 早速のレスありがとうございます 株など営業日がある場合どうすればよろしいでしょうか？ 祝日や年末年始などありますし 86 名前：１３２人目の素数さん mailto:sage [2011/06/17(金) 13:52:29.96 ] >>85 xtsパッケージのドキュメントを眺めてから再質問してください。 ttp://cran.r-project.org/web/packages/xts/ 日本語が必要なら、Rパッケージガイドブック pp.130-137 87 名前：１３２人目の素数さん mailto:sage [2011/06/17(金) 13:59:31.54 ] >>86 ありがとうございます。 パッケージガイドブックを買ってこようと思います。 これとR-tipsで頑張ります 88 名前：１３２人目の素数さん mailto:sage [2011/06/23(木) 16:49:30.62 ] RでVoigt関数(Gauss関数とLorentz関数のconvolution)でフィッティングするには、畳み込み積分を直接書くしかないのでしょうか。 89 名前：あんでぃはストーカー ◆AdkZFxa49I [2011/06/23(木) 21:21:17.46 ] あんでぃ 90 名前：１３２人目の素数さん mailto:sage [2011/06/24(金) 09:20:34.16 ] format="%Y/%m" と書くと 2011/06 と出力されるけど、幅を節約して 11/06 とするにはどうすれば? 91 名前：１３２人目の素数さん mailto:sage [2011/06/24(金) 09:22:44.15 ] あ、ぐぐったら一瞬で判ってしまった %y だ 92 名前：１３２人目の素数さん mailto:sage [2011/06/24(金) 10:43:36.69 ] >>90,91 そのあたりは共通だから覚えておいて損はないよ。 > format(Sys.Date(),"%Y") [1] "2011" > format(Sys.Date(),"%y") [1] "11" > format(Sys.Date(),"%EY") [1] "平成23年" > format(Sys.Date(),"%Ey") [1] "23" 93 名前：１３２人目の素数さん mailto:sage [2011/06/24(金) 12:32:27.69 ] おおー なんでEが平成なんだろう 94 名前：１３２人目の素数さん mailto:sage [2011/06/24(金) 13:01:48.18 ] >>93 マニュアルには、 %E Modifier: use alternative format と書いてあるよ。つまり日本語ローケルだと、日本語化した形式にするという意味。 >>92を書き込んだ後、MacOSXとWindowsで試してみたが、 どちらも和暦にならなかった。GNUだけかも。ヘルプには環境依存って書いてあるしね。 95 名前：あんでぃはストーカー ◆AdkZFxa49I [2011/06/24(金) 16:34:26.88 ] 分からない。 あんでぃ 96 名前：R入門 [2011/06/26(日) 00:02:26.05 ] これまで、Java、Ruby、Scala など使っているものですが、 Rについても勉強したいと思っています。 Rのサイトをみると、結構、「使い方」的なサイトが多く、 汎用言語の視点から書かれた、サイト、本が少ない(？)ような気がします。 たとえば、ほかの言語と比較しながら書かれたような 本、サイトなどあれば、教えていただきたいのですが。 ありますでしょうか。 97 名前：１３２人目の素数さん mailto:sage [2011/06/26(日) 00:19:39.46 ] ぐぐりにくい、というのがRの最大の弱点 98 名前：１３２人目の素数さん mailto:sage [2011/06/26(日) 01:03:33.08 ] >>96 Rの基礎とプログラミング技法 www.amazon.co.jp/dp/4431712186 d.hatena.ne.jp/sesejun/20090331/p2 Rと他の言語を比較した本は知りませんが、 汎用言語としてのR言語を勉強するならこの本がお勧めです。 99 名前：１３２人目の素数さん mailto:sage [2011/06/27(月) 09:37:34.48 ] >>96 そういえば、Rユーザ会で、韓国から来た演者が汎用言語の視点から、Rの問題点などを 発表していたな。興味がなかったので内容は覚えていない。 汎用言語の視点で書かれた書籍はないと思う。 しいてあげるなら、>>98が薦めるリゲスさんの本しかない。 100 名前：85 mailto:sage [2011/07/04(月) 10:41:00.37 ] Rパッケージガイドブックを購入 xtsをインストールしたんですが データフレームをxtsにコンバートできません。 > x.xts<-as.xts(x) 以下にエラー as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format Zooオブジェクトに経由させることもできません > x.xts<-as.xts(read.zoo(x)) 以下にエラー xts(coredata(x), order.by = order.by, frequency = frequency, : order.by requires an appropriate time-based object データフレームの中身は 日付 株価 20040101　　　 1500 20040102 1450 といった具合に２列のデータなのですが？ 101 名前：１３２人目の素数さん mailto:sage [2011/07/04(月) 11:05:02.37 ] >>100 ここに書くのではなくて、 それを書いた人に凸してみたら？ > library(xts) 要求されたパッケージ zoo をロード中です > dat <- data.frame(日付=c(20040101,20040102),株価=c(1500,1450)) > dat 日付 株価 1 20040101 1500 2 20040102 1450 > as.xts(dat) 以下にエラー as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format > rownames(dat) <- c("2004-01-01","2004-01-02") > as.xts(dat) 日付 株価 2004-01-01 20040101 1500 2004-01-02 20040102 1450 まぁ、参考になれば。 102 名前：１３２人目の素数さん [2011/07/14(木) 14:34:30.11 ] Rで複素誤差関数 w(x) = exp(-x^2) erfc(-ix) を使いたいのだけれども、これはどうやって定義すればよいでしょうか。 相補誤差関数 erfc を qnorm 関数から定義すると、 qnorm が実数に対してしか定義されていないためエラーになってしまいます。 103 名前：１３２人目の素数さん [2011/07/15(金) 02:10:37.99 ] >>96 ＞　汎用言語の視点から書かれた、サイト、本が少ない(？)ような気がします。 Rの本を書いているのは、大学の先生が多いからね。 そして、大学の先生ってのはコンピュータのことや、 プログラミングのことをわかってないから。 Rの使い方だけの説明になってしまうのは しかたがないと思いますよ。 104 名前：１３２人目の素数さん mailto:sage [2011/07/15(金) 11:23:45.87 ] ×　大学の先生 ○　情報科学が専門ではなく疫学やら統計屋が仕事の教員 105 名前：１３２人目の素数さん [2011/07/16(土) 00:45:00.68 ] >>104 ○　大学の先生 Ｘ　情報科学が専門ではなく疫学やら統計屋が仕事の教員 特に情報科学の先生や院生などは、知ったかの馬鹿ばっかりです。 国立情報学研究所を見ればよくわかるでしょ？ 税金を使ってＰＣで遊んでるだけ。 全員がコネ採用。 106 名前：１３２人目の素数さん [2011/07/16(土) 01:04:55.28 ] 情報研究所て何やってるんだ？ 唯のお飾りか？ 107 名前：１３２人目の素数さん [2011/07/16(土) 16:16:24.06 ] 国立情報学研究所は、仕事したくない院生や教授が 税金を無駄遣いしながら、馬鹿同士楽しく PCで遊んでいる組織ですね。 どいつもこいつも馬鹿ばっかり。 コネがないと絶対に採用されません。 公募する前から採用者は決まってます。 108 名前：１３２人目の素数さん mailto:sage [2011/07/16(土) 16:49:55.56 ] なんか嫌なことでもあった？ 109 名前：１３２人目の素数さん [2011/07/16(土) 20:02:28.47 ] 国立情報学研究所 統計数理研究所 この二つは廃止で結構です。税金の無駄です。 110 名前：１３２人目の素数さん [2011/07/16(土) 20:09:53.71 ] 疫学や感染症のシミュレーションをやっている連中なんざw バカばっかりやなw　あれで数学やっているつもりのもいるからw 頓珍漢なことばかりを言って来るw 111 名前：猫は嘘吐き ◆MuKUnGPXAY mailto:age [2011/07/16(土) 20:20:25.68 ] 猫 112 名前：１３２人目の素数さん [2011/07/16(土) 20:58:33.45 ] というか、研究者というものは情熱と実力が必要だと思う。 しかし、単なる就職先のひとつとなっており、 学歴とコネのみで判断されるから。 実力で研究者になれるところって、少ないと思う。 113 名前：１３２人目の素数さん mailto:sage [2011/07/16(土) 21:06:58.59 ] そこの研究員全員の給料の上限を年収３００万ぐらいにすれば他から文句はでないと思うよ。 114 名前：１３２人目の素数さん mailto:sage [2011/07/16(土) 21:19:33.03 ] 職を失ったのか？ 115 名前：猫は無職 ◆MuKUnGPXAY mailto:age [2011/07/17(日) 03:35:47.40 ] 職を失ったんは結構前やがな。 猫 116 名前：１３２人目の素数さん [2011/07/17(日) 13:55:02.20 ] >>113 最近の調査では勤労者の平均年収が ３７０万円ぐらいらしいんですよ。 ただし最頻値は２５０万円ぐらいですから、 一部の高給取りが平均を上げていると思います。 また、公務員の平均は５００万円ぐらいらしいですが、 公務員だけの平均年収は公表されていません。 研究と称してPCで遊んでいるだけですから、 年収１５０万円でも十分だと思います。 117 名前：猫の夢と痴漢の埋葬 ◆MuKUnGPXAY mailto:age [2011/07/17(日) 14:08:35.01 ] >>116 『研究と称してＰＣで遊んでる人』だけをそういう待遇にしたら良いと 思います。但しソレをどうやって判定スルのかは現実には難しく、なの で『虚偽院生と優秀院生の分離』と同様の困難が伴い兼ねませんね。 では質問ですが、虚偽院生に対してはどういう扱いが宜しいのでしょう かね？　やはり追放ですかね？ 猫 118 名前：１３２人目の素数さん [2011/07/17(日) 14:12:09.83 ] 国立情報研究所も統計数理研究所も、仕事らしい仕事をしているのはほんの一部なのではないかな？ 統計数理研究所では女の子の飲み物調査という論文もあったぞw 119 名前：１３２人目の素数さん [2011/07/17(日) 14:24:35.05 ] >>117 虚偽院生とは、論文を書かない院生のことですか？ 全く論文を書く気のない院生は追放でしかたないと思います。 少なくとも論文を何本か書いていて、研究費も助成金も１円も 受け取っておらず、アルバイトをして授業料を納めている奴は、 他の人の研究の邪魔をしないかぎりは、追放でなくて放置でいいと思います。 つまり、まあおれのことですけどね。 追放すべきは高い給料、研究費や助成金をたんまりもらっているのに 結果を出さない馬鹿な奴です。 これは院生だけではなく、助手も助教、準教授、教授も含まれます。 それらは能力と業績だけで判定されるべきだと思います。 その業績というのをどう判定するのかが難しいわけですけどね。 120 名前：猫の夢と痴漢の埋葬 ◆MuKUnGPXAY mailto:age [2011/07/17(日) 14:29:59.14 ] >>119 全くその通りだと思います。私が言う虚偽院生とは： ★★★『基礎体力も無ければ動機も無く、唯大学院に入学してしまったゾンビ』★★★ を指し示す言葉です。ですから恐らくは貴方の事ではないでしょう。放置 されていても自分から何かが出来る人は決してゾンビではアリマセンから。 猫 121 名前：１３２人目の素数さん [2011/07/17(日) 14:41:00.64 ] ＞全くその通りだと思います。 つまり>>119の発言にまったく同意しているという解釈でよいですか？ 122 名前：猫の夢と痴漢の埋葬 ◆MuKUnGPXAY mailto:age [2011/07/17(日) 14:42:52.30 ] >>121 はい。『私は>>119氏の発言に全く同意』です。 猫 123 名前：１３２人目の素数さん [2011/07/17(日) 14:51:47.93 ] やっぱCOLDPLAYいいよねぇ 眼醒めたらゾンビみたいだったが、観てたら生気がもどってきた この辺のクオリティの高さは邦楽にはない、流石だとおもう。 124 名前：１３２人目の素数さん mailto:sage [2011/07/17(日) 23:18:34.72 ] コネっていうと響きは悪いけど、要するにコミュニティーを広げる努力をしてきたってことじゃない？ 周囲と意見交換が出来るということは、非常に大きな能力であり財産だと思うよ そして研究室への貢献度という点で言えば 論文は書かないけどPCで遊べる程度の実力はあって、周囲と意見交換が出来る奴　＞　一人黙々と論文を書く奴 だと自分は思う 125 名前：１３２人目の素数さん mailto:sage [2011/07/17(日) 23:41:22.12 ] なら年収も１３０万でいいな 月１０万でボーナス（一ヶ月)も出て国家的身分は安泰 やりがいなんかまったくないと思うかもしれないが、言うことなしの職場じゃないか 126 名前：１３２人目の素数さん mailto:sage [2011/07/17(日) 23:45:41.25 ] >>116 それこそ、そこのなんとか研究所に問い合わせをして詳細資料（統計データ)請求すればいいんじゃないの？ ただ気になったんだけど、rlang使ってるなら「加重平均」ってなんのことかはしってるよね。 127 名前：１３２人目の素数さん [2011/07/18(月) 02:46:14.96 ] ↑馬鹿だ（笑 128 名前：猫は爺惨状の蚊軍 ◆MuKUnGPXAY mailto:age [2011/07/18(月) 03:33:15.34 ] >>124 その考え方は決して認められない。大学は馬鹿がママゴトして遊ぶ場所 では断じてない。当然に： ★★★『黙々と論文を書く奴＞論文は書かない奴』★★★ であり、PCで遊ぶかどうかや、周囲と意見交換が出来るかどうかなんて、 そんな事はどうでもヨロシ。 論文を書くかどうかだけが問題である。勘違いも甚だしい。ケシカラン。 猫 129 名前：１３２人目の素数さん [2011/07/18(月) 03:34:40.12 ] >>128 まずお前が実績を出してから批判しような 130 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 04:44:10.88 ] >>126 「左右非対称分布の平均値は分布を的確に示せない」 のは統計の基本常識です。 集計方法と指標の定義: さまざまな平均値/中央値/最頻値 web-tan.forum.impressrd.jp/e/2008/07/11/3428 131 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 05:55:50.19 ] R使うくらいならmatlab使った方がいいと思う 132 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 06:00:29.17 ] それはない 133 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 07:56:35.43 ] 猫さんが正論をいっとる 134 名前：１３２人目の素数さん [2011/07/18(月) 07:59:15.18 ] 統計なんてやっている奴はバカw 135 名前：１３２人目の素数さん [2011/07/18(月) 10:43:44.13 ] 統計は低脳のやる学問もどきである どこの大学でも統計の先生の論文は下らんものばかり 136 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 10:54:16.08 ] 量産可能で下らなくない論文とは例えばどういうのですか？ 137 名前：１３２人目の素数さん [2011/07/18(月) 11:20:01.81 ] 統計ですと、同じ様な調査結果を論文として何度も書く事が出来ます。 138 名前：１３２人目の素数さん [2011/07/18(月) 11:22:33.75 ] 女子大生の飲み物調査 女子高校生の飲み物調査 女子中学生の飲み物調査 女子小学生の飲み物調査 女子幼稚園生の飲み物調査 ２０代OLの飲み物調査 ３０代のOLの飲み物調査 以下、女と男にかえても論文が書けます　統計数理研究所の教授がこうした論文を書いています 139 名前：１３２人目の素数さん [2011/07/18(月) 11:24:51.54 ] かけられた電話に出るまでの時間の統計的調査１　一部上場企業の場合 同上　２　２部上場企業の場合 同上　３　ヘラクレス上場企業の場合 同上　４　町役場の場合 以下　延々と書けます。 140 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 11:30:15.87 ] 今渦中の経済産業省官僚のやることなんかはもっと姑息ですけどねw 141 名前：１３２人目の素数さん [2011/07/18(月) 15:21:16.30 ] グランドカノニカルは糞 142 名前：１３２人目の素数さん [2011/07/18(月) 18:25:28.32 ] 等高線の書き方について質問です。 等高線を書きたいのですが、情報量が少ないため、contour関数では書けません。現在与えられているのは、20個の点に対しての座標(x,y,z)です。2次元上の布置まではできているのですが、等高線を用いてz座標を表現することができません。 どの関数を扱えばいいかなど、わかる方がいらっしゃいましたら教えてください。 143 名前：１３２人目の素数さん mailto:sage [2011/07/18(月) 18:59:47.11 ] 手で書いたら？ 天気図も昔はそうしたよ。 144 名前：１３２人目の素数さん [2011/07/18(月) 22:51:46.35 ] >>142 つか、授業聞いてなかったの？ 全部説明したはずだけどね。 geoRで作成しなさいって言ったでしょ？ 単位あげないよ！ 145 名前：１３２人目の素数さん mailto:sage [2011/07/31(日) 19:04:21.64 ] データフレームの列名に日本語で文字化けして困ってます。 Rコマンダーでcsv経由。 RGUIはMSGhosicに設定。 コマンドでフォント指定以外に 対策あれば教えてください。 146 名前：１３２人目の素数さん mailto:sage [2011/08/01(月) 00:30:07.52 ] >>145 Rというよりそのcsvファイルがくせもの 解決策は2つ。 1. csvの文字コードを文字化けしないものに変更してからRに読み込む 2. Rの中でiconv()で文字コードを変換する 好きな方を選べ 147 名前：１３２人目の素数さん mailto:sage [2011/08/01(月) 05:59:12.09 ] >>146 了解です。 クリップボードを経由してみたら、ますます 酷くなりました。 特に列の変数名がボロボロです。 これはなにゆえでしょうか？ 148 名前：１３２人目の素数さん mailto:sage [2011/08/01(月) 11:35:11.11 ] >>147みたいな学生を指導する教員は大変だな。 全く別のアプローチを提示してみます。 > dat <- read.csv(file.choose()) > dat bake bakeda 1 32 re 2 242 ga ここで列名が文字化けしていたとします。 本当に化けている状態を示せないので仮に"bake"とか"bakeda"とかにしています。 > names(dat) <- c("日本語","大丈夫") > dat 日本語 大丈夫 1 32 re 2 242 ga というようにnames()で、列名を自由に再設定できます。 列名が、文字化けしていて、なおかつ文字コードを変更する能力がないなら、 上記のような方法もあります。 149 名前：１３２人目の素数さん mailto:sage [2011/08/08(月) 14:36:57.39 ] [S] 東大, 弁護士, Re, 150 名前：１３２人目の素数さん mailto:sage [2011/08/08(月) 14:38:15.29 ] [A] TS10,SBR,VFK10,TKK,VF1,LCCR,SINX, VF1M4,VF1L2,VF1H2,EMPC,MPE, 4231,4213,3331,3313,145,53A3,6236, EMPCB,EMJ,LP,CJ,F4,LC,DNA,RNA, SINT,JEL,23458,2348,DBT,GMO, AB,APLWJKSJ,PES,WE,CA,RR,ASL, EPH,ITU,261036,CBS6,1358,G1,AS3,M5, 151 名前：１３２人目の素数さん [2011/08/11(木) 22:25:16.52 ] 質問です。 20個の変数があり、各変数は0 から 0.2 まで 0.05 ずつの値をとることが可能 また、この変数の総和は 1 となる。 この時の組み合わせをすべて求めたいのですが、 どうすれば良いでしょうか? よろしくお願いいたします 152 名前：１３２人目の素数さん mailto:sage [2011/08/12(金) 01:53:48.61 ] >>151 > seq(from=0,to=0.2,by=0.05) [1] 0.00 0.05 0.10 0.15 0.20 > factorial(5+20-1)/factorial(20)/factorial(5-1) [1] 10626 この程度の組み合わせなら全て計算して合計が1のものを選んだ方が早いかも。 153 名前：152 mailto:sage [2011/08/12(金) 02:20:56.85 ] とりあえず、20変数じゃなくて6変数の場合を提示するのであとは自分で考えて > n <- seq(from=0,to=0.2,by=0.05) > m <- expand.grid(n,n,n,n,n,n) > m.s <- apply(m,1,sum) > a <- m[m.s==1,] > head(a) Var1 Var2 Var3 Var4 Var5 Var6 3125 0.20 0.20 0.20 0.20 0.20 0.00 5625 0.20 0.20 0.20 0.20 0.15 0.05 6125 0.20 0.20 0.20 0.15 0.20 0.05 6225 0.20 0.20 0.15 0.20 0.20 0.05 6245 0.20 0.15 0.20 0.20 0.20 0.05 6249 0.15 0.20 0.20 0.20 0.20 0.05

---

---

[ 続きを読む ] / [ 携帯版 ]

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