欠損 Y 応答変数適合モデル 比例ハザード私はいくつかの分析を 1 つのスクリプトに結合している生存分析を実行しようとしています。 私がやろうとしているのは、彼らが列をサバイバル プラットフォームに入力し、それらの列を取り出して比例ハザード プラットフォームで使用することです。 私のスクリプトは、Y 応答を除いてうまく機能します。Y 応答変数は見つかりません。 私が変わればY( As Column( Ycolname ) ) を Y("Time to Failure") のようにすると、うまくいきます。しかし、変数を使用しようとすると、機能しません。
</CODE></P><P>質問が 2 つあります。なぜ私の方法は Censor と X では機能するのに、Y では機能しないのですか。</P><P>2 番目の質問は、by 列または頻度列に入力しない場合、入力する列がないことを適合モデルにどのように伝えるかです。 As Column (variable) を使用すると、正しい列が見つからないと表示されます。</P><P> </P><P> </P><P> </P><PRE><CODE class=" language-jsl">Names Default To Here( 1 );
Clear Log();
dt = Current Data Table();
If( Not( Is Scriptable( dt ) ),
dtsel = Pick File( "Open File", {"JMP|jmp", "All files|*"} );
If( dtsel == "",
Stop(),
Try( dt = Open( dtsel ), Stop() )
);
);
dtname = Data Table( dt ) << Get Name;
validated = 0;
New Window( "Survival/Reliability",
<<MODAL></MODAL>
分析マックOSスクリプティングRe: 欠損 Y 応答変数適合モデル 比例ハザードご指導ありがとうございます、とても参考になりました！
私の問題は、最初の画面で頻度に入れられた列がない場合、
周波数名 =引数( 生存[ 4 ]、 1 );
これは空として返されます。
したがって、スクリプトを実行しようとした場合:
プラットフォーム = dt <<モデルをあてはめる(
頻度(桁( 周波数名 ) )、
沿って（桁（ 名前で ） ）、
Y(桁( Ycolname ) )、
検閲(桁（検閲名））、
効果(桁( エフェクト名 ) ),
個性（ 「比例ハザード」 )、
検閲コード( 「1」 )、
モデルの実行( 尤度比検定( 1 )、尤度信頼区間( 1 )、リスク比( 1 ) )
)
頻度の列がないと言うでしょうが、それは本当でした。
だから私がしたことは、8つのネストされたifループを配置して、どの列が空であるかを確認することでした。
必ずしも最良の方法ではないかもしれませんが、うまくいきました。
Re: 欠損 Y 応答変数適合モデル 比例ハザード問題を理解しているかどうかわかりません。 これが私が試して見たものです。
最初に Rats.jmp を変更してケースをモックアップし、コードをテストできるようにします。
列 By を追加し、データを 2 つのグループに均等に分割しました。
コードの最後の部分を次のように変更しました。
if (isempty(byname),
Platform = dt << Fit Model(
// Freq(),
// By(),
Y( Column( Ycolname ) ),
Censor( Column( Censorname ) ),
Effects( Column( Effectname ) ),
Personality( "Proportional Hazard" ),
Censor Code( "1" ),
Run Model( Likelihood Ratio Tests( 1 ), Likelihood Confidence Intervals( 1 ), Risk Ratios( 1 ) )
),
Platform = dt << Fit Model(
// Freq(),
By( column(byname) ),
Y( Column( Ycolname ) ),
Censor( Column( Censorname ) ),
Effects( Column( Effectname ) ),
Personality( "Proportional Hazard" ),
Censor Code( "1" ),
Run Model( Likelihood Ratio Tests( 1 ), Likelihood Confidence Intervals( 1 ), Risk Ratios( 1 ) )
)
);
コードを実行するために 2 つのパスを試しました。
まず、コードを実行し、次のように役割を入力します。
いいえ「によって」。コードの動作は以前と同じです。
次に、コードを実行し、次のように役割を入力します。
「By」を使用すると、最後にエラー メッセージが表示されます。
中間結果が異なると思いますが、コードは状況を正しく処理していません。
Re: 欠損 Y 応答変数適合モデル 比例ハザード私の2番目の質問について何か提案はありましたか？
バイ列または頻度列に入力しない場合、入力する列がないことを適合モデルに伝えるにはどうすればよいですか。 列 (変数) を使用すると、正しい列が見つからないと表示されます。
Re: 欠損 Y 応答変数適合モデル 比例ハザード私はこれを試してみましたが、うまくいきませんでしたが、コードを貼り付けると完全に機能しました。 理由はわかりませんが、うまくいきました。ありがとうございます。
Re: 欠損 Y 応答変数適合モデル 比例ハザードAs Column を Column に変更すると、問題を再現できません。 これが私が実行する完全なコードです。 最後の 3 つの As Column を Column に変更しました。
Names Default To Here( 1 );
Clear Log();
dt = Current Data Table();
If( Not( Is Scriptable( dt ) ),
dtsel = Pick File( "Open File", {"JMP|jmp", "All files|*"} );
If( dtsel == "",
Stop(),
Try( dt = Open( dtsel ), Stop() )
);
);
dtname = Data Table( dt ) << Get Name;
validated = 0;
New Window( "Survival/Reliability",
<<MODAL></MODAL>
これが私の手順と私が見たものです。 JMP14、JMP15、JMP16を確認しました。 あなたのコードが新しい名前リテラルをどのように使用しているかにより、あなたのバージョンはそれ以上前のものではないと思います。
- Rats.jmpを開く
- スクリプトを実行する
- ダイアログが表示されたら指定し、[OK] をクリックします。
- いくつかのレポートが表示されます。 そして、これは Fit Model 起動ダイアログが次のように表示されます:
ここで私が見ているものが見えない場合、私たちのどちらも気付いていない何かがあると疑うことしかできません. もしそうなら、あなたは連絡したいかもしれませんsupport@jmp.com助けを得るために。
Re: 欠損 Y 応答変数適合モデル 比例ハザードEval(PH) を一番下に移動すると、適合モデルが表示されなかったのと同じことがわかりました。そのため、移動した場所に移動しました。 Ratsでも同じエラーが発生しましたが、ファイルではこのエラーではありませんでした。理由はわかりません。 式を取り出してこのように実行すると、エラーがなくなり、前と同じ問題が発生します. Expr(PH) がスクリーンショットの下に貼り付けられていないと、他の列は表示されますが、Y は表示されません. このスクリプトは正しく、私の最新のものです。 何かアイデアがありましたら、お知らせください。 これは私にとって非常にイライラすることであり、その理由がわかりません!
Names Default To Here( 1 );
Clear Log();
dt = Current Data Table();
If( Not( Is Scriptable( dt ) ),
dtsel = Pick File( "Open File", {"JMP|jmp", "All files|*"} );
If( dtsel == "",
Stop(),
Try( dt = Open( dtsel ), Stop() )
);
);
dtname = Data Table( dt ) << Get Name;
validated = 0;
New Window( "Survival/Reliability",
<<MODAL></MODAL>Re: 欠損 Y 応答変数適合モデル 比例ハザードスクリプトを再確認していただけますか？ ここに表示されているダイアログにアクセスする方法がわかりません。
エラーメッセージが表示されました：
Eval(PH) をスクリプトの一番下に移動し、As Column を Column に変更すると、Fit Model ダイアログがまったく表示されません。
ところで、Rats.jmp を使用してテストします。 「days」を Y、「Censor」を Censor、「Group」を Grouping にします。
Re: 欠損 Y 応答変数適合モデル 比例ハザードColumn は As Column と同じことを行い、Y 以外の情報を入力します:
Re: 欠損 Y 応答変数適合モデル 比例ハザードAs Column の代わりに Column を使用すると、必要な結果が得られますか? どのような状況で、列ではなく列として使用する必要があるかは実際にはわかりません。 おそらく、この機能についてより詳しい知識を持っている私の同僚の 1 人がコメントできます。
