0

パッケージ coxme を使用しており、最適なものを選択するために、モデルから AIC を抽出したいと考えています。

1) これを直接行う方法がわかりませんでした。関数 coxme() のコードを変更しない限り、これは不可能だと思いますが、間違っていれば幸いです。

2)次のコマンドで関数のコードを調べました:

coxme:::print.coxme

ストックAICに変数を追加してコードを見たが、たとえば「coxme2」と呼び(最初に​​coxme2 <-を追加するだけ)、それを使用しようとすると(他の追加なしで)、エラーが発生します:

colnames<-(のエラー*tmp*、値 = c("NULL"、"Integrated"、"Fitted": 'dimnames' [2] の長さは配列の範囲と等しくありません

要約すると、coxme 関数はうまく機能しますが、そのコードをコピーして貼り付けるだけでは機能しません。この問題を解決するにはどうすればよいですか?

4

2 に答える 2

1

試しましたgetAnywhere(print.coxme)か?R で関数をコピーして変更するために必要なすべてのコードが提供されます。

編集

書くだけ

copyOfAIC <<- temp

以下の

dimnames(temp) <- list(c("Integrated loglik", " Penalized loglik"), 
                           c("Chisq", "df", "p", "AIC", "BIC"))

そして、探している値のコピーを取得します

于 2013-06-24T11:56:01.093 に答える