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

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

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

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

Q&A

0回答

518閲覧

相関係数最大化を最適化組み合わせで求めたい!

nk_bankruptcy

総合スコア1

Python

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

0グッド

0クリップ

投稿2020/11/30 11:21

csvファイルで系列x1~x5とyを読み込む。系列x1~x5から任意の系列をピックアップし、それら系列の合計から新系列を作成する。その新系列と系列yの相関係数が最大となるような組み合わせを求めたい。

import numpy as np
import pandas as pd
import cvxpy
df = pd.read_csv('teratail.csv')
df_x = df.drop('y', axis = 1)
df_y = df['y']

def correlation_calc(x, df_x, df_y):
cnt_rows = len(df_x)
goukei = np.zeros(len(df_x))

for i in range(cnt_rows): goukei[i] = np.dot(x, df_x.iloc[i,:]) corr_coef = np.corrcoef(goukei,df_y) result = corr_coef[0,1] return result

x = cvxpy.Variable(len(df_x), boolean = True)
objective = cvxpy.Maximize(correlation_calc(x, df_x, df_y))
prob = cvxpy.Problem(objective)
prob.solve(solver=cvxpy.ECOS_BB)
result = [round(ix,0) for ix in x.value]

エラーメッセージは下記の通りです。

TypeError Traceback (most recent call last)
TypeError: only size-1 arrays can be converted to Python scalars

The above exception was the direct cause of the following exception:

ValueError Traceback (most recent call last)
<ipython-input-4-df78a7bc155b> in <module>
1 x = cvxpy.Variable(len(df_x), boolean = True)
----> 2 objective = cvxpy.Maximize(correlation_calc(x, df_x, df_y))
3
4 prob = cvxpy.Problem(objective)
5 prob.solve(solver=cvxpy.ECOS_BB)

<ipython-input-3-a06e9252fe40> in correlation_calc(x, df_x, df_y)
4
5 for i in range(cnt_rows):
----> 6 goukei[i] = np.dot(x, df_x.iloc[i,:])
7
8 # 相関係数を求める

ValueError: setting an array element with a sequence.

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

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

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

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

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

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.35%

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

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

質問する

関連した質問