問題タブ [weibull]

R を使用して、行ごとに異なる「スケール係数」を持つワイブル確率の行で行列を埋めるのに問題があります。この時点で、各列は同じように埋められます。

私はいくつかの方法を試しましたが、以下に私がやったことを示します。私が見逃している基本的なものがあると思います！

26 の倍率のリストを取得するには:

行列を準備します:

行列の各行にワイブル関数を適用します (各 x 値 (風速) に対して)

また：

どちらの方法も機能しません。どんな助けでも非常に役に立ちます！

イベント履歴分析を実行するためのデータ セットを準備しました。イベント数を数えてみると、

走ると、

私は、

イベント数が event_info データ セットの数と異なる理由がわかりません。

残念ながら、データを投稿することはできません。

イベントがehaでどのようにカウントされるかについて誰かが光を当てることができますか

私は生存分析の初心者です。何も観測されなかった測定値をデータセットに含めることが重要かどうかを理解しようとしています。

次の架空のデータがあります: 3 人の患者、

1. 最初の患者は時間 12 で病気にかかっていることが観察され、
2. 2 人目は 13 時に 1 回、14 時に 1 回、病気がないことが 2 回観察されます。
3. 3 番目の患者は、時間 1 で病気にかかっていることが観察されます。

次の2つの例を試しました

と

結果は、2つの間に明確な違いがあることを示しています.患者が病気ではないという観察を含めることが重要である理由を誰かが説明できるかどうか疑問に思っていました. ありがとう！

SASで以前に作成したワイブル分布を使用して、生存分析を複製しようとしています-現在、ライセンスのないマシンで作業しているため、R（両方ともWindowsから）を使用しています。私の（右打ち切り）入力データは次のようになります。

Event と Censor はバイナリで、Score は約 150 から 450 の間で変化し、Time は 1 から 28 の間で変化します。入力データセットには約 30,000 行あります。

私がしようとすると：

次の警告メッセージが表示されます。

survreg.fit(X, Y, weights, offset, init = init, controlvals = control, : 反復を使い果たし、収束しませんでした。

そして出力なし。

このメッセージをオンラインで (およびこのサイトを介して) 検索しましたが、問題の原因を示すものを見つけることができませんでした。SASで同じデータを（procロジスティックと）liferegに通しても収束の問題はありませんでした。

この質問は 、多変量モデルをプロットしたいという点で、ここの質問に続きます。

パッケージのリリースノートsurvregの例から生存関数をプロットしたいとしましょう

モデルの場所

私が理解していることから、次の問題は問題に対処しているように見えますが、データセットを提供せず、1 つの変数のみを使用します。この例では、2 つの連続変数と 1 つの因子があります。

私の目標は、競合するリスク モデルのテストに使用できるデータ セットをシミュレートすることです。関数で簡単な例を試しているだけですsurvsim::crisk.simが、期待する結果にはなりません。

これらの数値は と同じであると予想されbeta0.evます。私が間違っているかもしれないことへのポインタ、または競合するリスクデータをシミュレートする方法に関するその他の提案。

補足: シミュレートされたデータのイベントが、リスクごとに異なるワイブル分布に従って発生することを望みます。データの階層とクラスターを指定できるようにしたいと考えています。打ち切りは、ワイブル分布またはベルヌーイ分布に従うことができます。

私は、相関データを含む独自のデータ セットに競合するリスク モデルを当てはめています。フィットは機能しますが、標準で奇妙な結果が得られます。私のモデルの要約のエラーとナイーブSE。誰かが同様の問題に直面していて、それを解決する方法を知っているかどうか疑問に思っていました.

データのワイブル パラメータ(つまり、形状とスケール)を特定したいと考えています。

この回答が提案したことは既に試しましたが、Python 3.4 を使用しています。

これにより、次の出力が得られます。

私のcsvファイルのデータは、ワイブル関数のy値ではなく、x入力値として読み取られたと思います。bin を使用して 2 番目の列 (または行) を追加すると、文字列値を float に変換できないというエラーが表示されます。

内部のデータをワイブル関数の y 値として使用するには、csv ファイルをどのように変更する必要がありますか?

私の問題は、この行を理解していないことにあると思います。

ここで何1, 1を表しているのですか？パラメータが負であってはなりません。

このデータをワイブル分布に当てはめようとしています:

私xとy変数は次のとおりです。

プロットは次のようになります。

私はこのようなものを探しています: 適合プロット

それにワイブル曲線を当てはめたい。次のように R で nls 関数を使用しています。

この関数は常に次のエラーをスローします。

そのため、最初にさまざまな開始値を試しましたが、成功しませんでした。開始値を「適切に」推測する方法がわかりません。次にSSweibull(x, Asym, Drop, lrc, pwr)、selfStart関数である関数を使用しました。現在、SSWeibull 関数は Asym、Drop、lrc、および pwr の値を想定していますが、これらの値が何であるかについての手がかりはありません。

誰かが私が進む方法を理解するのを手伝ってくれれば幸いです。

データの背景: bugzilla からいくつかのデータを取得しました。"y" 変数は特定の月に報告されたバグの数で、"x" 変数はリリース後の月数です。