質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.46%
Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

2670閲覧

pythonにおける基底変換の基底について

退会済みユーザー

退会済みユーザー

総合スコア0

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2020/09/25 07:03

編集2020/09/25 08:31

前提・実現したいこと

100次元の元データが300個あります.(300×100次元)
このデータについて主成分分析で導出した100次元の固有ベクトルが1本あります.

次に,元データとは別の100次元のデータ300個があります.
これらについて,先ほど導出した主成分軸(固有ベクトル)上でどの位置に投影されるかを確認したい(基底変換を行いデータをフィットさせたい)です.

リンク内容
上記のサイトを参考にしていますが,基底変換行列を求める部分で詰まってしまいました.

発生している問題・エラーメッセージ

計算後のデータが合わず進めることができません.
おそらく,元データの基底(E)の指定の方法がおかしいのだと考えていますが,どのように指定すればよいかわかりません.

(300, 100) (100,) numpy.linalg.LinAlgError: 1-dimensional array given. Array must be at least two-dimensional

該当のソースコード

# ファイル名を指定 X = pd.read_csv('hogehoge.csv', sep=",", header=None).values print(X.shape) #主成分軸取り出し U1 = pca.components_[0, :] print(U1.shape) #------------------------------------------------------------------------------- #テストデータの基底変換 #------------------------------------------------------------------------------- #元データの基底 E = np.eye(100) #基底変換行列計算 P = np.linalg.solve(E, U1) #------------------------------------------------------------------------------- # テストデータの主成分軸投影 #------------------------------------------------------------------------------- z = np.linalg.solve(P, X)

試したこと

P = np.linalg.solve(E, U.T) print(P.shape) >>(100,)

基底変換行列Pは正方行列になるはずなのですが,1行のベクトルになってしまいます.
基礎的な数学になってしまうのですが,よろしくお願いいたします.

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

最初の回答はほとんど関係ないようだったので
相変わらず出てる

numpy.linalg.LinAlgError: 1-dimensional array given. Array must be at least two-dimensional

については、reshapeで明示的にshapeを(100, 1)指定すると回避できます。
ただ、その場合でもsolveのところで正方行列じゃないという別のエラーが出ます。
これはlianlg.solveのアルゴリズムによるもので、ドキュメントでは正方行列以外についてはlinalg.lstsqを使うようにとなっているので、そっちを使ってみてください。


Python

1P = np.linalg.solve(E, U)

正しくは

Python

1P = np.linalg.solve(U1, E)

じゃないですか?

投稿2020/09/25 07:41

編集2020/09/25 10:31
yudedako67

総合スコア2047

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

退会済みユーザー

退会済みユーザー

2020/09/25 07:45

回答ありがとうございます。 UとEを入れ替えたところ, numpy.linalg.LinAlgError: 1-dimensional array given. Array must be at least two-dimensional というエラーが出てしまいます.
yudedako67

2020/09/25 08:20

変数名は'U'ではなく'U1'が正しいのではないかと思ったのですが、そこはそのままですか?
退会済みユーザー

退会済みユーザー

2020/09/25 08:30

はい,すみません.その通りです.どちらも100次元ベクトルです.コードは修正しておきます.
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.46%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問