トップページ⇒math
989コメント535KB

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

■ このスレッドは過去ログ倉庫に格納されています
0001132人目の素数さん2011/01/09(日) 17:31:49
R は統計計算とグラフィックスのための言語・環境です。
統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。

●関連サイト
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
0004132人目の素数さん2011/01/11(火) 11:56:09
>>3
RはShemeのような継続がないので全く同じ動作は再現できないと思います。
しかしベクトルの外積を使えば複数のambを組み合わせるのと似た事ができます。

独習 Scheme 三週間 Chapter 14 非決定性 (Schemeのamb)
http://www.sampou.org/scheme/t-y-scheme/t-y-scheme-Z-H-16.html
R言語
http://ja.wikipedia.org/wiki/R%E8%A8%80%E8%AA%9E
ベクトルに関する Tips 大全
http://www.okada.jp.org/RWiki/?%A5%D9%A5%AF%A5%C8%A5%EBTips%C2%E7%C1%B4
0005132人目の素数さん2011/01/11(火) 12:29:53
IncanterはClojureベースでR似の統計計算とグラフィックス環境です。
Rを単体で利用している人にはパッケージ不足な代用品にすぎませんが、
JavaアプリにはJVM系言語で書かれたIncanterの方が組み込みやすいです。

Incanter
http://incanter.org/
Data Sorcery with Clojure
http://data-sorcery.org/
naoya_t@Incanter
http://sites.google.com/site/incanterja/
0006132人目の素数さん2011/01/12(水) 09:20:46
どうせjavaならKawaで書いてくれればいいのに
0007132人目の素数さん2011/01/12(水) 12:40:29
>>6
Incasterのトップページにこういう意味の文が載っています。
「Incanterは動的型付けの関数型言語であるClojureのパワーと、
データアクセスと処理と視覚化のためのJVM上の豊富なライブラリを利用します」
IncasterはSchemeでは成立しないでしょう。

Clojureをつくったわけ by Rich Hickey
http://www.geidai.ac.jp/~marui/clojure/rationale/
>なぜ私は新しいプログラミング言語を書いたのだろうか。
>次のような特徴を持つ言語が見つからなかったから、という理由からだ。
> * 並行処理のために設計された
> * 定評のあるプラットフォームと相性がよい
> * 関数型プログラミングのための
> * Lisp。

> 1.2 関数型プログラミングは良いものだ
>  * Clojureは動的型付けをする関数型言語
>   o すべてのデータ構造がイミュータブルで永続的、再帰をサポート
>   o ヘテロジニアス(異種混合)なコレクション、戻り型
>   o 動的な多態性(ポリモーフィズム)
0008132人目の素数さん2011/01/12(水) 12:44:30
×Incaster
○Incanter
0009132人目の素数さん2011/01/19(水) 23:45:18
統計解析フリーソフト R 【第3章】
http://kamome.2ch.net/test/read.cgi/math/1224142396/527
> 527 :132人目の素数さん [↓] :2010/02/18(木) 17:43:21
> >>526
> ・スレッドのレス数が10以下
> ・レスが一週間書き込まれていない
> 二つの条件が満たされるとスレッドが沈むので
> そうならないようにしておけば沈まないよ

        ∧∧  ミ _ ドスッ
        (   ,,)┌─┴┴─┐
       /   つ.  保  守 │
     〜′ /´ └─┬┬─┘
      ∪ ∪      ││ _ε3
               ゛゛'゛'゛
0010132人目の素数さん2011/01/25(火) 22:44:23
sage
0011132人目の素数さん2011/01/26(水) 10:54:15
誰かRでbmpをimagematrixオブジェクトとして読み込む方法知らないか?
0012132人目の素数さん2011/01/26(水) 17:48:00
11だけど、EBImage使うことにした。
他にいい方法知ってる人いたらおせえて
0013132人目の素数さん2011/01/29(土) 03:04:25
read.matrixでひとつのファイルから複数の行列を読み込んだりすることは可能でしょうか?

初歩的な質問でしたらすいません、教えていただける方宜しくお願いします。
0014132人目の素数さん2011/02/06(日) 00:18:48
R初心者なのですが、小数の精度について解説されたサイト等ありませんでしょうか。
0015132人目の素数さん2011/02/06(日) 00:40:42
>>14
http://www.amazon.co.jp/dp/4901683500/
この本の14章あたりにあるよ
.Machine に浮動小数点演算関連定数のリストが収められてる
0016132人目の素数さん2011/02/13(日) 06:47:34
軸ラベルと凡例が日本語のグラフで、日本語のみに日本語フォントを使用してeps出力する方法をどなたかご存知ないでしょうか?
pdfなら、日本語関連を最後にまとめて、直前にpar(family=Japan1など)でいけるのですが・・・
0017132人目の素数さん2011/02/15(火) 22:46:52
質問です。

