R 删除虚线置信区间线并修改图例

如何从绘图中删除置信区间虚线

对于图例，为什么第2行是虚线而第1行是实线？有没有办法把两条实线都画出来？我不确定这段代码是做什么的-lty=1:2，cex=0.8

我已经用相应的命令附加了下面的一个数据集。谢谢

> km2 <- survfit(Surv(times_start, times_end, type="interval2")~ 1)
> summary(km2)
Call: survfit(formula = Surv(times_start, times_end, type = "interval2") ~ 
    1)

    time n.risk n.event survival std.err lower 95% CI upper 95% CI
    1.31  284.0    2.04    0.993 0.00743       0.9782        1.000
    2.06  282.0    1.03    0.989 0.00856       0.9724        1.000
    2.81  280.9    4.18    0.974 0.01770       0.9395        1.000
    3.56  276.7    2.10    0.967 0.01977       0.9278        1.000
    4.31  274.7    8.42    0.937 0.03557       0.8659        1.000
    5.06  266.2    2.11    0.930 0.03678       0.8556        1.000
    5.81  264.1   13.79    0.881 0.05659       0.7642        1.000
    6.38  250.3    4.24    0.867 0.05838       0.7440        1.000
    7.12  246.1    9.55    0.833 0.06455       0.6940        1.000
    7.88  236.5    7.43    0.807 0.06785       0.6577        0.990
    8.62  229.1    3.18    0.796 0.06866       0.6432        0.984
    9.38  225.9    9.55    0.762 0.07262       0.5963        0.974
   10.12  216.4    4.24    0.747 0.07355       0.5769        0.967
   10.88  212.1   18.03    0.683 0.08015       0.4883        0.957
   11.62  194.1    2.12    0.676 0.08017       0.4793        0.953
   12.56  192.0    8.68    0.645 0.08074       0.4415        0.944
   13.31  183.3    3.33    0.634 0.08069       0.4274        0.940
   14.06  180.0   13.67    0.586 0.08082       0.3690        0.929
   14.81  166.3    2.34    0.577 0.08058       0.3595        0.927
   15.56  164.0   10.33    0.541 0.07952       0.3176        0.921
   16.31  153.6    4.35    0.526 0.07886       0.3005        0.920
   17.06  149.3    3.32    0.514 0.07831       0.2875        0.919
   17.81  146.0    3.89    0.500 0.07761       0.2724        0.919
   18.56  142.1    2.25    0.492 0.07718       0.2638        0.919
   19.31  139.8    5.37    0.473 0.07605       0.2435        0.921
   20.06  134.5    3.36    0.462 0.07530       0.2309        0.923
   20.81  131.1    9.03    0.430 0.07291       0.1983        0.932
   21.56  122.1    5.54    0.410 0.07128       0.1790        0.941
   23.81  116.5   19.40    0.342 0.06341       0.1182        0.990
   24.56   97.1   24.24    0.257 0.05029       0.0575        1.000
 5013.31   72.9   72.88    0.000 0.00000           NA           NA

> plot(km, xmax=50, xlab="Time", ylab="Log cumulative hazard", col="blue", fun="cloglog")
> lines(km2, xlab="Time", ylab="Log cumulative hazard", col="red", fun="cloglog")
> legend("topleft", legend=c("Line 1", "Line 2"), col=c("red", "blue"), lty=1:2, cex=0.8)

>平方公里摘要（平方公里）
调用：survfit（公式=Surv（times\u start，times\u end，type=“interval2”）~
1)
时间风险事件存活标准误差下限95%CI上限95%CI
1.31  284.0    2.04    0.993 0.00743       0.9782        1.000
2.06  282.0    1.03    0.989 0.00856       0.9724        1.000
2.81  280.9    4.18    0.974 0.01770       0.9395        1.000
3.56  276.7    2.10    0.967 0.01977       0.9278        1.000
4.31  274.7    8.42    0.937 0.03557       0.8659        1.000
5.06  266.2    2.11    0.930 0.03678       0.8556        1.000
5.81  264.1   13.79    0.881 0.05659       0.7642        1.000
6.38  250.3    4.24    0.867 0.05838       0.7440        1.000
7.12  246.1    9.55    0.833 0.06455       0.6940        1.000
7.88  236.5    7.43    0.807 0.06785       0.6577        0.990
8.62  229.1    3.18    0.796 0.06866       0.6432        0.984
9.38  225.9    9.55    0.762 0.07262       0.5963        0.974
10.12  216.4    4.24    0.747 0.07355       0.5769        0.967
10.88  212.1   18.03    0.683 0.08015       0.4883        0.957
11.62  194.1    2.12    0.676 0.08017       0.4793        0.953
12.56  192.0    8.68    0.645 0.08074       0.4415        0.944
13.31  183.3    3.33    0.634 0.08069       0.4274        0.940
14.06  180.0   13.67    0.586 0.08082       0.3690        0.929
14.81  166.3    2.34    0.577 0.08058       0.3595        0.927
15.56  164.0   10.33    0.541 0.07952       0.3176        0.921
16.31  153.6    4.35    0.526 0.07886       0.3005        0.920
17.06  149.3    3.32    0.514 0.07831       0.2875        0.919
17.81  146.0    3.89    0.500 0.07761       0.2724        0.919
18.56  142.1    2.25    0.492 0.07718       0.2638        0.919
19.31  139.8    5.37    0.473 0.07605       0.2435        0.921
20.06  134.5    3.36    0.462 0.07530       0.2309        0.923
20.81  131.1    9.03    0.430 0.07291       0.1983        0.932
21.56  122.1    5.54    0.410 0.07128       0.1790        0.941
23.81  116.5   19.40    0.342 0.06341       0.1182        0.990
24.56   97.1   24.24    0.257 0.05029       0.0575        1.000
5013.31 72.9 72.88 0.000 0.00000不适用
>绘图（km，xmax=50，xlab=“时间”，ylab=“记录累积危险”，col=“蓝色”，fun=“cloglog”）
>线条（km2，xlab=“时间”，ylab=“记录累积危险”，col=“红色”，fun=“cloglog”）
>图例（“左上角”，图例=c（“第1行”，“第2行”），列=c（“红色”，“蓝色”），lty=1:2，cex=0.8）

我不确定您正在使用哪些数据，但我使用了

lung

（在执行

库（生存）

后应该可以使用这些数据）：

---

km2
？plot.survfit
和
线条。survfit
将非常有用，而且
？par
谢谢！它起作用了。我做了lty=1:1，使两者都成为固体。我尝试删除它们，但图例中没有显示线条。不过，非常感谢你的帮助。我有一个后续问题-我应该使用什么代码来查找该图的斜率？
km2 <- survfit(Surv(time, status) ~ 1, data = lung)

plot(km2, xmax=50, xlab="Time", 
     ylab="Log cumulative hazard", 
     col="blue", fun="cloglog", conf.int = FALSE)