間隔打ち切りデータを使用して生存分析を行っており、 incox パッケージの intcox() 関数を使用して Cox 回帰を実行しようとしています。私はすでに survfit() を使用して分析の一部を行っており、すべて正常に機能しました。
intcox を使用しようとすると、常に問題が発生します。
> intcox(Surv(tempo2,tempo1,type="interval2")~dados$sexo)
Error in copy.data[ord, ] :
object of type 'environment' is not subsettable
> intcox(Surv(tempo2,tempo1,type="interval2")~sexo, data=dados)
Error in if (any(derivs.wert$g1 <= 0)) { :
missing value where TRUE/FALSE needed
In addition: Warning messages:
1: In Surv(data$mix, lokal.cens) : Invalid status value, converted to NA
2: In coxph(formula, data) : X matrix deemed to be singular; variable 1
dados$sexo は 3 つの水準を持つ因子であり、dados は 156 個の観測値と 52 個の変数を持つリストです。Kaplan-Meier 分析を実行する場合や、間隔データなしで coxph() を使用する場合は問題ありません。他の変数を使用しても問題は同じようです。
R 3.0.1 を使用しています
**更新**
私は何も変更していませんが、エラーは次のようになります。
> intcox(Surv(tempo2,tempo1,type="interval2")~dados$sexo)
Error in intcox(Surv(tempo2, tempo1, type = "interval2") ~ dados$sexo) :
Invalid cens status
> intcox(Surv(tempo2,tempo1,type="interval2")~sexo, data=dados)
Error in intcox(Surv(tempo2, tempo1, type = "interval2") ~ sexo, data = dados) :
Invalid cens status