【R言語】統計解析フリーソフトＲ 第４章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/01/09(日) 17:31:49統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
●関連サイト
The R Project
http://www.r-project.org/
RjpWiki
http://www.okada.jp.org/RWiki/
リンク集
http://www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8
0047１３２人目の素数さん
2011/04/15(金) 20:26:39.00遅延評価Ｃで書くのすごい大変なのに
どうやってるんだろう
0048１３２人目の素数さん
2011/04/19(火) 11:28:58.25tar tf R-2.13.0.tar.gz | grep -oE '\.[a-z0-9]*$' | sort |uniq -c | sort -gr |head
599 .c
240 .h
194 .po
169 .mo
94 .in
87 .afm
52 .win
50 .save
45 .f
27 .gmo
RjpwikiのQ&A(初心者コース)に書き込んだmichiさんはこっちに来ないのかな。
0049１３２人目の素数さん
2011/04/27(水) 15:52:48.74Nattyまだこないね。
0050１３２人目の素数さん
2011/04/27(水) 21:40:58.88このソフトpropensity scoreを用いて解析をするときの参考になるサイトや書籍がありましたら教えてください。
0051１３２人目の素数さん
2011/04/27(水) 23:20:30.27おっ、同業者か。
nonrandomパッケージのVignettesを読まれてはいかが。
また、疫学のSIGのMLもあるので、参加されてはどうか。
ttps://stat.ethz.ch/mailman/listinfo/r-sig-epi
0052１３２人目の素数さん
2011/04/28(木) 13:47:23.33ありがとございます。
この休みにでもRのインストールからはじめてみます。
スレ違いかもしれませんが、propensity scoreを用いたマッチングでの解析を行うのに適したソフトは何でしょうか？
0053１３２人目の素数さん
2011/04/28(木) 13:50:34.33nonrandomパッケージのVignettesを読んでからもう一度質問してください。
ttp://cran.r-project.org/web/packages/nonrandom/vignettes/nonrandom.pdf
0054１３２人目の素数さん
2011/04/28(木) 17:38:07.68ありがとうございます。
0055１３２人目の素数さん
2011/04/30(土) 09:11:31.24ttp://cran.r-project.org/bin/linux/ubuntu/
0056１３２人目の素数さん
2011/04/30(土) 13:11:01.00Releaseのファイルサイズが0でまだ使えないね
$ w3m -dump http://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
今あるパッケージはこれだけ。
0057１３２人目の素数さん
2011/05/03(火) 20:04:04.55すでにデータAを学習したネットワークにさらに別のデータBを学習させるにはどうすればいいのでしょうか？先にAとBを結合して読ませるのは時間がかかるために困っています。
他のパッケージなら出来るという情報もお待ちしています。
0058１３２人目の素数さん
2011/05/05(木) 21:13:39.01MySQLに返しつつ、gnuplotで描画しようと思っています
RubyとかからMySQL、R、gnuplotを操作すべきか、
RからMySQLとgnuplotを操作すべきかで迷ってます
みなさんならどんな風にしますか？
0059１３２人目の素数さん
2011/05/06(金) 13:30:23.98このスレ的には、Rの中からMySQLにアクセスして、
Rで描画する方向が第一選択だよな。
Rのグラフィック能力とgnuplotってほぼ同等と思っていたけど、
Rを使っているのにgnuplotを加えて使うのは理由は？
>>57
反応がないから、Rjpwikiの方に書かれてはいかが。
私はnnetパッケージを使わないので分からない。
0061１３２人目の素数さん
2011/05/06(金) 23:00:08.59Rで描けるようになればR+MySQLなのでシンプルでいいですね。
ちょっとやってみます。ありがとうございました。
0062１３２人目の素数さん
2011/05/09(月) 23:39:55.72y_n<-60
n<-100
posterior<-((pi^y_n)*(1-pi)^(n-y_n))/beta(y_n+1,n-y_n+1)
の95％信頼区間を求めたいんですがその関数がよくわからないです誰かわかる人いらっしゃいますか？
0063１３２人目の素数さん
2011/05/10(火) 09:08:37.12Bayesian credible intervalは信頼区間じゃないよ
0064１３２人目の素数さん
2011/05/10(火) 15:28:06.310065１３２人目の素数さん
2011/05/14(土) 21:58:27.95だと上手く出力されないのはどうしてでしょう？
素人の質問ですみません
0066１３２人目の素数さん
2011/05/14(土) 22:42:04.29> 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が表示される。
0067１３２人目の素数さん
2011/05/15(日) 17:30:01.07for文の中ではprint()が必要です。
for (i in 1:5) {print(sprintf("%s", ticker[i,1]))}
普通は>>66の言うとおり、
as.character(ticker[1:5,1])だけどね。
0068１３２人目の素数さん
2011/05/25(水) 00:02:27.53それぞれの20個のデータの平均値と標準偏差を求めるプログラムstatsample(m,data)を作成する。
プログラムソースが分かる人は教えてください。
0069１３２人目の素数さん
2011/05/25(水) 09:29:38.65とりあえず、たたき台を示すので、データ数が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"
0070１３２人目の素数さん
2011/05/25(水) 19:42:32.88datamatrix<-matrix(data,20,length(data)/20)
apply(datamatrix,2,mean)
apply(datamatrix,2,sd)
とやるのが普通じゃないか。
データ数が20で割り切れない場合のwarningも表示されるし。
0071１３２人目の素数さん
2011/05/25(水) 19:45:58.57このシュミレーションを行うために、不良品率ｘ％でn個のサンプルの中に不良品を発生させるプログラムsampling(x,y)を以下の手順で作成する。
(1)不良品の出現確率ｘ/100とする。発生させた一様乱数がｘ/100より小さければ、不良品(1)、そうでなければ、良品(0)とする。
(2)不良品出現確率ｘ/100とした試行をn回繰り返し、不良品が出現した場合の数を求める。
その数がsampling(ｘ,n)の出力とする。
以下の手順をプログラムし、関数sampling(ｘ,n)を作る。
誰か教えてください。
0072１３２人目の素数さん
2011/05/25(水) 21:51:26.020073１３２人目の素数さん
2011/05/26(木) 00:30:46.840074１３２人目の素数さん
2011/05/26(木) 01:25:55.190075１３２人目の素数さん
2011/05/26(木) 09:18:07.75宿題？
ヒントだけ
> runif(10) < 5/100
[1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE
0076１３２人目の素数さん
2011/05/26(木) 22:53:09.940077１３２人目の素数さん
2011/05/31(火) 02:28:47.96見ようとしたら、
.Internal( sample....
となってて、中身が見えません。
他の言語で作成しているのだろうと思いましたけど、
sample関数のソースって、どこにあるんでしょうか？
0078７７
2011/05/31(火) 02:32:14.67CRANからソースを取得しました。
0079あんでぃは非存在 ◆AdkZFxa49I
2011/06/06(月) 14:26:34.020080１３２人目の素数さん
2011/06/07(火) 04:41:09.160081１３２人目の素数さん
2011/06/07(火) 22:07:44.630082１３２人目の素数さん
2011/06/07(火) 22:18:53.870083１３２人目の素数さん
2011/06/17(金) 13:02:41.63株価の日次データが数年分あるんですが
グラフのX軸の目盛りを年次にしたい場合はどうすればいいでしょう？
0084１３２人目の素数さん
2011/06/17(金) 13:23:39.86a <- ts(1:2000,start=c(2011,6),frequency=365)
plot(a)
0085１３２人目の素数さん
2011/06/17(金) 13:35:41.47早速のレスありがとうございます
株など営業日がある場合どうすればよろしいでしょうか？
祝日や年末年始などありますし
0086１３２人目の素数さん
2011/06/17(金) 13:52:29.96xtsパッケージのドキュメントを眺めてから再質問してください。
ttp://cran.r-project.org/web/packages/xts/
日本語が必要なら、Rパッケージガイドブック pp.130-137
0087１３２人目の素数さん
2011/06/17(金) 13:59:31.54ありがとうございます。
パッケージガイドブックを買ってこようと思います。
これとR-tipsで頑張ります
0088１３２人目の素数さん
2011/06/23(木) 16:49:30.620089あんでぃはストーカー ◆AdkZFxa49I
2011/06/23(木) 21:21:17.460090１３２人目の素数さん
2011/06/24(金) 09:20:34.16と書くと
2011/06
と出力されるけど、幅を節約して
11/06
とするにはどうすれば?
0091１３２人目の素数さん
2011/06/24(金) 09:22:44.15%y だ
0092１３２人目の素数さん
2011/06/24(金) 10:43:36.69そのあたりは共通だから覚えておいて損はないよ。
> format(Sys.Date(),"%Y")
[1] "2011"
> format(Sys.Date(),"%y")
[1] "11"
> format(Sys.Date(),"%EY")
[1] "平成23年"
> format(Sys.Date(),"%Ey")
[1] "23"
0093１３２人目の素数さん
2011/06/24(金) 12:32:27.69なんでEが平成なんだろう
0094１３２人目の素数さん
2011/06/24(金) 13:01:48.18マニュアルには、
%E Modifier: use alternative format
と書いてあるよ。つまり日本語ローケルだと、日本語化した形式にするという意味。
>>92を書き込んだ後、MacOSXとWindowsで試してみたが、
どちらも和暦にならなかった。GNUだけかも。ヘルプには環境依存って書いてあるしね。
0095あんでぃはストーカー ◆AdkZFxa49I
2011/06/24(金) 16:34:26.88あんでぃ
0096R入門
2011/06/26(日) 00:02:26.05Rについても勉強したいと思っています。
Rのサイトをみると、結構、「使い方」的なサイトが多く、
汎用言語の視点から書かれた、サイト、本が少ない(？)ような気がします。
たとえば、ほかの言語と比較しながら書かれたような
本、サイトなどあれば、教えていただきたいのですが。
ありますでしょうか。
0097１３２人目の素数さん
2011/06/26(日) 00:19:39.460098１３２人目の素数さん
2011/06/26(日) 01:03:33.08Rの基礎とプログラミング技法
http://www.amazon.co.jp/dp/4431712186
http://d.hatena.ne.jp/sesejun/20090331/p2
Rと他の言語を比較した本は知りませんが、
汎用言語としてのR言語を勉強するならこの本がお勧めです。
0099１３２人目の素数さん
2011/06/27(月) 09:37:34.48そういえば、Rユーザ会で、韓国から来た演者が汎用言語の視点から、Rの問題点などを
発表していたな。興味がなかったので内容は覚えていない。
汎用言語の視点で書かれた書籍はないと思う。
しいてあげるなら、>>98が薦めるリゲスさんの本しかない。
010085
2011/07/04(月) 10:41:00.37xtsをインストールしたんですが
データフレームを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
といった具合に２列のデータなのですが？
0101１３２人目の素数さん
2011/07/04(月) 11:05:02.37ここに書くのではなくて、
それを書いた人に凸してみたら？
> 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
まぁ、参考になれば。
0102１３２人目の素数さん
2011/07/14(木) 14:34:30.11相補誤差関数 erfc を qnorm 関数から定義すると、 qnorm が実数に対してしか定義されていないためエラーになってしまいます。
0103１３２人目の素数さん
2011/07/15(金) 02:10:37.99＞ 汎用言語の視点から書かれた、サイト、本が少ない(？)ような気がします。
Rの本を書いているのは、大学の先生が多いからね。
そして、大学の先生ってのはコンピュータのことや、
プログラミングのことをわかってないから。
Rの使い方だけの説明になってしまうのは
しかたがないと思いますよ。
0104１３２人目の素数さん
2011/07/15(金) 11:23:45.87○ 情報科学が専門ではなく疫学やら統計屋が仕事の教員
0105１３２人目の素数さん
2011/07/16(土) 00:45:00.68○ 大学の先生
Ｘ 情報科学が専門ではなく疫学やら統計屋が仕事の教員
特に情報科学の先生や院生などは、知ったかの馬鹿ばっかりです。
国立情報学研究所を見ればよくわかるでしょ？
税金を使ってＰＣで遊んでるだけ。
全員がコネ採用。
0106１３２人目の素数さん
2011/07/16(土) 01:04:55.28唯のお飾りか？
0107１３２人目の素数さん
2011/07/16(土) 16:16:24.06税金を無駄遣いしながら、馬鹿同士楽しく
PCで遊んでいる組織ですね。
どいつもこいつも馬鹿ばっかり。
コネがないと絶対に採用されません。
公募する前から採用者は決まってます。
0108１３２人目の素数さん
2011/07/16(土) 16:49:55.560109１３２人目の素数さん
2011/07/16(土) 20:02:28.47統計数理研究所
この二つは廃止で結構です。税金の無駄です。
0110１３２人目の素数さん
2011/07/16(土) 20:09:53.71バカばっかりやなw あれで数学やっているつもりのもいるからw
頓珍漢なことばかりを言って来るw
0111猫は嘘吐き ◆MuKUnGPXAY
2011/07/16(土) 20:20:25.680112１３２人目の素数さん
2011/07/16(土) 20:58:33.45しかし、単なる就職先のひとつとなっており、
学歴とコネのみで判断されるから。
実力で研究者になれるところって、少ないと思う。
0113１３２人目の素数さん
2011/07/16(土) 21:06:58.590114１３２人目の素数さん
2011/07/16(土) 21:19:33.030115猫は無職 ◆MuKUnGPXAY
2011/07/17(日) 03:35:47.40猫
0116１３２人目の素数さん
2011/07/17(日) 13:55:02.20最近の調査では勤労者の平均年収が
３７０万円ぐらいらしいんですよ。
ただし最頻値は２５０万円ぐらいですから、
一部の高給取りが平均を上げていると思います。
また、公務員の平均は５００万円ぐらいらしいですが、
公務員だけの平均年収は公表されていません。
研究と称してPCで遊んでいるだけですから、
年収１５０万円でも十分だと思います。
0117猫の夢と痴漢の埋葬 ◆MuKUnGPXAY
2011/07/17(日) 14:08:35.01『研究と称してＰＣで遊んでる人』だけをそういう待遇にしたら良いと
思います。但しソレをどうやって判定スルのかは現実には難しく、なの
で『虚偽院生と優秀院生の分離』と同様の困難が伴い兼ねませんね。
では質問ですが、虚偽院生に対してはどういう扱いが宜しいのでしょう
かね？ やはり追放ですかね？
猫
0118１３２人目の素数さん
2011/07/17(日) 14:12:09.83統計数理研究所では女の子の飲み物調査という論文もあったぞw
0119１３２人目の素数さん
2011/07/17(日) 14:24:35.05虚偽院生とは、論文を書かない院生のことですか？
全く論文を書く気のない院生は追放でしかたないと思います。
少なくとも論文を何本か書いていて、研究費も助成金も１円も
受け取っておらず、アルバイトをして授業料を納めている奴は、
他の人の研究の邪魔をしないかぎりは、追放でなくて放置でいいと思います。
つまり、まあおれのことですけどね。
追放すべきは高い給料、研究費や助成金をたんまりもらっているのに
結果を出さない馬鹿な奴です。
これは院生だけではなく、助手も助教、準教授、教授も含まれます。
それらは能力と業績だけで判定されるべきだと思います。
その業績というのをどう判定するのかが難しいわけですけどね。
0120猫の夢と痴漢の埋葬 ◆MuKUnGPXAY
2011/07/17(日) 14:29:59.14全くその通りだと思います。私が言う虚偽院生とは：
★★★『基礎体力も無ければ動機も無く、唯大学院に入学してしまったゾンビ』★★★
を指し示す言葉です。ですから恐らくは貴方の事ではないでしょう。放置
されていても自分から何かが出来る人は決してゾンビではアリマセンから。
猫
0121１３２人目の素数さん
2011/07/17(日) 14:41:00.64つまり>>119の発言にまったく同意しているという解釈でよいですか？
0123１３２人目の素数さん
2011/07/17(日) 14:51:47.93眼醒めたらゾンビみたいだったが、観てたら生気がもどってきた
この辺のクオリティの高さは邦楽にはない、流石だとおもう。
0124１３２人目の素数さん
2011/07/17(日) 23:18:34.72周囲と意見交換が出来るということは、非常に大きな能力であり財産だと思うよ
そして研究室への貢献度という点で言えば
論文は書かないけどPCで遊べる程度の実力はあって、周囲と意見交換が出来る奴 ＞ 一人黙々と論文を書く奴
だと自分は思う
0125１３２人目の素数さん
2011/07/17(日) 23:41:22.12月１０万でボーナス（一ヶ月)も出て国家的身分は安泰
やりがいなんかまったくないと思うかもしれないが、言うことなしの職場じゃないか
0126１３２人目の素数さん
2011/07/17(日) 23:45:41.25それこそ、そこのなんとか研究所に問い合わせをして詳細資料（統計データ)請求すればいいんじゃないの？
ただ気になったんだけど、rlang使ってるなら「加重平均」ってなんのことかはしってるよね。
0127１３２人目の素数さん
2011/07/18(月) 02:46:14.960128猫は爺惨状の蚊軍 ◆MuKUnGPXAY
2011/07/18(月) 03:33:15.34その考え方は決して認められない。大学は馬鹿がママゴトして遊ぶ場所
では断じてない。当然に：
★★★『黙々と論文を書く奴＞論文は書かない奴』★★★
であり、PCで遊ぶかどうかや、周囲と意見交換が出来るかどうかなんて、
そんな事はどうでもヨロシ。
論文を書くかどうかだけが問題である。勘違いも甚だしい。ケシカラン。
猫
0129１３２人目の素数さん
2011/07/18(月) 03:34:40.12まずお前が実績を出してから批判しような
0130１３２人目の素数さん
2011/07/18(月) 04:44:10.88「左右非対称分布の平均値は分布を的確に示せない」
のは統計の基本常識です。
集計方法と指標の定義: さまざまな平均値/中央値/最頻値
http://web-tan.forum.impressrd.jp/e/2008/07/11/3428
0131１３２人目の素数さん
2011/07/18(月) 05:55:50.190132１３２人目の素数さん
2011/07/18(月) 06:00:29.170133１３２人目の素数さん
2011/07/18(月) 07:56:35.430134１３２人目の素数さん
2011/07/18(月) 07:59:15.180135１３２人目の素数さん
2011/07/18(月) 10:43:44.13どこの大学でも統計の先生の論文は下らんものばかり
0136１３２人目の素数さん
2011/07/18(月) 10:54:16.080137１３２人目の素数さん
2011/07/18(月) 11:20:01.810138１３２人目の素数さん
2011/07/18(月) 11:22:33.75女子高校生の飲み物調査
女子中学生の飲み物調査
女子小学生の飲み物調査
女子幼稚園生の飲み物調査
２０代OLの飲み物調査
３０代のOLの飲み物調査
以下、女と男にかえても論文が書けます 統計数理研究所の教授がこうした論文を書いています
0139１３２人目の素数さん
2011/07/18(月) 11:24:51.54同上 ２ ２部上場企業の場合
同上 ３ ヘラクレス上場企業の場合
同上 ４ 町役場の場合
以下 延々と書けます。
0140１３２人目の素数さん
2011/07/18(月) 11:30:15.870141１３２人目の素数さん
2011/07/18(月) 15:21:16.300142１３２人目の素数さん
2011/07/18(月) 18:25:28.32等高線を書きたいのですが、情報量が少ないため、contour関数では書けません。現在与えられているのは、20個の点に対しての座標(x,y,z)です。2次元上の布置まではできているのですが、等高線を用いてz座標を表現することができません。
どの関数を扱えばいいかなど、わかる方がいらっしゃいましたら教えてください。
0143１３２人目の素数さん
2011/07/18(月) 18:59:47.110144１３２人目の素数さん
2011/07/18(月) 22:51:46.35つか、授業聞いてなかったの？
全部説明したはずだけどね。
geoRで作成しなさいって言ったでしょ？
単位あげないよ！
0145１３２人目の素数さん
2011/07/31(日) 19:04:21.64Rコマンダーでcsv経由。
RGUIはMSGhosicに設定。
コマンドでフォント指定以外に
対策あれば教えてください。
0146１３２人目の素数さん
2011/08/01(月) 00:30:07.52Rというよりそのcsvファイルがくせもの
解決策は2つ。
1. csvの文字コードを文字化けしないものに変更してからRに読み込む
2. Rの中でiconv()で文字コードを変換する
好きな方を選べ
0147１３２人目の素数さん
2011/08/01(月) 05:59:12.09了解です。
クリップボードを経由してみたら、ますます
酷くなりました。
特に列の変数名がボロボロです。
これはなにゆえでしょうか？
■ このスレッドは過去ログ倉庫に格納されています