cox-regression

0热度

1回答

我很困惑如何调整coxph函数中的变量。我知道可以使用strata（）进行分层，但是如何调整变量？线性模型中，可以通过以下（这里有一个我从这个link看到例子） fit.diamOnMachine <- lm(diameter˜machine) diam.adjusted <- residuals(fit.diamOnMachine) fit.diamadjmach <- lm(stren

3热度

1回答

Cox回归分层变量的参数估计和方差（地层/生存包）

我已经运行Cox回归，使用survival包来计算暴露A的死亡风险比。我发现age变量违反了比例风险假设（与cox.zph）并使用strata(age)分层进一步模型的年龄。我需要age变量的参数估计值，以及方差和协方差矩阵（以计算费率提升期间）......我不知道在哪里找到它们！我错过了什么，或者我误解了什么strata正在做什么？这里是一个重复的例子，使用从survival包lung数据

0热度

1回答

如何在R中保存ggcoxzph图？

我安装了Cox模型，然后使用cox.zph测试了比例风险假设。我使用ggcoxzph绘制了缩放的Schoenfeld残差，结果非常整齐，我想保存它。由于我的变量是一个因子，因此ggcoxzph会返回一个由4个图组成的网格，并且Global pvalue位于顶部。 > coxt %>% cox.zph %>% ggcoxzph %>% class [1] "ggcoxzph" "ggsurv"

2热度

1回答

JuMP.jl中的Cox比例危险

我正在第一次尝试Julia中的JuMP.jl，似乎无法解决错误。这是我的设置。使用DataFrames，DataFramesMeta，跳跃，Ipopt #time to event times = [143,164,188,189,190,192,206,209,213,216,220,227,230,234,246,265,304,216,244, 142,156,163,198,205,

1热度

1回答

JuMP中Cox比例风险的得分函数

我在将方程转换为代码方面经验不足。我被卡在将Juul中部分可能性的分数函数转换为代码以在JuMP中进行评估。 score function其中，当在0解决beta时，是最大值。我做了一个简单的小数据集。 Using DataFrames, DataFramesMeta, JuMP, Ipopt #build DataFrame times = [6,7,10,15,19,25] is_ce

0热度

1回答

从1月1日开始按年份划分生存数据集

我有以下生存数据集，我希望在每年的1月1日之前将时间间隔分开。例如，对于personid 1220，我会在1912-01-01,1913-01-01,1914-01-01,1915-01-01进行拆分。我试图使用survSplit，但他们只能做数字向量。如果有其他方法，请让我知道吗？在下面的数据集中，time = EndDate - StartDate。这是我到目前为止： test.ts <-

0热度

1回答

使用“pec”R软件包对肺数据集上的“coxph”功能进行预测

我想使用Coxph R函数为Cox比例风险模型获得从惩罚性R软件包中获得的肺数据集的预测。我有以下示例代码。 library(survival) library(pec) library(penalized) data("lung") data <- lung trainind <- sample(1:n,n*0.7) testind <- (1:n)[-trainind

0热度

1回答

使用类别作为新变量？

我自教过我几天后，我陷入了一个cox回归分析。我管理使用cut()函数到每个连续的变量我有分成2个分类组。现在我想知道是否可以在cox回归分析中单独合并这些类别。例如：我有3个变量（A，B，C），每个变量由2个类别组成：A-，A + & B-，B + & C-，C +。现在，如果我运行命名变量的coxph()，我只能得到“+”类别的结果（据我了解是因为“ - ”类别被用作参考组）。然而，由于C

1热度

2回答

R coxph（）错误：找不到对象'Ccoxmart'

我每次运行使用R中的生存包的Cox模型时都会收到以下错误。此错误在最近几天内出现。为了说明这个错误，我使用它在https://stat.ethz.ch/R-manual/R-devel/library/survival/html/coxph.html给出一个标准的例子命令： # Fit a stratified model, clustered on patients library(survi

0热度

1回答

Andersen-Gill计数过程中的无风险区间使用survival生存的Cox回归公式:: tmerge（）

我想使用tmerge()函数来转换数据集以用于Cox回归框架的Andersen-Gill扩展重复的事件。见Therneau的excellent vignette。我想说明的是，个人不受事件发生后30天重复的事件，也就是我想单独退出的风险暂时设置，例如，如果当个人是不是在事件发生风险，它被忽略。原始方法是迭代地添加所有事件，然后简单地将30添加到tstart变量。但是，这可能导致实例为tstar