問題タブ [splines]

私は与えられx = c(50, 37, 25, 0)ますy = c(30, 52, 65, 70)。これらは、4 つの点の x 座標と y 座標を含むベクトルです。4 つの点すべてを通る滑らかな曲線を描き、この関数の実際の方程式を取得する必要があります。

でこれを達成するにはどうすればよいRですか?

スプラインのドキュメントを読んでいましたが、この背後にある数学に精通していないため、かなり混乱しました。

下の図では、Matlab でspap2関数を使用して、ノイズの多いデータを平滑化しています。結果はとても良いです。Eigenライブラリは、この機能 Splinesをサポートしています。同様の結果を得るために、アイゲンで例を探しています。Matlabの場合、私は使用しましたspap2(4, 4, time, noisyY);

データはこの形式で提供されますtime noisyData

指定された時間に予測平均を取得したい経時データがあります。モデルには、2 つの項、それらの相互作用、および時間変数のスプライン項が含まれます。予測平均を取得しようとすると、「mm %*% fixef(m4) のエラー: 非適合引数」が表示されます。

私の問題を説明するために、lmer の睡眠データセットを使用しました。まず、データをインポートし、インタラクション用の変数「age」を作成します

次に、lmer モデルを実行します

最後に、予測平均を取得するために必要なデータと行列を作成します

エラーが発生するのはこの時点です: Error in mm %*% fixef(m4) : non-conformable arguments

予測を使用して平均を取得できます

これは正常に機能しますが、信頼区間を計算できるようにしたいので、適合する行列が必要です。

問題は lmer がエイリアスを作成することである可能性があることをどこかで読みました（その投稿が見つかりません）。このコメントは、同様のタスクに effect() を使用できないことに関して行われました。この問題を克服する方法がよくわかりませんでした。さらに、投稿が少し古かったことを思い出し、エイリアスの問題がもはや関係ないことを願っていました。

私が間違っている可能性があることについて誰かが提案をしている場合は、フィードバックをいただければ幸いです。ありがとう。

ggplot2を使用して、制限された 3 次スプラインを使用して近似を説明したいと思いますgeom_smooth()が、正しく機能していないようです。以下に短い例を示します。

これにより、次の画像が得られます: スプラインの比較 でgeom_smooth()正しい結果が得られない 理由がわかりません。明らかに回避策があります (上図のように) がgeom_smooth()、正しい結果を生成する方法はありますか?

I am having difficulty finding information in fitting splines using rjags (my motivation is to try to recreate a glm in jags to impute missing dependent values). Anyhow I can find very little info on this finding only this answer on Cross-validated: https://stats.stackexchange.com/questions/79973/how-to-analyze-this-data-using-rjags-or-any-other-way/80650#80650. However, I cannot understand the splines code there (and do not have the reputation there to ask a question!). For one thing I don't understand why that code loops over both S & G.

Therefore, I have made a toy linear model in jags:

My question is: how to I add a restricted cubic spline to the 'b' estimate in the linear model ?

R3 次スプラインを使用して補間したい一連の測定値があります。これらは区分的な多項式にすぎないため、後で補間関数を代数的に統合したいと思います。したがって、係数が必要です。これらを入手する方法はありますか？
呼び出しsplines::interpSpline(foo, bar)$coefても、実際の多項式係数が返されないようです。

不完全な楕円に似た 2D 座標雲を当てはめるつもりでした。R の「splines」パッケージのbs()関数を使用してフィッティング曲線 (多項式スプラインの B スプライン基底行列を生成するために使用される) を取得しようとしたとき、次のようになりました。

しかし、私が意図したのはこのようなものでした。

パラメータを bs() にフィードしてプロットするには、これに従いました

私が意図したことを達成する方法はありますか？

ありがとう！

「Linear Quantile mixed models」、Geraci M. and Bottai M. (2013) で公開されたいくつかの結果を複製しようとしています。データは以下からダウンロードできます: http://www.hsph.harvard.edu/fitzmaur/ala/tlc.txt

自然な 3 次スプライン (df=3、ブレーク ポイント 0、4、6) とランダム効果を使用して、線形分位数混合モデルを使用しています。を使用してスプラインを作成しました

分位点混合モデルは

.

次のように、スプラインと処理 (相互作用と主効果) をモデル化すると、なんとか出力を得ることができます。

しかし、以下のように T1 と処理項 (主効果と交互作用) を入力しようとすると、エラーが発生します。

私は R と分位点混合モデルにかなり慣れていないので、これが私のせいなのか、コーディング中なのか、知識不足なのかわかりません!