裏 RjpWiki さん プロフィール

  •  
裏 RjpWikiさん: 裏 RjpWiki
ハンドル名裏 RjpWiki さん
ブログタイトル裏 RjpWiki
ブログURLhttp://blog.goo.ne.jp/r-de-r
サイト紹介文文字通り,RjpWiki の裏を行きます
参加カテゴリー
更新頻度(1年)情報提供125回 / 365日(平均2.4回/週) - 参加 2014/11/25 07:58

裏 RjpWiki さんのブログ記事

  • Python で平均値? でも,R より遅い(その2)
  • コメントをいただいた> pythonは標準ではRで言うlistで配列を扱うので, 標準では速度的に不利です。それをRと同じようにベクトル(配列)で扱えるようにするのがnumpyですので,> numpy.mean(numpy.arange(1,10000000))> みたいにして比較するなら良いかと思います。> 上記のpythonの例だと, Rで言うなれば x にたような比較になってしまっています。> # つまりリスト処理のオーバーヘッドの計測をしてしまっている.なるほど,遅すぎ [続きを読む]
  • R か Python か?
  • ごく普通の数学問題を解くプログラムを R と Python3 で書いた場合,後者は前者の数倍遅いプログラムの可読性では,Python3 の「ブロック:インデントルール」で大変,読みにくいオブジエクト・オリエンティッドと関数の関係に一致性がなく,理解困難ということで,Python を使う意味を見いだせない [続きを読む]
  • 崩れないように積み上げて!
  • 崩れないように積み上げて!締め切りが 2017/11/21 10:00 AM なので,その 1 分後に投稿されるように予約直方体の箱を重ねて置くことを考えます。ただし、上の箱は下の箱よりも小さくないと、崩れてしまう可能性があります。そこで、大きな箱の上に小さな箱を置くことを考えます。ここで「小さい」とは縦と横の長さがともに短いことにします。つまり、中央に重ねて置くと、上から見たときに下の箱の輪郭が見えるようにします。また [続きを読む]
  • Python で平均値? でも,R より遅い
  • 1 から 10000000 までの整数の平均値Python3>>> import time; import statistics; start=time.time();statistics.mean
    (list(range(1,10000001)));print(time.time()-start)5000000.59.593008995056152 # 9 秒もかかる?>>> import time; impor
    t numpy; start=time.time();numpy.mean(list(range(1,10000001)));print(time.time()-start)5000000.51.1983740329742432 # numpy でも
    1 秒R> optons(digits=16)> system.time(print(mea [続きを読む]
  • Python で素因数分解 --- R の勝ちかな?
  • Python で素因数分解http://my-notes.hatenablog.com/entry/2017/11/19/022246結果の表示が特殊sympy がインストールされていない場合は,事前に一度だけ$ python3 -m pip install sympyその後,$ python3>>> import sympy as sym>>> sym.factorint(12878){2: 1, 47: 1, 137: 1}>>> sym.factorint(111111111111111111111111111111111111111111111111111111111111111){3: 2, 37: 1, 43: 1, 239: 1, 1933: 1, 4649: 1, 10837: 1, 233 [続きを読む]
  • マス目を回す
  • マス目を回す締め切りが 2017/11/17 10:00 AM なので,その 1 分後に投稿されるように予約【概要】4×4 のマス目があります。マス目のうちいくつかは穴になっています。abc-/d-ef/g-hi/opqr codeiqのような形で盤面の初期状態(abc-/d-ef/g-hi/opqr)と一連の操作(codeiq)を与えます。操作は一連の文字です。各文字は「その文字を中心として、隣接(斜め含む)しているマスの文字を時計回りに回す」という操作を示しています。ただし [続きを読む]
  • そこまで強弁する?
  • PIAAC データ解析https://oku.edu.mie-u.ac.jp/~okumura/stat/piaac.html> まったく相関がないように見える。ただ,強引に相関係数を求めて検定してみると :> -0.02 という非常に小さな負の相関があり,非常に小さいにもかかわらずp値はほぼ0で非常に有意である。このあたりが大量のデータを扱う上での落とし穴である。> この落とし穴をもうちょっと詳しく調べてみよう。例えば日本(392)だけに限って調べる。 :> やはり目で [続きを読む]
  • PIAACデータ解析
  • 奥村先生が,「PIAACデータ解析」を書いている。(途中みたいだが)https://oku.edu.mie-u.ac.jp/~okumura/stat/piaac.htmlそのきっかけはhttps://twitter.com/tmaita77/status/930020451899678720https://twitter.com/kohske/status/930218156508946432のようであるが,以下のようなプロットを作成すれば,実態が見えてくる。データは31カ国の wh: 週あたり仕事時間(連続量) と ll: 新しいことを学ぶのが好き(5段階)ただし,1 [続きを読む]
  • カウントゲームで先手が勝つのは何通り?
  • カウントゲームで先手が勝つのは何通り?締め切りが 2017/11/14 10:00 AM なので,その 1 分後に投稿されるように予約某SNSにおいて、チャットでAIと対戦できる「カウントゲーム」があります。ある数からスタートし、交互に最大3つまでの数をカウントダウンしていき、最後に「0」を言った方が負け、というゲームです。例えば、AとBが対戦し、Aが15からスタートして、以下のように進むとBの負けとなります。A「15, 14」B「13, 12」A [続きを読む]
  • js-STAR
  • js-STARversion 8.9.6j(β版) ちゃんとプログラム実行結果の検証をしているのか?みっともないぞ[直接確率計算2×2]で,どうやったら,「両側検定 : p=0.0000 ns (.10「カイ二乗検定の結果」で,4つのセルが全く同じ場合(例えば,10,10,10,10)でもカイ二乗値が0ではない( 0.1 になる)とは,なんたる惨状![直接確率計算2×2] 観測値1 観測値2-----------------------------------------群 [続きを読む]
  • スイッチを反転しても同じ数だけ点灯する?
  • スイッチを反転しても同じ数だけ点灯する?締め切りが 2017/11/07 10:00 AM なので,その 1 分後に投稿されるように予約設問 家庭に必ずある分電盤。その中にはブレーカーがあり、家庭内の電気スイッチやコンセントなどにつながっています。ここでは、1つのブレーカーが2つのスイッチにつながっているものとします。また、それぞれのスイッチに対して電球が1つずつ付いています。スイッチをONにしても、そのスイッチにつながるブレ [続きを読む]
  • 展開図上の反対側
  • 展開図上の反対側締め切りが 2018/0202 10:00 AM なので,その 1 分後に投稿されるように予約マス目を6つ指定します。6つのマス目を立方体の展開図として組み上げた場合に、最初の面の反対側がどの面になるのかを計算してください。【入出力】入力はglkmnqのようになっています。区切り文字なしで、面を示すアルファベット(図を参照)が6つ並んでいます。出力するのは、最初の面の反対側の面のアルファベットです。glkmnqに対応す [続きを読む]
  • 作詞支援ツールを作ろう
  • 作詞支援ツールを作ろう締め切りが 2017/11/04 10:00 AM なので,その 1 分後に投稿されるように予約設問 あなたが所属しているプロジェクトでは、様々な楽曲を扱っていましたが、ある日作詞家から作詞支援ツールの開発を依頼されました。それは、特定の言葉と同じ韻を踏む(母音が一致している)言葉を探すというもので、例えば「試合(SIAI)」と「気合(KIAI)」というように、 特にヒップホップでは欠かせない要素です。もちろん、 [続きを読む]
  • 隣の人と異なる仮装
  • 隣の人と異なる仮装締め切りが 2017/10/31 10:00 AM なので,その 1 分後に投稿されるように予約いよいよハロウィンの季節がやってきました。ハロウィンと言えば仮装ですね。ただ、せっかく仮装しても他の人と同じ衣装になるのは避けたいもの。そこで、全員が横一列に並んだときに、隣の人とは異なる仮装をすることにしました。m 人が仮装をしたとき、衣装の数がちょうど n 種類である「並び方」が何通りあるかを求めてください。 [続きを読む]
  • カーペット
  • カーペット締め切りが 2017/10/27 10:00 AM なので,その 1 分後に投稿されるように予約【概要】上図のように、正の整数が全て並んでいます。数をひとつ指定しますので、その数に上下左右に隣接しているマスの数を、昇順にコンマ区切りで出力して下さい。【入出力】入力は3のようになっています。ふつうに10進数です。出力は、入力で指定された数に上下左右に隣接しているマスの数を、昇順にコンマ区切りで出力して下さい。2,4,9の [続きを読む]
  • 「ループ・トラッキング」問題
  • 「ループ・トラッキング」問題締め切りが 2017/10/26 10:00 AM なので,その 1 分後に投稿されるように予約自然数 n に対し、関数 Fn(x) を次のように定義します(floor():床関数)。例えば n=10, x=1 のとき、F10(1) = floor(4×1×9÷10) = 3 です。さて、整数 k(0 ≦ k ≦ n)に対して、関数 Fn による変換を繰り返し行います。例えば n=10, k=1 のとき、F10 で 1 を変換すると 3 となり、さらに変換すると 8 となり [続きを読む]
  • urinomics !!!
  • > urinomics> 出典:『Wiktionary』 (2010/09/03 02:49 UTC 版)> 語源> urine + -omics> 名詞> urinomics (uncountable)> (physiology) The identification of the totality of the constituents of the urine of an organism 少し前に言っていた,「コイケノミクス」にしておけばよかったのにね(藁)せいぜい,かわいこぶって,「ゆりちゃのみくす」って。茶飲みを連想してほのぼのするか?しないなあ。に,してもだ。「な [続きを読む]
  • 正六角形ブロックの回転
  • 正六角形ブロックの回転締め切りが 2017/10/06 10:00 AM なので,その 1 分後に投稿されるように予約【概要】下図のような、正六角形マス目の集まりがあります。マス目のうちのいくつかにブロックが入っています。中心を指定するので、そこを中心にブロックを時計回りに 120度回したらどうなるのかを計算して下さい。【入出力】入力はa 000/0110/00000/0000/000のようになっています。空白の前は中心を示す記号です。a, b のいずれ [続きを読む]
  • ISBNのチェックディジットを計算して!
  • ISBNのチェックディジットを計算して!締め切りが 2017/10/03 10:00 AM なので,その 1 分後に投稿されるように予約設問 入力ミスなどを防ぐため、チェックディジットがよく使われます。代表的な例として、書籍を管理するときに使われるISBN(Wikipedia)があり、チェックディジットが使われています。ISBNには、かつて使われていた「ISBN-10」と現在主に使われている「ISBN-13」があります。最近出版されている書籍には両方が記載さ [続きを読む]
  • 「カウント・スリー」問題
  • 「カウント・スリー」問題締め切りが 2017/09/28 10:00 AM なので,その 1 分後に投稿されるように予約自然数を 1 から順に書き並べていきます。このとき、3 の数字が現れる回数を数えます。自然数 n に対し、ちょうど n 個目の 3 の数字が現れたときに書いている数を F(n) と定義します。例えば F(10)=35 です。下の通り、10 個目の 3 は、35 を書いているときに現れます。1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, [続きを読む]
  • パズルゲーム「2048」の組み合わせは何通り?
  • パズルゲーム「2048」の組み合わせは何通り?締め切りが 2017/09/26 10:00 AM なので,その 1 分後に投稿されるように予約設問 2048というパズルゲームがあります。(Wikipedia)iPhoneやAndroidなどのアプリだけでなく、Web上で楽しむこともできます。実際の動かし方は試してみるとよくわかるのですが、ここでは動かし方は問いませんので、以下の遊び方は理解しなくても問題ありません。以下、Wikipediaによる「遊び方」からの抜粋 [続きを読む]
  • 棒の長さを最小にするモビール
  • 棒の長さを最小にするモビール締め切りが 2017/09/19 10:00 AM なので,その 1 分後に投稿されるように予約設問 おしゃれなインテリアとして「モビール」があります。Wikipedia:モビールここでは、同じ大きさの n 個の飾りを、糸と棒を使ってバランスを取ることを考えます。ただし、棒の長さは整数で表され、棒の長さを整数で分割した位置にしか糸を吊るせないものとします。一般的には一つの棒の途中など複数箇所で飾りを吊るし [続きを読む]
  • wilcox.test と wilcox_test 前者だったら残念ですね!
  • http://my-notes.hatenablog.com/entry/2017/09/15/173257| My Notes| 統計学とかR(R言語)とかPython3とかプログラミングとかの覚え書きとか走り書きとか。 座右の銘にしたい: All work and no play makes Jack a dull boy.| 2017-09-15| R(R言語) ノンパラメトリック検定(独立サンプルの比較、独立した2群の中心位置の比較、Mann-Whitney (マン・ウイットニー) 検定(U検定))、wilcox.test() で,Mann-Whitney (マン・ウイットニー [続きを読む]
  • R の記事も多くなったけど,問題のあるものもある
  • > My Notes> 統計学とかR(R言語)とかPython3とかプログラミングとかの覚え書きとか走り書きとか。 座右の銘にしたい: All work and no play makes Jack a dull boy. > 2017-09-14 > R(R言語) ノンパラメトリック検定(適合度検定、カイ二乗検定、基準値との比較、一様性の検定)、chisq.test()http://my-notes.hatenablog.com/entry/2017/09/14/001050で,観測値が理論比に従うかの検定についての記事がある。=========== [続きを読む]
  • 円テーブルで席替え
  • 円テーブルで席替え締め切りが 2017/09/12 10:00 AM なので,その 1 分後に投稿されるように予約設問 先生1人と生徒 n-1 人の、合わせて n 人が丸いテーブルに座っています。途中で席替えをして、隣の人を変えることにしました。例えば、n = 5 のとき、最初に左図のように座っていると、右図のように席替えをすればすべての人の隣が変わることになります。このとき、先生は動かないものとします。標準入力から n が与えられたとき [続きを読む]