我想上运行的一组数据我已经收集了生存分析。 在该数据帧(M3)中,每一行是一个新的患者和每列是我已确定的突变。 我做了一个二进制数据表来表示每个患者是正的还是负的突变。 我可以运行为每列(突变)一survfit功能,但我通过他们有数百个,并希望循环。 我写了下面的代码,但不认为这是正确的(什么都没有输出)。
for (i in m3[,2:256]) {survdiff(Surv(m3$Overall.Survival, m3$Status) ~ i,
data = m3)}
有一次,我收集这些数据我想和每个突变(列)作为行表,并从该survfit对象作为列放的p值。
我不知道为什么我没有对任何输出循环,更要如何产生新的数据帧。 我相信我会被它子集化。