单变量coxph函数在多变量中的应用
单变量coxph函数是生存分析中的一种统计方法，用于估计因变量（通常是生存时间）与一个自变量之间的关系。它基于Cox比例风险模型，可以用来分析生存数据中的风险因素。
在多变量分析中，单变量coxph函数可以被扩展为多变量coxph函数，用于同时考虑多个自变量对生存时间的影响。多变量coxph函数可以通过引入多个自变量来建立一个更加全面的模型，以更准确地评估各个自变量对生存时间的影响。
多变量coxph函数的应用场景包括但不限于以下几个方面：
- 生存分析：用于评估不同因素对生存时间的影响，例如医学研究中评估不同治疗方法对患者生存时间的影响。
- 风险评估：用于评估不同因素对风险的影响，例如保险行业中评估不同因素对保险索赔风险的影响。
- 事件发生分析：用于评估不同因素对事件发生的影响，例如市场研究中评估不同因素对产品销售的影响。
腾讯云提供了一系列与生存分析相关的产品和服务，包括但不限于：
- 腾讯云机器学习平台（https://cloud.tencent.com/product/tiia）：提供了强大的机器学习算法和工具，可以用于构建生存分析模型。
- 腾讯云大数据平台（https://cloud.tencent.com/product/emr）：提供了分布式计算和存储能力，可以用于处理大规模的生存数据。
- 腾讯云人工智能平台（https://cloud.tencent.com/product/ai）：提供了丰富的人工智能算法和工具，可以用于辅助生存分析的模型构建和预测。
以上是对单变量coxph函数在多变量中的应用的简要介绍，希望能对您有所帮助。
相关·内容
我试图用coxph安全地使用purrr：：，以便捕获错误消息。我制作了一个安全的coxph版本，如下所示
library(survival)
library(purrr)
coxph_safe <- safely(coxph)
当我唯一的输入是公式和数据时，这是非常有效的，但是，如果我添加了另一个输入，例如子集或权重，我会得到以下错误消息：
simpleError (替换(子集)，数据，env)：..3在不正确的上下文中使用，不.向内看
当需要额外的输入时，有人知道如何安全地应用于coxph吗？我也得到了同样的错误使用安静，而不是安全，如果我做了一个安全版本的lm，并指定一个子集。我使
浏览 0提问于2019-07-31得票数 2
回答已采纳
当将ggforest()应用于coxph对象时，我得到以下错误消息 error in ggforest(res.cox3, data = Selection_cox) :
class(model) == "coxph" are not all TRUE res.cox3是coxph()的输出，其中包含tt术语strata，并且是类： > class(res.cox3)
[1] "coxph.penal" "coxph" 对于下面的虚拟数据，我得到了相同的消息： set.seed(132456)
'dummy survival
浏览 140提问于2020-01-24得票数 1
1回答
我试图用一个分类变量来计算危险比，但用一个参考因子来定义它们中的任何一个都没有意义。我想要进行的比较是每组单独与其他样本进行比较(即标记为A的样本与未标记为A的所有样本；标记为B的样本与未标记为B的所有样本)。默认情况下，coxph会将每个组与参考组进行比较(即标记为B的样本与标记为A的样本；标记为C的样本与标记为A的样本)。但是，如果我尝试将它们编码为伪变量，以这种方式强制它，最后一个总是给出NA结果。有没有办法做我正在尝试的事情，或者有没有办法修复我目前的方法？
示例：
coxph(formula = Surv(Time, Status) ~ class_A + class_B +
浏览 2提问于2021-04-23得票数 0
1回答
如何使用函数在R中运行一组模型
、、、、
library(survival)
justices <- read.csv("http://data.princeton.edu/pop509/justices2.csv")
PREDS = c("age", "year")
m = coxph(Surv(tenure, event == 1) ~ age + year, data = justices)
summary(m)
exp(coef(m)[1])
exp(confint(m,level=(1-0.05/1))[1,])
DVMOD <- function(PREDS
浏览 11提问于2020-05-16得票数 1
回答已采纳
1回答
我对如何调整coxph函数中的变量感到非常困惑。我知道可以用strata()进行分层，但是如何调整变量呢？
在线性模型中，可以通过以下方式调整变量
(这是我在这个中看到的一个示例)
fit.diamOnMachine <- lm(diameter˜machine)
diam.adjusted <- residuals(fit.diamOnMachine)
fit.diamadjmach <- lm(strength ˜ diam.adjusted + machine)
在coxph中，人们会做同样的事情吗？例如，我有1000个人的年龄、性别、BMI、到达事件的时间、结果和蛋白
浏览 0提问于2017-10-09得票数 1
2回答
我想要创建一个用户定义的函数来在一个变量列表上运行coxph。这个公式运行得很好：
summary(coxph(
Surv(var1, factor(var1)) ~ var3,
data = df_wts))
但是，如果我试图将它转换成函数，我就会得到一个“对象未找到”的警告：
cox_fun <- function(x, y) {
summary(coxph(
Surv(x, factor(y)) ~ var3,
data = df_wts))
已经试过了所有的东西，如(X)，{{}，[[]]，！，as.name()等等。
为什么会发生这种事？环境？谢谢!
浏览 4提问于2022-04-29得票数 0
回答已采纳
在R中，当比例检验(带coxph)表明违背Cox模型中的比例假设时，合并协变量与时间之间的相互作用项的最佳方法是什么？我知道你既可以用地层，也可以用时间项进行交互，我对后者很感兴趣。我还没有找到一个明确的解释，并举例说明如何在互联网上这样做。在最常见的使用Rossi数据集的例子中，Fox建议这样做，
coxph(formula = Surv(start, stop, arrest.time) ~ fin + age + age:stop + prio, data = Rossi.2)
建模和年龄之间是否有区别:停止和年龄:开始？公式必须使用这种格式吗？如果我使用两个参数格式的Surv，下面的内
浏览 1提问于2017-08-24得票数 6
1回答
我与R中的coxme包做了斗争，我想使用一个函数，比如survfit() --它通常用于coxph()模型--来绘制调整后的生存曲线，并在不同的参数值下找到中位生存期。
如果我使用coxph (无随机效应)来拟合模型，我可以做以下工作：
library(KMsurv)
data(burn)
my.surv <- with(burn, Surv(T1, D1))
cox_nr = coxph(my.surv ~ Z1 , data = burn)
survfit(cox_nr, newdata = data.frame(Z1 =1))
这提供了生存估计。但如果我把同样的模型和臀部搭配
浏览 2提问于2012-10-31得票数 10
我有一个由多个计算出来的类mids数据集。我使用with()函数来用coxph()函数估计m个不同的数据集。但是，我在自己的函数中使用with()函数时遇到了问题。
下面的代码是一个简单的示例，它在Surv(enter、exit、event)对象“enter”未找到时复制了错误: Error
list<-"X1+X2"
var.used<-formula(paste("Surv(enter,exit,event)~",list,sep=""))
with.coxph<-function(form,dataset){
with
浏览 3提问于2014-12-04得票数 0
1回答
下面的示例是通过添加“年龄”从文档示例中修改的
d.coxph <- (survfit(Surv(time, status) ~ sex+age, data = lung))
autoplot(d.coxph)
我将得到以下错误：
levels<-(*tmp*，value = if (nl == nL) as.character(labels) level paste0(标签)：因子级别41被复制)中的错误
输入一个帧号，或0以退出1:自动绘图(d.coxph)>levels<-中的错误(*tmp*，value = if (nl == nL) as.character
浏览 0提问于2018-05-01得票数 1
回答已采纳
1回答
有没有办法对多个输入运行coxph。在这里，我使用了输入hsa_let_7b_5p。
coxph(Surv(Time, Status)~ hsa_let_7b_5p, data=as.data.frame(test))
Call:
coxph(formula = Surv(Time, Status) ~ hsa_let_7b_5p, data = as.data.frame(test))
coef exp(coef) se(coef) z p
hsa_let_7b_5p 0.169 1.184 0.173
浏览 1提问于2016-05-04得票数 1
1回答
library(survival)
library(survminer)
library(dplyr)
ovarian=ovarian
ovarian$weighting = sample(1:100,26,replace=T)
fitWEIGHT <- coxph(Surv(futime, fustat) ~ age + rx,data=ovarian,weight=weighting)
fitNOWEIGHT <- coxph(Surv(futime, fustat) ~ age + rx,data=ovarian)
在上面的例子中，fitWEIGHT的R平方等于1。然而
浏览 0提问于2019-03-04得票数 0
我想在我自己的R包中使用R包survival。
我的问题是predict.coxph函数不是由survival导出的，这就是为什么我不能在DESCRIPTION文件中使用survival::predict.coxph在Imports:下列出生存。
现在，我试图将survival列在Depends:下，只使用predict(fittedmodel) (以及没有survival::的coxph和Surv )。但是我使用NOTEs获得了下面的R CMD check --as-cran
"*检查R代码中的依赖项.请注意依赖字段中未导入的包：‘生存’这些包需要从(命名空间文件中)导入，以便在
浏览 1提问于2021-03-26得票数 1
1回答
当使用变量交互作用时，我试图通过绘制cox回归来绘制调整后的生存曲线。
读取survfit.coxph页面https://stat.ethz.ch/R-manual/R-devel/library/survival/html/survfit.coxph.html时，我看到了参数"newdata“
newdata
a data frame with the same variable names as those that appear in the coxph formula. It is also valid to use a vector, if the data frame w
浏览 1提问于2015-12-01得票数 2
2回答
我想在函数中使用survfit()和basehaz()，但它们不起作用。你能看一下这个问题吗？谢谢你的帮助。以下代码导致错误：
library(survival)
n <- 50 # total sample size
nclust <- 5 # number of clusters
clusters <- rep(1:nclust,each=n/nclust)
beta0 <- c(1,2)
set.seed(13)
#generate phmm data set
Z <- cbind(Z1=sample(0:1,n,replace=TRUE),
浏览 3提问于2012-01-25得票数 4
我正在运行一个cox比例风险模型，该模型重复超过20,000次(给定不同基因组合的计算结果)。当我运行几百或数千次时，我没有得到任何错误。然而，根据绝对的概率，一些基因最终(我假设)完全分离，这阻止了收敛并产生了这样的警告。 Warning message:
In coxph.fit(X, Y, istrat, offset, init, control, weights = weights, ... :
Loglik converged before variable 1 ; coefficient may be infinite. 这会终止该函数。因此，我想做一个tryCatc
浏览 1879提问于2021-09-14得票数 0
1回答
我想知道是否有人在R中实现了带有Firth校正的cox回归的“预测”函数(包coxphf)。
使用coxphf的预测是否会以与coxph相同的方式进行计算？
谢谢!
下面是一个示例：
library(coxphf)
library(survival)
train <- data.frame(time=c(4,3,1,1,2,2,3),
status=c(1,1,1,0,1,1,0),
x=c(0,2,1,1,1,0,0))
test <- data.frame(x=c(0,2,1,1,1,0,0))
fit.coxph <-
浏览 0提问于2016-03-12得票数 0
我对我试图写的函数的行为感到困惑。我的例子来自于survival包，但我认为问题比这更普遍。基本上，下面的代码
library(survival)
data(bladder) ## this will load "bladder", "bladder1" and "bladder2"
mod_init <- coxph(Surv(start, stop, event) ~ rx + number, data = bladder2, method = "breslow")
survfit(mod_init)
会产生我感兴趣
浏览 0提问于2016-05-21得票数 5
回答已采纳
1回答
我想探索以下变量及其双向交互作用作为可能的预测因子:兄弟姐妹数量(nsibs)、断奶年龄(wmonth)、母亲年龄(mthage)、种族、贫困、出生体重(bweight)和母亲吸烟(smoke)。
我创建了我的Cox回归公式，但我不知道如何与预测因子形成双向交互：
coxph(Surv(wmonth,chldage1)~as.factor(nsibs)+mthage+race+poverty+bweight+smoke,data=pneumon)
final<-step(coxph(Surv(wmonth,chldage1)~(as.factor(nsibs)+mthage+race
浏览 1提问于2017-05-18得票数 0
1回答
如何用在R中实现的校正组预测方法绘制连续协变量(假设数值的第20和第80百分位数)的预测生存曲线？
例如,
library(survival)
library(survminer)
fit <- coxph( Surv(stop, event) ~ size + strata(rx), data = bladder )
ggadjustedcurves(fit, data=bladder, method = "conditional", strata=rx)
现在，这是有用的，因为我得到了两条生存曲线，它们是由rx (0或1)分层的，条件方法是对膀胱数据集的作用。然而
浏览 0提问于2019-09-02得票数 2
1回答
我正在建立一系列Cox回归模型，并根据新的数据从这些模型中得到预测。在某些情况下，我能够得到预期的事件数量，但在其他情况下则无法得到预期的事件数量。
例如，如果写出coxph调用中的公式，则计算预测。但是，如果公式存储在一个对象中，并且调用了该对象，我将得到一个错误。如果我试图在dplyr管道突变函数中创建预测(对于我正在编写的函数，这将是使预测正常工作的最理想的地方)，我也无法得到预测。
任何帮助都是非常感谢的！
谢谢,
丹尼尔
require(survival)
require(tidyverse)
n = 15
# creating tibble of tibbles.
results
浏览 0提问于2017-11-15得票数 0
回答已采纳
1回答
我试图在coxph包下使用survival函数。通常，它将被称为：
coxph(Surv(time,event) ~ age+gender+salary, data=THEDATA)。
但是，我在THEDATA中有多个列。我怎么能轻易给他们打电话呢？例如，我想在第8- 12列上建立基于变量的Cox模型。
coxph(Surv(THEDATA$time,THEDATA$event）~ THEDATA[,8] + THEDATA[,9] + THEDATA[,10] + THEDATA[,11] + THEDATA[,12])
如何使用更有效的代码来完成这项工作？
浏览 1提问于2014-08-14得票数 0
1回答
我正在尝试用ipython中的rpy2与R一起执行生存分析。
以下公式失败：
%Rpush df
surv_fit = surv.coxph("Surv(start, end, status) ~ x", data=df)
> RRuntimeError: Error: $ operator is invalid for atomic vectors
但是，如果我使用r魔法完成所有这些工作，它将会工作：
%Rpush df
%R surv.fit <- coxph(Surv(start, end, status) ~ x, data=df)
Surv对象有问题吗？
浏览 0提问于2014-11-01得票数 0
我正在做一个使用coxph()函数的多点触摸属性问题。它是一个包含大约100万个数据的大型数据集，但目前我运行的是大约100000个数据的子集。我已经从我的数据中删除了所有缺少的值。我收到一个错误
Error in if (any(infs)) warning(paste("Loglik converged before variable ",
:missing value where TRUE/FALSE needed
下面是Cox函数：
SurvObj <- Surv(Final_Data$NormalizedStartTime,Final_Data$Normal
浏览 0提问于2016-10-19得票数 2
我在竞争风险上下文中使用Cox比例风险(survival::coxph)模型--即多个事件类型，每个观察都有一个端点。我很难使用coxph.predict函数来显示给定一组协变量和后续时间的预期事件数。下面是一个在生存包中使用mgus2数据集的示例：
library(survival)
#Modify data so each subject transitions only once to a state.
crdata <- mgus2
crdata$etime <- pmin(crdata$ptime, crdata$futime)
crdata$event <- i
浏览 5提问于2022-06-14得票数 1
当我对survfit对象使用summary()时，我得到了一个很好的data.frame，其中包含
time n.risk n.event survival std.err lower 95% CI upper 95% CI
但是coxph对象上的summary()给出了不同的结果。到目前为止，我还没有在文档中找到答案。
谢谢!汤姆
浏览 0提问于2012-10-26得票数 0
回答已采纳
我想用Cox比例风险模型的coxph R函数从惩罚的R包中得到对肺数据集的预测。
我有下面的示例代码。
library(survival)
library(pec)
library(penalized)
data("lung")
data <- lung
trainind <- sample(1:n,n*0.7)
testind <- (1:n)[-trainind]
frm <- as.formula(paste("Surv(time, status)~",paste(names(data[,-c(2,3)]), col
浏览 7提问于2017-09-11得票数 0
回答已采纳
我在R的coxph()中使用了以下函数来拟合cox风险模型。我想报告正确的统计数据；但是，输出中没有95%的CI。
Surv(days, censor) ~ gender + age + treatment, data_1)
我只得到了以下几列。
coef exp(coef) se(coef) z p
浏览 24提问于2016-04-29得票数 3
回答已采纳
我目前正在努力解决R方面的以下问题：
我想用一个带有时变变量的Cox模型来通过提升来执行我的变量选择。iCoxBoost接受一个Surv对象作为公式，它允许开始停止标记.但是，当我使用它时，如下所示
iCoxBoost(Surv(tstart, tstop, status) ~ ., data=data)
我收到以下错误消息：
‘[.default’(响应，"time")中的.default错误:下标超出界限
另一方面，coxph(Surv(tstart, tstop, status) ~ ., data=data)和iCoxBoost(Surv(time, status)
浏览 4提问于2022-03-04得票数 0
我运行了一个coxph模型和一个脆弱模型，但现在我想更改连续变量(age)的危险比，以5个单位的增量表示，而不是1个单位。在R中有没有一个函数可以执行这样的任务？如果是这样，该函数是否也适用于脆弱模式？我使用了frailtypack包。
library('survival')
data(veteran)
cox <- coxph(Surv(time, status) ~ age, data = veteran)
summary(cox)
# Call:
# coxph(formula = Surv(time, status) ~ age, data = veter
浏览 0提问于2018-08-18得票数 2
我有一个模型对象，它显示硬盘中的大小比R中的大，经过一些搜索后，我设法得到了问题的原因，如下所示
format(object.size(JMFit1$model_info$coxph_components$TermsU), units='Mb')
[1] "0 Mb"
pryr::object_size(JMFit1$model_info$coxph_components$TermsU)
28.5 MB
但是JMFit1$model_info$coxph_components$TermsU返回
>JMFit1$model_info$coxph_compon
浏览 0提问于2018-11-22得票数 1
1回答
在我的生存分析任务中，我使用了cox比例模型来计算我的数据集的不同组中的协调性指数(c- index )值。我想知道如何将p值添加到我的c-index图中，以比较不同的组，使其看起来像这个图？
下面是我的代码：
surv <- with(group, Surv(group$survival, group$time))
# calculate survival
sum.surv_1 <- with(group, summary(coxph(surv ~ group$1)))
sum.surv.1_2 <- with(group, summary(coxph
浏览 0提问于2015-05-17得票数 5
1回答
在R？
、、
我用权重和时间自变量的时间依赖变量拟合了coxph模型，如下所示：
id start stop death weight age smoke
1 0 1 0 60 60 0
2 0 1 0 55 57 1
2 1 2 1 60 58 1
...
mod <- coxph(Surv(start, stop, death) ~ weight + age + factor(smoke) + cluster(id),
data=data)
我想绘制体重
浏览 0提问于2018-09-19得票数 0
运行coxph诊断以生成多变量coxph回归模型的dfBeta图。
require(survival)
require(survminer)
# Dummy data
set.seed(100)
permth_int <- seq(10)
mortstat <- c(rep(1,9),0)
age <- runif(10,50,70)
sex <- sample(seq(2),10,replace=T)
ill <- sample(seq(3),10,replace=T)
# Create survival obj
浏览 0提问于2018-03-05得票数 0
回答已采纳
我想为我的生存分析对象(suvfit/coxph)计算AIC值。当我试着去做的时候，它说：
> AIC(cox)
Error in UseMethod("logLik") :
no applicable method for 'logLik' applied to an object of class "coxph"
据我所知，这是一个软件限制。有人能帮我解决这个问题而不用手工计算AIC值吗？
浏览 3提问于2013-10-30得票数 10
回答已采纳
1回答
我使用pysurvival库来建立Cox比例风险模型。与其得到生存曲线，我更感兴趣的是得到点预测。在库中，函数predict_survival返回一个数组--类似于表示生存函数的预测，我假设它可以用来获得预期的值--但我只是找不到正确的方法。
下面我附了一个假例子。
# Initializing the simulation model
sim = SimulationModel( survival_distribution = 'log-logistic',
risk_type = 'linear',
浏览 6提问于2021-04-16得票数 0
回答已采纳
2回答
我现在在做教练。我的设置:我有一个参考(没有处理)，然后三个不同的处理(A，B，和C)。我也有A，B和C的相互作用，(例如，处理A和B的样品，或A和C，等等.)。我为这些治疗创建了虚拟变量，编码为1或2 (1 =接受治疗，2=没有接受治疗)。我使用as.factor()加载这些变量。
example:
A<-as.factor(Data$A)
我可以这样做，并得到一个结果，表明接受治疗B(又名B= 1)有利于寿命(coef为正)。这三者在某种程度上都具有重要意义：
> coxph1<-coxph(Surv(Lifespan,Status)~A+B+C
> summary
浏览 0提问于2018-06-28得票数 0
回答已采纳
2回答
在一项研究中，评估了一个大型中西部大学医学中心消毒操作规程的改变。这项研究的主要兴趣是比较两种身体清洁方法。first法是1983年1月至1984年6月期间专门使用的方法，包括一种常规的沐浴护理方法(先用10%聚维酮碘去除表面污染，然后用拨号皂定期洗澡)。从1984年6月至December1985研究期结束，最初使用4%葡萄糖酸氯己定进行身体清洁。84例患者接受新洗浴液洗必泰，70例患者作为对照组接受常规沐浴护理。数据集中包括一个协变量，用于测量燃烧的总表面积。数据是(刻录)。我想测试的是: 1-两组患者的生存功能是否有差异。2.两组的生存功能是否有差异，可根据烧伤总面积进行调整。
librar
浏览 2提问于2020-04-01得票数 1
回答已采纳
使用森林模型软件包的林地确实是多变量cox回归的好方法。但是，我在用原始变量和数据因素的名称替换最终表示标签(即变量:Age 2表示“年龄”；因子:0表示"<60"，1表示"≥60")时遇到了问题，≥对R编码的知识有限，但我尝试了expss包向变量和因素添加标签。然而，coxph()并不适用于标签，而是处理因素。
这是我的编码：
Cox比例模型：
mcox<-coxph(pblsurv~age2+sex1+origin,data = pbl)
使用森林模型包的森林地块：
print(forest_model(mcox))
最后一个图中的变量名是age
浏览 0提问于2019-06-15得票数 0
我正在使用生存软件包做一个cox回归(coxph)，我无法计算出我做错了什么。
我的代码是
prentice.fit <- coxph(Surv(age_ent, outage, diabetes) ~ a_menopause + cluster(id), data="mydata.cc", ties="breslow")
summary(prentice.fit)
但是我不断地得到一个错误:eval中的错误(predvars，data，env)：无效的‘env’类型‘字符’参数
我检查了所有变量的类别。age_ent和中断是数字，a_menopause和
浏览 2提问于2015-11-27得票数 0
回答已采纳
2回答
在 answer中，它说调试是默认的cmake构建配置。
但我有一个不同的观察结果：
在我的CMakeLists.txt中，我根据当前的构建配置选择库的调试和发布版本。
target_link_libraries(MyApp debug Widgets_d)
target_link_libraries(MyApp optimized Widgets)
似乎当我在没有指定-DCMAKE_BUILD_TYPE标志的情况下调用cmake时，使用的是Widgets而不是Widgets_d (当我删除Widget并尝试构建时，会抱怨lib不在那里)。因此，这意味着默认情况下构建配置是优化的，而不是调试的
浏览 0提问于2014-11-23得票数 8
回答已采纳
set<-c(1,1,1,2,2,2)
gender<-c(1,0,1,0,1,0)
smoke<-c(1,1,0,0, 1,0)
case_control<-c(1,0,0,1,0,0)
data<-data.frame(set, gender, smoke, case_control)
data$gender<-factor(data$gender, levels=c(0,1), labels=c("female", "male"))
data$smoke<-factor(data$smoke,
浏览 2提问于2020-03-14得票数 1
首先，我从这个链接中了解到，使用" function“函数可能会做我想要的事情。然而，我还没能从以现在的方式去思考它，到让它在我的情况下真正发挥作用(甚至真正知道从哪里开始)。我是R的初学者，所以如果这是一个真正的“新”问题，我会提前道歉。我的数据是由事件变量(肿瘤复发)和时间变量(随访时间/复发时间)以及复发危险因素(t期、肿瘤大小、dx年龄等)组成的数据框架。一些危险因素是明确的，有些是连续的。我一直在手工进行单变量分析，就像这个例子univariateageatdx<-coxph(survobj~agedx)一样，然后收集数据。对于多个因素来说，这变得非常繁琐，并且对于几个
浏览 0提问于2015-09-24得票数 2
回答已采纳
在this tutorial之后，我正在尝试将R中的多变量Cox回归分析应用于我的数据集。特别是，我正在尝试应用以下函数coxph() install.packages(c("survival", "survminer"));
library("survival");
library("survminer");
data("lung");
res.cox <- coxph(Surv(time, status) ~ age + sex + ph.ecog, data = lung)
summary(res
浏览 49提问于2020-10-08得票数 1
回答已采纳
1回答
我想在考克斯生存模型之后创建一个森林地块。但是，我想只显示图中的一些协变量吗？有人知道这是否可能吗？也许用ggforest2？谢谢
library(survival)
library(survminer)
model <- coxph(Surv(time, status) ~ sex + rx + adhere,
data = colon )
ggforest(model)
colon <- within(colon, {
sex <- factor(sex, labels = c("female", "male&
浏览 0提问于2019-10-11得票数 0
1回答
一些日子以来，我一直在教自己R，我一直在做cox回归分析。
我使用cut()函数将每个连续变量划分为两个分类组。现在，我想知道是否可以将这些类别分别合并到cox回归分析中。
例如:我有3个变量(A，B，C)，每个变量由2类组成: A-，A+ & B-，B+ & C-，C+。现在，如果我使用命名变量来运行coxph()，我只会得到"+“类别的结果(据我所知，这是因为"-”类别被用作引用组)。然而，由于C-似乎对生存有负面影响，我更感兴趣的是收到该类别的结果。
另外，我想知道是否有一种方法，我可以将每个类别定义为一个新的组/变量，并将它们单独组合起来，以了解每个类别
浏览 4提问于2017-04-21得票数 0
回答已采纳
当我尝试在R中创建一个类，包含一个应该接受glm类对象或coxph类对象的插槽时，在包构建过程中我得到了一个错误，那就是Error in setClassUnion(name = "glmORcoxph", members = c("glm", "coxph")) : the member classes must be defined: not true of "coxph"。我试着把@import，@importClassesFrom，@importFrom Roxygen指令放到我的源代码中，但是没有成功。这是我用来创建Do
浏览 3提问于2018-01-26得票数 0
回答已采纳
我有一个数据集，其中有3个组接受不同媒体的曝光。一组被暴露在三种媒体中的一种。因此，我的coxph模型是分层的：
# My treatment variable is loaded as a factor.
fullModel <- coxph(Surv(time, status) ~ strata(treatment), data = d)
当我尝试做模型诊断时，我会得到以下错误：
test.assump <- cox.zph(fullModel)
Error in cox.zph(fullModel) :
there are no score residuals for a
浏览 5提问于2022-06-24得票数 1
回答已采纳
2回答
我使用sklearn-learn在python中构建了一个随机森林，并对所有分类变量应用了“一热”编码。问:如果我将"one hot“应用到我的DV，我是将它的所有虚拟列都作为DV应用，还是应该以不同的方式处理DV？
浏览 4提问于2018-12-03得票数 3
1回答
我想使用R包smcfcs和命令smcfcs。
CRAN:
例子如下。
imps <- smcfcs(ex_linquad, smtype="lm", smformula="y~z+x+xsq",
method=c("","","norm","x^2",""))
imps <- smcfcs(ex_coxquad, smtype="coxph", smformula="Surv(t,d)~z+x+xsq",
method=c("&#
浏览 1提问于2016-10-25得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
