Python3 – 2次元配列からの要素の取り出し

概要

リスト、ndarrayの2次元配列について、いろいろな要素の取り出し方を整理する。

リストについては、以下の配列を準備する。

ndarrayについては、以下の配列を準備する。

各要素の10の位が行を、1の位が列を表している。

特定要素の取り出し

リストの場合は[行番号][列番号]で取り出す。

ndarrayの場合もリストと同じ。

ただしndarrayにはもう一つの方法があって、[行番号, 列番号]でも指定できる。

行の取り出し

リストの場合[行番号]で該当する行を取り出せる。

スライスを使って複数行を取り出すこともできる。

ndarrayの場合もリストと同じ。

スライスの場合も同じ。

ndarray用の指定の場合、列番号を省略。

行をスライス指定する場合。

列の取り出し

ndarrayでの列の取り出し方は以下の通り。ただし、取り出された列は行方向の配列に変更される点に注意。

スライスを使って複数列を取り出せる。

スライスを使うと、単一の列を列ベクトルとして取り出せる。

reshape()で行ベクトルを列ベクトルに変換してもよい。

 

なお、リスト、ndarrayとも、以下の指定方法では意図した列の抽出にならず、行の抽出になってしまう。

内包表記による列指定

リスト、ndarrayとも、内包表記を使って列を取り出せる。ただし結果は行ベクトルとなる。

複数列の取り出しは以下の通り。

リストの場合、表現が少し違うが、[行番号][列番号]で要素が取り出せる。

 

ndarray – 要素・行・列の挿入・追加

1次元配列

要素の挿入

ndarrayの1次元配列への要素の挿入にはnumpy.insert()関数を使う。第1引数に対象となる配列、第2引数に挿入位置、第3引数に挿入する値を指定。

指定された位置に値が挿入され、それ以降の要素は一つずつ後ろにずれる。

要素の追加

配列の最後に要素を追加するには、以下の方法がある。

  • numpy.append()
  • numpy.insert()
  • numpy.hstack()

numpy.append()

第1引数に配列、第2に引数に追加する要素を指定。メソッドではなく、モジュール関数。1次元配列に1次元配列を追加すると、単に要素が追加される。

numpy.insert()

挿入位置として最後の要素の次の位置を指定する。この値は配列のサイズと等しい。

numpy.hstack()

注意点として、引数として指定するのはタプルで、そのタプルの要素をつなげた結果が返される点で、追加される要素を配列としなければならない。

2次元配列

numpy.insert()で単純に位置指定だけすると、1次元配列として扱われてしまう。

行の挿入・追加

1次元配列を行として挿入するには、numpy.insert()で第4引数にaxis=0を指定する。

配列の最後に行を追加するには、numpy.insert()で最終行の次の位置=行数を指定する。

shapeプロパティーで各次元のサイズの配列が得られるので、行数として第0要素の値を使う。

numpy.vstack()を使ってもよい。vstack()は複数の配列を縦方向につなげていく関数で、つなげる関数群を要素とするタプルを引数とする点に注意。

列の挿入・追加

numpy.insert()で列を挿入する場合、axis=1を指定する。列の追加の場合でも、追加する配列は行形式の1次元配列で指定。

numpy.insert()で配列に列を追加する場合、指定位置はshapeの第1要素の値。

numpy.hstack()を使ってもよい。hstack()の注意点は以下の通り。

  • 引数をタプルとすること
  • 追加する配列は列形式とする必要があり、1次元の配列を列形式に変換すること

 

Python3 – リストの要素・行・列の抽出

概要

リストのスライスについてまとめる。リスト要素にアクセスするために[]で要素を指定する、その方法。

1次元

要素の指定方法

まず標準的な方法。

[開始値]:[終了値]を指定して、複数の要素をリストとして取り出せる。開始値・終了値のいずれかを省略すると、それぞれ先頭・最後尾に対応し、両方とも省略すると”全ての要素”を意味する。

Pythonのリストでは、要素位置は0から始まり、”終了値”は”その1つ手前の要素まで”を意味する。

特に開始値と終了値を省略した[:]の表記は、すべての要素を示す。

[開始値]<[終了値]でなければならない(終了値が最終要素の位置+1 であることに注意)。この条件に満たないときは、エラーにならず空のリストが返される。

終端からの指定

要素位置を負の数で指定した場合、最後尾の要素からの位置指定になる。-1は最後尾の要素の位置で、-2、‐3はその1つ前、2つ前の要素。

