問題タブ [r-grid]

次の相関関係を最大化したいと思います。

ここで、qhat は長さ 90 の時系列です。mu_1,mu_2 ,mu_3 , andmu_4 are essentially the parameters I would like to run the grid search over;I is the 4 X 4 identity matrix;A is a 4 X 4 matrix of coefficients; hereinv D`represents the matrix inverse operator;は、次元 4 X 90 (4 つの変数の時系列) の行列です。同様の問題を検索して見つけ、ニーズに合わせてソリューションをハックしようとしました。私はRで次のことを試しました（冗長なコードをお詫びします。私が試みていることが明確になるように、すべてを分割したかったのです）：

オブジェクトのparam_spaceサイズは 14641 X 4 になります。as.vector上記の関数を使用して、数値以外のオブジェクトを使用して行列演算を実行する場合に R で発生するエラーを回避しました。t(varset)上記の一般化された式の D は、実際のコードでは here に置き換えられています。オブジェクトの寸法t(varset)は 4 X 90 です。

役に立つかもしれないし、役に立たないかもしれない情報は次のとおりです。単にコピーして貼り付けると

コンソールに入力すると、上で書いたコードを実行したときと同じエラーが表示されます。すなわち、

さらに: 警告メッセージ:

そのため、グリッド検索を不適切に設定したと考えています。どんな助けでも大歓迎です。これは実際には私の最初の投稿なので、このリソースにも感謝しています。うるさいわけではありませんが、R でのループを回避しようとしています。

-クリス

アメリカとカナダの格子セルに色をつけたいです。私の目標は、この質問と非常に似ています: R Plot Filled Longitude-Latitude Grid Cells on Map ただし、その質問は米国のみを扱っており、カナダを追加する方法がわかりません。

ここにあるコードを変更することで、米国とカナダの地図を描くことができました: https://groups.google.com/forum/#!topic/ggplot2/KAKhoE0GO4U

ここにあるグリッド セルの偽の属性データを作成するコードを次に示します。

以下はggplot2、以前のスタック オーバーフローの投稿で使用されたコードで、Lower 48 のマップに属性のグリッドを重ね合わせたものです。

ggplot2 つのステートメントを組み合わせて、偽の属性データのグリッドを米国とカナダの地図に重ねるにはどうすればよいですか? アドバイスありがとうございます。

ASCII TOMS グリッド形式のファイルを R に読み込もうとしています。R で開くような方法で読み込むことができました。ただし、iy は線形行列として開きます。ファイルの内容の概要は、次の場所にあります。

データ セットのサンプルは、次の場所からダウンロードできます。

このデータ セットは 2006 年 1 月のもので、使用する必要のあるものがかなりあるため、簡単にアクセスできるように名前を変更しました。私はそれを使用してそれを読みました:

各緯度のデータが 1 行になるように、これを R に読み込むにはどうすればよいですか? これは、適切なデータ構造を構築するのに役立つと思います。注: 試してみると、簡単に理解できました:
1 行がタイトルを示すような構造になっていることを意味します。だから私の問題は、次の「lat = ...」の前にこれらすべての要素を1行に読み込むことです。

さらに、これを使用して一連のファイルをループしようとしました:

boxplotヒストグラムの下に配置したい。私はすでにこれを行う方法を考え出しましたが、boxplotと ヒストグラムは同じサイズであり、 をスリムにしたいと考えていboxplotます。幅を小さくすると、端までのスペースは同じままです。しかし、全体の幅を減らしたい。

これが私がこれまでに持っているものです：

編集:
grid.arrange() を次のように使用する場合:

比率は私が望んでいたとおりですが、軸が再び整列するように上記の最初の例を調整する方法がわかりません。

誰？

grid.curve 関数でプロットされた曲線の色と太さを変更するにはどうすればよいですか? 関数パラメーターに引数を追加しようとしましたが、サポートされていません。どうにかしてデフォルトの par() を変更できますか?

データの回帰平面をプロットしたい:

構造（list（l = c（96.4155、76.803、71.5615、68.193、65.6975、74.627、67.82、64.26、62.06、60.35、68.284、63.7、61.04、59.05、57.569、57.569、57.569） 58.27, 56.42, 54.91, 53.65), a = c(-0.8375, -20.47, -22.1875, -22.6125, -22.5845, -2.2415, -12.91, -16.16, -17.9, -18.77, -2.4, -10.28, - 13.64、-15.33、-16.37、-2.4495、-8.66、-11.74、-13.46、-14.55、-2.4645、-7.58、-10.38、-12.06、-13.09)、b = c(1.437、-24.5915、-29.3275) , -31.892, -33.524, -5.151, -15.22, -19.61, -22.65, -24.8, -6.239, -13.24, -17.11, -19.64, -21.55, -6.7845, -12.18, -15.45, -17.68, - 19.47, -7.124, -11.53, -14.38, -16.41, -17.89)), .Names = c("L", "a", "b"), row.names = c(NA, 25L), クラス = "data.frame")

hereの例を試しました。そして、それが私の試みです：

3d プロットは問題ないように見えますが、グリッドでは次のエラーが発生します。

私のコードのどこに問題がありますか?私の問題を解決する別の方法はありますか? 助けてくれてありがとう！

アップデート：

plane3d() で試してみましたが、次のエラーが発生しました:

しかし、他の解決策を探したところ、私の意図が間違っていたことがわかりました。このようなものを探しています...

データ
の意図

...しかし、それを達成する方法、またはRに関する私の実際の知識でさえ可能かどうかはわかりません。誰か何か考えはありますか?

アップデート2:

私はこのコードで scatter3d() 関数を試しました:

しかし、このエラーが発生しました:

自由度を調整するには choose.k を使用する必要があることは既にわかりましたが、gam 関数については常に説明されており、scatter3d 関数でそれを使用する方法がわかりません。何か案が？

graphics低レベルのパッケージ機能を使用して、幅の 2 倍の高さの「交差した」長方形を描画するには、次のように呼び出します。

この図には優れた機能があります。縦横比が保持されるため、プロット ウィンドウのサイズを変更しても、同じ高さと幅の比率が得られます。

gridグラフィックスで同等の結果を得るにはどうすればよいですか?

私はRを初めて使用しますが、グリッドの背後にあるプロットを作成する必要があります.グリッドは1000ステップで0.004から0.08の間隔を持つ必要があります. 間隔のシーケンスは既にありますが、グリッドと軸を正しく取得する方法がわかりません。データをポイントとしてきれいにplot(data)取得し、線を使用してきれいな線を取得します。私は試しgrid(NA,1000,lwd=2)ましたが、これはプロットのグリッド線しか得られません。axis(1, xaxp = seq(0.004, 0.08,by = (0.08-0.004)/1000), las=2)これを試してみると、グラフィカルパラメーターxaxpの長さが間違っているというエラーが表示されるため、途方に暮れています。