統計解析R たぶんpart3くらい (587レス)
統計解析R たぶんpart3くらい http://mevius.5ch.net/test/read.cgi/tech/1340339592/
上
下
前
次
1-
新
通常表示
512バイト分割
レス栞
1: デフォルトの名無しさん [sage] 2012/06/22(金) 13:33:12.40 あっても良さそうなんですが,需要ないの? 過去スレのうちのひとつ 統計解析R http://toro.2ch.net/test/read.cgi/tech/1285597198/ http://mevius.5ch.net/test/read.cgi/tech/1340339592/1
461: デフォルトの名無しさん [] 2023/01/19(木) 11:03:41.55 ID:KKbwsgWQ RのS5クラス使って継承する場合、コンストラクタ中で継承するクラスのコンストラクタを呼ぶことはできないのかな? つまりクラスAとクラスBを定義して、クラスBがクラスAを継承している場合、クラスBのコンストラクタ中でクラスAのコンストラクタを呼ぶ方法 ちょっとググったらS4クラスは出来そうなんだが、S5クラスのは解決策が見当たらんかった 別途自前でクラスAのコンストラクタ本体を定義して、クラスBのコンストラクタからそれを呼ぶしかないか?! http://mevius.5ch.net/test/read.cgi/tech/1340339592/461
462: デフォルトの名無しさん [] 2023/01/20(金) 19:06:37.10 ID:tvsG8hGF R6はどうだろう http://mevius.5ch.net/test/read.cgi/tech/1340339592/462
463: デフォルトの名無しさん [] 2023/02/04(土) 08:53:58.96 ID:QWz+y8KS dplyrを1.1.0にしてみた .byはまあまあ便利だけどtally()には使えなかった http://mevius.5ch.net/test/read.cgi/tech/1340339592/463
464: デフォルトの名無しさん [] 2023/02/12(日) 21:15:52.26 ID:dcVKzcVl ナンノコッチャと思って調べたら [A Heisenbug lurking in async Python | Hacker News](https://news.ycombinator.com/item?id=34754276) 専門用語だったでござる。 [Heisenbug - Wikipedia](https://en.wikipedia.org/wiki/Heisenbug) http://mevius.5ch.net/test/read.cgi/tech/1340339592/464
465: デフォルトの名無しさん [] 2023/02/14(火) 18:38:12.15 ID:dFdG4M+K ヒヤリハットみたいな駄洒落はどこでもあるのね http://mevius.5ch.net/test/read.cgi/tech/1340339592/465
466: デフォルトの名無しさん [] 2023/02/16(木) 19:25:35.83 ID:76X7tEQ5 こうした言葉を思いつくのは才能なんだろうね。 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 () }); ``` アレレ?ということでデバッグ文を入れてみる。 ``` {r a_29042, dependson = ""} principle = with (new.env (), { plus_debug = \(x) { cat ("Schr?dinger says ", x, "\n", sep = "") \(y) x + y } x = 1L plus_1 = plus_debug (x) x = 1000L plus_1 (2) |> print () }); ``` デバッグ文の有無で挙動が変わる。`uncertainly`の挙動は、バグではなく、 言語仕様だと思った方が良いと思う。 `uncertainly`の挙動を防ぐには関数`force`を使えば良い。 * [How to not fall into R's 'lazy evaluation trap](https://stackoverflow.com/questions/29084193) http://mevius.5ch.net/test/read.cgi/tech/1340339592/466
467: デフォルトの名無しさん [] 2023/02/16(木) 19:27:35.51 ID:76X7tEQ5 タイトルに"遅延評価"が入っているが、遅延評価自体の問題ではなく、 Rの実装に原因がある。Haskellも遅延評価だが、`uncertainly`のような 挙動が起きれば、Haskellはその存在意義を失う。 ``` ghci import Data.Functor.Identity (Identity (..)) :{ uncertainly :: (Monad f) => f Int uncertainly = do x <- pure 1 plus_1 <- pure (x +) x <- pure 1000 pure $ plus_1 2 :} runIdentity uncertainly flip ($) () uncertainly maybe 0 id uncertainly ``` <pre> 3 3 3 </pre> 言語仕様レベルでのトレードオフかもしれない。 * 局所最適化と大域最悪化 * 整合性とコンパイル時間 * 等など http://mevius.5ch.net/test/read.cgi/tech/1340339592/467
468: デフォルトの名無しさん [] 2023/02/19(日) 11:12:54.40 ID:CwNaXfgw magrittrの例を少し改変 # fns first <- function(x){ message("first") invisible(x) } second <- function(x){ message("second") invisible(x) } # lazy NULL %>% first() %>% second() # eager NULL %!>% first() %!>% second() http://mevius.5ch.net/test/read.cgi/tech/1340339592/468
469: デフォルトの名無しさん [] 2023/02/19(日) 23:57:13.24 ID:5hVGJ6gZ ビックリパイプは知らなかった。"抜かりなし"だね。ただし、ビックリパイプは ハイゼンバグとは関係ないかな。 ハイゼンバグから離れて、宇宙大戦争について書いてみたい。 `lapply`等のループ系の関数は、最終的にC関数`R_forceAndCall`を呼び出す。 * [apply.c](https://github.com/wch/r-source/blob/trunk/src/main/apply.c) Rは、"フォースと共にあらんことを"などと呑気なことを言っている場合ではなく、 フォースと共にあらねば死んでしまう。 前回のコード`uncertainly`は変数`x`の使い回しが敗因だが、ループ系の関数 では変数を使い回すしかない。そのために、ループ系の関数ではフォースが必須 になっている。98パーセントぐらいの使用例では、フォース抜きでも動作する (当社調べ)が、残りの2パーセントでコケる。僅か2パーセントでも、 プログラミングをギャンブルにしないためには、穴を塞いでおく必要がある。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/469
470: デフォルトの名無しさん [] 2023/02/23(木) 21:58:25.40 ID:OvjjRMrU do.call()ならどうだろう http://mevius.5ch.net/test/read.cgi/tech/1340339592/470
471: デフォルトの名無しさん [] 2023/02/26(日) 00:14:58.50 ID:Qcp1BikG 上のコード`uncertainly`の挙動を"プロミスの罠"と書くことにする。 ここでの"プロミス"は、JSの"プロミス"ではなくて、Rでの"変数"の 実装方法を指す。 * [6 Functions | Advanced R](https://adv-r.hadley.nz/functions.html) この記事には次の一節がある。 > You cannot manipulate promises with R code. Promises are like a **quantum state:**: ... 多分、ここでの"量子状態"はハイゼンバグと同じ現象を指しているんだと思う。 観測すると、状態が変化してしまう。 関数`do.call`自体はプロミスの罠と関係しないと思う。 ``` {r a_28697, dependson = ""} do.call (`+`, list (1L, 2L)) ``` Pythonだと次のコードに対応する。 ``` {python a_10786, dependson = ""} (lambda x, y: x + y) (* range (1, 3)) ``` 殆どのプログラミング言語で、関数の引数リストは [一級市民](https://en.wikipedia.org/wiki/First-class_citizen) でないように思う。一級市民でない代わりに、一級市民の配列からの変換が 用意されている。Pythonでは`*`という関数がビルトインで用意されている。 Rでは引数リストを直接作れるかもしれないが、`do.call`で配列から 変換するのが一般的だと思う。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/471
472: デフォルトの名無しさん [] 2023/02/26(日) 23:03:36.96 ID:Qcp1BikG 書けるかな? http://mevius.5ch.net/test/read.cgi/tech/1340339592/472
473: デフォルトの名無しさん [] 2023/02/26(日) 23:09:42.38 ID:Qcp1BikG 自分の知る限り、プロミスの罠にハマるのは次のパターンに限られる。 ~~~ {.r} x = "hello" g = f (x) x = "world" g () ~~~ "関数を返す関数"`f`に、"変数"`x`を代入した時にプロミスの罠が可能性が 出てくる。Rの評価戦略はHaskellと同じ [コールバイニード](https://en.wikipedia.org/wiki/Evaluation_strategy) に分類されている。上の例では、次の場合にプロミスの罠が発生する。 1. 関数`f (x)`の中で引数`x`がニードにならず、 1. 返り値の関数`g`に渡される。 次の例はプロミスの罠が発生する。 ``` {r a_20914, dependson = ""} rude = with (new.env (), { const = \(x) \(...) x x = "hello" g = const (identity (x)) x = "world" g () |> print () }); ``` 関数への代入はニードでないので(by definition)、関数適用の連鎖 `const (identity (x))`の中にはニードがない。そのために、最終的に`g ()`が コールされるまで引数`x`は評価されない。コード`g ()`がコールされた時点で 引数`x`の定義を探しに行くので、プロミスの罠にハマる。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/473
474: デフォルトの名無しさん [] 2023/02/26(日) 23:11:33.81 ID:Qcp1BikG コールバイニードを実現するためには、コンパイラーが内部的に次のように 書き換える必要があると思う。 ``` {r a_1451, dependson = ""} polite = with (new.env (), { const = \(x) \(...) x x_1 = "hello" g = const (identity (x_1)) x_2 = "world" g () |> print () }); ``` [静的単一代入](https://en.wikipedia.org/wiki/Static_single-assignment_form) と呼ばれる操作と同じだと思う。この書き換えはループ処理では必須になる。 コード`rude`はアカン奴として切り捨てることができるかもしれないが、 ループ処理では変数の多重定義が避けられない。Rでは、何らかの理由で、 この書き換えが省略されているために、プロミスの罠が発生するのかな?と 思っている。 コード`rude`からコード`polite`への書き換えを上とは異なる形で行ってみる。 Rではできないので、Pythonを使う。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/474
475: デフォルトの名無しさん [] 2023/02/26(日) 23:12:32.78 ID:Qcp1BikG ``` {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" return g () def rhs (): return (lambda x: (lambda g: (lambda x: g ()) ("world")) (const (id (x)))) ("hello") assert lhs () == rhs () none (); ``` 関数`lhs`から関数`rhs`への書き換え規則は次のようになる。[rule]{#rule} ~~~ {.python} y = f (x) z = g (x, y) rest (x, y, z) == (lambda y: z = g (x, y) rest (x, y, z) ) (f (x)) == (lambda y: (lambda z: rest (x, y, z)) (g (x, y))) (f (x)) ~~~ http://mevius.5ch.net/test/read.cgi/tech/1340339592/475
476: デフォルトの名無しさん [] 2023/02/26(日) 23:25:12.65 ID:Qcp1BikG この書き換え規則を適用した結果、関数`rhs`では、静的単一代入への変換は、 [ラムダ計算](https://en.wikipedia.org/wiki/Lambda_calculus)のアルファ変換 に置き換わる。この意味では、静的単一代入とアルファ変換は同じことになる。 実際のコードでは、分岐、ループ、副作用が入ってくるので、こんな単純は 話では済まないと思うが、極度に単純化すると、静的単一代入とアルファ変換は、 単に、方言の違いということになる。 Rの場合は、書き換え規則[rule](#rule)は等価な変換にならない可能性がある。 特に、プロミスの罠が現れた場合は、挙動が変わると思う。 余談だが、書き換え規則[rule](#rule)をモナドに拡張したものがHaskellの [do-記法](https://en.wikibooks.org/wiki/Haskell/do_notation)になっている。 ~~~ {.python} y <- f (x) z <- g (x, y) rest (x, y, z) == kleisli-extension (lambda y: z <- g (x, y) rest (x, y, z) ) (f (x)) ~~~ モナドが恒等関数の場合に[rule](#rule)に一致するので、do-記法を [rule](#rule)の拡張と言って差し支えないと思う。 おしまい http://mevius.5ch.net/test/read.cgi/tech/1340339592/476
477: デフォルトの名無しさん [] 2023/03/01(水) 23:45:23.13 ID:X+BX12eb [トランスフォーマー](https://en.wikipedia.org/wiki/Transformer_(machine_learning_model)) についての記事: オーサーヒアーかな? * [Coinductive guide to inductive transformer heads](https://news.ycombinator.com/item?id=34970877) 中を読んでいないが、アテンション機構を [ホップ代数](https://en.wikipedia.org/wiki/Hopf_algebra)で説明、もしくは、 ホップ代数を使ってアテンション機構の類似物を導出するという話のようだ。 この手の救世主的な話はがっかりすることが多いので、過大な期待は禁物かも。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/477
478: デフォルトの名無しさん [] 2023/03/02(木) 07:01:07.09 ID:FFiKbXDz g = do.call(const, list(x)) ではどうだろう http://mevius.5ch.net/test/read.cgi/tech/1340339592/478
479: デフォルトの名無しさん [sage] 2023/03/02(木) 07:26:39.75 ID:rDVggB7z RのS5クラス定義(setRefClass)ってsave、loadでシリアライズ、デシリアライズ出来ないのか デシリアライズされたメンバ関数にアクセスできない http://mevius.5ch.net/test/read.cgi/tech/1340339592/479
480: デフォルトの名無しさん [] 2023/03/04(土) 20:31:08.80 ID:0X36s/jF ## Rにおけるオブジェクト指向への取り組み 趣旨は新規フレームワークR7の紹介だが、既存のフレームワークの俯瞰図 としても使いやすいと思う。 * [What is R7? A New OOP System for R](https://www.jumpingrivers.com/blog/r7-oop-object-oriented-programming-r/) ## オブジェクト指向とは? 次の記事を取り上げてみる。 * [Object Oriented Programming Features of Rust](https://doc.rust-lang.org/book/ch17-00-oop.html) 内容については賛否両論だと思う。しかし、箇条書きできるということが、 オブジェクト指向が単一の概念でないことを示しているように思う。 ## S3はオブジェクト指向か? 比較のために、S3のキーワード`default`を使わないで書く。 ``` {r a_3770, dependson = ""} none = with (new.env (), { to_string = \(x, ...) UseMethod ("to_string", x) as_string = \(x, ...) UseMethod ("as_string", x) to_string.Base = \(x, ...) as_string (x, ...) as_string.Base = \(x, ...) to_string (x, ...) to_string.Derived = \(x, ...) paste0 ("hello ", class (x)) as_string.Derived = as_string.Base structure (identity, class = "Derived") |> as_string () |> print () }); ``` `Base`クラスは、関数`to_string`か関数`as_string`の"どちらか1つを実装せえ"と言っている。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/480
481: デフォルトの名無しさん [] 2023/03/04(土) 20:34:04.91 ID:0X36s/jF 似た内容を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で書いてみる。 ``` haskell class Base_ a where to_string :: a -> String to_string = as_string as_string :: a -> String as_string = to_string instance Base_ () where to_string :: () -> String to_string = ("hello " ++) . show as_string () ``` <pre> "hello ()" </pre> http://mevius.5ch.net/test/read.cgi/tech/1340339592/481
482: デフォルトの名無しさん [] 2023/03/04(土) 20:37:50.63 ID:0X36s/jF ワンセットの関数群をインターフェースと書く。 1. Pythonのクラスでは、 * 構造体の定義とインターフェースの定義と実装が同時に行われる。 1. HaskellのクラスやRのS3では、 * 構造体とインターフェースは別々に定義される。 * インターフェースの実装は既存の"構造体"に対して行われる。 HaskellやRの"構造体"には、通常の構造体に加えて、関数も含まれる。 構造体とインターフェースが別々に定義される点で、RのS3はHaskellのクラスに 近いように思う。 おしまい http://mevius.5ch.net/test/read.cgi/tech/1340339592/482
483: デフォルトの名無しさん [] 2023/03/26(日) 10:03:46.90 ID:z5V/p30e join_by()が便利で感動した 気象データのマージがすごくシンプルになった http://mevius.5ch.net/test/read.cgi/tech/1340339592/483
484: デフォルトの名無しさん [] 2023/03/26(日) 22:59:22.57 ID:Dr7EPpsh 使用イメージがわかない。どのあたりが便利ですか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/484
485: デフォルトの名無しさん [] 2023/03/27(月) 22:57:14.70 ID:huaNY8vp 任意の地点・期間における測定値についての風向風速などの気象の影響をみるため近傍アメダスの気象データをマージする際に、便利になったと実感しました 具体的にはマージの際にbetween(時刻,)の条件が追加できるようになったのが大きいです 処理が重いかコードが長いかだった上のマージ処理が、join_by(地点名, between(time, start_time, end_time))で済むようになりました http://mevius.5ch.net/test/read.cgi/tech/1340339592/485
486: デフォルトの名無しさん [] 2023/03/30(木) 00:47:43.72 ID:7ArmqfaI [非等価ジョイン](https://en.wikipedia.org/wiki/Relational_algebra)は 最近になって関数`dplyr::xxx_join`に導入されたらしい。 * Rで非等価結合 (2) * [dplyr 1.1.0: Joins](https://www.tidyverse.org/blog/2023/01/dplyr-1-1-0-joins/) [ラッキー](https://www.youtube.com/watch?v=xBJMr1v5Zuw)。 関数`vctrs::vec_locate_matches`には、以前から限定された形の非等価ジョイン のオプションがあったので、実用的な形にするのに時間がかかったのかもしれない。 リスト内包表記が使いたいのでPythonで書く。 ``` {python a_32372, dependson = ""} #| cache: false # def slow_matches (pred, zipper): #{ def go (a, b): #{ a = pandas.DataFrame ({"" : a}).groupby ("").groups b = pandas.DataFrame ({"" : b}).groupby ("").groups key = ((j, k) for j in a.keys () for k in b.keys () if pred (j, k)) val = (zipper (a [j], b [k]) for (j, k) in key) val = zip (* itertools.chain (* val)) return (numpy.fromiter (val, dtype = numpy.int32) for val in val) #} return go #} http://mevius.5ch.net/test/read.cgi/tech/1340339592/486
487: デフォルトの名無しさん [] 2023/03/30(木) 00:49:00.50 ID:7ArmqfaI 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)) go (lambda a, b: abs (a - b) <= 1, itertools.product) go (lambda a, b: abs (a - b) <= 1, lambda a, b: itertools.islice (zip (a, b), 1)) #} none (); ``` `dplyr::xxx_join`の数々のオプションは理解できていないが、`NA`絡み以外は、 コールバック`pred`と`zipper`でかなりカバーしてるんじゃないかと思う。 ジョインの操作自体は明朗会計だが、実行時性能がボッタクリ価格に転嫁される というプログラミング社会が抱える構造的な問題が反映されている。 おしまい http://mevius.5ch.net/test/read.cgi/tech/1340339592/487
488: デフォルトの名無しさん [] 2023/03/30(木) 18:24:14.72 ID:Yqa5DxWU 一応Rスレだしコード一式を載せるならRにしてほしいな http://mevius.5ch.net/test/read.cgi/tech/1340339592/488
489: デフォルトの名無しさん [] 2023/03/31(金) 01:01:10.23 ID:3FvfgxM7 正直、書き方を思い出すためだけに、話題にかこつけて、わざと異なる プログラミング言語で書くこともある。しかし、今回は純粋に、処理の道筋を 単刀直入な形で表せるプログラミング言語を選択した。 [線形回帰](https://en.wikipedia.org/wiki/Linear_regression)を例にとる。 Rで線形回帰はとても書きやすい。しかし、線型回帰をRのイディオムとして 覚えるよりも、線形代数として覚えた方が長い目で見た時に利益になると思う。 同じように、等価ジョインは [プルバック](https://en.wikipedia.org/wiki/Pullback_(category_theory)#Sets) に対応する。等価ジョインを、Rのイディオムとして覚えるよりも、 単純な`for`ループの形で覚えた方が長い目で見た時に利益になると思う。 リスト内包記法は`for`ループを簡潔に書くための記法になっている。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/489
490: デフォルトの名無しさん [] 2023/04/15(土) 20:49:08.02 ID:Mf5VNEJ8 プルバックの話の続きをアップした。 * [JSFiddle](https://jsfiddle.net/ytkhdpcq/show) プルバックのカリー化の話になっている。プルバックの話はこれでお終い。 洗練された話は次の記事にある。 * [locally cartesian closed category in nLab](https://ncatlab.org/nlab/show/locally+cartesian+closed+category) 表計算の立場で見ると、"何が問題なのかわからない"問題かもしれない。 ある規則を満たす表でサマリーするとプルバックのカリー化になる。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/490
491: デフォルトの名無しさん [] 2023/04/24(月) 19:43:15.59 ID:Jj998bTN 4.3.0入れてみた パイプがまたひとつ便利になっていい感じ Windowsでパスが長くても通るようになったのも嬉しい http://mevius.5ch.net/test/read.cgi/tech/1340339592/491
492: デフォルトの名無しさん [] 2023/04/26(水) 12:00:07.89 ID:GTOaZ7cP どなたか教えてくだされ 対照実験の医療統計をしたいと思って 傾向スコアマッチングにしようと思ってます 肺がん手術を実施した人たち200症例のうち10人だけ 手術後に合併症を起こしていました この場合 「合併症を起こしたグループ」・・・10人 「合併症を起こさなかったグループ」・・・190人になり 患者の背景で傾向スコアマッチングをキャリパー幅0.2でペア分けすると 各群8人ずつの比較になってしまいます。 これで優位差を出す統計するには数が少なすぎるのであまりにも乱暴な統計になりますか? それとも元々のサンプルが200人から選りすぐりのマッチングさせた8人ずつを選んでいるので8人ずつと少なくても問題ないですか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/492
493: デフォルトの名無しさん [] 2023/04/26(水) 18:48:06.52 ID:ZF8W3LJH 共変量のバイアスがより少ないサンプルが得られますが、サイズが8人の二群比較なのは変わらないんじゃないかと思います http://mevius.5ch.net/test/read.cgi/tech/1340339592/493
494: デフォルトの名無しさん [] 2023/04/26(水) 23:58:34.61 ID:GTOaZ7cP >>493 つまりサイズが8人ずつの比較では 優位差を言うには弱いと思いますか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/494
495: デフォルトの名無しさん [] 2023/04/27(木) 08:26:17.11 ID:vIGDv5il >>494 検定はサンプルサイズ込みの判定なので8人であっても二群に有意な差があったと言って差し支えありません 有意水準が甘ければ弱いとは感じますが 一方で、検定は標本のバイアスが除去できていることを保証してくれません http://mevius.5ch.net/test/read.cgi/tech/1340339592/495
496: デフォルトの名無しさん [] 2023/05/13(土) 14:17:33.54 ID:LcdFqRfa スレッドの質をガクッと下げるが、Rの非標準評価で遊んでみた。 * [JSFiddle](https://jsfiddle.net/fo1qjt87/show) http://mevius.5ch.net/test/read.cgi/tech/1340339592/496
497: デフォルトの名無しさん [] 2023/05/16(火) 14:12:25.04 ID:T6FZDpAX 4.3.0のWindows版はdir()にバグあるね リストがまともに取得できない 4.2.3に戻した http://mevius.5ch.net/test/read.cgi/tech/1340339592/497
498: デフォルトの名無しさん [] 2023/05/19(金) 19:12:10.55 ID:LBQCZv0F EZR使ってる人いる? http://mevius.5ch.net/test/read.cgi/tech/1340339592/498
499: デフォルトの名無しさん [] 2023/05/19(金) 23:40:20.00 ID:AtidLON9 EZRまだ更新続いてて驚いた すごいな http://mevius.5ch.net/test/read.cgi/tech/1340339592/499
500: デフォルトの名無しさん [sage] 2023/05/21(日) 15:55:45.05 ID:Ifo/fxoU 10年以上前からあるスレでようやく500ゲット! http://mevius.5ch.net/test/read.cgi/tech/1340339592/500
501: デフォルトの名無しさん [] 2023/05/21(日) 20:10:58.14 ID:F9CBOLLZ 750に達するのにもう10年を要するだろう http://mevius.5ch.net/test/read.cgi/tech/1340339592/501
502: デフォルトの名無しさん [sage] 2023/05/22(月) 00:40:13.35 ID:Cnh2SEE8 tibbleを利用してデータ探索しているんだけど、2つの変数で条件つけて散布図つくるにはどんな手順踏んだらいいのか… http://mevius.5ch.net/test/read.cgi/tech/1340339592/502
503: デフォルトの名無しさん [] 2023/05/22(月) 06:06:37.46 ID:wj8Gj0ut 変数1をプロットの色にして変数2をサイズにするとかじゃだめなの? http://mevius.5ch.net/test/read.cgi/tech/1340339592/503
504: デフォルトの名無しさん [sage] 2023/05/22(月) 09:44:06.92 ID:/JtWuqEq それも一つの方法ですね。ありがとう ただ、目的以外の条件をプロットしたくないのです。 tibbleの構成としては、観測が4時点、変数がたとえば6個(うち1つは個体識別番号で、変数としてはA-Eの5種類あるとしましょう)のような感じで、 散布図のx軸に「観測時点1における全個体の変数A」 y軸に「観測時点4における全個体の変数E」 をプロットしたいのです。 個体識別番号をキューにして、異なる時点の異なる種類の変数を散布図にするとでもいいましょうか… http://mevius.5ch.net/test/read.cgi/tech/1340339592/504
505: デフォルトの名無しさん [] 2023/05/22(月) 11:53:24.61 ID:wj8Gj0ut その要件だと以下のような感じでいけると思います 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) ) df |> pivot_wider(id_cols = id, names_from = time, values_from = c(a, e)) |> ggplot(aes(a_1, e_4)) + geom_point() http://mevius.5ch.net/test/read.cgi/tech/1340339592/505
506: デフォルトの名無しさん [sage] 2023/05/22(月) 14:40:05.73 ID:7X6vWo+R ありがとうございます! それをヒントに試してみます。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/506
507: デフォルトの名無しさん [sage] 2023/05/24(水) 23:00:43.06 ID:1oF+cvFs >>505 これ、もしかして途中ではないですか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/507
508: デフォルトの名無しさん [] 2023/05/25(木) 05:49:23.72 ID:NLUeSyH/ >>507 要件の理解が違っていなければ散布図作成まで行きますよ time1でのaの値とtime4でのeの値をid毎にプロットしています http://mevius.5ch.net/test/read.cgi/tech/1340339592/508
509: デフォルトの名無しさん [sage] 2023/05/25(木) 07:49:04.33 ID:DfoFl+x4 >>508 失礼しました こちらのビューワの問題で、ブラウザでみたらスクリプトの最後まで表示されました💦 http://mevius.5ch.net/test/read.cgi/tech/1340339592/509
510: デフォルトの名無しさん [sage] 2023/05/26(金) 13:34:57.26 ID:gPviguuO >>508 教えてもらったスクリプトでドンピシャでした、ありがとうございます。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/510
511: デフォルトの名無しさん [] 2023/05/27(土) 00:00:44.30 ID:qqVzW/zi どなたか統計素人の私に教えてください 医療で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検定で優位差があるかどうかを調べる。 ↑Q1. この認識であってますか? Q2. この正規分布を計るのはペアを作った後のデータ(n=60)で、正規分布かどうか?等分散かどうか?を見ますか? それともマッチング前のデータ(n=200)で正規分布かどうか?等分散かどうか?を確認するべきですか? Q3. アミラーゼじゃなくて他の連続変数の項目(BMIだったり、血圧だったり、白血球数だったり)の優位差を調べる場合も 全て上の流れで一つ一つの項目ごとで正規分布か?等分散か?など確認していって適宜、該当する算出法を項目ごとに採用して優位差を測る必要がありますか? もし良ければ教えていただけるとありがたいです http://mevius.5ch.net/test/read.cgi/tech/1340339592/511
512: デフォルトの名無しさん [] 2023/05/27(土) 19:37:35.09 ID:H+hysMph 素人だと思ううちは傾向スコアには手を出さないほうがいいような気がする http://mevius.5ch.net/test/read.cgi/tech/1340339592/512
513: デフォルトの名無しさん [] 2023/05/27(土) 22:19:54.02 ID:Z7qvu7kx >>512 上記の場合はどういうかんじでやるのが 最適解になるのかだけでも教えてもらえないでしょうか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/513
514: デフォルトの名無しさん [sage] 2023/06/02(金) 21:25:38.69 ID:AkfCPeTx writexlパッケージで保存すると時刻のタイムゾーンは反映されないのね http://mevius.5ch.net/test/read.cgi/tech/1340339592/514
515: デフォルトの名無しさん [] 2023/06/13(火) 17:38:19.62 ID:J5T/voP2 ggraggedパッケージがなかなか便利 facet_grid()を詰めて並べられる http://mevius.5ch.net/test/read.cgi/tech/1340339592/515
516: デフォルトの名無しさん [] 2023/06/18(日) 13:00:54.90 ID:13gXw1hf R4.3.1にした dir()関連が直ってるようでひと安心 標準パイプの新機能もようやく使える http://mevius.5ch.net/test/read.cgi/tech/1340339592/516
517: デフォルトの名無しさん [sage] 2023/06/18(日) 16:20:13.18 ID:ECqif9s2 しっぽりしてえな http://mevius.5ch.net/test/read.cgi/tech/1340339592/517
518: デフォルトの名無しさん [sage] 2023/06/19(月) 07:17:50.83 ID:lFY/Y6gc ggplot2のカラーパレットのデフォルトが変わったりしました? 棒グラフの色の割り当てが変わった気がする… http://mevius.5ch.net/test/read.cgi/tech/1340339592/518
519: デフォルトの名無しさん [] 2023/06/19(月) 19:42:23.18 ID:CYQnxNgF >>518 変わってないはず ただggplot2のデフォルトはグループ数に応じて関数で色を割り当ててるのでグループ数が異なると違うパレットのように感じるとは思います http://mevius.5ch.net/test/read.cgi/tech/1340339592/519
520: デフォルトの名無しさん [sage] 2023/06/19(月) 20:01:30.62 ID:b0yPHKbc 株価データを取得したいがためにRを弄り始めたけど中々難しいね ググってコピペしたコードを一つ一つパッケージと関数の挙動を調べているけど やりたいこと全てができるようになるまでは遠そうだ http://mevius.5ch.net/test/read.cgi/tech/1340339592/520
521: デフォルトの名無しさん [sage] 2023/06/19(月) 22:06:40.20 ID:J5tIl656 >>520 データをどこから持ってくるかによって随分難易度は違う。 ファイルとしてダウンロードできるところもあるし、 スクレイピングが必要なところもある。場合によってはせれにうむとか。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/521
522: デフォルトの名無しさん [] 2023/06/22(木) 23:31:13.20 ID:eZc4H/Dz 企業が求めるエンジニア言語ランキング 2位「R」 https://www.itmedia.co.jp/news/spv/2306/20/news183.html 3位以下を大きく引き離しててすごいけど… Rユーザの少なさを反映しただけな気もする http://mevius.5ch.net/test/read.cgi/tech/1340339592/522
523: デフォルトの名無しさん [] 2023/06/25(日) 11:51:31.98 ID:JqE+uKYL Rの求人なんてみたことない… http://mevius.5ch.net/test/read.cgi/tech/1340339592/523
524: デフォルトの名無しさん [] 2023/06/25(日) 15:23:19.17 ID:XZytU2Vy >求人数のトップ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の需要があったみたいだ 不思議 http://mevius.5ch.net/test/read.cgi/tech/1340339592/524
525: デフォルトの名無しさん [sage] 2023/06/26(月) 18:52:56.83 ID:wRy8Kzvb 結局企業で欲しいのは統計処理だから、統計処理ならRだろうっていう短絡的な 求人が多かったってことじゃない? http://mevius.5ch.net/test/read.cgi/tech/1340339592/525
526: デフォルトの名無しさん [] 2023/06/28(水) 18:59:42.22 ID:z2KrS2jG いくつか求人サイト見てみたけどRのみなんて求人はほぼ皆無でほとんどがPythonとの併記だった これ集計ミスかもね 同じ会社の別の期間のランキング見るとRは10位以内に入ってないし http://mevius.5ch.net/test/read.cgi/tech/1340339592/526
527: デフォルトの名無しさん [sage] 2023/06/28(水) 19:11:13.31 ID:pwcxWBkk 金が有ーる(R)といい なんちって http://mevius.5ch.net/test/read.cgi/tech/1340339592/527
528: デフォルトの名無しさん [] 2023/06/29(木) 08:36:21.11 ID:szIurvDh 業∞!!!! 字∞!!!!! 苗∞!!!!!! 米∞!!!!!!! 畑∞!!!!!!!! 灰∞!!!!!!!!! 糞∞!!!!!!!!!! 机∞!!!!!!!!!!! http://mevius.5ch.net/test/read.cgi/tech/1340339592/528
529: デフォルトの名無しさん [sage] 2023/07/23(日) 03:41:51.31 ID:SQtd3M8i 結構何年か前にRは結構古臭いし、言語仕様もわかりづらいし衰退するのでは、と書いているのを 読んだことがあるけど、逆に今だと他の商用統計パッケージが廃れる勢いなのにRって立場が盤石な 気がする なんで? http://mevius.5ch.net/test/read.cgi/tech/1340339592/529
530: デフォルトの名無しさん [] 2023/07/24(月) 10:46:27.45 ID:uQbqJHBb AIブームでPythonに全部もってかれると思ったけどそうでもなかったね 先発優位で統計屋や大学教授が昔からRユーザーってのもあるのかも? ユーザー視点からだとtidyverseの存在が大きいと思う 505のように地味に手間がかかりそうな処理が関数ひとつでひょいと処理できるとかすごい http://mevius.5ch.net/test/read.cgi/tech/1340339592/530
531: デフォルトの名無しさん [sage] 2023/07/24(月) 12:50:03.13 ID:IaCKgEWP ど〜も〜!ご存じオナニーマン チン太郎でーす! http://mevius.5ch.net/test/read.cgi/tech/1340339592/531
532: デフォルトの名無しさん [sage] 2023/07/25(火) 14:51:45.41 ID:iTChcdyR >>529 今でもC/C++が盤石なのと同じ理由な気がする http://mevius.5ch.net/test/read.cgi/tech/1340339592/532
533: デフォルトの名無しさん [sage] 2023/07/25(火) 22:12:29.13 ID:KSU0u1F8 とくに大学の研究室とかは、同じような機能のプログラムを別言語に コンバートする作業は雑用でしかないので誰も積極的にはやらない。 RでできるのならこのままRでいいじゃん、もしくはフォートランでもw みたいな感じ やったところでなにか新しい知見が得られるわけもなし。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/533
534: デフォルトの名無しさん [sage] 2023/07/26(水) 14:36:34.35 ID:/vrstjVu 昔の医学分野の統計処理はSASがデフォだったけど COVID-19で疫学の人(西浦先生とか)がRガンガン使ってるの見て疫学以外でもRが広がってる 学会行くとggplotで書いたグラフいっぱい見る http://mevius.5ch.net/test/read.cgi/tech/1340339592/534
535: デフォルトの名無しさん [] 2023/07/30(日) 15:01:13.36 ID:/zgOfttX ggplot2のデフォルトテーマは微妙なので変えているけどいつも悩む ggthemesをもとにそこから装飾することが多い http://mevius.5ch.net/test/read.cgi/tech/1340339592/535
536: デフォルトの名無しさん [sage] 2023/07/31(月) 13:12:54.07 ID:sgBBFIN2 >>533 むしろ結果が違ってたらωωω http://mevius.5ch.net/test/read.cgi/tech/1340339592/536
537: デフォルトの名無しさん [] 2023/08/08(火) 23:42:45.27 ID:D9D8OCal Rしばらく使ってないから忘れそう http://mevius.5ch.net/test/read.cgi/tech/1340339592/537
538: デフォルトの名無しさん [] 2023/08/23(水) 14:24:09.90 ID:LaIh/7Co ついにExcelでPython使えるようになるらしい マイクロソフトはRのベンダー買収してMRANとかやってたのに時代の流れかな http://mevius.5ch.net/test/read.cgi/tech/1340339592/538
539: デフォルトの名無しさん [sage] 2023/08/24(木) 07:11:35.29 ID:UfAeCzV0 Pythonでexcel弄るのは 15年以上前から出来てるのを今更 excel上のpythonで実行する方が遅い可能性まである http://mevius.5ch.net/test/read.cgi/tech/1340339592/539
540: デフォルトの名無しさん [] 2023/08/25(金) 13:03:41.78 ID:8Z9WZ/AE PythonからExcelをいじることとExcelでPythonが動くことは全く違う http://mevius.5ch.net/test/read.cgi/tech/1340339592/540
541: デフォルトの名無しさん [] 2023/08/25(金) 13:58:38.54 ID:5+gJach+ くりゃうどω http://mevius.5ch.net/test/read.cgi/tech/1340339592/541
542: デフォルトの名無しさん [sage] 2023/08/25(金) 14:03:48.50 ID:HPTKrBcu pythonは何故かPCでJupyter notebookが不具合で使えなくてな 導入時点でコケるとその言語をやる気失くす…Rはその点、導入が簡単で良かった http://mevius.5ch.net/test/read.cgi/tech/1340339592/542
543: デフォルトの名無しさん [sage] 2023/08/27(日) 01:44:38.85 ID:VsjIkj4j Excelも最近はTidydata志向で、パワークエリやパワーピボットで再現性のある加工や分析が できるようになってきているけど、使い始めはなかなかいいかな、と思ったがちょっと凝ると やっぱり複雑化してメンテできないVBAみたいになっちゃうなw ま、RやPythonだからメンテ性がいいというわけでもない気はするけどw http://mevius.5ch.net/test/read.cgi/tech/1340339592/543
544: デフォルトの名無しさん [] 2023/08/27(日) 20:22:33.62 ID:RTNKd3te Excelマクロ禁止してきた企業はPython機能も禁止しそう Rを社員が使うケースなんてのはもうレア度SSRなんだろうね http://mevius.5ch.net/test/read.cgi/tech/1340339592/544
545: デフォルトの名無しさん [sage] 2023/08/28(月) 13:50:39.79 ID:62AYCckx Rで学ぶなんとか分析とかいう本いっぱい出ているし、結構個人でやる気のあるやつは覚えるんじゃね? R人口増えている印象。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/545
546: デフォルトの名無しさん [] 2023/08/28(月) 13:57:24.90 ID:3YABDYXJ 昨日の統数研のワークショップ見て、大学や企業によってはライブラリのインストールに制限あったりでデフォルトで解析できるのがRが多言語に置き換わらない理由の一つなのかもなと思った 今後も「選択と集中」続けるどころか懸賞金型推進するようだから商用に金払うのも辛くなるとこも出てきそうね http://mevius.5ch.net/test/read.cgi/tech/1340339592/546
547: デフォルトの名無しさん [] 2023/08/29(火) 18:08:36.20 ID:vrxlGHM5 数が増えてることと数が少ないこととは違う これ混同する人は結構多い http://mevius.5ch.net/test/read.cgi/tech/1340339592/547
548: デフォルトの名無しさん [sage] 2023/08/30(水) 15:28:57.25 ID:zHni/Vbc 当社比200% 元が少なかっただけやろ みたいな http://mevius.5ch.net/test/read.cgi/tech/1340339592/548
549: デフォルトの名無しさん [] 2023/09/01(金) 21:33:15.97 ID:T0/12uyc 自分は今の部署にいる間だけは統計解析するけど異動したらもう解析なんてしないからRのモチベは低いな 人事の巡り合わせでいつかまた今の部署に戻るかもしれないけど http://mevius.5ch.net/test/read.cgi/tech/1340339592/549
550: デフォルトの名無しさん [sage] 2023/09/02(土) 10:50:19.16 ID:aOHtbbKy 隙あらばどうでもいい自分語り http://mevius.5ch.net/test/read.cgi/tech/1340339592/550
551: デフォルトの名無しさん [] 2023/09/02(土) 14:49:02.94 ID:DQO4DB1k 過疎スレでその煽りは馬鹿じゃねえの http://mevius.5ch.net/test/read.cgi/tech/1340339592/551
552: デフォルトの名無しさん [sage] 2023/09/02(土) 15:24:55.18 ID:mCX3wjBN 本人だろ age てるのは http://mevius.5ch.net/test/read.cgi/tech/1340339592/552
553: デフォルトの名無しさん [] 2023/09/02(土) 15:27:52.78 ID:DQO4DB1k ageを最近知ったみたいに強調してる馬鹿ww http://mevius.5ch.net/test/read.cgi/tech/1340339592/553
554: デフォルトの名無しさん [] 2023/11/26(日) 13:47:18.77 ID:TX3Gb+56 [プッシュアウト](https://en.wikipedia.org/wiki/Pushout_%28category_theory%29) について書いてみた。 [JSFiddle](https://jsfiddle.net/zLpqd9mj/show) 集合の圏でのプッシュアウトは [等価ジョイン](https://en.wikipedia.org/wiki/Join_%28SQL%29#Equi-join) の圏論的双対にあたる。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/554
555: デフォルトの名無しさん [] 2023/11/26(日) 19:19:27.89 ID:fOcIWwuE おーまんちょー! http://mevius.5ch.net/test/read.cgi/tech/1340339592/555
556: デフォルトの名無しさん [] 2023/12/02(土) 19:24:46.52 ID:UAUVJeJE 前回のおーまんちょーでは一部のコードで途中から行番号が表示されるようになっていた。 多分、JSFiddleのエディターで`tidy`というボタンを押してしまったからだと思う。 実は、`prettydoc::html_pretty`で行番号が表示されないのを直したいと思っていた。 しかし、どこから手を付けたら良いかわからなかった。それが、ハプニングのおかげで 手の入れ場所がわかった。多分、HTMLとしては不整合な状況になっているが、 意図した通りに表示されるようになった。 [ありがとう、おーまんちょー](https://jsfiddle.net/105nohbq/show) http://mevius.5ch.net/test/read.cgi/tech/1340339592/556
557: デフォルトの名無しさん [sage] 2023/12/15(金) 00:55:03.28 ID:KmwBi9fF group_by()はmap()と組み合わせて使えないんだね http://mevius.5ch.net/test/read.cgi/tech/1340339592/557
558: デフォルトの名無しさん [] 2023/12/17(日) 12:11:13.68 ID:j39jK9SN 一般的な意味では、そんなことないと思うよ。 状況によってはgroup_byとmapを組み合わせると上手く行くことがあると思う。 ぐるぐる図のくだりでグダグダなところがあったので、 [修正した](https://jsfiddle.net/6z7u3d09/show)。 アップした後にタイポを見つけたが、覆水盆に返らず。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/558
559: デフォルトの名無しさん [sage] 2023/12/17(日) 23:02:11.11 ID:/RmbQ2uq 多変数の時系列データをクラスタリングするpackageをご存知であれば教えていただけますか…? http://mevius.5ch.net/test/read.cgi/tech/1340339592/559
560: デフォルトの名無しさん [] 2023/12/18(月) 00:40:20.74 ID:+7hNKB64 知らない。 人がいれば、意味のない返事はしない方が良いかもしれないが、過疎だからね。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/560
561: デフォルトの名無しさん [sage] 2023/12/18(月) 23:51:29.23 ID:t5S7JJTy あなたは地球上で何日過ごしてますか? という問いに答える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 age=12*(0:9) + j-i age=age[18<=age][1:5] age AD=paste0(this_year-age,"-",tag) days=\(x) as.numeric(Sys.Date() - as.Date(x)) rbind(age,days=sapply(AD,days)) あってますか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/561
562: デフォルトの名無しさん [] 2023/12/19(火) 00:23:50.21 ID:jHCQnODd [Eto'o](https://en.wikipedia.org/wiki/Samuel_Eto%27o)はもう引退しました。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/562
563: デフォルトの名無しさん [] 2023/12/19(火) 00:58:14.13 ID:jHCQnODd マジで聞いているんなら、テストツールの使い方を覚えた方が建設的だと思う。 個人的には、Haskellからの流れで [quickcheck](https://github.com/RevolutionAnalytics/quickcheck) しか経験したことがないが、他にもあると思う。 コードが合っているかどうかわからないが、パット見問題なさそうに見える。 問題があったとしても、デバッグで十分対応できるように見える。 [リーナスの法則](https://en.wikipedia.org/wiki/Linus%27s_law): あたなが十分な暇人ならばバグはいずれ塞がれる。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/563
564: デフォルトの名無しさん [sage] 2023/12/19(火) 08:05:39.46 ID:5FW9R46b >>563 サンクス http://mevius.5ch.net/test/read.cgi/tech/1340339592/564
565: デフォルトの名無しさん [] 2023/12/20(水) 23:21:30.67 ID:H48HUSxu 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? ChatGPT: Certainly! To test the provided R code, you can follow these steps: 1. Understand the Code: * Make sure you understand the logic of the code and how it calculates age based on the Eto and birthdate information. 1. Run the Code: * Copy and paste the code into an R script or R console. * Ensure that there are no syntax errors or issues when running the code. 1. Test Cases: * Define several test cases with different combinations of Eto, birthdate (tag), and current date. * For example, you can use different values for eto and tag, and verify the calculated age against your expectations. 1. Expected Results: * Manually calculate the expected age for each test case based on the provided Eto information and current date. * Compare the manually calculated age with the age calculated by the R code. 1. Edge Cases: * Test the code with edge cases, such as birthdates close to the boundary of age ranges, to ensure the code handles them correctly. http://mevius.5ch.net/test/read.cgi/tech/1340339592/565
566: デフォルトの名無しさん [] 2023/12/21(木) 22:26:55.69 ID:Jt1+wTI1 出禁になった。書けるかな? てなわけで、ちゃっとじーぴーてーの残りの部分は割愛する。 意地悪クソ野郎と親切じーぴーてーにも共通することろがあって、 最後の締めはやっぱりリーナスの法則になる。言い方が全然違うけど。 親切じーぴーてーが例に挙げたテストコードがテストの基本形になると思う。 くいっくちぇっくのようなテストツールでは、テストケースをランダムに発生させるが、 そうしたテストツールが使えない場合も多々ある。テストツールが使えない場合は、 親切じーぴーてーが示したようなベタな方法で行くしかない。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/566
567: デフォルトの名無しさん [sage] 2023/12/21(木) 22:41:24.08 ID:6SaY/e5W >>566 書き込めるかどうか判らぬままに 黙々とその低能感溢れる文章を書いている姿を想像すると泣けてくる 書き込めてよかったな http://mevius.5ch.net/test/read.cgi/tech/1340339592/567
568: デフォルトの名無しさん [] 2023/12/22(金) 00:02:15.06 ID://WOG+3q 親切じーぴーてーとの会話を見返すと、関数を細分化して、(ナウ、干支、月日) という三つ組を引数とする関数を作るのもありかな?という気がしてきた。 このように関数を細分化すると、テストケースが作りやすくなる。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/568
569: デフォルトの名無しさん [sage] 2023/12/22(金) 15:14:15.66 ID:BVAe/0CV かけるかな http://mevius.5ch.net/test/read.cgi/tech/1340339592/569
570: デフォルトの名無しさん [sage] 2024/01/04(木) 21:13:27.82 ID:r+ggAGcN クラスタリングとか次元削減に関する教科書的な本ってありますか? http://mevius.5ch.net/test/read.cgi/tech/1340339592/570
571: デフォルトの名無しさん [] 2024/01/05(金) 00:08:11.56 ID:Cbq0Hg0X かけるかな 統計の人でも機械学習の人でもないので、返事を書くのは心苦しいものがあるが、 過疎なのでしょうがない。統計処理と機械学習の区別がつかないぐらいの素人だよ。 ありきたりだけど: * [朱鷺の杜Wiki](http://ibisforest.org/index.php?FrontPage) グーグルとアヒルで検索してみたけど、この検索については英語の方が情報が取りやすいと思った。 英語でも日本語でも本屋の広告が多くなるが、英語の方が俺が推薦図書はこれだ! みたいなのが多いような気がした。検索で見つけた本だけど非売品と商品の両方がある。 * [Modern Dimension Reduction](https://arxiv.org/abs/2103.06885) * [Modern Dimension Reduction](https://www.cambridge.org/core/elements/abs/modern-dimension-reduction/D1506D1D7D51CC4EF375E4633D54F377) このようなケースでは立ち読みしてから買うこともできる。時々、非売品と商品の 両方が用意されていることがあるので、注意した方が良いかもしれない。 殊、次元削減で言えば、線形代数みないなもので、ありとあらゆる場所で使われるので、 俺の次元削減にジャストフィットするものを探すのは難しいかもしれない。 ちなみに、専門的には、学校の教室で研究会を開くのが統計で、リゾート地で 国際会議を開くのが機械学習らしい。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/571
572: デフォルトの名無しさん [sage] 2024/05/02(木) 14:41:31.88 ID:9naKJa9o 統計解析向けのプログラミング言語「R」の実装に欠陥、任意コード実行のおそれ https://forest.watch.impress.co.jp/docs/news/1588724.html http://mevius.5ch.net/test/read.cgi/tech/1340339592/572
573: デフォルトの名無しさん [] 2024/07/06(土) 23:18:34.28 ID:dmyvQpwi 書けるかな? [Rのリサイクル規則](https://jsfiddle.net/c617uwp9) 中にも書いたが、polyfill.ioはシリアスな案件なので、外だししておく。 [polyfill.io問題のQuarto・R Markdownへの影響](https://zenn.dev/yutannihilation/articles/252c7795433429) http://mevius.5ch.net/test/read.cgi/tech/1340339592/573
574: デフォルトの名無しさん [sage] 2024/07/17(水) 11:23:25.89 ID:jPCMTQ1M ubuntu24.04対応のRstudio来ねーな・・・ http://mevius.5ch.net/test/read.cgi/tech/1340339592/574
575: デフォルトの名無しさん [sage] 2024/08/20(火) 00:42:46.77 ID:c8o+ygTe こねーなw 難航してるのかな。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/575
576: デフォルトの名無しさん [] 2024/09/10(火) 01:06:21.09 ID:SFqXDqGn 参考になるかな? * [RStudio fails to run on Ubuntu 24.04 - Issue #14336 - rstudio/rstudio - GitHub](https://github.com/rstudio/rstudio/issues/14336) 他人事なら正論が吐ける。Dockerで試すのが吉。 実際には面倒なのでなかなかDockerで試してからとはいかないが、 他人事なら何でも言える。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/576
577: デフォルトの名無しさん [] 2024/09/24(火) 00:03:09.15 ID:B2G5MWgw Hacker-Newsネタ 1. [React for R | Hacker News](https://news.ycombinator.com/item?id=41554193) 2. [Why wordfreq will not be updated | Hacker News](https://news.ycombinator.com/item?id=41578483) 3. [Data Science Agent and Code Transformation | Hacker News](https://news.ycombinator.com/item?id=41622080) 4. [Terence Tao on O1 | Hacker News](https://news.ycombinator.com/item?id=41540902) 5. [Brainfuck Enterprise Solutions | Hacker News](https://news.ycombinator.com/item?id=41620198) 6. [What happened to the Japanese PC platforms? | Hacker News](https://news.ycombinator.com/item?id=41612984) 一年に一本くらいRについての投稿があるような気がする。 2番目は、個人的には、考えさせられる内容になっている。 3番目と4番目は夢のある話で、5番目は楽しい話になっている。 6番目はCPU業界栄枯盛衰の物語としても読むことができる。 Intelの運命や如何に。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/577
578: デフォルトの名無しさん [sage] 2024/09/30(月) 19:17:21.97 ID:87Zg0Lej Rstudio ubuntu24版キタ http://mevius.5ch.net/test/read.cgi/tech/1340339592/578
579: デフォルトの名無しさん [sage] 2024/09/30(月) 19:32:04.73 ID:87Zg0Lej とりあえず騒がれてた起動の不具合は解消されてるみたいだね http://mevius.5ch.net/test/read.cgi/tech/1340339592/579
580: デフォルトの名無しさん [] 2024/10/10(木) 01:08:17.58 ID:Q3cGnE1w * [Nobel Prize in Physics awarded to John Hopfield and Geoffrey Hinton | Hacker News](https://news.ycombinator.com/item?id=41775463) ノーベル物理学賞が機械学習系のお二方に決まった。 業績から言って文句なしのお二方だが、それが物理学賞?という疑問が湧く。 話に聞いていた物理学賞の選考基準は「物理現象の新発見」だった。 これが理由でarxivの作者が物理学賞を受賞しないと思っていた。 理数系の研究におけるarxivの貢献は過去10年分の物理学賞と化学賞と フィールズ賞とチューリング賞を全て合わせたより大きいかもしれない。 個人的には次の2点が引っかかる。 * 何故、このタイミングで? * お二方とも高名である。 こんなことを考えているうちに化学賞も似たようなことになった。 おいくら兆円?が頭をよぎる。 AI関連の投資への呼び水として利用されている疑念が拭えない。 もちろん、何事も変わらなければ消えてしまうし、変わりすぎても消えてしまう。 物理学や化学の定義も時代と共に変わっていくものだと思う。 ノーベル賞の選考基準も時代と共に変わっていくものだろう。 しかし、今回のノーベル賞は偶然の一致の度合いが大きすぎる気がする。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/580
581: デフォルトの名無しさん [sage] 2024/10/10(木) 12:14:46.24 ID:bwO1aVtU 数学症も無いし CS賞部門も増設すべきだよな http://mevius.5ch.net/test/read.cgi/tech/1340339592/581
582: デフォルトの名無しさん [sage] 2024/10/11(金) 14:00:28.70 ID:W/ZPwHBO 放送大学の秋山先生の講座ってどうなの? http://mevius.5ch.net/test/read.cgi/tech/1340339592/582
583: デフォルトの名無しさん [] 2025/01/21(火) 21:45:31.47 ID:Q1AkQbXB Rには関係ないが、資料として置いておく。 * [PCAST Reports (2020-2025)](https://www.erichorvitz.com/PCAST_reports_2020-2025.htm) * [Documents & Reports | PCAST | The White House](https://web.archive.org/web/20250116064747/https://www.whitehouse.gov/pcast/documents-reports/) 中身は全く読んでいないが、次の討論会を見る限り、政治色はないと思う。 * [PCAST: Discussion of Artificial Intelligence (AI) Enabling Science and AI Impacts on Society - YouTube](https://www.youtube.com/watch?v=gZb7Yr4C8po) http://mevius.5ch.net/test/read.cgi/tech/1340339592/583
584: デフォルトの名無しさん [] 2025/02/04(火) 00:46:14.41 ID:V6kAFtDi この話題は政治色抜きでは語れないと思うが、それを抜きにして、 アメリカ政府の発行したデータを使っている人がいれば、注意が必要になると思う。 * [Backing Up US Federal Databases | The n-Category Cafe](https://golem.ph.utexas.edu/category/2025/02/backing_up_us_federal_database.html) * [CDC data are disappearing | Hacker News](https://news.ycombinator.com/item?id=42897696) http://mevius.5ch.net/test/read.cgi/tech/1340339592/584
585: デフォルトの名無しさん [] 2025/02/15(土) 21:36:58.32 ID:jpL9Lr/5 次の動画の音声だけ聞いていたところ、只者でない気配を感じた。訳者の腕もあると思う。 * [GPT解説2 アテンションの仕組み (Attention, Transformer) | Chapter6, 深層学習 - YouTube](https://www.youtube.com/watch?v=j3_VgCt18fA) そこで、動画を見ると、丁寧なアニメが動いているじゃありませんか。 手練が作っていると思うが、それでも相当制作時間がかかっていると思う。 おまけだが、声質やゆったりさ加減がiPadの広島弁の解説に似ているように感じた。 あんたあ中国新聞読むじゃろ。あんたが数学に合わせる必要はないんよ。 数学の方があんたのええように合わせるけえ。だといいね。 外のホームページでは、カテゴリ分けされているので、外の方が番組を選びやすいかもしれない。 * [3Blue1Brown](https://www.3blue1brown.com/#lessons) http://mevius.5ch.net/test/read.cgi/tech/1340339592/585
586: デフォルトの名無しさん [] 2025/03/10(月) 01:40:13.55 ID:trsagf4Z 0. [I'll be back](https://www.youtube.com/watch?v=kw41Yyfyh_0) 1. [Public health data disappeared. RestoredCDC.org is bringing it back | Hacker News](https://news.ycombinator.com/item?id=43248610) CDC was back. 続報 2. [Introduction to Stochastic Calculus | Hacker News](https://news.ycombinator.com/item?id=43160779) 伊藤さん is ubiquitously back. 3. [Differentiable Logic Cellular Automata | Hacker News](https://news.ycombinator.com/item?id=43286161) 車輪の大発明 will be back again and again. [量子解析学](https://en.wikipedia.org/wiki/Quantum_calculus) の言い出しっぺはオイラーやヤコビらしいが、これは [Brzozowski微分](https://en.wikipedia.org/wiki/Brzozowski_derivative) と通常の微分を結ぶ。読んでいないが、おそらく、記事では、 何らかの活性化関数の極限として論理式が再現するようにしているんじゃないかと思う。 量子力学は、ポアソン括弧の構造を保つように、ニュートン力学を変形する。 プランク定数がゼロの極限でニュートン力学になる。 量子解析学は、ライプニッツ則の構造を保つように、微分を差分に書き直す。 変形パラメータqがゼロの極限でBrzozowski微分になる。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/586
587: デフォルトの名無しさん [] 2025/05/22(木) 01:15:55.60 ID:+HZBQ7+B RStudioを更新したらreticulateが動作しなくなった。 * [Release reticulate 1.41.0 - rstudio/reticulate - GitHub](https://github.com/rstudio/reticulate/releases/tag/v1.41.0) 今までどうやって動作させてたんだ?記憶にございません。 というわけで、一からやり直した。ワイの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)でいける。 似た話題で、 * [Show HN: Rv, a Package Manager for R | Hacker News](https://news.ycombinator.com/item?id=44006735) http://mevius.5ch.net/test/read.cgi/tech/1340339592/587
メモ帳
(0/65535文字)
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.022s