要素位置を負数で指定する場合も、[開始値]<[終了値]というルールは同じ。

先頭要素、最後尾要素の指定

1次元のリストの先頭要素、最後尾要素の指定方法は以下の通り。

2次元

シンプルな2次元配列のリストの場合。

行の抽出

第1引数だけを指定すると、2次元配列の場合は指定した行を抽出。スライスで指定した場合は対応する行が抽出される。

特定の行のスライスは以下のように指定。

列の抽出

第1引数で行のスライスを指定したうえで列方向の要素を取り出そうとi意図して、以下のように指定するとエラー。

さらに第1引数を[:][n]のように指定すると、n行目が取り出されてしまう。

リストの内包表記で、取り出す行を変数にするとうまくいく。

 

Python3 – リストのコピー

変数の代入

リストなどのミュータブルオブジェクトのコピーには注意を要する。

単に変数を代入した場合、元のオブジェクトへの変更は同じ参照の全ての変数に影響する。

copy()関数・copy()メソッド

copyモジュールのcopy()関数を使うと、リストの各要素のインスタンスが新たに生成され、その内容がコピーされる。このため、コピー元への変更の影響はない。

deepcopy()関数による完全な複製

2次元配列など、ミュータブルオブジェクトの要素がミュータブルな場合は、copy()だけでは各要素への参照が変わらないため、下位の要素の変更が波及する。このような上位要素のみに対するコピーをシャローコピーと呼ぶ。

ミュータブルオブジェクトを要素に持つミュータブルオブジェクトの完全な複製を得るには、copy.deepcopy()関数を使う。

リストのコピーの別方法

copy()関数はcopyモジュールをインポートする必要があるが、以下の方法はcopyモジュールがなくても同じ結果が得られる。

copy()メソッド

リストのcopy()メソッドは、copy.copy()関数と同じく各要素の複製を生成する。

[:]

リストのスライスは元のリストの要素を生成し、新たなリストを返す。スライスを全要素[:]とすることで、元のリストの複製が得られる。

ただし、copy()メソッドもスライスによる方法もシャローコピーなので、多次元配列などのディープコピーはできない。

 

matplotlib.pyplot – グラフ

グラフの描画

pyplot

figure, axes

グラフの種類

3次元グラフ

tips

 

matplotlib.pyplot – 基本的なグラフの描画方法

概要

pyplotのみを使った基本的なグラフの描画方法。一般的にはmatplotlib.pyplotpltというエイリアスでインポートする。

この例では、pyplotの関数のみを用いて1つのグラフに2つの曲線を描画し、タイトルや軸ラベル、凡例などを追加している。

また、軸目盛の間隔を設定したり、グリッドを描かせたりしている。

実行結果は以下の通り。

ライブラリーのインポート

グラフ描画のためのパッケージ、matplotlob.pyplotpltとしてインポートしている。また、関数描画のためにnumpynpとしてインポート。

描画する関数

ここではpyplotplot()で描画するため、2つの三角関数を準備している。

xの値は共通で、numpy.linspace()で−π~+πの間で等間隔にデフォルトの50個の点を準備。endpoint=Trueで終端の点も含めている。

そのxの値に対して、sin、cosの関数の値をy1y2に保存。numpy.linspace()から得られるxや、それに対してNumpyのユニバーサル関数で計算されるy1y2ndarray

各種設定

タイトルと軸ラベル

pyplot.title()の引数でグラフのタイトル、pyplot.xlabel()pyplot.ylabel()の引数でx軸とy軸のラベルを設定している。

軸の範囲

pyplot.xlim()pyplot.ylim()で、それぞれの軸の下限値と上限値を設定している。x軸については−π~+πとしている。

軸目盛

pyplot.xticks()pyplot.yticks()の引数で、軸目盛に表示する値を設定している。引数は、表示させたい軸の値をリストとして与える。

x軸については−π、−π/2、0、π/2、πの値を直接与え、y軸についてはnumpy.linspace()関数で、-1~+1の間で等間隔な5点(-1、-0.5、0、0.5、1)を計算している。

グリッド(補助線)

グラフの描画エリア内に、軸目盛に対応した細い補助線を描くため、pyplot.grid()の引数をTrueで指定している。

グラフの描画

pyplot.plot()の引数にxy1及びy2を指定してグラフを描いている。2つのplot()関数を書くことで、関数は重なって描かれる。

それぞれのグラフに対してlabel引数を指定していて、それらは凡例表示に使われる。

凡例の表示

