【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
0032１３２人目の素数さん
2011/03/08(火) 19:58:24.45模試の大学合否判定って、A判定が80％以上、E判定が20％未満だな。
97.5％以上とか2.5％以下とかにしないのはなぜなんだろう。
0033１３２人目の素数さん
2011/03/08(火) 20:05:19.440034１３２人目の素数さん
2011/03/08(火) 20:27:35.03受験者数と合格者数を考えてみれば
003532
2011/03/08(火) 21:28:41.29p値が5％以下じゃないとモデルとして不適切などという手法ではないのでしょうね。
0036１３２人目の素数さん
2011/03/11(金) 11:14:42.550037１３２人目の素数さん
2011/03/16(水) 18:58:17.39筑波のCRANは停電というより地震による直接被害で、電気はきていますが、
漏水がひどい建物にあるので念のため落としています。 -- 岡田 2011-03-16 (水) 10:27:28
0038１３２人目の素数さん
2011/03/17(木) 18:18:07.470039１３２人目の素数さん
2011/03/23(水) 14:38:52.94連想リストのことなの？
0040１３２人目の素数さん
2011/03/24(木) 18:04:13.33?formula
>>22
Infになるなら、そもそもモデルが不適切
>>24
?class
>>31
R環境の中からLANGを変更できるよ
>>36
誰も統計を取っていないが、昔はEmacs以外の選択肢が事実上なかった
>>39
全然違う
全レスって嫌われるのかなぁ
0041１３２人目の素数さん
2011/03/24(木) 18:46:14.980042１３２人目の素数さん
2011/04/11(月) 08:54:06.06kernlab の中の SVM 関数 ksvmを使っているのですが、
例えばkernlab の中spamを分類すると、決定関数にsign関数で
符号化されて「spam」か「nonspam」の2値で結果が出力されます。
これを2値に符号化される前の決定関数で出力する方法はあるのでしょうか。
どなたかご存知の方はいらっしゃるでしょうか。
0043１３２人目の素数さん
2011/04/12(火) 16:56:34.83今の所特に日本語周りの不具合は無し
というわけで不要になったTinn-Rをアンインスト
NppToRは窓配置フリーが便利なのでとりあえず残しておく
0044１３２人目の素数さん
2011/04/12(火) 17:16:53.99Windows版ではそうかもしれないが、MacOSXとUbuntuでは
完全に対応していない。
004544
2011/04/13(水) 09:55:41.26MacOSXで、問題があった日本語変数名も使えるようになった。
> Sys.getlocale()
[1] "C/en_US.UTF-8/C/C/C/C"
こんな感じだけど。
Ubuntu amd64の方も、エディタ部分で日本語が怪しかったけど、
これも直った感じ。
0046１３２人目の素数さん
2011/04/13(水) 17:17:16.150047１３２人目の素数さん
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.19■ このスレッドは過去ログ倉庫に格納されています