データを区分線形モデルに適合させようとしています。セグメント化されたパッケージは、係数に制約を適用する方法が見つからない限り、本当にうまく機能します。線形セグメントの傾きをある範囲の値に制限する必要があります(たとえば、セグメント1では0〜0.1、セグメント2では> 0.5)。ここに制約のない例があります。ありがとう。 library(segmented)
#generate data for
私は時間の経過とともに従属変数である行列(5,000 x 5,000)を持ち、時間の経過とともに同じ形式でいくつかの行列を持ちます。両方のマトリックスには時々NAが含まれているので、これらはna.excludeで処理する必要があります。 私は私の問題を説明するためにいくつかのサンプルデータを作っ: y <- matrix(rnorm(25000),5000,5000)
x1 <- matrix(