単回帰についてですが、データフレームの一列目に被説明変数が格納されており、
残りの列に説明変数が格納されています。

各説明変数を用いて単回帰を行いたいのですが、いちいち変数を指定するのが面倒です。
apply関数を利用しようと思ったのですがうまくいきません。

何か良い方法がありますでしょうか?
0018142011/02/16(水) 21:17:20
>>15
規制で遅くなりました。ありがとうございます。
0019132人目の素数さん2011/02/18(金) 00:46:45
>>17
データフレームが縦横逆になっているということですか?
それなら関数t()を使えば通常のデーター形式になると思いますが。例えば
#irisデータを表示
iris
#縦と横が逆
t(iris)
#もう一度やると元に戻る
t(t(iris))
0020132人目の素数さん2011/02/18(金) 16:26:08
アウトライア分類と変数選択が同時にできる
線形回帰法ないのでしょうか
0021132人目の素数さん2011/02/22(火) 14:38:21.11
チルダ記号ってどういう意味なんでしょうか
関数(変数1 〜 変数2)
なんて使い方です
0022132人目の素数さん2011/02/25(金) 18:20:32.50
質問します。
glmを使い,説明変数の数を変えて複数のモデルを作りました。
複数モデルの中から最もあてはまりの良いものを知りたいので
AICを比較しようと考えたのですが,AICがInfとなり,表示してくれませんでした。
使っているソフトはRです。
AICが非常に小さいくて(おそらく2以下?)表示されないと考えているのですが,
これでは比較できません。
対処法をご存じの方がいらっしゃいましたら,教えていただきたく存じます。
0023132人目の素数さん2011/02/26(土) 08:00:23.96
最近MacにRcmdrをインストールしても文字が汚くならないというとのサイトを見つけたのですが、
ttp://ito-hi.blog.so-net.ne.jp/2011-02-06
「パッケージの依存関係を解決」にチェックしてインストールすると
虹色ぐるぐるでインストールされません

何が問題だと考えられますか?

環境は
iMac
Mac OS X 10.6.6
Intel Core 2 Duo
2.66 GHz

R 2.12.1
0024132人目の素数さん2011/02/28(月) 14:43:10.70
よく変数名に
foo.bar
と . でくぎられた名前があるのですが
Rにもc++みたいなclassがあるのでしょうか
0025132人目の素数さん2011/03/01(火) 17:49:13.03
>>20
lmをransacで繰り返し呼んだら?
0026132人目の素数さん2011/03/02(水) 03:00:19.18
>>23
>何が問題だと考えられますか?

自己解決力もないのにマイナーな環境で動かす冒険心。
0027132人目の素数さん2011/03/02(水) 03:29:20.55
>>26
確かにその冒険心は問題だわねw
0028132人目の素数さん2011/03/03(木) 17:19:55.02
cranに
RANSAC
ないのは、そのぐらい自分で作れってことなの?

