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

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

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

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

**デフォルトの名無しさん** · 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()使ってるってことなのか

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

>>200
うちのマシーンの環境がめちゃくちゃの可能性もありますが、
新規インストール直後のコンソールで保存したワークスペースをそれで読み込んでも
同じ乱数が再現されるので
保存した時のその後の乱数生成パターンも復元されると解釈しちゃいました。

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

はぇ～
そういや最近出たブルーバックスのRの本が売れてるみたいだね

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

>>201
環境の自動復旧は無効にしておいた方がいいよ。ググれば無効にする方法が見かるハズ…

>>202
あの本、ブログか何かで宣伝してたよ。センセーショナルな見出しだったな。

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

この記事だ
ttps://lab-on.jp/article/23

**デフォルトの名無しさん** · 2018/03/20(火) 10:02:37.61

その本を買ってＰＣで真似してます
Ｒコマンダーのお陰で敷居は下がりました
Ｒの壁は
出力データの読み方、意味するところの理解
質のよいデータはあるので頑張る

**デフォルトの名無しさん** · 2018/03/20(火) 11:50:41.25

敷居が高いと思ってる人のためのセミナーが明日、白金台の明学で開催されるよ

http://www.socialpsychology.jp/seminar/seminar2.html

**デフォルトの名無しさん** · 2018/03/20(火) 12:11:32.56

共立出版のRで学ぶデータサイエンスシリーズの15巻とか16巻は
ずっと欠けてるけどいつ出るのですか？

**デフォルトの名無しさん** · 2018/03/21(水) 13:46:02.45

>>204
jupyter notebook と pandas 最強だな

**デフォルトの名無しさん** · 2018/03/27(火) 23:23:26.50

エクセルソルバーのGRGみたいな非線形最適化制約有りをやりたいんだけど、Ｒでもできる？

**デフォルトの名無しさん** · 2018/03/29(木) 21:47:07.36

>>205です
読み終えた、とても勉強になりました
困難さはRではなく統計学でしたが
これは実践するしかないですね
気になったとこは頭に入れておきます
P112の3次元散布図ですが
Microsoft R Open version3.4(64bit)では
作図が表示されませんでした

**デフォルトの名無しさん** · 2018/03/30(金) 15:16:43.56

おめ

**デフォルトの名無しさん** · 2018/05/17(木) 10:20:38.86

Rをかじっても
やはり、考え方がわかってないとだめだね
ということで、
「入門はじめての分散分析と多重比較」を読んでる

**デフォルトの名無しさん** · 2018/05/17(木) 17:50:52.94

Rだけとりあえず使えても、統計学がわかるわけじゃないからねｗ
しかし慣れるとデータ扱うだけでも便利だなぁ。
でもちょっと不安なのが、変数が隠せないから自分で書いた変数とパッケージ内の関数が
偶然だぶったりしないもんなのかな・・・ってところ。
どうなん？

**デフォルトの名無しさん** · 2018/05/17(木) 18:14:17.63

そのあたりは「名前空間」で調べて見るといいよ

**デフォルトの名無しさん** · 2018/05/17(木) 21:17:22.96

ふむ。
ありがとう。

**デフォルトの名無しさん** · 2018/05/19(土) 07:34:08.00

ダブっても問題無いんじゃない？
パッケージの変数を上書きしてしまう理由として考えられるのは
1. 知らずに上書きしてしまう 2. 知っててあえて上書きする
だけど、どちらの場合も元の変数を使わないだけの気がする
安心するには実際にパッケージを壊してみるのも手かも
理由は異なるけど変数スコープの問題でwithを使ってる

