我已经尝试对数据进行分类,我正在尝试对以下示例数据进行生存分析.n是每组的单位数,时间,失败指标的组合.
> df <- structure(list(group = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("", "A", "B"), class = "factor"), t = c(0L, 1L, 2L, 3L, 1L, 2L, 3L, 0L, 1L, 2L, 3L, 1L, 2L, 3L), failure = c(0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L), n = c(40000L, 30000L, 20000L, 10000L, 5L, 4L, 3L, 20000L, 15000L, 14000L, 11000L, 10L, 6L, 4L)), .Names = …我在R中有两个不同的字符向量,我想将它们组合起来用于列名:
groups <- c("Group A", "Group B")
label <- c("Time","Min","Mean","Max")
当我尝试使用粘贴时,我得到结果:
> paste(groups,label)
[1] "Group A Time" "Group B Min" "Group A Mean" "Group B Max"
是否有一个简单的功能或设置可以粘贴这些以获得以下输出？
[1] "Group A Time" "Group A Min" "Group A Mean" "Group A Max" "Group B Time"
[6] "Group B Min" "Group B Mean" "Group B Max"
如果我有一个值,并且我想要的附加项列表乘以该值:
n = 10
terms = [1,2,3,4]
是否可以使用列表推导来执行以下操作:
n *= (term for term in terms) #not working...
或者是唯一的方法:
n *= reduce(lambda x,y: x*y, terms)
这是在Python 2.6.2上.谢谢!
我是R的新手,正在尝试绘制survfit生存曲线.
在玩弄survfit物体时,我发现我得到了2个不同的图表:
library(survival)
#example survfit object
mysurvfit <- survfit(Surv(time, status)~1, data=aml)
#default survfit plot, survival curve with upper & lower conf intervals
plot(mysurvfit, mark.time=FALSE, conf.int=TRUE)
#create another curve by accessing surv, upper, lower
#(I'd expect this to produce the same as above, but it doesn't)
lines(mysurvfit$surv, col="blue",lty=1)
lines(mysurvfit$upper, col="blue",lty=2)
lines(mysurvfit$lower, col="blue",lty=2)
为什么这些曲线不同？我对幸存物体的遗漏是什么？