統計解析R たぶんpart3くらい (587ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞
ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡
212: 2018/05/17(木)10:20:38.86 ID:X8D7/hwr(1) AAS
Rをかじっても
やはり、考え方がわかってないとだめだね
ということで、
「入門はじめての分散分析と多重比較」を読んでる
266: 2019/04/28(日)22:29:54.86 ID:Ce5AO+sY(1) AAS
* [March: "Top 40" New CRAN Packages](外部ﾘﾝｸ:www.r-bloggers.com
[lenses](外部ﾘﾝｸ:cran.r-project.org
というライブラリが入っている
ゲット/セットのペアからスタートする
[nlab](外部ﾘﾝｸ:ncatlab.org
ではプットと書かれているがセットと同じ
[イントロ](外部ﾘﾝｸ:cfhammill.github.io
も簡潔な説明だと思う
「てか、`iris $ Sepal.Length [3]`という書き方で何も困っていないし」
と思うかもしれないけどそれは正常な感覚だと思う
[この例](外部ﾘﾝｸ:ncatlab.org
はもっとビミョーな気分になるかもしれない
「ここで小麦粉を$1/3$カップ入れます」の$1/3$という書き方は
[随伴](外部ﾘﾝｸ:ncatlab.orgを表している
正の自然数からスタートした場合、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)`の
モナドになる
[この例](外部ﾘﾝｸ:en.wikipedia.org
はRの`list (list (0 : 0), list (0 : 1), ..., list (0 : n))`という
リストがつくる圏からスタートする
UMAPの土台に使われているらしい
* [UMAP](外部ﾘﾝｸ:cran.r-project.org
* [How UMAP Works](外部ﾘﾝｸ:umap-learn.readthedocs.io
300: 2019/09/13(金)20:18:38.86 ID:wTJZs2gr(3/4) AAS
Rの関数`ecdf`は次のように定義されている。 $$ \begin{split}
\ecdf(x) &:= \ecdf_\xi(x) := \frac{1}{n} \sum_{j=1}^n \isa{x \ge \xi_j}, \\
\isa{\operatorname{expr}} &:= \begin{cases}
1, & \text{ iff } \operatorname{expr} = \mathtt{true}, \\
0, & \text{ otherwise}. \end{cases} \end{split} \eqtag{eq:ecdf}
$$ ここで、$\xi_1<\cdots<\xi_n\in\bbR$を観測された値とする。
これを次のように積分したものを$\ebdf$と書くことにする。
$$ \ebdf(x) := \ebdf_\xi(x) := \int_{y=-\infty}^x \ecdf(x)
= \frac{1}{n} \sum_{j=1}^n \relu(x - \xi_j). \eqtag{eq:ebdf} $$
`ebdf_1`はこの式を、`ebdf_0, ebdf_huge, ebdf_2`は式変形したものを実装している。
``` {r}
with (dog_data, {
plot (ecdf (xi)); x = seq (min (xi) - 1, max (xi) + 1, len = 1e+3);
bdf = ebdf_2 (xi); plot (x, bdf (x), type = 'l', main = 'ebdf');
points (xi, bdf (xi)); lines (x, prelu (x - mean (xi)), lty = lty_dash); });
```
コードを見ると、`ebdf_2`は、
[ReLU](外部ﾘﾝｸ:en.wikipedia.orgを
[活性化関数](外部ﾘﾝｸ:en.wikipedia.orgとする
[ResNet](外部ﾘﾝｸ:en.wikipedia.orgに
なっていることがわかる。さらに、`ebdf_2`を微分すると、ReLUとシグモイドが
混在したResNetになる。
``` {r}
ecdf_2 = function (xi) { n = length (xi); xi = sort (xi);
qi = cumsum (xi); function (x) {
purrr::reduce (.init = cbind (n * x - qi [n], n), .x = n : 1, .f = function (out, j) { f = out [, 1]; d_f = out [, 2]; g = (f + qi [j]) / j - xi [j];
cbind (f - plure (g), d_f - (g < 0)); }) / n; }; };
```
続く
337: 2021/03/07(日)00:01:20.86 ID:9nPgQwYX(1) AAS
マイクロソフトがMicrosoft365Rというパッケージを公開してる
365(Office365)をRから操作できるらしい
371: 2021/05/09(日)15:07:00.86 ID:YBysisGe(1) AAS
>>370
デフォルト環境で0.29が0.28と表示されるわけないわ
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.024s