my = list (`%>%` = purrr::`%>%`, json = jsonlite::toJSON, id = function (x, ...) x);
my = with (my, {
id (1 : 10) %>% json %>% print; # my環境の変数を使う
rocal = 123; rocal %>% print; # rocalは外から見えない
my $ environment = function (...) 'hello world'; # 勧められないけど無問題
my $ try = tryCatch; # あえて上書き
my;
});
with (my, {
try (typeof (rocal), error = function (ex) 'rocal is local') %>% print;
environment (try) %>% print; # hello world
base::environment (try) %>% print; # my$tryの環境はbase
base::environment (try) = base::environment (); # my$tryの環境をmyに変更
tryCatch (try (stop ('try'), error = function (ex) 'catch') %>% print
, error = id) %>% print; # my$tryが壊れた my$environmentが悪さ
base::try (stop ('try')) %>% print; # base::tryは無傷
base::environment (try) = parent.env (base::environment ()); # my$tryの修理
try (stop ('try'), error = function (ex) 'catch') %>% print; # my$tryの復活
});

てな感じ
多くのサンプルでwithにデータフレームを代入してるけどリストでもオケー
小さな計算をいくつもする時にお手軽だと思う

**デフォルトの名無しさん** · 2018/05/19(土) 09:59:21.70

今から学ぶならJuliaにせい
Rよりずっと速い、書きやすい

**デフォルトの名無しさん** · 2018/05/19(土) 11:12:53.58

うむ

**デフォルトの名無しさん** · 2018/05/19(土) 13:40:20.91

まあ、rocalとか書いちゃうのが

**デフォルトの名無しさん** · 2018/05/23(水) 19:12:36.38

僕の知り合いの知り合いができたパソコン一台でお金持ちになれるやり方
役に立つかもしれません
グーグルで検索するといいかも『ネットで稼ぐ方法　モニアレフヌノ』

6O1Q6

**デフォルトの名無しさん** · 2018/05/24(木) 10:36:25.46

6O1Q6

**デフォルトの名無しさん** · 2018/06/03(日) 23:20:33.44

>>136 あるよ

``` {python}
def joking (a, b):
while True: a, b = b, a + b; yield b;
```

``` {r}
with (new.env (), {
obj = reticulate::py $ joking (0, 1);
sapply (1 : 10, function (...) reticulate::iter_next (obj));
});
```

もっとマシというより楽しい解答はこれかな　読んでないけど

