問題タブ [piecewise]

したがって、MATLAB で heaviside 関数と ezplot 関数を使用して区分関数である f(t) をプロットしようとしています。現在、私はMATLABにまったく慣れていません。このエラーが発生する理由を誰かが知っていれば、それは役に立ちます。

f = '12+(-2t+8)*ヘヴィサイド(t-2)+(2t-12)*ヘヴィサイド(t-6)'

f = 12+(-2t+8)*ヘビシド(t-2)+(2t-12)*ヘビシド(t-6)

ezplot(f)

inlineeval の使用エラー (15 行目) インライン式のエラー ==> 12+(-2t+8).*heaviside(t-2)+(2t-12).*heaviside(t-6) エラー: 予期しない MATLAB 式です。

inline/feval のエラー (34 行目) INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);

ezplotfeval のエラー (52 行目) z = feval(f,x(1));

ezplot>ezplot1 のエラー (469 行目) [y, f, loopflag] = ezplotfeval(f, x);

ezplot のエラー (145 行目) [hp, cax] = ezplot1(cax, f{1}, vars, labels, args{:});

私はR次のような非線形モデルを扱ってきました。

そして、ベルヌーイ変数がAlpha1*timeand 切片に影響するかどうかを並べます。そのようなベルヌーイ変数は次のようになります。

持つ

がこのモデルに適合することがわかりましたが、このカテゴリ変数の効果がなければ、パッケージは非線形モデルの項ではなく変数に対する効果を推定できることnlsもわかりました。nlme私の質問は次のとおりです。これらのパッケージのどれが私を助けることができますか? そして、このベルヌーイ変数をコードに追加するにはどうすればよいでしょうか? よろしくお願いします。

次のようなプロットがある場合

プロット内にラベルを追加する人はいますか?

たとえば、x=1 の位置に「磁場に入る」などのラベルを付けたい場合、何か方法はありますか?

sympy を使用して区分関数を作成しました。

0 から 2XXの間の任意の値で計算するように機能します(対象範囲のみ)。x

それをプロットすると、負の数を分数で累乗できないというエラーが返されます。これxは、分数の累乗が適用されるのは が 1 より大きい場合のみであるため、決して起こらないはずです。

どういうわけかフルレンジで両方のピースを処理していますか？何かご意見は？

次のような間隔で定義された関数があるとします。

で微分すると

我々が得る

一方、私は取得したい

そのような結果を得る方法はありますか？

この例では、 segmented.lm (またはその他の関数) が明らかなブレークポイントを見つけられるように助けていただければ幸いです。

次のエラーで返されます。

solve.default(crossprod(x1), crossprod(x1, y1)) のエラー: システムは計算上特異です: 逆条件数 = 1.51417e-20

K を変更すると:

別のエラーが発生します：

segmented.lm(lm(y ~ x, data = data), seg.Z = ~x, psi = NA, control = seg.control(K = 1)) のエラー: 間隔内の 1 つのデータのみ: ブレークポイント(s ) 境界または近すぎる

R（バージョン3.1.1）のオフセットを使用して、ポアソン回帰でセグメント化された関数を使用しようとしています。私のセグメント化されたパッケージ i バージョン 0.4-0.0.

以前の投稿で、再現可能な例が示されました。

ただし、例を再現できません。次のエラーが表示されます。

他の誰かが同じ問題を抱えていますか？それとももっと良い - 説明？

ありがとうございました

/キラ

異なる時間間隔で異なるハザード率を持つ区分的指数分布から乱数を生成したいと思います。

私が理解したのは、標準指数分布の無記憶性を適用できるということです。

次のコードがこの目的に適しているかどうかは誰にもわかりますか? 特に、ベクトル 'S' には、対応するハザード率 (ベクトル 'lambda' に含まれる) に関連付けられた区間の上限が含まれます。

質問が明確であることを願っています。

私は長い間解決策を探していましたが、それでもうまくいきませんでした。私はすべての有益な提案にとても感謝しています!

Rで区分線形関数を定義しようとしています。これは、異なる間隔で異なる線形式を持つ関数を定義したいということです。これらの間隔は、xWerte と yWerte で指定されたデータ ポイント (およびこれらのデータ ポイントの数) に依存するため、切断ポイントの数とさまざまな関数ルールは修正されません (そのため、このがらくたを for ループではなく for ループで試していますifelse 構造...)

私がこれまで試みてきたのは、カット (区間の切断点) と値 = これらの区間のそれぞれの関数規則を定義することです。私が行ったことを理解できるようにコードを書き留めますが、この部分は機能しているように見えますが、私の質問は私のコードのこの部分には関係ありません...

これはすべて、私がやりたいことをしているようです。 さて、ここに私の問題があります：

これは私を示しています

これは vals[length(vals)]=vals[6] です。しかし、searchInt(0.8) = 5 なので、vals[5]=0.63806 と表示されるはずです。

異なる関数ルールを vals に追加することに成功したことを知っています。

Rが見せてくれる

したがって、さまざまな間隔の正しい関数ルールは存在しますが、fx によって適切に呼び出すことはできません。R は、以前は vals に異なる関数規則があったことを「忘れた」ように思えます。searchInt(x) の値に関係なく、vals に追加された最後の関数ルールを呼び出すだけです。

これはローカル変数とグローバル変数の問題である可能性があることがわかりましたが、修正方法がわかりません:(

読んでくれて、私を助けてくれてありがとう！

PS: 間違いをお詫びします。私はネイティブ スピーカーではありませんが、最善を尽くしました ;)