統計解析R たぶんpart3くらい (587レス)
統計解析R たぶんpart3くらい http://mevius.5ch.net/test/read.cgi/tech/1340339592/
上
下
前次
1-
新
通常表示
512バイト分割
レス栞
抽出解除
レス栞
253: デフォルトの名無しさん [] 2019/02/22(金) 22:52:15.13 ID:+pFzFvdM ``` {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; } ``` http://mevius.5ch.net/test/read.cgi/tech/1340339592/253
402: デフォルトの名無しさん [sage] 2021/07/03(土) 00:33:16.13 ID:gqk6Zgdk 四捨五入や浮動小数点の話を蒸し返す感じになるけど… 1.275-1と0.275を第二位でround()してみる 前者は0.27で後者は0.28になる x1 <- c(1.275-1, 0.275) round(x1, 2) [1] 0.27 0.28 sprintf()で眺めてみるとなんとなく理由がわかる気がする sprintf("%.24f", x1) [1] "0.274999999999999911182158" "0.275000000000000022204460" 0.265でも試してみると両者とも0.26になる x2 <- c(1.265-1, 0.265) round(x2, 2) [1] 0.26 0.26 sprintf()で眺めてみると… sprintf("%.24f", x2) [1] "0.264999999999999902300374" "0.265000000000000013322676" 偶数丸めで切り下げになるケースだけ誤差の範囲が広いのかな?? http://mevius.5ch.net/test/read.cgi/tech/1340339592/402
424: デフォルトの名無しさん [sage] 2021/10/03(日) 22:36:37.13 ID:A1Dyon2q 適当に日本語のrmdからread_csv試してみたけど動作もメッセージも特に問題ないみたい 自分が条件を勘違いしてるだけかもだけど http://mevius.5ch.net/test/read.cgi/tech/1340339592/424
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
531: デフォルトの名無しさん [sage] 2023/07/24(月) 12:50:03.13 ID:IaCKgEWP ど〜も〜!ご存じオナニーマン チン太郎でーす! http://mevius.5ch.net/test/read.cgi/tech/1340339592/531
533: デフォルトの名無しさん [sage] 2023/07/25(火) 22:12:29.13 ID:KSU0u1F8 とくに大学の研究室とかは、同じような機能のプログラムを別言語に コンバートする作業は雑用でしかないので誰も積極的にはやらない。 RでできるのならこのままRでいいじゃん、もしくはフォートランでもw みたいな感じ やったところでなにか新しい知見が得られるわけもなし。 http://mevius.5ch.net/test/read.cgi/tech/1340339592/533
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
メモ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.022s