lm繰り返し呼べばいいだけなんだけど
0029132人目の素数さん2011/03/06(日) 11:25:08.31
Mac の Carbon Emacs (23.2) で ESS 使い始めたんだけど、
R を呼び出すとメッセージが英語になる。
Terminal.app から R を呼んだ時は日本語表記なんだけど、
ESS 内でも日本語メッセージ使わせるためには何か設定が必要?
0030132人目の素数さん2011/03/06(日) 17:28:12.08
環境変数にLANG=Cでも設定してるんだろうか?
0031292011/03/07(月) 10:57:20.97
~/.bash.rc では LANG=ja_JP.UTF-8 にしていて、
init.el では (set-language-environment 'Japanese) になってる。

ためしに、Terminal から明示的に LANG=ja_JP.UTF-8 emacs & で Emacs を
起動するとメッセージが日本語になった。
デフォルトで Emacs が LANG=C で起動されている??
0032132人目の素数さん2011/03/08(火) 19:58:24.45
Rそのものの話題ではないけど、
模試の大学合否判定って、A判定が80%以上、E判定が20%未満だな。
97.5%以上とか2.5%以下とかにしないのはなぜなんだろう。
0033132人目の素数さん2011/03/08(火) 20:05:19.44
高校生を絶望の谷へ突き落とすような事はしてはいけないのさ。
0034132人目の素数さん2011/03/08(火) 20:27:35.03
>>32
受験者数と合格者数を考えてみれば
0035322011/03/08(火) 21:28:41.29
統計解析の本を見ていて、回帰分析と判別分析の違いかなと思えてきた。
p値が5%以下じゃないとモデルとして不適切などという手法ではないのでしょうね。
0036132人目の素数さん2011/03/11(金) 11:14:42.55
emacsからR使ってる人が多いんでしょうか?
0037132人目の素数さん2011/03/16(水) 18:58:17.39

筑波のCRANは停電というより地震による直接被害で、電気はきていますが、
漏水がひどい建物にあるので念のため落としています。 -- 岡田 2011-03-16 (水) 10:27:28
0038132人目の素数さん2011/03/17(木) 18:18:07.47
床がだいぶ乾いたので、筑波CRANミラー復旧させてみました。 -- 岡田 2011-03-17 (木) 13:17:57
0039132人目の素数さん2011/03/23(水) 14:38:52.94
Rのリストって
連想リストのことなの?
0040132人目の素数さん2011/03/24(木) 18:04:13.33
>>21
?formula

>>22
Infになるなら、そもそもモデルが不適切

>>24
?class

>>31
R環境の中からLANGを変更できるよ

>>36
誰も統計を取っていないが、昔はEmacs以外の選択肢が事実上なかった

>>39
全然違う

全レスって嫌われるのかなぁ
0041132人目の素数さん2011/03/24(木) 18:46:14.98
ステキです^^
0042132人目の素数さん2011/04/11(月) 08:54:06.06
Rでsuppot vector machineの分類の問題が解きたいので、
kernlab の中の SVM 関数 ksvmを使っているのですが、
例えばkernlab の中spamを分類すると、決定関数にsign関数で
符号化されて「spam」か「nonspam」の2値で結果が出力されます。
これを2値に符号化される前の決定関数で出力する方法はあるのでしょうか。
どなたかご存知の方はいらっしゃるでしょうか。
0043132人目の素数さん2011/04/12(火) 16:56:34.83
(祝) RStudio日本語対応!!

今の所特に日本語周りの不具合は無し
というわけで不要になったTinn-Rをアンインスト
NppToRは窓配置フリーが便利なのでとりあえず残しておく
0044132人目の素数さん2011/04/12(火) 17:16:53.99
>>43
Windows版ではそうかもしれないが、MacOSXとUbuntuでは
完全に対応していない。
0045442011/04/13(水) 09:55:41.26
すまん。0.92.44から0.93.75に上げたら、
MacOSXで、問題があった日本語変数名も使えるようになった。
> Sys.getlocale()
[1] "C/en_US.UTF-8/C/C/C/C"
こんな感じだけど。

Ubuntu amd64の方も、エディタ部分で日本語が怪しかったけど、
これも直った感じ。

0046132人目の素数さん2011/04/13(水) 17:17:16.15
線形回帰で変数選択とロバスト推定の同時が出来る方法ないですか
0047132人目の素数さん2011/04/15(金) 20:26:39.00
Rって中身はCなんですか?
遅延評価Cで書くのすごい大変なのに
どうやってるんだろう
0048132人目の素数さん2011/04/19(火) 11:28:58.25
>>47
tar 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さんはこっちに来ないのかな。
0049132人目の素数さん2011/04/27(水) 15:52:48.74
ttp://cran.r-project.org/bin/linux/ubuntu/
Nattyまだこないね。
0050132人目の素数さん2011/04/27(水) 21:40:58.88
言語初心者ですので、教えていただきたいのですが
このソフトpropensity scoreを用いて解析をするときの参考になるサイトや書籍がありましたら教えてください。
0051132人目の素数さん2011/04/27(水) 23:20:30.27
>>50
おっ、同業者か。
nonrandomパッケージのVignettesを読まれてはいかが。
また、疫学のSIGのMLもあるので、参加されてはどうか。
ttps://stat.ethz.ch/mailman/listinfo/r-sig-epi
0052132人目の素数さん2011/04/28(木) 13:47:23.33
>>51
ありがとございます。
この休みにでもRのインストールからはじめてみます。
スレ違いかもしれませんが、propensity scoreを用いたマッチングでの解析を行うのに適したソフトは何でしょうか?
0053132人目の素数さん2011/04/28(木) 13:50:34.33
>>52
nonrandomパッケージのVignettesを読んでからもう一度質問してください。
ttp://cran.r-project.org/web/packages/nonrandom/vignettes/nonrandom.pdf
0054132人目の素数さん2011/04/28(木) 17:38:07.68
>>53
ありがとうございます。
0055132人目の素数さん2011/04/30(土) 09:11:31.24
nattyきたー
ttp://cran.r-project.org/bin/linux/ubuntu/
0056132人目の素数さん2011/04/30(土) 13:11:01.00
>>55
Releaseのファイルサイズが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
今あるパッケージはこれだけ。
0057132人目の素数さん2011/05/03(火) 20:04:04.55
ニューラルネットのパッケージnnetについて質問させてください。
すでにデータAを学習したネットワークにさらに別のデータBを学習させるにはどうすればいいのでしょうか?先にAとBを結合して読ませるのは時間がかかるために困っています。
他のパッケージなら出来るという情報もお待ちしています。
0058132人目の素数さん2011/05/05(木) 21:13:39.01
MySQLからダウンロードしたデータをRで計算して
MySQLに返しつつ、gnuplotで描画しようと思っています

RubyとかからMySQL、R、gnuplotを操作すべきか、
RからMySQLとgnuplotを操作すべきかで迷ってます

みなさんならどんな風にしますか?
0059132人目の素数さん2011/05/06(金) 13:30:23.98
>>58
このスレ的には、Rの中からMySQLにアクセスして、
Rで描画する方向が第一選択だよな。

Rのグラフィック能力とgnuplotってほぼ同等と思っていたけど、
Rを使っているのにgnuplotを加えて使うのは理由は?

>>57
反応がないから、Rjpwikiの方に書かれてはいかが。
私はnnetパッケージを使わないので分からない。
0060572011/05/06(金) 22:06:00.16
>>59
そうしてみます。ありがとうございました。
0061132人目の素数さん2011/05/06(金) 23:00:08.59
gnuplotの方が慣れているというだけです。
Rで描けるようになればR+MySQLなのでシンプルでいいですね。

ちょっとやってみます。ありがとうございました。
0062132人目の素数さん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%信頼区間を求めたいんですがその関数がよくわからないです誰かわかる人いらっしゃいますか?
0063132人目の素数さん2011/05/10(火) 09:08:37.12
>>62
Bayesian credible intervalは信頼区間じゃないよ
0064132人目の素数さん2011/05/10(火) 15:28:06.31
>>63 だね. bayequentistかな
0065132人目の素数さん2011/05/14(土) 21:58:27.95
for (i in 1:5) {sprintf("%s", ticker[i,1])}

だと上手く出力されないのはどうしてでしょう?
素人の質問ですみません
0066132人目の素数さん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が表示される。
0067132人目の素数さん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])だけどね。
0068132人目の素数さん2011/05/25(水) 00:02:27.53
データ配列の最初から20個ずつをまとめ、
それぞれの20個のデータの平均値と標準偏差を求めるプログラムstatsample(m,data)を作成する。

