統計解析R たぶんpart3くらい (587ﾚｽ)
上下前次1-新
1(1): 2012/06/22(金)13:33 AAS
あっても良さそうなんですが,需要ないの？
過去スレのうちのひとつ
統計解析R
2chｽﾚ:tech
461: 2023/01/19(木)11:03 ID:KKbwsgWQ(1) AAS
RのS5クラス使って継承する場合、コンストラクタ中で継承するクラスのコンストラクタを呼ぶことはできないのかな？
つまりクラスAとクラスBを定義して、クラスBがクラスAを継承している場合、クラスBのコンストラクタ中でクラスAのコンストラクタを呼ぶ方法
ちょっとググったらS4クラスは出来そうなんだが、S5クラスのは解決策が見当たらんかった
別途自前でクラスAのコンストラクタ本体を定義して、クラスBのコンストラクタからそれを呼ぶしかないか？！
462: 2023/01/20(金)19:06 ID:tvsG8hGF(1) AAS
R6はどうだろう
463: 2023/02/04(土)08:53 ID:QWz+y8KS(1) AAS
dplyrを1.1.0にしてみた
.byはまあまあ便利だけどtally()には使えなかった
464: 2023/02/12(日)21:15 ID:dcVKzcVl(1) AAS
ナンノコッチャと思って調べたら
[A Heisenbug lurking in async Python | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
専門用語だったでござる。
[Heisenbug - Wikipedia](外部ﾘﾝｸ:en.wikipedia.org
465: 2023/02/14(火)18:38 ID:dFdG4M+K(1) AAS
ヒヤリハットみたいな駄洒落はどこでもあるのね
466: 2023/02/16(木)19:25 ID:76X7tEQ5(1/2) AAS
こうした言葉を思いつくのは才能なんだろうね。
Rだと簡単にハイゼンバグの例を作れた。
``` {r a_10850, dependson = ""}
uncertainly = with (new.env (), {
plus = \(x) \(y) x + y
x = 1L
plus_1 = plus (x)
x = 1000L
plus_1 (2) |> print ()
});
省18
467: 2023/02/16(木)19:27 ID:76X7tEQ5(2/2) AAS
タイトルに"遅延評価"が入っているが、遅延評価自体の問題ではなく、
Rの実装に原因がある。Haskellも遅延評価だが、`uncertainly`のような
挙動が起きれば、Haskellはその存在意義を失う。
``` ghci
import Data.Functor.Identity (Identity (..))
:{
uncertainly :: (Monad f) => f Int
uncertainly = do
x <- pure 1
plus_1 <- pure (x +)
省16
468: 2023/02/19(日)11:12 ID:CwNaXfgw(1) AAS
magrittrの例を少し改変
# fns
first <- function(x){
message("first")
invisible(x)
}
second <- function(x){
message("second")
invisible(x)
}
省4
469: 2023/02/19(日)23:57 ID:5hVGJ6gZ(1) AAS
ビックリパイプは知らなかった。"抜かりなし"だね。ただし、ビックリパイプは
ハイゼンバグとは関係ないかな。
ハイゼンバグから離れて、宇宙大戦争について書いてみたい。
`lapply`等のループ系の関数は、最終的にC関数`R_forceAndCall`を呼び出す。
* [apply.c](外部ﾘﾝｸ:github.com
Rは、"フォースと共にあらんことを"などと呑気なことを言っている場合ではなく、
フォースと共にあらねば死んでしまう。
前回のコード`uncertainly`は変数`x`の使い回しが敗因だが、ループ系の関数
では変数を使い回すしかない。そのために、ループ系の関数ではフォースが必須
になっている。98パーセントぐらいの使用例では、フォース抜きでも動作する
省2
470: 2023/02/23(木)21:58 ID:OvjjRMrU(1) AAS
do.call()ならどうだろう
471: 2023/02/26(日)00:14 ID:Qcp1BikG(1/6) AAS
上のコード`uncertainly`の挙動を"プロミスの罠"と書くことにする。
ここでの"プロミス"は、JSの"プロミス"ではなくて、Rでの"変数"の
実装方法を指す。
* [6 Functions | Advanced R](外部ﾘﾝｸ:adv-r.hadley.nz
この記事には次の一節がある。
> You cannot manipulate promises with R code. Promises are like a **quantum state:**: ...
多分、ここでの"量子状態"はハイゼンバグと同じ現象を指しているんだと思う。
観測すると、状態が変化してしまう。
省14
472: 2023/02/26(日)23:03 ID:Qcp1BikG(2/6) AAS
書けるかな？
473: 2023/02/26(日)23:09 ID:Qcp1BikG(3/6) AAS
自分の知る限り、プロミスの罠にハマるのは次のパターンに限られる。
~~~ {.r}
x = "hello"
g = f (x)
x = "world"
g ()
~~~
"関数を返す関数"`f`に、"変数"`x`を代入した時にプロミスの罠が可能性が
出てくる。Rの評価戦略はHaskellと同じ
[コールバイニード](外部ﾘﾝｸ:en.wikipedia.org
省17
474: 2023/02/26(日)23:11 ID:Qcp1BikG(4/6) AAS
コールバイニードを実現するためには、コンパイラーが内部的に次のように
書き換える必要があると思う。
``` {r a_1451, dependson = ""}
polite = with (new.env (), {
const = \(x) \(...) x
x_1 = "hello"
g = const (identity (x_1))
x_2 = "world"
g () |> print ()
});
省9
475: 2023/02/26(日)23:12 ID:Qcp1BikG(5/6) AAS
``` {python a_14213, dependson = ""}
##| cache: false
#
def none ():
id = lambda a: a
const = lambda x: lambda * y: x
def lhs ():
x = "hello"
g = const (id (x))
x = "world"
省19
476: 2023/02/26(日)23:25 ID:Qcp1BikG(6/6) AAS
この書き換え規則を適用した結果、関数`rhs`では、静的単一代入への変換は、
[ラムダ計算](外部ﾘﾝｸ:en.wikipedia.orgのアルファ変換
に置き換わる。この意味では、静的単一代入とアルファ変換は同じことになる。
実際のコードでは、分岐、ループ、副作用が入ってくるので、こんな単純は
話では済まないと思うが、極度に単純化すると、静的単一代入とアルファ変換は、
単に、方言の違いということになる。
Rの場合は、書き換え規則[rule](#rule)は等価な変換にならない可能性がある。
特に、プロミスの罠が現れた場合は、挙動が変わると思う。
余談だが、書き換え規則[rule](#rule)をモナドに拡張したものがHaskellの
[do-記法](外部ﾘﾝｸ:en.wikibooks.orgになっている。
省13
477: 2023/03/01(水)23:45 ID:X+BX12eb(1) AAS
[トランスフォーマー](外部ﾘﾝｸ:en.wikipedia.org
についての記事: オーサーヒアーかな？
* [Coinductive guide to inductive transformer heads](外部ﾘﾝｸ:news.ycombinator.com
中を読んでいないが、アテンション機構を
[ホップ代数](外部ﾘﾝｸ:en.wikipedia.orgで説明、もしくは、
ホップ代数を使ってアテンション機構の類似物を導出するという話のようだ。
この手の救世主的な話はがっかりすることが多いので、過大な期待は禁物かも。
478: 2023/03/02(木)07:01 ID:FFiKbXDz(1) AAS
g = do.call(const, list(x)) ではどうだろう
479: 2023/03/02(木)07:26 ID:rDVggB7z(1) AAS
RのS5クラス定義（setRefClass）ってsave、loadでシリアライズ、デシリアライズ出来ないのか
デシリアライズされたメンバ関数にアクセスできない
480: 2023/03/04(土)20:31 ID:0X36s/jF(1/3) AAS
## Rにおけるオブジェクト指向への取り組み
趣旨は新規フレームワークR7の紹介だが、既存のフレームワークの俯瞰図
としても使いやすいと思う。
* [What is R7? A New OOP System for R](外部ﾘﾝｸ:www.jumpingrivers.com
## オブジェクト指向とは？
次の記事を取り上げてみる。
* [Object Oriented Programming Features of Rust](外部ﾘﾝｸ:doc.rust-lang.org
内容については賛否両論だと思う。しかし、箇条書きできるということが、
オブジェクト指向が単一の概念でないことを示しているように思う。
省14
481: 2023/03/04(土)20:34 ID:0X36s/jF(2/3) AAS
似た内容をPythonで書いてみる。
``` {python a_8560, dependson = ""}
class Base:
def to_string (self): return self.as_string ()
def as_string (self): return self.to_string ()
class Derived (Base):
def to_string (self): return "hello " + type (self).__name__
Derived ().as_string ()
```
Haskellで書いてみる。
省14
482: 2023/03/04(土)20:37 ID:0X36s/jF(3/3) AAS
ワンセットの関数群をインターフェースと書く。
1. Pythonのクラスでは、
* 構造体の定義とインターフェースの定義と実装が同時に行われる。
1. HaskellのクラスやRのS3では、
* 構造体とインターフェースは別々に定義される。
* インターフェースの実装は既存の"構造体"に対して行われる。
HaskellやRの"構造体"には、通常の構造体に加えて、関数も含まれる。
構造体とインターフェースが別々に定義される点で、RのS3はHaskellのクラスに
近いように思う。
おしまい
483: 2023/03/26(日)10:03 ID:z5V/p30e(1) AAS
join_by()が便利で感動した
気象データのマージがすごくシンプルになった
484: 2023/03/26(日)22:59 ID:Dr7EPpsh(1) AAS
使用イメージがわかない。どのあたりが便利ですか？
485: 2023/03/27(月)22:57 ID:huaNY8vp(1) AAS
任意の地点・期間における測定値についての風向風速などの気象の影響をみるため近傍アメダスの気象データをマージする際に、便利になったと実感しました
具体的にはマージの際にbetween(時刻,)の条件が追加できるようになったのが大きいです
処理が重いかコードが長いかだった上のマージ処理が、join_by(地点名, between(time, start_time, end_time))で済むようになりました
486: 2023/03/30(木)00:47 ID:7ArmqfaI(1/2) AAS
[非等価ジョイン](外部ﾘﾝｸ:en.wikipedia.orgは
最近になって関数`dplyr::xxx_join`に導入されたらしい。
* Rで非等価結合 (2)
* [dplyr 1.1.0: Joins](外部ﾘﾝｸ:www.tidyverse.org
[ラッキー](動画ﾘﾝｸ[YouTube]。
関数`vctrs::vec_locate_matches`には、以前から限定された形の非等価ジョイン
のオプションがあったので、実用的な形にするのに時間がかかったのかもしれない。
リスト内包表記が使いたいのでPythonで書く。
``` {python a_32372, dependson = ""}
#| cache: false
省12
487: 2023/03/30(木)00:49 ID:7ArmqfaI(2/2) AAS
def none (): #{
a = numpy.random.choice (8, 2 * 8, replace = True)
b = a + 4
def go (pred, zipper): #{
(j, k) = slow_matches (pred, zipper) (a, b)
out = pandas.DataFrame ({"a" : a [j], "b" : b [k]})
print (out)
#}
go (lambda a, b: a == b, itertools.product)
go (lambda a, b: a == b, lambda a, b: itertools.islice (zip (a, b), 1))
省10
488: 2023/03/30(木)18:24 ID:Yqa5DxWU(1) AAS
一応Rスレだしコード一式を載せるならRにしてほしいな
489: 2023/03/31(金)01:01 ID:3FvfgxM7(1) AAS
正直、書き方を思い出すためだけに、話題にかこつけて、わざと異なる
プログラミング言語で書くこともある。しかし、今回は純粋に、処理の道筋を
単刀直入な形で表せるプログラミング言語を選択した。
[線形回帰](外部ﾘﾝｸ:en.wikipedia.orgを例にとる。
Rで線形回帰はとても書きやすい。しかし、線型回帰をRのイディオムとして
覚えるよりも、線形代数として覚えた方が長い目で見た時に利益になると思う。
同じように、等価ジョインは
[プルバック](外部ﾘﾝｸ:en.wikipedia.org
に対応する。等価ジョインを、Rのイディオムとして覚えるよりも、
単純な`for`ループの形で覚えた方が長い目で見た時に利益になると思う。
省1
490: 2023/04/15(土)20:49 ID:Mf5VNEJ8(1) AAS
プルバックの話の続きをアップした。
* [JSFiddle](外部ﾘﾝｸ:jsfiddle.net
プルバックのカリー化の話になっている。プルバックの話はこれでお終い。
洗練された話は次の記事にある。
* [locally cartesian closed category in nLab](外部ﾘﾝｸ:ncatlab.org
表計算の立場で見ると、"何が問題なのかわからない"問題かもしれない。
ある規則を満たす表でサマリーするとプルバックのカリー化になる。
491: 2023/04/24(月)19:43 ID:Jj998bTN(1) AAS
4.3.0入れてみた
パイプがまたひとつ便利になっていい感じ
Windowsでパスが長くても通るようになったのも嬉しい
492: 2023/04/26(水)12:00 ID:GTOaZ7cP(1/2) AAS
どなたか教えてくだされ
対照実験の医療統計をしたいと思って
傾向スコアマッチングにしようと思ってます
肺がん手術を実施した人たち200症例のうち10人だけ
手術後に合併症を起こしていました
この場合
「合併症を起こしたグループ」・・・10人
「合併症を起こさなかったグループ」・・・190人になり
患者の背景で傾向スコアマッチングをキャリパー幅0.2でペア分けすると
省3
493(1): 2023/04/26(水)18:48 ID:ZF8W3LJH(1) AAS
共変量のバイアスがより少ないサンプルが得られますが、サイズが8人の二群比較なのは変わらないんじゃないかと思います
494(1): 2023/04/26(水)23:58 ID:GTOaZ7cP(2/2) AAS
>>493
つまりサイズが8人ずつの比較では
優位差を言うには弱いと思いますか？
495: 2023/04/27(木)08:26 ID:vIGDv5il(1) AAS
>>494
検定はサンプルサイズ込みの判定なので8人であっても二群に有意な差があったと言って差し支えありません
有意水準が甘ければ弱いとは感じますが
一方で、検定は標本のバイアスが除去できていることを保証してくれません
496: 2023/05/13(土)14:17 ID:LcdFqRfa(1) AAS
スレッドの質をガクッと下げるが、Rの非標準評価で遊んでみた。
* [JSFiddle](外部ﾘﾝｸ:jsfiddle.net
497: 2023/05/16(火)14:12 ID:T6FZDpAX(1) AAS
4.3.0のWindows版はdir()にバグあるね
リストがまともに取得できない
4.2.3に戻した
498: 2023/05/19(金)19:12 ID:LBQCZv0F(1) AAS
EZR使ってる人いる？
499: 2023/05/19(金)23:40 ID:AtidLON9(1) AAS
EZRまだ更新続いてて驚いた
すごいな
500: 2023/05/21(日)15:55 ID:Ifo/fxoU(1) AAS
10年以上前からあるスレでようやく500ゲット！
501: 2023/05/21(日)20:10 ID:F9CBOLLZ(1) AAS
750に達するのにもう10年を要するだろう
502: 2023/05/22(月)00:40 ID:Cnh2SEE8(1) AAS
tibbleを利用してデータ探索しているんだけど、2つの変数で条件つけて散布図つくるにはどんな手順踏んだらいいのか…
503: 2023/05/22(月)06:06 ID:wj8Gj0ut(1/2) AAS
変数1をプロットの色にして変数2をサイズにするとかじゃだめなの？
504: 2023/05/22(月)09:44 ID:/JtWuqEq(1) AAS
それも一つの方法ですね。ありがとう
ただ、目的以外の条件をプロットしたくないのです。
tibbleの構成としては、観測が4時点、変数がたとえば6個(うち1つは個体識別番号で、変数としてはA-Eの5種類あるとしましょう)のような感じで、
散布図のx軸に「観測時点1における全個体の変数A」
y軸に「観測時点4における全個体の変数E」
をプロットしたいのです。
個体識別番号をキューにして、異なる時点の異なる種類の変数を散布図にするとでもいいましょうか…
505(1): 2023/05/22(月)11:53 ID:wj8Gj0ut(2/2) AAS
その要件だと以下のような感じでいけると思います
library(tidyverse)
df <- tibble(
time = rep(1:4, each = 5),
id = rep(1:5, 4),
a = rnorm(20),
b = rnorm(20),
c = rnorm(20),
d = rnorm(20),
e = rnorm(20)
省4
506: 2023/05/22(月)14:40 ID:7X6vWo+R(1) AAS
ありがとうございます！
それをヒントに試してみます。
507(1): 2023/05/24(水)23:00 ID:1oF+cvFs(1) AAS
>>505
これ、もしかして途中ではないですか？
508(2): 2023/05/25(木)05:49 ID:NLUeSyH/(1) AAS
>>507
要件の理解が違っていなければ散布図作成まで行きますよ
time1でのaの値とtime4でのeの値をid毎にプロットしています
509: 2023/05/25(木)07:49 ID:DfoFl+x4(1) AAS
>>508
失礼しました
こちらのビューワの問題で、ブラウザでみたらスクリプトの最後まで表示されました💦
510: 2023/05/26(金)13:34 ID:gPviguuO(1) AAS
>>508
教えてもらったスクリプトでドンピシャでした、ありがとうございます。
511: 2023/05/27(土)00:00 ID:qqVzW/zi(1) AAS
どなたか統計素人の私に教えてください
医療でA群100名とB群100名をいろんな観点から
比較するとします
背景を揃えるために傾向スコア分析で
キャリパー0.2で取って30ペア(合計60人)を抽出しました
例えば
①アミラーゼの値とかを統計にかけたいときは、その60人のアミラーゼのデータをまずは正規分布かどうかをShapiro-wilk検定で測って、0.05を上回っていたら「正規分布」と見なして、今度はこれら60人のアミラーゼのデータが等分散かどうかを確認するためにf検定を行なって等分散だったらt検定(Student's t-test)、違ったらWelch's t-testで優位差があるかどうかを調べる。
もしShapiro-wilk検定が0.05未満で非正規分布だったときはMann-whitney's U検定で優位差があるかどうかを調べる。
省6
512(1): 2023/05/27(土)19:37 ID:H+hysMph(1) AAS
素人だと思ううちは傾向スコアには手を出さないほうがいいような気がする
513: 2023/05/27(土)22:19 ID:Z7qvu7kx(1) AAS
>>512
上記の場合はどういうかんじでやるのが
最適解になるのかだけでも教えてもらえないでしょうか？
514: 2023/06/02(金)21:25 ID:AkfCPeTx(1) AAS
writexlパッケージで保存すると時刻のタイムゾーンは反映されないのね
515: 2023/06/13(火)17:38 ID:J5T/voP2(1) AAS
ggraggedパッケージがなかなか便利
facet_grid()を詰めて並べられる
516: 2023/06/18(日)13:00 ID:13gXw1hf(1) AAS
R4.3.1にした
dir()関連が直ってるようでひと安心
標準パイプの新機能もようやく使える
517: 2023/06/18(日)16:20 ID:ECqif9s2(1) AAS
しっぽりしてえな
518(1): 2023/06/19(月)07:17 ID:lFY/Y6gc(1) AAS
ggplot2のカラーパレットのデフォルトが変わったりしました？
棒グラフの色の割り当てが変わった気がする…
519: 2023/06/19(月)19:42 ID:CYQnxNgF(1) AAS
>>518
変わってないはず
ただggplot2のデフォルトはグループ数に応じて関数で色を割り当ててるのでグループ数が異なると違うパレットのように感じるとは思います
520(1): 2023/06/19(月)20:01 ID:b0yPHKbc(1) AAS
株価データを取得したいがためにRを弄り始めたけど中々難しいね
ググってコピペしたコードを一つ一つパッケージと関数の挙動を調べているけど
やりたいこと全てができるようになるまでは遠そうだ
521: 2023/06/19(月)22:06 ID:J5tIl656(1) AAS
>>520
データをどこから持ってくるかによって随分難易度は違う。
ファイルとしてダウンロードできるところもあるし、
スクレイピングが必要なところもある。場合によってはせれにうむとか。
522: 2023/06/22(木)23:31 ID:eZc4H/Dz(1) AAS
企業が求めるエンジニア言語ランキング 2位「R」
外部ﾘﾝｸ[html]:www.itmedia.co.jp
3位以下を大きく引き離しててすごいけど…
Rユーザの少なさを反映しただけな気もする
523: 2023/06/25(日)11:51 ID:JqE+uKYL(1) AAS
Rの求人なんてみたことない…
524: 2023/06/25(日)15:23 ID:XZytU2Vy(1) AAS
>求人数のトップ3は1位「C言語」（39万2559件）、2位「R言語」（30万3759件）、3位「Go」（8万3922件）だった。
>このランキングは、SalesNow DBに1月1日から5月31日までに載った114万5294件の有効求人から算出した。
5ヶ月で30万件もあったらしい
Pythonが4万件弱なところを見るとデータ解析用言語を列記したわけじゃなく純粋にRの需要があったみたいだ
不思議
525: 2023/06/26(月)18:52 ID:wRy8Kzvb(1) AAS
結局企業で欲しいのは統計処理だから、統計処理ならRだろうっていう短絡的な
求人が多かったってことじゃない？
526: 2023/06/28(水)18:59 ID:z2KrS2jG(1) AAS
いくつか求人サイト見てみたけどRのみなんて求人はほぼ皆無でほとんどがPythonとの併記だった
これ集計ミスかもね
同じ会社の別の期間のランキング見るとRは10位以内に入ってないし
527: 2023/06/28(水)19:11 ID:pwcxWBkk(1) AAS
金が有ーる(R)といい
なんちって
528: 2023/06/29(木)08:36 ID:szIurvDh(1) AAS
業∞！！！！
字∞！！！！！
苗∞！！！！！！
米∞！！！！！！！
畑∞！！！！！！！！
灰∞！！！！！！！！！
糞∞！！！！！！！！！！
机∞！！！！！！！！！！！
529(1): 2023/07/23(日)03:41 ID:SQtd3M8i(1) AAS
結構何年か前にRは結構古臭いし、言語仕様もわかりづらいし衰退するのでは、と書いているのを
読んだことがあるけど、逆に今だと他の商用統計パッケージが廃れる勢いなのにRって立場が盤石な
気がする
なんで？
530: 2023/07/24(月)10:46 ID:uQbqJHBb(1) AAS
AIブームでPythonに全部もってかれると思ったけどそうでもなかったね
先発優位で統計屋や大学教授が昔からRユーザーってのもあるのかも？
ユーザー視点からだとtidyverseの存在が大きいと思う
505のように地味に手間がかかりそうな処理が関数ひとつでひょいと処理できるとかすごい
531: 2023/07/24(月)12:50 ID:IaCKgEWP(1) AAS
ど〜も〜！ご存じオナニーマン チン太郎でーす！
532: 2023/07/25(火)14:51 ID:iTChcdyR(1) AAS
>>529
今でもC/C++が盤石なのと同じ理由な気がする
533(1): 2023/07/25(火)22:12 ID:KSU0u1F8(1) AAS
とくに大学の研究室とかは、同じような機能のプログラムを別言語に
コンバートする作業は雑用でしかないので誰も積極的にはやらない。
RでできるのならこのままRでいいじゃん、もしくはフォートランでもｗ
みたいな感じ
やったところでなにか新しい知見が得られるわけもなし。
534: 2023/07/26(水)14:36 ID:/vrstjVu(1) AAS
昔の医学分野の統計処理はSASがデフォだったけど
COVID-19で疫学の人（西浦先生とか）がRガンガン使ってるの見て疫学以外でもRが広がってる
学会行くとggplotで書いたグラフいっぱい見る
535: 2023/07/30(日)15:01 ID:/zgOfttX(1) AAS
ggplot2のデフォルトテーマは微妙なので変えているけどいつも悩む
ggthemesをもとにそこから装飾することが多い
536: 2023/07/31(月)13:12 ID:sgBBFIN2(1) AAS
>>533
むしろ結果が違ってたらωωω
537: 2023/08/08(火)23:42 ID:D9D8OCal(1) AAS
Rしばらく使ってないから忘れそう
538: 2023/08/23(水)14:24 ID:LaIh/7Co(1) AAS
ついにExcelでPython使えるようになるらしい
マイクロソフトはRのベンダー買収してMRANとかやってたのに時代の流れかな
539: 2023/08/24(木)07:11 ID:UfAeCzV0(1) AAS
Pythonでexcel弄るのは
15年以上前から出来てるのを今更
excel上のpythonで実行する方が遅い可能性まである
540: 2023/08/25(金)13:03 ID:8Z9WZ/AE(1) AAS
PythonからExcelをいじることとExcelでPythonが動くことは全く違う
541: 2023/08/25(金)13:58 ID:5+gJach+(1) AAS
くりゃうどω
542: 2023/08/25(金)14:03 ID:HPTKrBcu(1) AAS
pythonは何故かPCでJupyter notebookが不具合で使えなくてな
導入時点でコケるとその言語をやる気失くす…Rはその点、導入が簡単で良かった
543: 2023/08/27(日)01:44 ID:VsjIkj4j(1) AAS
Excelも最近はTidydata志向で、パワークエリやパワーピボットで再現性のある加工や分析が
できるようになってきているけど、使い始めはなかなかいいかな、と思ったがちょっと凝ると
やっぱり複雑化してメンテできないVBAみたいになっちゃうなｗ
ま、RやPythonだからメンテ性がいいというわけでもない気はするけどｗ
544: 2023/08/27(日)20:22 ID:RTNKd3te(1) AAS
Excelマクロ禁止してきた企業はPython機能も禁止しそう
Rを社員が使うケースなんてのはもうレア度SSRなんだろうね
545: 2023/08/28(月)13:50 ID:62AYCckx(1) AAS
Rで学ぶなんとか分析とかいう本いっぱい出ているし、結構個人でやる気のあるやつは覚えるんじゃね？
R人口増えている印象。
546: 2023/08/28(月)13:57 ID:3YABDYXJ(1) AAS
昨日の統数研のワークショップ見て、大学や企業によってはライブラリのインストールに制限あったりでデフォルトで解析できるのがRが多言語に置き換わらない理由の一つなのかもなと思った
今後も「選択と集中」続けるどころか懸賞金型推進するようだから商用に金払うのも辛くなるとこも出てきそうね
547: 2023/08/29(火)18:08 ID:vrxlGHM5(1) AAS
数が増えてることと数が少ないこととは違う
これ混同する人は結構多い
548: 2023/08/30(水)15:28 ID:zHni/Vbc(1) AAS
当社比200%
元が少なかっただけやろ
みたいな
549: 2023/09/01(金)21:33 ID:T0/12uyc(1) AAS
自分は今の部署にいる間だけは統計解析するけど異動したらもう解析なんてしないからRのモチベは低いな
人事の巡り合わせでいつかまた今の部署に戻るかもしれないけど
550: 2023/09/02(土)10:50 ID:aOHtbbKy(1) AAS
隙あらばどうでもいい自分語り
551: 2023/09/02(土)14:49 ID:DQO4DB1k(1/2) AAS
過疎スレでその煽りは馬鹿じゃねえの
552: 2023/09/02(土)15:24 ID:mCX3wjBN(1) AAS
本人だろ age てるのは
553: 2023/09/02(土)15:27 ID:DQO4DB1k(2/2) AAS
ageを最近知ったみたいに強調してる馬鹿ww
554: 2023/11/26(日)13:47 ID:TX3Gb+56(1) AAS
[プッシュアウト](外部ﾘﾝｸ:en.wikipedia.org
について書いてみた。
[JSFiddle](外部ﾘﾝｸ:jsfiddle.net
集合の圏でのプッシュアウトは
[等価ジョイン](外部ﾘﾝｸ:en.wikipedia.org
の圏論的双対にあたる。
555: 2023/11/26(日)19:19 ID:fOcIWwuE(1) AAS
おーまんちょー！
556: 2023/12/02(土)19:24 ID:UAUVJeJE(1) AAS
前回のおーまんちょーでは一部のコードで途中から行番号が表示されるようになっていた。
多分、JSFiddleのエディターで`tidy`というボタンを押してしまったからだと思う。
実は、`prettydoc::html_pretty`で行番号が表示されないのを直したいと思っていた。
しかし、どこから手を付けたら良いかわからなかった。それが、ハプニングのおかげで
手の入れ場所がわかった。多分、HTMLとしては不整合な状況になっているが、
意図した通りに表示されるようになった。
[ありがとう、おーまんちょー](外部ﾘﾝｸ:jsfiddle.net
557: 2023/12/15(金)00:55 ID:KmwBi9fF(1) AAS
group_by()はmap()と組み合わせて使えないんだね
558: 2023/12/17(日)12:11 ID:j39jK9SN(1) AAS
一般的な意味では、そんなことないと思うよ。
状況によってはgroup_byとmapを組み合わせると上手く行くことがあると思う。
ぐるぐる図のくだりでグダグダなところがあったので、
[修正した](外部ﾘﾝｸ:jsfiddle.net。
アップした後にタイポを見つけたが、覆水盆に返らず。
559: 2023/12/17(日)23:02 ID:/RmbQ2uq(1) AAS
多変数の時系列データをクラスタリングするpackageをご存知であれば教えていただけますか…？
560: 2023/12/18(月)00:40 ID:+7hNKB64(1) AAS
知らない。
人がいれば、意味のない返事はしない方が良いかもしれないが、過疎だからね。
561: 2023/12/18(月)23:51 ID:t5S7JJTy(1) AAS
あなたは地球上で何日過ごしてますか？
という問いに答えるRのスクリプト
干支と誕生日（月日）から今日が生まれて何日めか計算
eto=("ひつじ")
tag="10-17"
Eto=c("子(ね＝ねずみ)","丑(うし)","寅(とら)","卯(う＝うさぎ)","辰(たつ＝龍(りゅう))","巳(み＝へび)","午(うま)","未(ひつじ)","申(さる)","酉(とり＝にわとり)","戌(いぬ)","亥(い＝いのしし)")
# Eto
i=grep(eto,Eto)
this_year=as.numeric(strsplit(as.character(Sys.Date()),'-')[[1]][1])
j=this_year%%12 - 3
省7
562: 2023/12/19(火)00:23 ID:jHCQnODd(1/2) AAS
[Eto'o](外部ﾘﾝｸ:en.wikipedia.orgはもう引退しました。
563(1): 2023/12/19(火)00:58 ID:jHCQnODd(2/2) AAS
マジで聞いているんなら、テストツールの使い方を覚えた方が建設的だと思う。
個人的には、Haskellからの流れで
[quickcheck](外部ﾘﾝｸ:github.com
しか経験したことがないが、他にもあると思う。
コードが合っているかどうかわからないが、パット見問題なさそうに見える。
問題があったとしても、デバッグで十分対応できるように見える。
[リーナスの法則](外部ﾘﾝｸ:en.wikipedia.org
あたなが十分な暇人ならばバグはいずれ塞がれる。
564: 2023/12/19(火)08:05 ID:5FW9R46b(1) AAS
>>563
サンクス
565: 2023/12/20(水)23:21 ID:H48HUSxu(1) AAS
ChatGPTはなんていい奴なんだろう。
クソ野郎: I'd like to test the following code:
``` r
拝借
```
* The code is written in the programming language R.
* "子", "丑", ..., are the Chinese / Japanese zodiac called Eto in Japanese.
This code have to calculate the age for given birth date in Eto, e.g.
`eto = "ひつじ", tag = date = "10-17"`. Are there any advice to test this code?
省14
566(1): 2023/12/21(木)22:26 ID:Jt1+wTI1(1) AAS
出禁になった。書けるかな？
てなわけで、ちゃっとじーぴーてーの残りの部分は割愛する。
意地悪クソ野郎と親切じーぴーてーにも共通することろがあって、
最後の締めはやっぱりリーナスの法則になる。言い方が全然違うけど。
親切じーぴーてーが例に挙げたテストコードがテストの基本形になると思う。
くいっくちぇっくのようなテストツールでは、テストケースをランダムに発生させるが、
そうしたテストツールが使えない場合も多々ある。テストツールが使えない場合は、
親切じーぴーてーが示したようなベタな方法で行くしかない。
567: 2023/12/21(木)22:41 ID:6SaY/e5W(1) AAS
>>566
書き込めるかどうか判らぬままに
黙々とその低能感溢れる文章を書いている姿を想像すると泣けてくる
書き込めてよかったな
568: 2023/12/22(金)00:02 ID://WOG+3q(1) AAS
親切じーぴーてーとの会話を見返すと、関数を細分化して、(ナウ、干支、月日)
という三つ組を引数とする関数を作るのもありかな？という気がしてきた。
このように関数を細分化すると、テストケースが作りやすくなる。
569: 2023/12/22(金)15:14 ID:BVAe/0CV(1) AAS
かけるかな
570: 2024/01/04(木)21:13 ID:r+ggAGcN(1) AAS
クラスタリングとか次元削減に関する教科書的な本ってありますか？
571: 2024/01/05(金)00:08 ID:Cbq0Hg0X(1) AAS
かけるかな
統計の人でも機械学習の人でもないので、返事を書くのは心苦しいものがあるが、
過疎なのでしょうがない。統計処理と機械学習の区別がつかないぐらいの素人だよ。
ありきたりだけど：
* [朱鷺の杜Wiki](外部ﾘﾝｸ[php]:ibisforest.org
グーグルとアヒルで検索してみたけど、この検索については英語の方が情報が取りやすいと思った。
英語でも日本語でも本屋の広告が多くなるが、英語の方が俺が推薦図書はこれだ！
みたいなのが多いような気がした。検索で見つけた本だけど非売品と商品の両方がある。
* [Modern Dimension Reduction](外部ﾘﾝｸ:arxiv.org
* [Modern Dimension Reduction](外部ﾘﾝｸ:www.cambridge.org
省6
572: 2024/05/02(木)14:41 ID:9naKJa9o(1) AAS
統計解析向けのプログラミング言語「R」の実装に欠陥、任意コード実行のおそれ
外部ﾘﾝｸ[html]:forest.watch.impress.co.jp
573: 2024/07/06(土)23:18 ID:dmyvQpwi(1) AAS
書けるかな？
[Rのリサイクル規則](外部ﾘﾝｸ:jsfiddle.net
中にも書いたが、polyfill.ioはシリアスな案件なので、外だししておく。
[polyfill.io問題のQuarto・R Markdownへの影響](外部ﾘﾝｸ:zenn.dev
574: 2024/07/17(水)11:23 ID:jPCMTQ1M(1) AAS
ubuntu24.04対応のRstudio来ねーな・・・
575: 2024/08/20(火)00:42 ID:c8o+ygTe(1) AAS
こねーなｗ
難航してるのかな。
576: 2024/09/10(火)01:06 ID:SFqXDqGn(1) AAS
参考になるかな？
* [RStudio fails to run on Ubuntu 24.04 - Issue #14336 - rstudio/rstudio - GitHub](外部ﾘﾝｸ:github.com
他人事なら正論が吐ける。Dockerで試すのが吉。
実際には面倒なのでなかなかDockerで試してからとはいかないが、
他人事なら何でも言える。
577: 2024/09/24(火)00:03 ID:B2G5MWgw(1) AAS
Hacker-Newsネタ
1. [React for R | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
2. [Why wordfreq will not be updated | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
3. [Data Science Agent and Code Transformation | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
4. [Terence Tao on O1 | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
5. [Brainfuck Enterprise Solutions | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
6. [What happened to the Japanese PC platforms? | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
一年に一本くらいRについての投稿があるような気がする。
2番目は、個人的には、考えさせられる内容になっている。
3番目と4番目は夢のある話で、5番目は楽しい話になっている。
省2
578: 2024/09/30(月)19:17 ID:87Zg0Lej(1/2) AAS
Rstudio ubuntu24版ｷﾀ
579: 2024/09/30(月)19:32 ID:87Zg0Lej(2/2) AAS
とりあえず騒がれてた起動の不具合は解消されてるみたいだね
580: 2024/10/10(木)01:08 ID:Q3cGnE1w(1) AAS
* [Nobel Prize in Physics awarded to John Hopfield and Geoffrey Hinton | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
ノーベル物理学賞が機械学習系のお二方に決まった。
業績から言って文句なしのお二方だが、それが物理学賞？という疑問が湧く。
話に聞いていた物理学賞の選考基準は「物理現象の新発見」だった。
これが理由でarxivの作者が物理学賞を受賞しないと思っていた。
理数系の研究におけるarxivの貢献は過去10年分の物理学賞と化学賞と
フィールズ賞とチューリング賞を全て合わせたより大きいかもしれない。
個人的には次の2点が引っかかる。
* 何故、このタイミングで？
* お二方とも高名である。
省7
581: 2024/10/10(木)12:14 ID:bwO1aVtU(1) AAS
数学症も無いし
CS賞部門も増設すべきだよな
582: 2024/10/11(金)14:00 ID:W/ZPwHBO(1) AAS
放送大学の秋山先生の講座ってどうなの？
583: 01/21(火)21:45 ID:Q1AkQbXB(1) AAS
Rには関係ないが、資料として置いておく。
* [PCAST Reports (2020-2025)](外部ﾘﾝｸ:www.erichorvitz.com
* [Documents & Reports | PCAST | The White House](外部ﾘﾝｸ:web.archive.org
中身は全く読んでいないが、次の討論会を見る限り、政治色はないと思う。
* [PCAST: Discussion of Artificial Intelligence (AI) Enabling Science and AI Impacts on Society - YouTube](動画ﾘﾝｸ[YouTube]
584: 02/04(火)00:46 ID:V6kAFtDi(1) AAS
この話題は政治色抜きでは語れないと思うが、それを抜きにして、
アメリカ政府の発行したデータを使っている人がいれば、注意が必要になると思う。
* [Backing Up US Federal Databases | The n-Category Cafe](外部ﾘﾝｸ:golem.ph.utexas.edu
* [CDC data are disappearing | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
585: 02/15(土)21:36 ID:jpL9Lr/5(1) AAS
次の動画の音声だけ聞いていたところ、只者でない気配を感じた。訳者の腕もあると思う。
* [GPT解説2 アテンションの仕組み (Attention, Transformer) | Chapter6, 深層学習 - YouTube](動画ﾘﾝｸ[YouTube]
そこで、動画を見ると、丁寧なアニメが動いているじゃありませんか。
手練が作っていると思うが、それでも相当制作時間がかかっていると思う。
おまけだが、声質やゆったりさ加減がiPadの広島弁の解説に似ているように感じた。
あんたあ中国新聞読むじゃろ。あんたが数学に合わせる必要はないんよ。
数学の方があんたのええように合わせるけえ。だといいね。
外のホームページでは、カテゴリ分けされているので、外の方が番組を選びやすいかもしれない。
省1
586: 03/10(月)01:40 ID:trsagf4Z(1) AAS
0. [I'll be back](動画ﾘﾝｸ[YouTube]
1. [Public health data disappeared. RestoredCDC.org is bringing it back | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
CDC was back.
続報
2. [Introduction to Stochastic Calculus | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
伊藤さん is ubiquitously back.
3. [Differentiable Logic Cellular Automata | Hacker News](外部ﾘﾝｸ:news.ycombinator.com
車輪の大発明 will be back again and again.
[量子解析学](外部ﾘﾝｸ:en.wikipedia.org
の言い出しっぺはオイラーやヤコビらしいが、これは
省7
587: 05/22(木)01:15 ID:+HZBQ7+B(1) AAS
RStudioを更新したらreticulateが動作しなくなった。
* [Release reticulate 1.41.0 - rstudio/reticulate - GitHub](外部ﾘﾝｸ:github.com
今までどうやって動作させてたんだ？記憶にございません。
というわけで、一からやり直した。ワイのPython環境:
* システムにPython via sudo apt
* ユーザーにもPythonパッケージ via pip3
ワイの処方箋:
* in .Rprofile: Sys.setenv (RETICULATE_PYTHON = "/usr/bin/python3")
Pythonのvenvを使う場合は、Sys.setenv (RETICULATE_PYTHON_ENV = path-to-venv)でいける。
似た話題で、
省1
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 1.203s*