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

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

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

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

Q&A

解決済

1回答

631閲覧

繰り返し処理を用いて、行列を1000個ほど作りたい場合どうすればよいでしょうか?

hino0303

総合スコア1

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

0グッド

0クリップ

投稿2021/02/09 10:35

前提・実現したいこと

以下の行列を1000個ほど作りたい場合、どうすればよいでしょうか?

X <- matrix(c(2,-1,1,-1,1,0,1,0,1),3,3)%%matrix(rnorm(35,0,1),3,5);X

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

エラーは出ませんが、コンソール画面に何も表示されません

該当のソースコード

R

1for(i in 1:k) { 2X <- matrix(c(2,-1,1,-1,1,0,1,0,1),3,3)%*%matrix(rnorm(3*5,0,1),3,5);X 3}

試したこと

実際に作った1000個のXに対して、分散共分散行列を求めてその第1固有値λを用いて経験分布を作りヒストグラムにしたい。

補足情報(FW/ツールのバージョンなど)

ここにより詳細な情報を記載してください。

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

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

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

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

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

guest

回答1

0

ベストアンサー

オブジェクトX自体はちゃんと作成されているはずです。
ただ、繰り返しのたびに以前のXを上書きしているので、最終的に一つのmatrixができるだけです。k個のmatrixを作りたいなら、リストにしましょう。

X <- list() for(i in 1:k){ X[[i]] <- matrix(c(2,-1,1,-1,1,0,1,0,1),3,3)%*%matrix(rnorm(3*5,0,1),3,5); }

コンソール画面に何も表示されません

Rのコマンドラインは、素のオブジェクトを入力した場合のみはそれをすぐにコンソールに表示してくれますが、{}で括った中にオブジェクト名を書き込んでもそれを表示してくれたりはしません。コードの中ではprintなどを使いましょう。

投稿2021/02/09 12:47

編集2021/02/09 12:47
KojiDoi

総合スコア13692

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.35%

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

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

質問する

関連した質問