对不同的数据帧重复使用Surv函数
是指在统计学中，使用Surv函数对不同的数据帧进行重复使用。Surv函数是生存分析中的一个重要函数，用于计算生存时间和事件发生的状态。它通常用于分析生存数据，如医学研究中的患者生存时间、产品的寿命等。
Surv函数的输入通常是两个参数，第一个参数是生存时间（或观察时间），第二个参数是事件发生的状态（如死亡、失效等）。Surv函数的输出是一个Surv对象，可以用于构建生存分析模型，如Kaplan-Meier曲线、Cox比例风险模型等。
对于不同的数据帧，重复使用Surv函数可以进行多种生存分析的比较和统计。通过对不同数据帧的生存时间和事件状态进行Surv函数的重复使用，可以比较不同组之间的生存曲线、计算生存率、估计风险比等。
在云计算领域中，对不同的数据帧重复使用Surv函数可以用于分析云服务的可靠性和可用性。通过收集云服务的运行时间和故障事件，可以使用Surv函数对不同的数据帧进行生存分析，评估云服务的可靠性，并提供相应的改进措施。
腾讯云提供了多个与云计算相关的产品，如云服务器、云数据库、云存储等。这些产品可以帮助用户构建和管理云计算环境，提供稳定可靠的云服务。具体的产品介绍和链接地址可以参考腾讯云官方网站：https://cloud.tencent.com/
相关·内容
经过一些挖掘(1，2，3)，如果我理解正确的话，似乎有一些关于在函数中使用公式导致作用域问题的帖子。有些人建议使用环境、赋值或<<-来解决这个问题，但我被如何使用它们弄糊涂了(也不明白为什么首先会有问题)。 让我们试试这个玩具代码： library(survival)
library(survminer)
set.seed(1)
give_p_val <- function() {
df <- data.frame('OS' = ovarian[, 'futime'], 'Survival_event' = ova
浏览 33提问于2021-09-15得票数 1
回答已采纳
我正在运行一些Surv()函数，有一件事我不喜欢，也不理解，那就是为什么这个函数不接受"data=“参数。这很烦人，因为我想在相同的数据帧上执行相同的Surv()函数，但每次都按不同的标准进行过滤。
例如，我的数据框名为"ikt“，我希望通过"plan _TYPE2==‘LD’”进行过滤，并使用一个strata变量“PLAN2”。我尝试了以下方法，但不起作用：
library(survival)
library(dplyr)
ikt<-data.frame(organ_yrs=(seq(1,20)),
organ_status=rep(c(
浏览 0提问于2018-08-02得票数 0
我有一个映射的森林地块的数据框架，其中所有的树干都有X，Y坐标，胸径(厘米)和存活(0,1) (在这里找到，名为"MFP14_surv_forSO.csv"：)。
我创建了一个空间点模式，然后是一个平滑的大茎死亡内核(im对象)，作为树冠“空隙”(或光线穿透量)的代表。我认为大茎是直径>9厘米的茎，而'0‘成活表明茎死了。
我想在所有的树干上插入这个“空隙”的度量。我已经知道模式中有几个点位于窗口之外，所以我将这些点排除在使用inside.owin()对数据框架进行子集的分析中。
surv14 <- read.csv("MFP14_surv_forS
浏览 3提问于2016-11-14得票数 0
回答已采纳
要为两条生存曲线之间的差异生成p值，我需要执行一系列乏味的步骤，这些步骤需要与起始数据帧的一列中表示的三个值中的两个值的所有组合重复，比方说组A、B和C： # Load packages
library(survival)
library(survminer)
# Generate practice dataset
surv_time <- sample(100:200, size = 50, replace = TRUE)
surv_group <- sample(c("A", "B", "C"), size = 50, rep
浏览 70提问于2021-11-07得票数 1
回答已采纳
我有两个独立的数据集，一个是5421个，另一个是1000个。我想要做的是使用外部数据集(main_dat，n=1000)验证从主数据集(n=5421，test_dat)获得的Cox模型。但是，我在Hmisc包中使用rcorr.cens()并在rms中使用val.surv获得了一个错误消息。以下是我一直在做的事情：
library(rms)
surv.obj=with(main_dat,Surv(survival,surv_cens)) ## to use with rcorr.cens
phmodel=cph(surv.obj~sex+age+treatment, x=TRUE, y=TRU
浏览 5提问于2015-02-06得票数 0
我在surv_rate中有必要列的数据。这就是我正在运行的，它基于以下示例：
fit2 <- survfit(Surv(YearsToStatus, VitalStatus) ~ PR_statusIndex + HER2_statusIndex + ER_statusIndex,
data = surv_rate)
ggsurv <- ggsurvplot(fit2, data = surv_rate, ggtheme = theme_bw)
我得到了这个错误：
Error: Don't know how to add ggtheme t
浏览 0提问于2018-07-12得票数 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
回答已采纳
2回答
我正在尝试对6500条记录的数据帧运行RFSRC，其中包含59个变量：
rfsrc_test <- rfsrc(Surv(TIME, DIED) ~ ., data=test, nsplit=10, na.action = "na.impute")
当我在1500条记录上运行它时，它似乎可以工作，但在整个数据集上崩溃。
它在没有任何特定错误的情况下崩溃了--有时它会给出“异常处理错误”。
你有什么想法怎么调试这个吗？我在数据库中浏览了一下奇怪的行，但一无所获。
浏览 0提问于2018-11-29得票数 2
1回答
我知道如何绘制生存曲线(由survival软件包)由autoplot从ggfortify。这里我使用ovarian数据帧，也来自生存包。
library(survival)
library(ggfortify)
autoplot(survfit(Surv(ovarian$futime, ovarian$fustat) ~ ovarian$resid.ds))
可以只提供列名的survfit函数，并分别指定数据框架，如下所示：
autoplot(survfit(Surv(futime, fustat) ~ resid.ds, data = ovarian))
这个很好用。我想要做的是在开始时使
浏览 1提问于2018-04-14得票数 0
回答已采纳
我有一个生存数据，但我不确定psm和cph之间有什么不同。我如何确定模型？
不同的模型将建立不同的名义图，但我不确定我应该使用哪种模型？
library(rms)
f2 <- psm(Surv(follow_time_5y, DEATH_5y) ~ age+ID_SEX+MH_CCI_total_score,
data =sci_20190505, dist='lognormal')
f2 <- cph(Surv(follow_time_5y, DEATH_5y) ~ age+ID_SEX+MH_CCI_total_score,
浏览 3提问于2019-05-16得票数 1
回答已采纳
在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回答
当我在顶层对象上分派时，所有函数都导出到名称空间。一切都如我所料。
myfun <- function(obj) {
UseMethod("myfun",obj)
}
#' @export
myfun.this <- function()
命名空间
S3Method(myfun,this)
export(myfun)
但是，当我试图对我的对象中的一个元素进行调度时.
#' Subsets recoded data by a condition
#'
#' This function filters a list, data
浏览 1提问于2016-02-17得票数 0
由于某些原因，我无法将我的图例的标题从"surv_status“(变量名)更新为"Survival Status”，尽管在代码块中所有不同的可能位置都添加了+labs(fill = "Survival status")。
有什么建议吗？谢谢!
ggplot(data_train, aes(x = axil_nodes, fill = as.factor(surv_status))) +
geom_bar(position = "fill") +
labs(fill = "Survival status") +
the
浏览 0提问于2021-03-24得票数 0
1回答
假设我有以下数据帧
library(survival)
library(multcomp)
data(cml)
cml$group<-sample(1:5, 507, replace=T)
plot(survfit(Surv(time=cml$time, cml$status)~cml$group))
(survdiff(Surv(time=cml$time, cml$status)~cml$group))
如何执行多个比较测试，例如group0与所有其他组的比较？还是每一组人都在一起？
有没有一种很好的方法来绘制这些多重比较(例如在aov()中的plot.TukeyHSD()中
浏览 2提问于2012-06-24得票数 1
1回答
我一直在准备生存分析和考克斯回归的R。然而，我的线经理是一个Stata用户，并希望输出显示方式类似于Stata会显示它。
# Stata code
. strate
. stsum, by (GROUP)
stsum会输出每组的风险时间和发病率，我不知道如何用R实现这一点。
数据看起来大致是这样的(由于它是在安全的环境中，所以我无法了解它)：
PERS GROUP INJURY FOLLOWUP
111 1 0 2190
222 2 1 45
333 1 1 560
444 2 0 1200
到目前为止，我一直在
浏览 3提问于2015-01-22得票数 0
回答已采纳
2回答
我正在尝试编写一个函数，它只需要几行代码，并允许我输入一个变量。下面的代码使用Surv函数(Survival包)创建一个对象。第二行接受有问题的变量，在本例中是列为Variable_X的列，并输出数据，然后可以使用ggsurvplot可视化这些数据。输出是Kaplan-Meier生存曲线。我想做的是有一个函数，这样我就可以输入f(Variable_X)，并对我从数据中选择的任何一列的输出KM曲线进行可视化。我希望f(y)输出KM，就像我把y放在~Variable_X当前所在的位置一样。我是R的新手，对函数的工作方式也很陌生，我已经尝试过下面的代码，但它显然不能工作。我正在通过datacamp和
浏览 3提问于2018-06-29得票数 0
1回答
我正在尝试对一种特定癌症类型中的数百个基因进行生存分析。我有两个数据帧(m2和m3)。m2包括样本ID以及总体存活率(样本已经存活了多长时间)和状态(如果样本是活着还是死了)的列。在m3中，我有一个列表示样本ID，第2:256列是不同的基因。如果每个样本都有该基因的突变，则表示为1，如果没有，则表示为0。我正在尝试确定哪些基因在比较它们在生存中的作用时具有统计学意义。我正在尝试运行一个for循环来运行这个survdiff函数并生成p值，但是总是得到一个错误。
for (x in 2:ncol(m3)) {survdiff(Surv(m2$Overall.Survival, m2$Status)
浏览 0提问于2018-02-11得票数 0
我有唯一识别的个人的数据，在群体中，这些数据被多次调查。有时在第一次调查后会出现新的个体，我想按组找到这些新的个体。例如， group <- rep(1:2, each=9)
surv <- c(1,1,2,2,2,3,3,3,3,1,1,1,2,2,2,2,3,3)
ind <- c("a","b","a","b","c","a","b","c","d", "a","b","c"
浏览 8提问于2021-02-11得票数 0
回答已采纳
我有一些极端的异常值抛弃了我的回归模型，我使用If-Then-Else语句删除了它们。然而，SAS完全消除了这些数据点，并在剩余的数据点中发现了新的异常值。有没有一种方法可以从分析中删除异常值，而不会在混合中投入更多？
我计算了Q3 + 1.5 * IQR，并按如下方式使用该值：
Data lungcancer; input trt surv age sex @@;
/* create a new variable diff */
diff = surv - 365;
/* create a new categorical variable resp */
If diff > 0 th
浏览 3提问于2016-12-03得票数 0
1回答
如何组合包含survival::Surv对象的两个数据文件，以便生成的数据帧中的字段与源数据帧具有相同的类？
我发现使用rbind会导致Surv对象被转换为矩阵。例如，我按照以下方式创建了df1：
library(survival)
df1 <- data.frame(obs = c('A','B','C','D','E')
, lo = c(10,20,30,40,50)
, hi = c(30,30,30,40,50))
df1$co
浏览 3提问于2020-07-11得票数 1
回答已采纳
3回答
建立了一个应用Cox回归模型对数据进行检验的函数，建立了基于协变量的生存函数，并对每次观测的生存概率进行了预测。
下面的示例使用了肺数据集，并且运行得很好。然而，应用到我自己的数据处理时间是冗长的。对于n= 60000，我只是在一个小时后停止了它，因为它对我打算使用这个程序是不实际的。
从代码结构来看，有明显的方法可以加快速度吗？
require(dplyr, survival, pec)
cox_model <- coxph(Surv(time, status) ~ sex, data = lung)
surv_preds <- function(model, query)
浏览 4提问于2017-11-12得票数 1
回答已采纳
1回答
我很难使用coxph进行预测。
我希望学习样本内数据的Cox PH模型，然后使用样本外数据的参数，如下所示：
# learn IS params
model.PH <- coxph(Surv(days.IS, outcome.IS) ~ predictor.IS)
# apply IS params to OOS data to make predictions
predictions.raw <- predict(model.PH, newdata = predictor.OS)
# binarise predictions
predictions.OS <- rep(
浏览 0提问于2014-09-02得票数 0
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
我用R写了下面的函数
Plots <- function (database, variable) {
fit <- survfit(Surv(database$Months) ~ variable, data = database)
ggsurvplot(fit, data = database, pval = TRUE)
surv_pvalue(fit, data = database)
}
当我像这样调用它时(DB，DB$variable)，我得到了以下错误：
Error in eval(inp, data, env) : object 'data
浏览 2提问于2019-11-22得票数 0
我有以下两个数据帧：
df1：
datetime actual forecast previous
0 2015-01-09 13:30:00 -0.2 0.2 0.2
1 2015-02-06 13:30:00 0.5 0.3 -0.2
2 2015-03-06 13:30:00 0.1 0.2 0.5
3 2015-04-03 12:30:00 0.3 0.2 0.1
浏览 0提问于2016-12-26得票数 0
我的目标是使用survfit对象，以90%的置信度水平估计中位数的上限和下限的中位数生存率。 churn_dat <-read_csv("https://raw.githubusercontent.com/square/pysurvival/master/pysurvival/datasets/churn.csv")
churn_dat <- churn_dat %>% filter(months_active > 0) #create a function of the dataframe by sizes
boot <- function(
浏览 71提问于2021-09-07得票数 1
我正在尝试为我需要多次执行的任务编写一个函数(在多个推定数据集上运行cox比例风险函数)。但是，当我将必要的对象传递给我的用户定义函数时，它会给出一个错误，指出找不到该对象。我认为这是因为对象是在cch()函数中使用"data=“参数指定的数据帧中定义的。有人能帮我吗？
示例数据：
my.list<-list(my.df1 <- data.frame(my.id = 1:100, my.time = rlnorm(100),
my.event= c(rbinom(50,1,0.2),rep(1,50)), my.det=rbinom(100,1,0.5),
sub= c
浏览 0提问于2014-03-25得票数 1
回答已采纳
1回答
R-适用于列表和相应名称的绘图。
、、、、
我有一个名为dataframes的列表：
library(survival)
library(survminer)
surv.days<- runif(n = 50, min = 0, max = 500)
censor<- sample(c(0,1), 50, replace=TRUE)
survdata<- data.frame(surv.days, censor)
survlist<- list(survdata, survdata)
names(survlist)<- c("name1", "name2")
rm(surv
浏览 3提问于2017-11-28得票数 1
回答已采纳
我遗漏了什么？
我有一个值的下拉列表。我认为一切都很好，除了当我将结果发布到另一个表时，我似乎无法从第一个表中获得要发布的gameID，它只是作为空白通过。我是否错过了与团队和gameID的关联？
用于获取下拉列表数据的查询：
$sql = "select s.*, (DATE_ADD(NOW(), INTERVAL " . SERVER_TIMEZONE_OFFSET . " HOUR) > gameTimeEastern or DATE_ADD(NOW(), INTERVAL " . SERVER_TIMEZONE_OFFSET . "
浏览 1提问于2017-10-20得票数 0
1回答
我有一个数据帧df，如下所示： df <- structure(list(times = c(724L, 1624L, 1569L, 2532L, 1271L, 2442L,
757L, 848L, 3675L, 1229L, 1582L, 1257L, 1270L, 555L, 357L, 1133L,
633L), Samples = structure(c(1L, 10L, 11L, 12L, 13L, 14L, 15L,
16L, 17L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L), .Label = c("Sample1",
&#
浏览 11提问于2019-04-08得票数 1
回答已采纳
我用下面的脚本来绘制曲线。
我认为中位生存线是一个很好的工具。然而，中位生存线是画成一条虚线的黑线，这是图形压倒性的。我能改变生存线的颜色或不透明度以减少surv.median.line函数的图形输出吗？
如果没有，我是否可以手动添加一个垂直/水平的中位生存线，在其中我可以改变颜色或不透明度？
j <- ggsurvplot(
fit,
data = p,
#fun="cumhaz",
risk.table = "abs_pct", #risk.table.col="strata"
浏览 0提问于2018-10-21得票数 1
回答已采纳
1回答
在rms中使用cph()生成模型后，调用model$y将返回一个生存对象。有没有一个函数可以“撤销”生存对象并返回一个数据帧？
我希望能够使用生存对象作为我正在编写的函数的参数，但我还需要响应数据。我试图避免使用数据作为参数，并在函数内部创建模型。
下面提供了一个最小的工作示例：
library(rms)
# generate data
time <- c(82, 73, 89, 79, 72, 87, 103, 83, 100, 79)
event <- c(0, 0, 1, 0, 1, 0, 0, 0, 1, 1)
covar <- c(15, 11, 11, 20, 1
浏览 0提问于2016-02-22得票数 2
1回答
我正在尝试对这些数据运行一个t.test (以“测试”男性是否真的比女性喝得更多)。我不完全确定我在做什么，哈哈-我想我已经知道如何做t.test了，但是现在出现了一个错误，说：
"Error in complete.cases(x, y) : not all arguments have the same length".
x是男性，y是女性-这些都是正确的x和y值吗？
我的完整代码：
surv <- read.csv("classsurvey.csv")
library(ggplot2)
ggplot(surv, aes(y=UnitsAlcohol,
浏览 2提问于2017-01-13得票数 1
我有一个由多个计算出来的类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
我有一个很简单的活动，像这样：
public class SurvivalActivity extends Activity {
private static final String KEY_LAYOUT_ID = "SurvivalLayoutId";
int mLayoutId;
private static final int[] mSurvivalLayouts = {
R.layout.survival1,
R.layout.survival2,
R.layout.survival3,
浏览 5提问于2012-03-02得票数 3
回答已采纳
我试图绘制一个生存图，并遇到了一个问题，当我试图将我的越权函数转换为一个主函数时，我可以为不同的数据集调用它。当我运行代码时
fit<- survfit(Surv(time, status) ~ sex, data = lung)
allsurv <- function(fit){
ggsurvplot(
fit,
pval = TRUE,
pval.coord = c(200, 0.10),
conf.int = TRUE,
xlab = "Days",
ggtheme = theme_light(),
surv.median.line = "hv&
浏览 0提问于2018-06-11得票数 4
回答已采纳
我在使用oracle数据库的Asp.net C# 3层工作。我为Read和Delete创建了一个函数，它可以工作，但是对于insert和update，我得到了将字符串转换为日期时间的错误。
以下是业务层代码：
public void InsertSurvey(string SURV_NAME, DateTime SURV_STARTDATE, DateTime S_ENDDATE)
{
try
{
string str = "select max(SURV_ID) as max_SURV_ID from
浏览 0提问于2012-10-01得票数 0
我一直在尝试用应用族和for循环进行批量单变量cox分析，但我不能编辑有用和有效的代码。求求你，请帮帮我！
我准备了两种数据。数据A，"data_train_list“是包含许多数据帧的列表。在我的代码中，每个数据帧中的每个列元素都将作为独立变量。数据B，"os_train_list“也是包括具有数据帧的数据的列表。并且该列表中的每个数据帧都将用作响应变量。
这是我的代码，
for (i in 1:length(data_train_list)) {comb_resBM[[i]]=apply(data_train_list[[i]],2,FUN = function(x){
浏览 1提问于2019-07-18得票数 1
回答已采纳
我试图用coxph安全地使用purrr：：，以便捕获错误消息。我制作了一个安全的coxph版本，如下所示
library(survival)
library(purrr)
coxph_safe <- safely(coxph)
当我唯一的输入是公式和数据时，这是非常有效的，但是，如果我添加了另一个输入，例如子集或权重，我会得到以下错误消息：
simpleError (替换(子集)，数据，env)：..3在不正确的上下文中使用，不.向内看
当需要额外的输入时，有人知道如何安全地应用于coxph吗？我也得到了同样的错误使用安静，而不是安全，如果我做了一个安全版本的lm，并指定一个子集。我使
浏览 0提问于2019-07-31得票数 2
回答已采纳
1回答
我需要从我的模型中提取出随机效应的后验估计值和间隔。
为了说明起见，与我正在使用的数据集类似的数据集将是基R中的ChickWeight数据集。
我提取固定效应的后验估计值和间隔的方法如下：
#load package
library(lme4)
#model
m.surv<-lmer(weight ~ Time + Diet + (1|Chick), data=ChickWeight)
#load packages
library(MCMCglmm)
library(arm)
#set up for fixed effects
sm.surv<-sim(m.surv)
smf
浏览 5提问于2019-11-24得票数 2
回答已采纳
我希望从coxme模型中，在给定的时间步骤中，为我的解释变量的给定值，获得事件概率的预测。我的数据结构如下：
# Generate data
set.seed(123)
mydata <- data.frame(Site = as.factor(sample(c("SiteA", "SiteB", "SiteC"), 1000, replace = TRUE)),
Block = as.factor(sample(c("A", "B", "C", &
浏览 4提问于2020-02-29得票数 2
2回答
在R中使用RMS包预测生存？
、、、、
我正在使用RMS包中的函数survest来生成生存概率。我希望能够获取我的数据的一个子集，并通过survest传递它。我已经开发了一个这样做的for循环。这将运行并输出每组预测值的生存概率。
for (i in 1:nrow(df)){
row <- df[i,]
print(row)
surv=survest(fit, row, times=365)
print(surv)
}
我的第一个问题是，是否有一种方法可以使用survest来预测中位生存期，而不是必须指定特定的时间范围，或者是否有更好的函数可以使用？
其次，我希望能够仅使用我的cox模型的五个预测因子中的四
浏览 0提问于2018-03-27得票数 1
这是回答的问题的更简单的变体。
# Create simple Weibull survival fit using library(survival)
surmo<-survreg( Surv(validtimes, status)~1, dist="weibull")
# Getting Kaplan-Meier
fKM<-survfit( Surv(validtimes, status)~1)
# Plot Kaplan-Meier
plot(fKM,xlab="Time,Days",conf.int=TRUE,mark.tim
浏览 0提问于2019-03-04得票数 1
回答已采纳
我试图在R中绘制一个Kaplan-Meier生存图，但我遇到了一些麻烦。 我是R的新手，所以请原谅我糟糕的代码。 library(survival)
data_time = c(0.19,0.75,0.27,0.26,0.22,0.91,0.21,0.091,0.19,0.37,0.093,0.92,0.046,0.93,042)
data_event = c(1,1,1,1,0,0,1,1,0,0,0,1,1,1,0)
surv_object = Surv(time = data_time, event = data_event)
survfit(surv_object) 这当然会给我一个
浏览 18提问于2019-12-11得票数 1
回答已采纳
我正在按照指令()创建一个参数化的生存学习者来使用MLR。我的密码在下面。
当我试图使MakeLearner(id = "AFT"，"surv.parametric")时，我会得到一个错误dist丢失，即使我已经在代码中指定dist默认值为"weibull“，也不会设置默认值。
makeRLearner.surv.parametric = function() {
makeRLearnerSurv(
cl = "surv.parametric",
package = "survival",
pa
浏览 3提问于2021-12-27得票数 0
1回答
我试图遵循一个教程，在R()的.The计算机，我使用的工作，没有USB端口或互联网连接-它只有R与少数几个库安装。我的工作计算机有“生存，护林员，ggplot2和dplyr”。然而，它并没有“加强”。我试图找出如何在没有“ggfortify”的情况下绘制教程中的图形。下面是我使用的代码：
#load libraries
library(survival)
library(ranger)
library(ggplot2)
library(dplyr)
#load data
data(veteran)
head(veteran)
#Part 1 :
浏览 6提问于2020-11-27得票数 0
回答已采纳
请考虑以下几点：
要创建生存曲线，可以使用survfit包的survival函数。
我的目标是编写一个函数(除其他外)创建这样的曲线，但该函数应该与不同的data.frames一起工作，其列名也不同。此外，分组变量将依赖于各自的数据集。
我设法将不同的data.frame名称传递给函数，但是为survfit和Surv函数提供列名对我来说不起作用。
任何帮助都是非常感谢的。
在我看来，这是一个不同于简单地将data.frame列名传递给函数的问题，如这里所讨论的：
# required libraries
library(survival)
library(flexsurv)
#### Exa
浏览 0提问于2018-07-05得票数 0
回答已采纳
2回答
我使用以下命令绘制了8个不同组的数据帧的存活概率：
fit2<-Surv((time=t2$uptimeDay,event=t2$solved,type='right')~t2$cluster)
plot(fit2,conf.int=F,xlim=c(0, 250),mark.time=c(1,50,100,200),mark=c(1,3,4,2,5,7,6,8,9,10),lwd=1,cex=0.7,lty = 1:11,xlab='Time(days)',ylab='Survival Probability')
这里的集群是一个介于1
浏览 2提问于2015-08-19得票数 0
1回答
我有一个生存kapplan-Meier曲线，我想要外推不同的模型曲线(例如，威布尔，Gompertz，...)使用flexsurv包。我已经成功地做了外推，但我没有找到解决方案来创建外推图的矩阵。
library(survival)
library(flexsurv)
kmsurvival <- survfit (Surv(time, status) ~ 1, data=lung)
summary(kmsurvival)
plot(kmsurvival, xlab="Time", ylab="Survival probability")
Gompertz&
浏览 8提问于2018-08-24得票数 1
回答已采纳
1回答
我需要从我的模型中提取每个固定效果的后验SE估计。 为了便于说明，与我正在使用的数据集类似的数据集是基本R中的ChickWeight数据集。 我提取固定效果的后验估计和间隔的方法如下所示： #load package
library(lme4)
#model
m.surv<-lmer(weight ~ Time + Diet + (1|Chick), data=ChickWeight)
#load packages
library(MCMCglmm)
library(arm)
#set up for fixed effects
sm.surv<-sim(m.surv)
sm
浏览 23提问于2020-12-16得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
