統計解析R たぶんpart3くらい

**デフォルトの名無しさん** · 2012/06/22(金) 13:33:12.40

あっても良さそうなんですが,需要ないの？

過去スレのうちのひとつ
統計解析R
http://toro.2ch.net/test/read.cgi/tech/1285597198/

**デフォルトの名無しさん** · 2015/06/20(土) 23:30:00.14

bash スクリプト内でRのスクリプトを起動して統計処理させてその結果を
1か0かで戻り値($?)として返すにはどうすれば良いでしょうか？

**デフォルトの名無しさん** · 2015/06/21(日) 00:47:25.51

シェルスクリプト test.sh の内容

#!/bin/bash
r --vanilla --slave < test.R

# コマンドの直後に、終了コードを見る
echo $?

**デフォルトの名無しさん** · 2015/06/21(日) 01:10:40.72

ありがとうございます。
test.R側で戻り値を返す際の指定の仕方が判らなかったのですが、
単に戻り値が入った変数を最終行に書く？

**デフォルトの名無しさん** · 2015/06/21(日) 01:21:55.21

Rで、returnを使えるかどうか、よく知らないけど、

test.R(Rスクリプト)側に、return 0か1を書けば、
シェルスクリプト側に、返り値を返せないのかい？

もし返り値を返せるのなら、その直後に、
echo $? で確認できる

この話題を続けるのなら、名前欄に100と入れておくれ

**100** · 2015/06/21(日) 19:38:19.06

手元のPCにインストールしてやってみたのですが、
test.Rで
nantoka <- 123
nantoka
とやって>>101のようにして $?を確認しても0ですね
return nantoka
と書くと今度は「エラー：予想外のシンボルです」と出て弾かれます
return はtest.R内定義された関数からの戻りにしか使えなさそうです。
やっぱtest.Rを呼び出した側にステータス値を返すのは出来ないんですかねぇ。

**デフォルトの名無しさん** · 2015/06/22(月) 12:27:26.12

>>104
ちょっとまて。
0で正解だろ。
Rの処理が失敗すれば1に、成功すれば0だから、それで正しいだろ。
何がしたいんだ。

------- ここから -----
#!/bin/bash

cat <<EOF > test.R
nantoka <- 123
hoge
EOF

R --vanilla --slave < test.R
echo $?
------- ここまで ------

これならちゃんと期待通り1が返るけど。

**100** · 2015/06/22(月) 23:35:39.56

＞0で正解だろ。
＞Rの処理が失敗すれば1に、成功すれば0だから、それで正しいだろ。

質問の仕方が悪くてすみません。具体的にやりたい事を書くと、mean()で計算した結果が、
ある範囲に入っていたらbashスクリプト側にスタータス値として 1を返し、範囲外だったら0
を返すをやりたいんですが、どう書いたら良いんでしょう。。。

**100** · 2015/06/22(月) 23:57:01.52

Rスクリプト側でステータス値を標準出力に印字して
bashスクリプト側で変数に代入する事で出来ました。
>>105さんのレスがヒントになりました。感謝！
ありがとうございました。

**デフォルトの名無しさん** · 2015/08/02(日) 12:32:42.11

RのpackageがGPL3のとき、それを使ったcodeもGPL3にしなければいけないのでしょうか？

packageがGPL2のときは、どうなるのでしょうか？ライブラリを動的にリンクすると解釈して、自分のコードはGPLにしなくていいのでしょうか？

**デフォルトの名無しさん** · 2015/08/28(金) 17:44:03.14

Rが嫌いすぎる。

1オリジンが気に食わない。
インデックス指定で、やたら+1とか-1を書かないといけない。バグの元凶。

名前空間の概念が皆無なのが気に食わない。デフォルトインストール状態でグローバルな名前空間に、
小文字二文字でしらっと置いてある識別子がただのサンプルデータ(例: lh)だったりする。

文法が気に食わない。「a <- b < -1」って書くべきところ、間違ってマイナスの前のスペースを
入れ忘れると、警告も無く別の意味になる。
設計したやつ頭悪すぎ。

できれば使いたくない、抹殺したいんだけど、仕事の分野で好んで使う人が多いんだよなぁ。

