問題タブ [loess]

100 列の数値を含むデータ フレームがあり、各列は最初の列の異なる循環シフトです。これらの列に 1 つずつ LOESS フィルターを適用する必要があります。私のコンテキストでは、共変量は単純で、 index だけ1, 2, 3, ..., <number of rows>です。

新しいデータ フレームで平滑化された値を取得するにはどうすればよいですか? ありがとうございました！

ローカル回帰モデリングでグローバルに適合した項の係数を抽出する方法はありますか?

関数 loess におけるグローバルに適合した用語の役割を誤解しているかもしれませんが、私が望んでいるのは次のとおりです。

x_1フィッティングは完璧に行われますが、グローバルにフィッティングされた項の推定値(上記の例では 0.25 に近い値) を抽出するにはどうすればよいですか?

手動で作成したブートストラップを使用して、黄土適合の信頼区間を構築しようとしました。しかし、私が得た信頼区間は、デフォルトで ggplot が返すものよりもはるかに狭かったです。不一致の理由を知っている人はいますか？

ここに私のdata.frameがあります:

そして、ブートストラップを構築しようとしたスクリプトは次のとおりです。

次に、ggplot のデフォルト CI (灰色の領域) と、ブートストラップされた CI の上限と下限 (ポイントとして) をプロットしました。

残念ながら、ご覧のとおり、私の CI は非常に狭いように見えます。

この不一致の理由を知っている人はいますか? 私は何を間違えましたか？

PS: 利用可能なブートストラップ機能があるにもかかわらず、なぜこれを手動で書いているのか疑問に思っている場合に備えて。これは、次のステップで黄土適合の一次導関数の CI を作成するために行います。したがって、これが機能する場合、各適合を数値的に微分し、97.5% と 2.5% のパーセンタイルを等しく取り、黄土適合の一次導関数の CI を取得します。

編集:ローランドが提案する分位数関数を使用しています。しかし、ブート機能の使い方がわかりませんでした。これは私が試したものです：

aosmith が見つからないことに気付いたので、return コマンドも追加しました。

aosmith が正しく指摘したように、関数は何にも割り当てていなかったため、実際には使用されませんでした。なんて初歩的なミス。また、return 関数を別の位置に配置しました。ここでは、まだあまり一致していないことがわかりますが、はるかに優れており、許容できると思います.

ありがとう@aosmith

R で LOESS 回帰を実行していますが、いくつかの小さなデータ セットで警告が表示されました。

警告メッセージ:

1: simpleLoess(y, x, w, span, degree = degree, parametric = parametric, : -2703.9 で使用される疑似逆数

2: simpleLoess(y, x, w, span, degree = degree, parametric = parametric, : 近傍半径 796.09

3: simpleLoess(y, x, w, span, degree = degree, parametric = parametric, : 逆条件数 0

4: simpleLoess(y, x, w, span, degree = degree, parametric = parametric, では、他にも近傍特異点があります。 6.1623e+005

これらのエラーについては、こちらの別の投稿で説明しています: R での黄土エラーの理解。

これらの警告は、LOESS 回帰用に設定されたスパンに関連しているようです。許容可能な平滑化スパンのパラメーターが 0.3 から 0.6 の間である他のデータ セットで行われた同様の方法論を適用しようとしています。場合によっては、これらの問題を回避するためにスパンを調整できますが、他のデータ セットでは、エラー/警告を回避するために許容レベルを超えてスパンを増やす必要がありました。

これらの警告が具体的に何を意味するのか、また、これが回帰を使用できる状況になるのかどうかについて興味がありますが、これらの警告が発生したこと、または回帰が完全に無効であることに注意する必要があります。

問題のあるデータ セットの例を次に示します。

そして、ここに改行なしのデータがあります

私が使用しているコードは次のとおりです。

ご協力ありがとうございます。追加情報が必要な場合はお知らせください。

編集：

ggplot2関数と を使用して、次のグラフィックを作成しますgeom_smooth()。私はコードを使用しました：

私のデータには、変数 x (値が 9 つしかありません) と y (メトリック) が含まれています。loess次に、コードで計算されたメソッドの投影点を追加します。

次に、これらの投影点をプロットに追加します。

青い点がggplotの青い黄土線上にないのはなぜですか? -methodに使用されるコードは、Rloessの標準の -function とは異なりますか?loess

情報: 8,000 を超える観測を含む元のデータには、疑似逆警告はありませんが、問題は同じです。

サンプル画像

関数を使用して、2 つの変数とlowessの間の回帰を当てはめています。の新しい値での適合値を知りたいです。たとえば、次の例で適合値を見つけるにはどうすればよいですか。それができることは知っていますが、誰かのプロットを再現したいのですが、彼は.xyxx=2.5loesslowess

次のデータセットがあります。

rate1 は時間依存関数で、数日かけて値 40 に近づく傾向があります。

rate1 が 40 になる日付を知りたいです。

私が必要とすることを行う 1 つの可能性は、rate1 を密集して設定し、目標値 (40) を減算することです。ポイントがゼロに近づくほど、目標値に近づきます。また、小さな「ノイズ」バッファを含めて、ゼロの範囲内のすべてのポイントを受け入れます。

Rで40に等しいrate1の時間値を見積もるにはどうすればよいですか? また、私が使用しているものよりもRに適切な方法があるかどうかを提案してください。すべての提案を歓迎します。