pyplot.legend()で凡例を設定している。

凡例にはpyplot.plot()で指定したラベルが使われるため、pyplot.legend()pyplot.plot()の後に書かなければならない。

また、凡例の表示位置をloc引数で指定していて、この場合は左上に凡例が表示される。凡例の位置は、upper/lowerleft/rightをスペースで繋いだ文字列で指定。

グラフの表示

pyplot.show()でグラフを表示させている。

 

Python3 – pandas

概要

pandasはPythonのライブラリで、Rのdata frameのようなデータフレームをDataFrameクラスとして提供している。

その由来はpanel-data-sらしいが、公式ドキュメントにはそういう記述は見当たらない。

DataFrame

 

pandas.DataFrame

情報・内容の取得

データの概観

基本的なパラメーターの取得

DataFrameの生成

要素の操作

列の操作

行の操作

DataFrameで直接行を指定して参照

DataFrameで直接行番号を指定する場合はスライスで指定する。1行取り出す場合でもdf[col:col+1]のようにスライスにしなければならない。

maskによる行の抽出

DataFrameの特定列に対する条件式は、その条件に応じたTrue/Falseを要素とするSeriesオブジェクトを返す。

そのSeriesオブジェクトをDataFrameの引数とすると、Trueに対応する行だけが抽出される。

これをまとめて次のようにも書ける。

抽出結果の特定の列を参照可能。結果はSeriesオブジェクト。

queryによる行の抽出

DataFrameのquery()メソッドでも行を抽出できる。

抽出結果から特定の列を参照可能。結果はSeriesオブジェクト。

行の追加

1行追加する場合

1行追加する場合、locで新たな行インデックスを指定して内容を与える。

複数行追加する場合

複数行(=配列)を追加する場合、その配列をDataFrameオブジェクトとしてからappend()メソッドを使う。ただし単純に追加すると、新たな行・列として拡張されてしまう。

これを回避するためには、追加する行の列名をcolumnsで明示する。ただしこの段階では、追加された行のインデックスが連続していない。

インデックスを連続させるためには、ignore_index=Trueを指定する。

DataFrameの結合

joinメソッド~横方向の結合

join()メソッドは元のDataFrameに引数のDataFrameを結合した結果を返す。結合は行のインデックスに着目して実行され、インデックスが同じデータは同じ行の列方向に追加される。すなわち、同じ行数で異なる列のDataFrameを結合するのに適している。

破壊的な処理ではなく、元のDataFrameはそのままで新しいDataFrameが生成される。

この例では2つのDataFrameの行が一致しており、列が異なるので、横方向に結合されている。同じ列名が存在する場合はValueErrorとなる。

3つ以上のDataFrameを結合したい場合(引数に2つ以上のDataFrameを指定したい場合)は、それら複数のDataFrameをタプルかリストでまとめて指定する。

合計の計算と合計欄

以下の配列を使う。

sum()メソッドはSeriesオブジェクトを返す。

合計欄としてDataFrameに追加する場合、行/列いずれかの合計を計算して追加し、その上で他方の合計を計算して追加すると、総計まで計算できる。

行の追加はlocプロパティーを通じて、列の追加は列名を指定して行う。

より詳しい合計や率の算出方法はこちら

ソート

 

以下のデータを使用

sort_values()~列名を指定してソート

デフォルトは昇順。元のオブジェクトは変更されず、ソート後のDataFrameオブジェクトが生成されて返される。

ascendingパラメーター

ascending=Falseを指定すると降順でソート。

reset_index()~インデックスの降り直し

ソート後にインデックスを振りなおす場合、reset_index()メソッドを使う。このメソッドも元のオブジェクトを変更しない。

降りなおした後のインデックスで元のインデックスを上書きする場合は、drop=Trueを指定。

複数列によるソート

複数列を指定してソートする場合、ソートの優先順にリストで指定。ascendingパラメーターに同じ要素数のTrue/Falseリストを与えて、列ごとの昇順/降順を指定可能。

axisの方向

各種演算・処理の引数でaxisを設定するときの、行方向・列方向がまぎらわしいので整理。

DataFrameの設定

列のインデックス設定

set_index()メソッドにより、第1引数keysで指定した列をインデックスに設定される。

表示行数・列数の設定

表示行数・列数の省略

行数・列数が多いDataFrameを表示させると、途中を省略して表示する。デフォルトで行数は60行を超えると省略され、列数は画面幅に応じて調整される。以下のコードは、pandasのオプションでdisplay.max_rowsdisplay.max_columnsも表示させている。

