Pythonでは基本的に、Global Interpreter Lockといって、一つのプロセスでしか動かせない、という制約がある。これは、複数スレッドを立ち上げてPythonコードを実行しようとしても、CPUでの計算を…

自分用のメモです。回転同士の距離がどういうものか、といったことはここに日本語で詳しく書いてあります(とても良いページ!)。 from pyquaternion import Quaternion import numpy…

はじめに 3次元回転(金谷健一著)を読んで、あるコスト関数を最小化するための3次元回転の最適化について、リー代数の方法、の理論的なところを少なくとも一部理解したので、実際に実装してみようというのが今回の目的。 リー代数と…

Scipyによる最適化計算 Pythonで最適化計算を行う場合、Scipyのoptimizationパッケージを使うことになると思う。これには、いろいろな最適化アルゴリズム(Nelder-Mead, Powell, BF…

自分が後で見返す用のまとめです。 このコースはスタンフォード大学のAndrew Ng先生が公開している機械学習の入門コースです。機械学習について勉強するにあたり、まず見ておくといいと思います。無料で視聴することができます…

数値の代入 数値を変数に代入するにはハイフンと大小記号を使って矢印を作ります。以下ではxに2を代入しています。xの中身を表示するにはxと入力してEnterを押下します。 x <- 2 plot(xy) ベクトルの代…

1説明変数の線形回帰 xとyの関係が定数a、bを用いて、y=ax+bというように表現されると予想できるとき、線形回帰を行います。このページで使用するデータとして、まず、以下の仮想的なデータを代入します。散布図を書くと、ど…

データの用意 このページでは、乱数を使って仮想的なデータを作成します。正規分布に従う乱数はrnorm関数で生成できます。rnorm関数は第一引数に生成する個数、第二引数に平均、第三引数に分散を指定します。 x1 <…

非線形回帰 xとyの間に非線形の関係がある場合にも、Rはパラメーターの推定を行うことができます。このページで使用するデータとして、ロジスティック成長曲線に従う誤差を伴ったデータを生成します。ロジスティック成長曲線の解は、…

以下ではPythonで塩基配列を操作するコードを紹介します。説明は最小限に留めてありますが、C++がわかる人ならばPython初心者でも構文が似ているので意味はわかると思います。   1.DNA配列中の各塩基の…