R语言中生存分析的“~1”代表什么意思

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近在研究生存分析，发现R语言中做生存分析时，类似这种kmsurvival<-survfit(Surv(time,time2,event)~1)，右边的括号要写“~1”才能运行，否则不能运行，我知道如果是“~x”的意思，但是不明白“~1”是什么意思？求大神告知（R语言小白，毕业论文要用空间生存模型，才来研究的）

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：生存分析 R语言 Survival SUR KMS

你好，想请教一个问题，一篇论文中对变量的定义和数据来源的说明中，发现有些变量是时间序列数据，而有些变量是截面数据，论文的目的是用Weibull hazard model做实证分析。

比如 variable1  它所描述的数据为1970-2015的时间序列数据
     variable2  它所描述的数据为2000-2010年的平均值
     variable3  ......
     variable4  ......

     实在想不通它的数据结构是什么样子的？

是不是以下这种形式呢？先谢谢啦。

year   variable1               varible2    variable3   variable4  ...
1970        数值                 缺失        ...         ...
1971        数值                 缺失
1972        数值                 缺失
1973        数值                 缺失
1974        数值                 缺失
1975        数值                 缺失
1976        数值                 缺失
1977        数值                 缺失
1978        数值                 缺失
…         …                  …
2000        数值        2000-2015年变量2的平均值
2001        数值        2000-2015年变量2的平均值
2002        数值        2000-2015年变量2的平均值
2003        数值        2000-2015年变量2的平均值
2004        数值        2000-2015年变量2的平均值
2005        数值        2000-2015年变量2的平均值
2006        数值        2000-2015年变量2的平均值
2007        数值        2000-2015年变量2的平均值
2008        数值        2000-2015年变量2的平均值
2009        数值        2000-2015年变量2的平均值
2010        数值        2000-2015年变量2的平均值
2011        数值        2000-2015年变量2的平均值
2012        数值        2000-2015年变量2的平均值
2013        数值        2000-2015年变量2的平均值
2014        数值        2000-2015年变量2的平均值
2015        数值        2000-2015年变量2的平均值

我也不太了解，你参考的这篇文章所在期刊权威吗？我还没有做这方面的实证，如果这样运行没问题，而论文中又没说明的话，可能是这种吧
，你可以发邮件尝试咨询下作者~

~在R里面用于构建公式对象，~左边为因变量，右边为自变量。

马甲1号 发表于 2016-6-16 22:58
~在R里面用于构建公式对象，~左边为因变量，右边为自变量。

可是右边是1代表什么意思呢，我后面看帮助文档，猜想，1可能代表只是做一个生存曲线，并未对任何自变量做拟合，不知道对不对

冰小二 发表于 2016-6-17 10:28
可是右边是1代表什么意思呢，我后面看帮助文档，猜想，1可能代表只是做一个生存曲线，并未对任何自变量做 ...

~1是表示模型只含有一个常数截距，无自变量，~0表示无截距亦无自变量，不过在survfit里面~1好像和~0没区别

事件发生记为1，删失为0，则~1.
若事件记为2，删失为0，则~2

在Surv函数中，默认censor数据1为样本在考察期间死亡。因此survfit中公式为~1，意思是对死亡数据进行拟合计算

波浪号右边代表strata变量，为1说明只有一条生存曲线