この場合、行数に61行をしているので行表示は省略され、列数はコンソールの幅に応じて省略されている(コンソールの幅を変えて実行すると、その幅に応じて納めようと省略される)。

表示行数の設定

表示させる最大行数を変更するには、pandas.set_option()display.max_rows属性を設定する。

表示列数の設定

表示させる最大列数を設定するには、pandas.set_option()display.max_columns属性を設定する。以下の例では20列すべてを表示させていて、コンソールの幅に収まらない列は表単位で改行されている。

ファイルの読み込み

CSVファイル

DataFrameにCSVファイルの読み込みには、pandas.read_csv()を使う。

pandas.read_csv(path)
pathは読み込むCSVのパス

ヘッダー付きCSVファイル

たとえば、以下のようなヘッダー付きのCSVファイル(ex_basic_w_header.csv)を考える。

このファイルが実行ファイルと同じディレクトリにある場合、以下のように__file__を使って読み込むとよい

データには自動的にインデックス番号が付加される。

ヘッダー無しCSVファイル

以下のようなヘッダーがないファイルの場合。

そのまま読み込んでしまうと1行目がヘッダーと解釈されてしまう。

引数にheader=Noneを指定すると1行目からデータとして読み込み、列名には番号が振られる。

引数としてnames=(列名1, 列名2, ...)を指定すると、ヘッダー無しとして読み込まれ、列名がセットされる。

行のスキップ

以下のように、先頭から何行か読み飛ばしたいファイルの時。

引数にskiprowsでスキップする行数を指定する。このとき、最初の行はヘッダーとして解釈される。

ヘッダー行の指定

ヘッダー行を含むファイルでヘッダーの開始行を指定したいときはheaderで開始位置を指定する。指定したヘッダーより前の行は無視される。

インデックス列の指定

以下のファイルは通し番号付きで、番号には列名がない(1行目の1列目が空)。

これをそのままデータフレームに読み込むと、1列目もデータと解釈されてインデックスが自動生成される。

引数にindex_col=0を指定すると、その列がインデックスとして認識される。

index_colを1以上の値にすると、それ以前の列は読み込まれない。ただし指定した列に列名が定義されていると、それもデータとみなされてしまう。

日本語

日本語などマルチバイトの文字を含む場合、対応する文字コードを指定する。

One-hot-encoding

DataFrameのget_dummies()メソッドは、属性データをone-hotの形にエンコードしてくれる。

Tips – 注意点など

配列で生成しない方がいい

コード上で簡単なDataFrameを生成するとき、配列で数値と文字列を混在させるとトラブルになる

DataFrameの行追加の速さ

DataFrameに行を追加する場合、基本的にはリストに変換してリストの状態で行を追加し、それをDataFrameに変換するのが最も速い。いくつかの方法による速さの違いはこちら

行単位/列単位の合計・率の計算

DataFrameオブジェクトの行単位/列単位の合計や、それらを使った率の計算方法

 

ndarray – 配列による要素指定

“Pythonではじめる機械学習”中、27番のコードで気になったので調べた結果。

Scitkit-learnを使ってk-最近傍法の予測結果を保存し、その結果を使って予測されたアヤメの種類を表示させている。

ここでpredictioniris_dataset['target_names']のいずれの型もnumpy.ndarrayであり、ndarrayの引数に整数ではなくてndarrayを使っている。

通常のリストでこれをやると、”リストのインデックスにリストを使ってはダメ”とエラーになる。

ところがndarrayの場合はその引数にndarrayを渡すことができて、結果はndarrayとして帰ってくる。

さらに、引数に渡す配列が複数の要素を持つ場合、各要素を引数とした場合に対応する元の配列の要素が配列として返される。

なお、元の配列がndarrayであれば、引数に普通のリストを渡しても同じ結果が得られる。

 

Python3 – os/os.path

主な関数

os.getcwd()
カレントディレクトリー。プログラム実行直後は実行中のディレクトリーを指している。実行ファイルのディレクトリー情報を得るには__file__を使うとよい
os.path.dirname(dir)
dirのディレクトリーのフルパス。
os.listdir(dir)
dirにあるファイル・ディレクトリーのリスト。
os.chdir(dir)
dirに移動。存在しない場合は’NotADirectoryError’をスローする。

実行例

実行結果

 

ファイル/ディレクトリーの判定

ファイルの判定はos.path.isfile(target)。ディレクトリーの判定はos.path.isdir(target)

実行結果