プログラムソースが分かる人は教えてください。
0069132人目の素数さん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"
0070132人目の素数さん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も表示されるし。
0071132人目の素数さん2011/05/25(水) 19:45:58.57
農家の果物の良品、不良品の検査や工場での製品の良品、不良品のチェックなどをシュミレーションにより実験する場合を考える。
このシュミレーションを行うために、不良品率x%でn個のサンプルの中に不良品を発生させるプログラムsampling(x,y)を以下の手順で作成する。
(1)不良品の出現確率x/100とする。発生させた一様乱数がx/100より小さければ、不良品(1)、そうでなければ、良品(0)とする。
(2)不良品出現確率x/100とした試行をn回繰り返し、不良品が出現した場合の数を求める。
 その数がsampling(x,n)の出力とする。
以下の手順をプログラムし、関数sampling(x,n)を作る。
誰か教えてください。
0072132人目の素数さん2011/05/25(水) 21:51:26.02
#include "/dev/console"
0073132人目の素数さん2011/05/26(木) 00:30:46.84
#
0074132人目の素数さん2011/05/26(木) 01:25:55.19
>>69
0075132人目の素数さん2011/05/26(木) 09:18:07.75
>>71
宿題?

ヒントだけ
> runif(10) < 5/100
[1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE
0076132人目の素数さん2011/05/26(木) 22:53:09.94
rbinom(1,n,x/100)
0077132人目の素数さん2011/05/31(火) 02:28:47.96
sample関数について、中でどうやってドローしているんだろうと、
見ようとしたら、
.Internal( sample....
となってて、中身が見えません。
他の言語で作成しているのだろうと思いましたけど、
sample関数のソースって、どこにあるんでしょうか?
0078772011/05/31(火) 02:32:14.67
自己解決しました。
CRANからソースを取得しました。
0079あんでぃは非存在 ◆AdkZFxa49I 2011/06/06(月) 14:26:34.02
あんでぃ
0080132人目の素数さん2011/06/07(火) 04:41:09.16
乱数のところ、ソースみてもわからん。
0081132人目の素数さん2011/06/07(火) 22:07:44.63
乱数だからな
0082132人目の素数さん2011/06/07(火) 22:18:53.87
メルセンヌツイスターとかいうやつか
0083132人目の素数さん2011/06/17(金) 13:02:41.63
すいません、初心者です
株価の日次データが数年分あるんですが
グラフのX軸の目盛りを年次にしたい場合はどうすればいいでしょう?
0084132人目の素数さん2011/06/17(金) 13:23:39.86
>>83
a <- ts(1:2000,start=c(2011,6),frequency=365)
plot(a)
0085132人目の素数さん2011/06/17(金) 13:35:41.47
>>84
早速のレスありがとうございます
株など営業日がある場合どうすればよろしいでしょうか?
祝日や年末年始などありますし
0086132人目の素数さん2011/06/17(金) 13:52:29.96
>>85
xtsパッケージのドキュメントを眺めてから再質問してください。
ttp://cran.r-project.org/web/packages/xts/
日本語が必要なら、Rパッケージガイドブック pp.130-137
0087132人目の素数さん2011/06/17(金) 13:59:31.54
>>86
ありがとうございます。
パッケージガイドブックを買ってこようと思います。
これとR-tipsで頑張ります
0088132人目の素数さん2011/06/23(木) 16:49:30.62
RでVoigt関数(Gauss関数とLorentz関数のconvolution)でフィッティングするには、畳み込み積分を直接書くしかないのでしょうか。
0089あんでぃはストーカー ◆AdkZFxa49I 2011/06/23(木) 21:21:17.46
あんでぃ
0090132人目の素数さん2011/06/24(金) 09:20:34.16
format="%Y/%m"
と書くと
2011/06
と出力されるけど、幅を節約して
11/06
とするにはどうすれば?
0091132人目の素数さん2011/06/24(金) 09:22:44.15
あ、ぐぐったら一瞬で判ってしまった
%y だ
0092132人目の素数さん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"

0093132人目の素数さん2011/06/24(金) 12:32:27.69
おおー
なんでEが平成なんだろう
0094132人目の素数さん2011/06/24(金) 13:01:48.18
>>93
マニュアルには、
%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.05
これまで、Java、Ruby、Scala など使っているものですが、
Rについても勉強したいと思っています。

Rのサイトをみると、結構、「使い方」的なサイトが多く、
汎用言語の視点から書かれた、サイト、本が少ない(?)ような気がします。

たとえば、ほかの言語と比較しながら書かれたような
本、サイトなどあれば、教えていただきたいのですが。
ありますでしょうか。


0097132人目の素数さん2011/06/26(日) 00:19:39.46
ぐぐりにくい、というのがRの最大の弱点
0098132人目の素数さん2011/06/26(日) 01:03:33.08
>>96
Rの基礎とプログラミング技法
http://www.amazon.co.jp/dp/4431712186
http://d.hatena.ne.jp/sesejun/20090331/p2

Rと他の言語を比較した本は知りませんが、
汎用言語としてのR言語を勉強するならこの本がお勧めです。
0099132人目の素数さん2011/06/27(月) 09:37:34.48
>>96
そういえば、Rユーザ会で、韓国から来た演者が汎用言語の視点から、Rの問題点などを
発表していたな。興味がなかったので内容は覚えていない。

汎用言語の視点で書かれた書籍はないと思う。
しいてあげるなら、>>98が薦めるリゲスさんの本しかない。
0100852011/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

といった具合に2列のデータなのですが?
0101132人目の素数さん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

まぁ、参考になれば。
0102132人目の素数さん2011/07/14(木) 14:34:30.11
Rで複素誤差関数 w(x) = exp(-x^2) erfc(-ix) を使いたいのだけれども、これはどうやって定義すればよいでしょうか。
相補誤差関数 erfc を qnorm 関数から定義すると、 qnorm が実数に対してしか定義されていないためエラーになってしまいます。
0103132人目の素数さん2011/07/15(金) 02:10:37.99
>>96
> 汎用言語の視点から書かれた、サイト、本が少ない(?)ような気がします。

Rの本を書いているのは、大学の先生が多いからね。
そして、大学の先生ってのはコンピュータのことや、
プログラミングのことをわかってないから。
Rの使い方だけの説明になってしまうのは
しかたがないと思いますよ。
■ このスレッドは過去ログ倉庫に格納されています