KMggplot2插件似乎只适用于包示例

Question

我正在尝试为Rcmdr使用KMggplot2插件。它似乎只适用于构建它的数据集dataKm。 当我尝试使用其他数据集（如肺）时，我没有得到任何情节 - 只是错误信息＆＃34;参数列的数量不匹配。这是每个情节尝试的代码。 这是我尝试使用肺部数据时的代码，即使没有尝试在风险列表中设置数字。

     library(survival, pos=17)
data(lung, package="survival") 
sapply(c("ggplot2", "grid"), require, character.only = TRUE)
Loading required package: ggplot2
Loading required package: grid
ggplot2    grid 
   TRUE    TRUE 
.df <- data.frame(x = lung$time, y = lung$status, z = factor("At risk"))
.df <- .df[do.call(order, .df[, c("z", "x"), drop = FALSE]), , drop = FALSE]
.fit <- survival::survfit(survival::Surv(time = x, event = y, type = "right") ~ z,  .df)
.fit <- data.frame(x = .fit$time, y = .fit$surv, nrisk = .fit$n.risk, nevent =    .fit$n.event, ncensor= .fit$n.censor, upper = .fit$upper, lower = .fit$lower)
.df <- unique(.df)
.df <- .fit <- data.frame(.fit, .df[, c("z"), drop = FALSE])
Error in data.frame(.fit, .df[, c("z"), drop = FALSE]) : arguments imply differing number of rows: 186, 199
.df <- .fit <- rbind(unique(data.frame(x = 0, y = 1, nrisk = NA, nevent = NA,     ncensor = NA, upper = 1, lower = 1, .df[, c("z"), drop = FALSE])), .fit)
Error in rbind(deparse.level, ...) : numbers of columns of arguments do not match
.cens <- subset(.fit, ncensor == 1)
.plot <- ggplot(data = .fit, aes(x = x, y = y, colour = z)) + geom_step(data =     subset(.fit, !is.na(upper)), aes(y = upper), size = 1, lty = 2, alpha = 0.5,    show_guide = FALSE, na.rm = FALSE) +
 geom_step(data = subset(.fit, !is.na(lower)), +   aes(y = lower), size = 1, lty = 2, alpha = 0.5, show_guide = FALSE, na.rm = FALSE)   + 
geom_step(size =1.5)+         geom_linerange(data = .cens, aes(x = x, ymin = y, ymax =  y + 0.02), size = 1.5) +     
scale_x_continuous(breaks = seq(0, 900, by = 300), limits  = c(0, 900)) +     scale_y_continuous(limits = c(0, 1), expand = c(0.01, 0)) + 
 scale_colour_brewer(palette = "Set1") + xlab("Time from entry") + 
ylab("Proportion of survival") + 
theme_gray(base_size = 14, base_family = "serif") 
Error in +geom_step(size = 1.5) : invalid argument to unary operator+
theme(legend.position = "none")
Error in inherits(x, "theme") : argument "e2" is missing, with no default
 print(.plot)
Error in eval(expr, envir, enclos) : object 'z' not found

Answer 1

我刚收到开发人员的回复，他们表示将在月底上传更新： ＆＃34;我们发现错误是由绑定数据处理引起的。可以使用以下代码。＆＃34;

    sapply(c("ggplot2", "grid"), require, character.only = TRUE)
.df <- data.frame(x = lung$time, y = lung$status, z = factor("At risk"))
.df <- .df[do.call(order, .df[, c("z", "x"), drop = FALSE]), , drop = FALSE]
.fit <- survival::survfit(survival::Surv(time = x, event = y, type = "right") ~ z, .df)
.fit <- data.frame(x = .fit$time, y = .fit$surv, nrisk = .fit$n.risk,
               nevent = .fit$n.event, ncensor= .fit$n.censor, upper = .fit$upper, lower = .fit$lower)
.df <- .df[!duplicated(.df$x), ]
.df <- .fit <- data.frame(.fit, .df[, c("z"), drop = FALSE])
.df <- .fit <- rbind(unique(data.frame(x = 0, y = 1, nrisk = NA, nevent = NA, ncensor = NA, upper = 1, lower = 1, .df[, c("z"), drop = FALSE])), .fit)
.cens <- subset(.fit, ncensor == 1) .plot <- ggplot(data = .fit, aes(x = x, y = y, colour = z)) + 
geom_step(data = subset(.fit, !is.na(upper)), aes(y = upper), size = 1, lty = 2, alpha = 0.5, show_guide = FALSE, na.rm = FALSE) +
geom_step(data = subset(.fit, !is.na(lower)), aes(y = lower), size = 1, lty = 2, alpha = 0.5, show_guide = FALSE, na.rm = FALSE) + 
geom_step(size = 1.5) +
geom_linerange(data = .cens, aes(x = x, ymin = y, ymax = y + 0.02), size = 1.5) +
scale_x_continuous(breaks = seq(0, 900, by = 300), limits = c(0, 900)) +
scale_y_continuous(limits = c(0, 1), expand = c(0.01, 0)) + 
scale_colour_brewer(palette = "Set1") + 
xlab("Time from entry") + 
ylab("Proportion of survival") + 
theme_gray(base_size = 14, base_family = "serif") + 
theme(legend.position = "none")
print(.plot)