我想用R计算出Cox回归。 我喜欢计算国家的生存概率。我对一个国家取决于多个变量的独裁统治之前的生存时间(或该国家成为民主国家的时间)感兴趣。
在下面,您可以从数据中看到一个示例。
ccode cname year gle_cgdpc p_polity2 time1 time2 dummy
4 AFG 1946 NA -10 1946 1947 0
4 AFG 1947 NA -10 1947 1948 0
4 AFG 1948 NA -10 1948 1949 0
4 AFG 1949 NA -10 1949 1950 0
4 AFG 1950 130.82 -10 1950 1951 0
4 AFG 1951 141.98 -10 1951 1952 0
4 AFG 1952 144.22 5 1952 1953 1
4 AFG 1953 145.22 -8 1953 1954 0
4 AFG 1954 145.22 7 1954 1955 1
4 AFG 1955 147.22 -10 1955 1956 0
5 ALB 1946 122.52 -9 1946 1947 0
5 ALB 1947 111.42 5 1947 1948 1
5 ALB 1948 113.29 -1 1948 1947 0
我创建了一个虚拟变量,当一个国家达到特定的民主指数时,该变量的值为1。
到目前为止,数据集还可以吗? 我现在可以只使用变量吗?
coxph(Surv(time1,time2,dummy) ~ var,method = "efron",data = data)
还是我完全错了?
非常感谢您的帮助!