**デフォルトの名無しさん** · 2015/08/28(金) 18:52:05.79

>>109
はははは、ウケる！

> デフォルトインストール状態でグローバルな名前空間に、
> 小文字二文字でしらっと置いてある識別子がただのサンプルデータ(例: lh)だったりする。
これは君の勘違いだから、よく確認して。

**デフォルトの名無しさん** · 2015/08/29(土) 19:00:44.80

そこらのLinuxディストリのリポジトリからインストールすると、
datasetsパッケージから読まれた状態にならない？

**デフォルトの名無しさん** · 2015/08/31(月) 02:17:53.94

googlabilityの低さだけでRは唾棄するに値する。

**デフォルトの名無しさん** · 2015/08/31(月) 02:41:18.19

**デフォルトの名無しさん** · 2015/08/31(月) 13:01:31.84

>>111
そうだね。
$ R --vanilla
[snip]
> search()
[1] ".GlobalEnv" "package:stats" "package:graphics"
[4] "package:grDevices" "package:utils" "package:datasets"
[7] "package:methods" "Autoloads" "package:base"
でも、
> ls(envir = .GlobalEnv)
character(0)
だよ。

**デフォルトの名無しさん** · 2015/08/31(月) 13:41:28.94

なお、起動時にライブラリが自動読み込みされるのが嫌なら、
~/.Rprofileなどで、
options(defaultPackages="")
とでもしておけばよいのでは。

**デフォルトの名無しさん** · 2015/08/31(月) 15:23:52.28

なにバッドノウハウひけらかしてドヤ顔してんのｗｗ
Rが低能向けゴミ言語であることの証左にしかならんね。

**デフォルトの名無しさん** · 2015/09/01(火) 01:31:57.27

じゃあフリーでLinuxでもWindowsでも使える
R以上の統計解析言語教えてくれ

**デフォルトの名無しさん** · 2015/09/01(火) 10:36:16.50

>>116
コンピュータ言語だと思うから、そのような感想が出てくるのでは？
統計ソフトウェアのマクロ機能と認識したらいい。

**デフォルトの名無しさん** · 2015/10/23(金) 10:35:38.68

すみません、質問です。Rの教科書をめくりながら、length(grep("^[0-9]{8}", code))
というのを理解しようとしているのですが、"^[0-9]{8}"というのはどういうことですか？

**デフォルトの名無しさん** · 2015/10/23(金) 21:46:06.88

どう見ても正規表現だらうが
r以前の問題

**デフォルトの名無しさん** · 2015/10/24(土) 00:07:21.80

^[0-9]{8}

^ は行頭。
[0-9]は、0～9のうちのどれか、つまり数字。
{8}は、直前の文字が、8回繰り返す

つまり、行頭に、8桁の数字がある行にマッチする

**デフォルトの名無しさん** · 2015/10/24(土) 10:04:48.60

ありがとうございます、正規表現でぐぐって勉強しました。

**デフォルトの名無しさん** · 2015/11/21(土) 00:24:32.92

striptimeで、時刻型に変換したデータを横軸にトレンドグラフをplotしています
このグラフに縦線を入れたいのですが、可能でしょうか？

date = c(20151101235959 , 20151103235959)
value = c(10 , 20)
plotdate = striptime(date , "%Y%m%d%H%M%S")
plot(plotdate , value , type="l")

ここまでは可能でした。

x軸に縦線を入れたいのですが、ablineでの描き方がわかりません
どなたかご教示いただけないでしょうか

**デフォルトの名無しさん** · 2015/11/24(火) 22:19:32.58

>>123
縦線って
> for(i in seq(from=plotdate[1], to=plotdate[2], length.out=10)){abline(v=i)}
こんな感じのものを求めているの？

**デフォルトの名無しさん** · 2015/12/18(金) 19:52:52.36