* [Python-like generators in R ・ GitHub](https://gist.github.com/klmr/d10623a0b4c7e1e9a6523eebee4913d1)

で、最終的な答えは「ない」だと思う

技術的な問題というより文化的なところに理由がある気がする
pythonにパイプ演算子がないのと同じ理由
実用的なyieldを実装しようとすると文法の修正が必要になると思う
言語仕様の複雑化という投資が必要になるので
それなりの需要が見込めないと導入には踏み切れないかもね

フィボナッチ数列ではyieldの有り難みは見えないけど
自動微分とは比較にならないほど汎用性の高い有り難い機能だと思う
だけど自動微分と似ているところがあって
どんな複雑なニューラルネットも泣けば自前で微分できるし
どんな複雑な関数のyieldも泣けば自前でイテレーターを作れる

**デフォルトの名無しさん** · 2018/07/04(水) 22:35:50.56

Y4M

**デフォルトの名無しさん** · 2018/07/05(木) 12:07:37.96

fib = lambda n: int(((1+sqrt(5))/2)**n/sqrt(5) + .5)

**デフォルトの名無しさん** · 2018/07/05(木) 13:10:49.78

>>217
初めてききました。

調べたら6年前につくられてるようですが
書籍はほとんどないようですね。
なんで普及してないのでしょうか。

**デフォルトの名無しさん** · 2018/07/05(木) 13:11:04.45

>>217
初めてききました。

調べたら6年前につくられてるようですが
書籍はほとんどないようですね。
なんで普及してないのでしょうか。

**デフォルトの名無しさん** · 2018/07/05(木) 19:49:27.60

Juliaのこと？
本屋では単行本は１冊しか見かけないな。
あと、「データサイエンティスト養成読本 R活用編」というムック本の中に若干の記事があるくらいかな。
他にある？

たぶん書籍が少ないのと、蓄積されたノウハウとか他人の作った関数とか少ないから
まだみんな食いつかないんじゃないかな。
自分もなんか良さげではあると思いつつRで済ましてるしｗ

**デフォルトの名無しさん** · 2018/08/19(日) 22:19:39.16

>>218

http://blogimg.goo.ne.jp/user_image/45/bc/b6502c7b8c0d62a058342de43d82770b.jpg

**デフォルトの名無しさん** · 2018/08/20(月) 13:18:00.41

JuliaはFORTRANっぽくて嫌い

**デフォルトの名無しさん** · 2018/10/13(土) 00:57:54.92

盛り上がっていたので貼っとく
読んでないけど
Faster R with FastR | Hacker News
https://news.ycombinator.com/item?id=18193557
FastRを使う人は少ないかもしれないけど
中のおしゃべりは楽しめるかも

**デフォルトの名無しさん** · 2018/10/21(日) 20:53:08.73

個人的にはRotaさんの名前がここで出てきたことに驚き
こういうところがHNの面白さかも
Deriving the Normal Distribution | Hacker News
https://news.ycombinator.com/item?id=18261892

正規分布からの連想で
コメントにあるインタラクティブなデモに感心
A tutorial on Principal Component Analysis | Hacker News
https://news.ycombinator.com/item?id=18256048

ブログの最後に書いてあるけど
現状はインタラクティブを実現するには
鬼プログラミングが必要なのかもしれない

インタラクティブからの連想で
Distill ― Latest articles about machine learning
https://distill.pub/

**デフォルトの名無しさん** · 2018/10/21(日) 21:25:09.40

ggmap使えなくなってた
Googleにクレカ登録が必須なのか

**デフォルトの名無しさん** · 2018/10/21(日) 21:51:09.06

ggmapってなんだか知らないけどこの辺の話かな？

Change in Google Geocoding API billing Issue #227 dkahle/ggmap GitHub
https://github.com/dkahle/ggmap/issues/227

お買い上げありがとう御座います

**デフォルトの名無しさん** · 2018/11/17(土) 16:28:05.14

>>39
20年後になりそう

**デフォルトの名無しさん** · 2018/11/18(日) 12:04:28.90

もう5年経ったか
機械学習データマイニングω流行ってるのに1/4スレも消費していない
次の3/4スレは15年じゃまだ余るだろう

**デフォルトの名無しさん** · 2018/11/23(金) 11:50:10.19

数学板のRスレのほうが賑わってるのがRらしい

**デフォルトの名無しさん** · 2018/11/23(金) 12:00:16.91

過疎ってるな

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

数学板って基地外が多いので嫌いｗ

**デフォルトの名無しさん** · 2018/11/23(金) 12:29:41.52

基地外かどうかは知らんがほぼ全員コミュ障

**デフォルトの名無しさん** · 2018/11/23(金) 12:44:03.46

そしてリアルの理学部数学科の雰囲気に非常によく似ている。

**デフォルトの名無しさん** · 2018/11/23(金) 19:08:15.85

Rのマニアルを見たら関数のパラメータがやから多くてびびった

**デフォルトの名無しさん** · 2018/12/02(日) 00:14:30.40

Rは難しいね
=と<-が違うとか、forループ使うなとか

**デフォルトの名無しさん** · 2018/12/02(日) 00:17:55.74

Rは簡単やろ？forは禁止されてなく使っても構わない。
けどベクトル演算に持ち込んだほうが圧倒的に速い処理があるってだけの話。

**デフォルトの名無しさん** · 2018/12/02(日) 00:57:29.97

プログラミングの基礎がないから、データフレームをnestしてmutateからのmap_dblという教科書的な処理すら少したつと忘れてしまうわ

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

それは難しいやつや

**デフォルトの名無しさん** · 2018/12/03(月) 00:40:44.81

コンパイラー以外では余り語られないかもしれないパーサーの話
The new pqR parser, and R’s “else” problem
https://www.r-bloggers.com/the-new-pqr-parser-and-rs-else-problem/
elseはC言語もしくはもっと古い言語の理論的バグと言ってもいいかも
Dangling else - Wikipedia
https://en.wikipedia.org/wiki/Dangling_else
yacc/bisonの問題にエラーメッセージが理解不能なことがある
Yacc is dead
https://arxiv.org/abs/1010.5023
yaccそのものというよりbisonとの間で情報の欠落が起きるらしい
Yacc is Not Dead (2010) | Hacker News
https://news.ycombinator.com/item?id=8782218
パーサーを自動生成するのは人間がプログラミングするには複雑過ぎるからだけど
複雑 in 理解不能 out
あれ？何処かで見たような

おまけ：現在のニューラルネットの興隆は彼との対話が起点らしい
Geoff Hinton Facts
http://yann.lecun.com/ex/fun/

**デフォルトの名無しさん** · 2018/12/14(金) 21:02:39.08

>232 斜め読みできる分量ではないけどオーサムらしい
Free online book: Geocomputation with R, a book on geographic data analysis, visualization and modeling : rstats
https://www.reddit.com/r/rstats/comments/a5tmht/free_online_book_geocomputation_with_r_a_book_on/
Geocomputation with R - the afterword | R-bloggers
https://www.r-bloggers.com/geocomputation-with-r-the-afterword/
各トピック毎に使えるソフトウェアの紹介があるので乗り換えの参考になるかも
中に書いてあるようにギットハブにもリソースがある
GitHub - Robinlovelace/geocompr: Open source book: Geocomputation with R
https://github.com/Robinlovelace/geocompr#geocomputation-with-r
Geocomputation with R
https://geocompr.github.io/

**デフォルトの名無しさん** · 2018/12/15(土) 11:21:36.85

>>247
ありがとう
こんな良さげな資料が公開されてたとは
でもggmapのように地図(日本語)に図形やテキストをプロットした画像を出力する手法は載ってないみたいで残念

**デフォルトの名無しさん** · 2018/12/16(日) 00:23:14.49

こういう話ではないと思うけど

``` {r}
leaflet::leaflet () %>% leaflet::addTiles () %>% leaflet::addMarkers (
lng = 174.768, lat = -36.852, popup = 'The birthplace of R'
) %>% leaflet::addLabelOnlyMarkers (
lng = 174.768 + 0.01, lat = -36.852, label = '此処は何処ですか？'
, labelOptions = leaflet::labelOptions (
noHide = T, textOnly = T, opacity = 1, style = list (
'font-size' = '2ex', 'color' = 'red'
)
)
) %>% leaflet::addRectangles (
lng1 = 174.768 - 0.02, lat1 = -36.852 - 0.02
, lng2 = 174.768 + 0.02, lat2 = -36.852 + 0.02
)
```

**デフォルトの名無しさん** · 2018/12/16(日) 18:28:48.34

leafletではさすがに希望には合わないかな
仕事で自分のクレカ使うのも嫌だし、こういう有料化は困る

**デフォルトの名無しさん** · 2019/02/02(土) 18:39:54.09

タイトル勝ち - envのハッシュテーブルとしての使い方
Hash Me If You Can | R-bloggers
https://www.r-bloggers.com/hash-me-if-you-can/
分散の占有率は考えたことがなかった
Principal Component Analysis (PCA) 101, using R : rstats
https://www.reddit.com/r/rstats/comments/akytig/principal_component_analysis_pca_101_using_r/

PCAネタでOjaの学習則 - ご本人による解説
Oja learning rule - Scholarpedia
http://www.scholarpedia.org/article/Oja_learning_rule
スティーフェル多様体上の勾配
The Geometry of Algorithms with Orthogonality Constraints
https://arxiv.org/abs/physics/9806030
動機と出発点は異なるが同一の微分方程式を導き出してる

対称行列の大きい方の固有ベクトルを求める問題は次のようにも書ける
given M: (n,n)対称行列, to be found X: (n,k)行列
argmax <X, M X> subject to <X, X> = 1, where <A, B> := tr (A^T, B)
subject to ...の部分がスティーフェル多様体の定義になっている
特に、k=1の場合は(n-1)次元単位球面になる
身近にある多様体の例になっていると思う

**デフォルトの名無しさん** · 2019/02/02(土) 19:02:21.25

間違え　subject to X^T X = (k,k)単位行列

**デフォルトの名無しさん** · 2019/02/22(金) 22:52:15.13

``` {Rcpp}
#include <Rcpp.h> // 「rcpp 参照渡し」で検索すると出てくる話題
// [[Rcpp::plugins(cpp14)]]
// [[Rcpp::export]]
SEXP unsafe_negate (SEXP out) { // コピーなしを確実にするために面倒だがSEXP
switch (TYPEOF (out)) { // 思いつく残りのキーワード: Rtools on Windows,
case REALSXP: { // knitr::all_rcpp_labels, RCPP_MODULE, Rcpp/dispatch.h.
Rcpp::NumericVector a (out); // RCPP_MODULE = boost::python
std::transform (a.begin (), a.end (), a.begin (), std::negate <double> ());
return out; // https://wiki.python.org/moin/boost.python/HowTo
} break; default: { // R 3.5以上で動くかわからない
throw std::runtime_error ("unsupported type");
} break; // https://purrple.cat/blog/2018/10/14/altrep-and-cpp/
} // デビアン系だけかもしれないけど、Rstudioの環境下ではC++が超お手軽
}
// [[Rcpp::export]]
Rcpp::NumericVector safe_negate (Rcpp::NumericVector inn) {
Rcpp::NumericVector out = Rcpp::no_init (inn.size ());
std::transform (inn.begin (), inn.end (), out.begin (), std::negate <double> ());
return out;
}
```

**デフォルトの名無しさん** · 2019/02/22(金) 22:55:11.29

残りを貼り忘れた

``` {r}
a = 1.0;
b = safe_negate (a);
cat (a, '->', b); # 1 -> - 1
a = 1.0;
b = unsafe_negate (a);
cat (a, '->', b); # - 1 -> - 1
```

**デフォルトの名無しさん** · 2019/02/24(日) 02:20:59.21

統計使ってオプションの自動売買したいんですけど
これ使えばできるようになりますかね？

**デフォルトの名無しさん** · 2019/02/24(日) 13:00:49.75

もう30年近く前からいっぱい使われてるよ。
それより統計知識と業務知識が先にないとあかんよ。

**デフォルトの名無しさん** · 2019/03/16(土) 22:01:31.43

最近ハマったRの文法
``` {r}
tryCatch ({
f = function (x) {
x + 1;
} + 1;
cat (f (1));
}, error = function (ex) {
cat (ex);
});
```
落とし穴というより他のプログラミング言語からの固定観念に縛られていた
jsと対比してみる
<pre id = 'dump'></pre>
``` {js}
$ (window).on ('load', function () {
try {
const f = function (x) {
x + 1;
} + 1;
$ ('#dump').html (f (1));
} catch (ex) {
$ ('#dump').html (ex);
}
});
```

**デフォルトの名無しさん** · 2019/03/21(木) 11:09:31.67

仕事で急遽Rを勉強している者です。
rvestを用いたスクレイピングについて解る方
いらっしゃいましたらご教授願います。

①下記ページを参考に rvest html_sessionで「次へ」のリンクを辿ってURLを抜き出し
http://estrellita.hatenablog.com/entry/2015/11/11/084310
②そのURLをリストに追加して
③read_html をかけようとしてるんですが
④no applicable method for 'xml_find_all' applied to an object of class "list" 　とエラーが出てしまいます。

無知で申し訳ないのですが、原因と解決法わかりますでしょうか？

**デフォルトの名無しさん** · 2019/03/21(木) 23:21:26.63

わかる人が来るまでのつなぎ - ３０年以内に来ればラッキーだと思うけど
* ['rvest::html' is deprecated, but rvest::read_html doesn't exist. Issue #191 tidyverse/rvest GitHub](https://github.com/tidyverse/rvest/issues/191)
* [rvest package | R Documentation](https://www.rdocumentation.org/packages/rvest/versions/0.3.2)
``` {r}
with (list (`%>%` = purrr::`%>%`, size = length, null = NULL), {
home = 'https://stackoverflow.com';
depth = 0;
done = list ();
todo = list ('/questions/28863775/scraping-linked-html-webpages-by-looping-the-rvestfollow-link-function');
while (0 < size (todo) & depth < 2) {
depth = depth + 1;
done = c (done, todo);
todo = purrr::reduce (.init = null, .x = todo, .f = function (out, path) {
url = paste0 (home, path);
text = xml2::read_html (url);
nodes = rvest::html_nodes (text, css = '.related a.question-hyperlink');
purrr::reduce (.init = out, .x = nodes, function (out, node) {
path = rvest::html_attr (node, 'href');
if (path %in% done) {
cat ('what a small world:', path, '\n');
out;
} else {
cat ('i am going to stalk:', rvest::html_text (node), '\n');
c (out, path);
}
});
});
}
});
```

**デフォルトの名無しさん** · 2019/03/21(木) 23:54:04.09

no applicable method for 'xml_find_all' applied to an object of class "list"

あちこちの変数をデバッグすれば？
まず、エラーの場所を特定するべき！

スクレイピングには、Ruby, Nokogiri, Selenium WebDriver が良い

**デフォルトの名無しさん** · 2019/03/22(金) 12:52:49.39

なんで、こういうツールってDOMでのスクレイピングじゃないの?
覚えるの面倒なんだ。

**デフォルトの名無しさん** · 2019/03/22(金) 22:45:06.06

バグ修正
if (path %in% done | path %in% out) {
cat ('what a small world:', path, '\n');
out;
} else {
cat ('i am going to stalk:', rvest::html_text (node), '\n');
c (out, path);
}

**デフォルトの名無しさん** · 2019/03/24(日) 22:07:13.75

楽しいライブラリ

* [Lego Mosaics Using R | Hacker News](https://news.ycombinator.com/item?id=19469142)

ライブラリの階層としては
[rgl](https://github.com/cran/rgl)/
[rayshader](https://github.com/tylermorganwall/rayshader)/
[brickr](https://github.com/ryantimpe/brickr)/
という感じかな？
Rのライブラリというと変化球勝負というイメージを持っているのだけど
rglは豪速球な気がする

**デフォルトの名無しさん** · 2019/03/26(火) 19:24:08.41

誰でも頭が良くなる、プログラムが書けるようになる方法が発見される 81137
https://you-can-program.hatenablog.jp

**デフォルトの名無しさん** · 2019/04/13(土) 23:21:10.50

* [TF-IDF in a nutshell](https://www.reddit.com/r/LanguageTechnology/comments/bb5bcr/tfidf_in_a_nutshell/)

[TF-IDF](https://en.wikipedia.org/wiki/Tf%E2%80%93idf)
が覚えられないので
[PMI](https://en.wikipedia.org/wiki/Pointwise_mutual_information)
と関連付けてみる

\newcommand{\nwd}[2]{\sharp\left\{{#1}\to{#2}\right\}}
$X$を単語の有限集合、$Y$を文書の有限集合とする
データを単語から文書への
[二部グラフ](https://en.wikipedia.org/wiki/Bipartite_graph)として見て
辺の統計を考える

| 記号 | 日本語の記号 |
|:--|:----------------------------------|
| $\nwd{x}{y}$ | 単語$x$の文書$y$での出現頻度 |
| $\nwd{x}{*} := \sum_{y\in Y}\nwd{x}{y}$ | 単語$x$の全文書での出現頻度 |
| $\nwd{*}{y} := \sum_{x\in X}\nwd{x}{y}$ | 文書$y$の長さ |
| $\nwd{*}{*} := \sum_{x\in X,\; y\in Y}\nwd{x}{y}$ | 全文書の長さ |

$$
\frac{P(x,y)}{P(x,*)P(*,y)}
:= \frac{\cfrac{\nwd{x}{y}}{\nwd{*}{*}}}
{\cfrac{\nwd{x}{*}}{\nwd{*}{*}}\cfrac{\nwd{*}{y}}{\nwd{*}{*}}}
= \underbrace{\frac{\nwd{x}{y}}{\nwd{*}{y}}}_{\approx\mathtt{TF}}
\underbrace{\frac{\nwd{*}{*}}{\nwd{x}{*}}}_{\approx\mathtt{IDF}}.
$$

**デフォルトの名無しさん** · 2019/04/28(日) 22:29:54.86

* [March: "Top 40" New CRAN Packages](https://www.r-bloggers.com/march-2019-top-40-new-cran-packages/)

[lenses](https://cran.r-project.org/web/packages/lenses/index.html)
というライブラリが入っている
ゲット/セットのペアからスタートする
[nlab](https://ncatlab.org/nlab/show/lens+%28in+computer+science%29)
ではプットと書かれているがセットと同じ
[イントロ](https://cfhammill.github.io/lenses/index.html)
も簡潔な説明だと思う
「てか、`iris $ Sepal.Length [3]`という書き方で何も困っていないし」
と思うかもしれないけどそれは正常な感覚だと思う

[この例](https://ncatlab.org/nlab/show/Grothendieck+group)
はもっとビミョーな気分になるかもしれない
「ここで小麦粉を$1/3$カップ入れます」の$1/3$という書き方は
[随伴](https://ncatlab.org/nlab/show/adjoint+functor)を表している
正の自然数からスタートした場合、Haskell風の書き方をして
直積に同値関係`Bunsu (a, b) == Bunsu (c, d) = a * d == b * c`を定義すると
単位射が`return a = Bunsu (a, 1)`で
積射が`join (Bunsu (Bunsu (a, b), Bunsu (c, d))) = Bunsu (a * d, b * c)`の
モナドになる

[この例](https://en.wikipedia.org/wiki/Simplicial_set)
はRの`list (list (0 : 0), list (0 : 1), ..., list (0 : n))`という
リストがつくる圏からスタートする
UMAPの土台に使われているらしい

* [UMAP](https://cran.r-project.org/web/packages/umap/vignettes/umap.html)
* [How UMAP Works](https://umap-learn.readthedocs.io/en/latest/how_umap_works.html)

**デフォルトの名無しさん** · 2019/05/12(日) 09:01:18.88

Rで正規表現を使って文字列から数値だけ抽出するにはどう書きますか？

**デフォルトの名無しさん** · 2019/05/12(日) 12:36:17.67

Ruby では、\d で数字を、+ で1文字以上、to_i で整数型に変換する

p ary = "1a23bc04".scan( /\d+/ ).map( &:to_i )
#=> [1, 23, 4]

ただし、これでは、負数を処理できない。
負数を処理するには、- を、? で、0か1文字

p ary = "-1a23bc-04".scan( /-?\d+/ ).map( &:to_i )
#=> [-1, 23, -4]

小数点や、e 表記は、もっと難しい。
これらは正規表現じゃなく、ライブラリを探すべき！

**デフォルトの名無しさん** · 2019/05/12(日) 19:02:18.45

どんな文字列か分からんけど単に数字を取り出すなら
stringr::str_extract_allで[[:digits:]+]を指定すりゃいいのでは？
パターンによってはreadr::parse_numberでいけるけど

**デフォルトの名無しさん** · 2019/05/13(月) 08:08:42.07

>>269
ありがとうございます。無事目的を果たすことが出来ました。

**デフォルトの名無しさん** · 2019/05/13(月) 11:34:06.46

>>270
[[:digit:]+]だったね。いつも間違える。

最終的にどうやったか書いてくれると同じように困ってる人が助かるよ