すみません、どうしてもわからなくて質問です。
a<-1:10
b<-1:10
for(i in 1:10){
for(j in i+1:10){
if(a[i]!=b[j]) print(

**デフォルトの名無しさん** · 2015/12/18(金) 19:54:57.43

すみません、どうしてもわからなくて質問です。
a<-1:10
b<-1:10
for(i in 1:10){
for(j in i+1:10){
if(a[i]!=b[j]) print("hello")　}}
これを実行したらif文の所にエラーが出るんです。どうしたら解決できますか。教えてください、お願いします。

**デフォルトの名無しさん** · 2015/12/20(日) 11:31:46.28

「i+1:10」のところの演算子で、+よりも:が優先されるから。
「(i+1):10」とすればよい。

**デフォルトの名無しさん** · 2015/12/22(火) 17:32:24.26

ありがとうございます。聞かないと分かりませんでした。

**デフォルトの名無しさん** · 2016/02/15(月) 22:29:22.41

R言語徹底入門、初版は何箇所か誤訳があるので購入は様子見。

**デフォルトの名無しさん** · 2016/02/16(火) 23:39:43.76

増刷がかかるか
誤訳が修正されるか

と考えると買っちまうかな、「徹底入門」じゃなくて「徹底解説」。

**デフォルトの名無しさん** · 2016/02/17(水) 03:10:57.33

増刷かかるほど売れないと思う。
でも仮に誤りがあったとしても得られるものは少なくないと思う。
岩波のWEBサイトで数ページPDFで読める。

R言語徹底解説
http://www.kyoritsu-pub.co.jp/bookdetail/9784320123939

**デフォルトの名無しさん** · 2016/02/17(水) 16:47:27.88

岩波w

**デフォルトの名無しさん** · 2016/02/19(金) 15:53:33.09

第二刷まだー？

**デフォルトの名無しさん** · 2016/02/19(金) 17:22:28.25

そろそろ第三刷

**デフォルトの名無しさん** · 2016/02/19(金) 17:25:50.06

>>134
いろいろ直ってる？

**デフォルトの名無しさん** · 2016/02/20(土) 11:42:28.63

Pythonのyieldに相当するものってないの?

**デフォルトの名無しさん** · 2016/02/26(金) 21:22:23.14

書泉ブックタワーで一刷が山積みになってたから
各種の誤りが訂正された二刷が出るのは当分先になるのかな？

**デフォルトの名無しさん** · 2016/02/28(日) 16:19:00.69

買った。第一刷。ハードカバーとは知らなんだ。

**デフォルトの名無しさん** · 2016/02/29(月) 22:28:35.19

技術書なんて誤訳だらけだろ。校正してるのは文系なんだから。
ほとんどがマイナー、専門すぎてサイトに修正一覧出て初版で終わりだよ。
買いそびれたら最後、ボッタクリ価格でアマゾンで買うことになる。
手に入らないと途端に、誤植もなく、素晴らしい書籍とかﾃﾞﾀﾗﾒなレビューが書かれて価格が釣り上がる。

**デフォルトの名無しさん** · 2016/02/29(月) 22:45:23.04

なんかあったの

**デフォルトの名無しさん** · 2016/02/29(月) 22:52:15.59

過去に散々酷い目にあったよ。一番酷いのはPrograming C++の糞訳本だったかな。
あと確率統計の本で題名詐欺の本もあったよ。まじみんな気をつけたほうがいいよ。アマゾンのレビューだけは。

**デフォルトの名無しさん** · 2016/03/01(火) 13:16:03.01

こんなとこで文句言ってるなら原著で読めって。

**デフォルトの名無しさん** · 2016/03/01(火) 13:59:13.93

>>142
今日もあっちこっちで馬鹿レスしてるなおまえってw

**デフォルトの名無しさん** · 2016/03/08(火) 03:16:56.45

>>141
最近のアマゾンのレビューはほんと嘘だらけ。

**デフォルトの名無しさん** · 2016/04/15(金) 16:50:29.31

R-Sourceが404な気がする

**デフォルトの名無しさん** · 2016/04/19(火) 23:37:39.26

もしかしてR-Tipsなくなってる？

**デフォルトの名無しさん** · 2016/06/15(水) 18:56:31.95

R言語徹底解説
八重洲ブックセンターで第三刷があった。
平積みの方は第一刷だったけどｗ

**デフォルトの名無しさん** · 2016/07/23(土) 19:10:41.11

Rでは２次元座標に円を描く事はできないんですか？

**デフォルトの名無しさん** · 2016/07/23(土) 19:46:19.77

できるよ
でも今なら Julia がおすすめでっせ

**デフォルトの名無しさん** · 2016/07/24(日) 11:38:41.23

どうやってできるんですか？

**デフォルトの名無しさん** · 2016/12/01(木) 20:05:10.85

過疎しているときに申し訳なくお聞きしたいのですが、
N(10,0.5^2)の,7.5以上の上側確率をRで一発で求めるコマンドはないでしょうか？

**デフォルトの名無しさん** · 2016/12/01(木) 20:14:41.90

すいません自己解決しました
pnorm(-(7.5-10.0)/0.5)でした
スレ汚しすいませんでした

**デフォルトの名無しさん** · 2016/12/01(木) 20:26:09.26

おめ

**デフォルトの名無しさん** · 2017/05/05(金) 07:53:03.19

Rを使っている人ならわかる人もいると思うので
GNU科学技術計算ライブラリ　リファレンス・マニュアル
https://ja.wikipedia.org/wiki/GNU_Scientific_Library
のどの関数を使って、どの関数を数値積分すればよいのか、
教えてほしい。
γ関数　https://ja.wikipedia.org/wiki/%E3%82%AC%E3%83%B3%E3%83%9E%E9%96%A2%E6%95%B0　を
使って　ControlChartConstant2Double = gsl_sf_gamma(InSizeDouble / 2) * Sqr(2 / (InSizeDouble - 1)) * gsl_sf_gammainv((InSizeDouble - 1) / 2)
とすると、管理図用定数 c4 が得られる。
d2, d3 を求めるための
http://excelshogikan.com/qc/qc11/controlchartconstant.html
φなんたらが、GSLのどの関数なのか、数値積分のどの関数を使うのか、
がわからない。

積分は、わからなければ台形公式でエィヤーとやってしまうけど。
n20位までは数値表があるのだが、nをちょっと大きくして、24時間ことのn24とかやろうとすると、のっていない。
１回計算してしまえば、数値表を使って参照するだけなのでRで処理で切ると思う。

**デフォルトの名無しさん** · 2017/05/05(金) 10:57:47.64

日本語を勉強する

**デフォルトの名無しさん** · 2017/05/05(金) 11:12:35.50

それな

**デフォルトの名無しさん** · 2017/07/20(木) 19:18:28.80

Sはどこへ行ったか

**デフォルトの名無しさん** · 2017/08/15(火) 19:07:36.81

羽鳥の本は英語ならたいていpdfで無料で読めるから、文句いう人は英語で読めばよいと思うの

**デフォルトの名無しさん** · 2017/09/05(火) 00:29:28.39

使い始めて1ヶ月
未だにエクセルより便利な場面に出会いませぬ
回帰分析程度じゃ意味ないかな

**デフォルトの名無しさん** · 2017/09/05(火) 04:24:09.07

したい事と選んだ手段がズレてるのでは?

**デフォルトの名無しさん** · 2017/09/05(火) 10:33:35.89

まずは、セルでチェックできる程度の規模のデータか、あるいはもっと大きいのかによるんじゃないかな。
そして作業として簡単なものなら、多少表計算ソフトに歩があると思う。
でもこっちに慣れるとぜんぶこっちでやっちゃうけどねｗ
欠点としてはライブラリの秩序のなさかなｗ

**デフォルトの名無しさん** · 2017/09/05(火) 12:46:08.24

RMarkdownを使い始めると分析とレポート作成が同時にできるので便利さを感じると思うけどな

**デフォルトの名無しさん** · 2017/09/06(水) 12:00:48.68

Pandas

**デフォルトの名無しさん** · 2017/09/06(水) 17:41:06.32

>>162
うむ。あれはいいものだ。

**デフォルトの名無しさん** · 2017/09/07(木) 08:59:11.93

jupyter

**デフォルトの名無しさん** · 2017/09/09(土) 10:59:21.18

なるほどと思ったブログがあったから貼っとくよ
yhat | R for Excel Users
http://shotwell.ca/blog/post/r_for_excel_users
ある特定の環境での成功体験が新たな環境への適応の障害となるのは
エクセルからRに限ったことではないと思う

**デフォルトの名無しさん** · 2018/02/09(金) 21:23:26.98

Rが盛り上がらなくて悲しい

**デフォルトの名無しさん** · 2018/02/10(土) 22:42:50.82

大盛り上がりと言ってもいいんじゃない？
https://insights.stackoverflow.com/trends?tags=r

**デフォルトの名無しさん** · 2018/02/11(日) 00:58:01.63

上限100%の上部90%以上カットしたらそうなるわ

**デフォルトの名無しさん** · 2018/02/11(日) 11:54:36.80

>>167
それは仕方ないと思うよ
統計解析手法開発用には向いているソフトだけど
統計解析を利用した処理アプリで使うのに向いている
ソフトウエアというわけじゃなく
無償なんで大学でも採用して広まっているという話だから。

**デフォルトの名無しさん** · 2018/02/11(日) 12:28:34.22

放送大学の講座は親切
https://vod.ouj.ac.jp/view/ouj/#/navi/vod?ca=496

**デフォルトの名無しさん** · 2018/02/11(日) 13:02:45.83

>>157
Rになった
S-plusやってたTIBCOもRの別実装始めたくらい

**デフォルトの名無しさん** · 2018/02/11(日) 13:37:06.23

そりゃ当然の流れやな。

**デフォルトの名無しさん** · 2018/02/11(日) 17:47:44.11

>>168
右肩上がりなのはとても良いこと
純粋なプログラミング言語じゃないから増えてるだけで御の字なのかねえ

**デフォルトの名無しさん** · 2018/02/14(水) 19:42:52.67

データをこねくり回してグラフ書くにはexcelより融通が利いてしかも早い

ような気がする

**デフォルトの名無しさん** · 2018/02/14(水) 22:01:38.62

>>175
そりゃそうでしょうね
Excelはこねくり回さないで使うためのツールだから

**デフォルトの名無しさん** · 2018/02/14(水) 23:37:10.00

グラフはggplot2が強力だよね
facetとかはExcelだと大変

デフォルトのテーマが少し古くさいけど

**デフォルトの名無しさん** · 2018/02/14(水) 23:51:00.08

R使い出してから色々と試すのが簡単なんで面白いわ

**デフォルトの名無しさん** · 2018/02/15(木) 01:24:08.19

機械学習とかビッグデータとかPythonと何かの組み合わせに圧倒されててRとかどマイナーになってしもうた

**デフォルトの名無しさん** · 2018/02/15(木) 06:34:31.40

適材適所でしょうね。
元々設計自体がデカいメモリーがあったあらと言う前提で
設計者にとっては作りが楽だけど、大容量で使う利用者が
対処しないといけないという安易なつくりなんですよね。

その代り少量データを対象とした統計手法研究などには
楽なんで。

**デフォルトの名無しさん** · 2018/02/15(木) 12:27:10.58

Rにはなんといってもtidyverseがあるからなあ
dplyrとかpurrrとか簡単で強力だから解析やるならR

まあ適材適所で、ベイズならstanだし、機械学習や文字ならPythonが優勢、事務処理ならExcelVBAってとこ

**デフォルトの名無しさん** · 2018/02/19(月) 19:54:51.25

>>171
統計って回帰分析から急に難解になる
それにじっと耐える生徒さんは偉い
R初心者の私にはとっても役立つ講義だ

**デフォルトの名無しさん** · 2018/02/22(木) 18:00:01.85

何かいいエディタないかな
Rに最初からついてるエディタ使ってるけど不便

**デフォルトの名無しさん** · 2018/02/22(木) 19:41:30.41

RStudioのエディタじゃだめかい？

**デフォルトの名無しさん** · 2018/02/23(金) 17:28:14.29

>>182
紫布さんかわいい

**デフォルトの名無しさん** · 2018/02/23(金) 20:13:41.57

>>184
RStudioは2,3回インストールしてみたけど良く分からないし重いしですぐ消してしまってたんだよね
慣れると便利なの？

**デフォルトの名無しさん** · 2018/02/23(金) 21:18:49.28

RStudioのエディタは入力補完機能が凄い便利。パッケージの関数だけでなく自分で作った変数や関数も補完してくれる。
関数のヘルプもキー一つで表示できるし、ノートブック機能を使えばエディタ内で実行結果を表示してくれる。

その他、パッケージ管理とかプロジェクト管理とか便利だと思う。

**デフォルトの名無しさん** · 2018/02/25(日) 16:50:33.30

具体的な話を聞くと便利そうな気がしてくる
やっぱりRStudioがベストなのかもね
しばらく使ってみようかな、ありがとう

**デフォルトの名無しさん** · 2018/02/26(月) 12:53:30.71

補完機能なしでプログラミングするのって、武器を持たずに福岡の街を歩くようなものだよ。

**デフォルトの名無しさん** · 2018/02/26(月) 18:57:28.32

こえぇぇぇ

**デフォルトの名無しさん** · 2018/03/03(土) 09:23:09.96

初心者には
Rコマンダーは便利だが
基本が理解するのが課題

**デフォルトの名無しさん** · 2018/03/14(水) 02:41:29.02

R version 3.4.3 でR consoleを起動した直後に乱数を発生させると
set.seedで指定していないのに前回起動時と同じ乱数が出てしまいます。たとえば
runif(5)だと
[1] 0.3108217 0.4556531 0.2194121 0.6496887 0.2677915
が毎回出てきます。
そういのものなのでしょうか？

**デフォルトの名無しさん** · 2018/03/14(水) 02:55:44.49

>>192
R x64 3.4.1（windows10）ではそれを再現できなかった
そういうものではないと思う

**デフォルトの名無しさん** · 2018/03/14(水) 13:07:49.78

>>193　ありがとうございます。そうなんですか。自分の環境はwin7sp1 64bit で
Rは5年前にインストールしその後何度か上書きインストールしたまに使っててたんですが
何らかの理由でRの中の設定が変更されたか破壊されてるんですかね
原因がわかったらまた書きます。

**デフォルトの名無しさん** · 2018/03/14(水) 14:55:10.28

>>192
R3.4.3@Win10だけどそういう状態にはならないねえ。RStudioとかVSCodeとかでも何度かやってみたけど全て違う値。
事前にhtmlwidgetとか使ってない？

**デフォルトの名無しさん** · 2018/03/14(水) 16:28:52.83

>>192です
試しに別のフォルダに新規インストールして何度か起動を繰り返して確認してみたら
正常に毎回違う乱数がでました。
なのでもう新規インストールで実行することにしました。
でも、どこに問題があるのか特定しないとまた起きそうなので少しづつパッケージを削除しつつ確認中・・・
>>195
ありがとうございます。
htmlwidgetsって何？って感じなんですが入ってました。
いつ何のために入れたのかもわからないのですが(笑)

**デフォルトの名無しさん** · 2018/03/14(水) 17:15:01.96

>>192です。
パッケージの問題ではなくて、パッケージの作者の方々疑ってごめんなさい

どうも>>195さんの指摘どおり事前になにかが読み込まれてること気づきました
起動直後にカーソルの↑を押すとなぜか
昔コンソールで打ち込んだコマンドが未だに出てきちゃうんです
たとえば、
> source('~/R/test.r')　

**デフォルトの名無しさん** · 2018/03/14(水) 17:57:53.40

.Rhistoryがおかしくなっちゃったのかな？まあ、原因が分かってなにより？

**デフォルトの名無しさん** · 2018/03/14(水) 18:20:45.09

>>192です。
原因がわかりました。
[Previously saved workspace restored]と出てることに今気づきまして
以前マイドキュメントに保存していたワークスペースが自動で読み込まれるためでした。
でもなぜか、"test.RData"の名前で保存したものは自動で読み込まれず
".RData"のように拡張子のみのファイル名が存在するとそれが読み込まれるようです。

新規インストールのRではがあっても読み込まれないのですが、
試しにC:\Program Files\R\R-3.4.3\bin\x64のフォルダ内にコピペしたら起きるようなので
変な名前やいけない場所には保存しては駄目ってことですね

ありがとうございました。

**デフォルトの名無しさん** · 2018/03/14(水) 18:21:46.77

.Rhistoryでset.seed()使ってるってことなのか