SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)

生存分析

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第1张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第2张图片

cox可以用在金融等行业
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第3张图片


SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第4张图片

lifereg过程

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第5张图片

response是因变量,censor是标识是否是截尾数据,list中是终值的标记是哪些

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第6张图片
数据
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第7张图片

censor是标识截尾值,当t<0的时候,censor等于1
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第8张图片

右终检值就是知道生存了x天以上,但具体多少不知道
左终检值就是知道生存了x天以下,但是具体多少不知道
区间终检值是知道在一个区间内

下面的结果可以看到只有第一个和最后一个相关性比较大
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第9张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第10张图片

简化以后:
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第11张图片
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第12张图片

在测试集上的表现

lifetest过程

生存分布函数sdf
累计分布函数,1-生存分布函数
概率密度函数是生存分布函数的导数
风险函数是概率密度函数除以生存分布函数
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第13张图片
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第14张图片
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第15张图片
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第16张图片

注意这里对数据的读入,自己写了一个变量C,那么读入的时候
fail和withdraw读入两个数,那么输出两行:
C和fail;C和withdraw
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第17张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第18张图片

Cox比例风险模型

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第19张图片

phreg过程

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第20张图片

时间序列

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第21张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第22张图片

因变量取决于之前某一些时间的值,被称为自回归
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第23张图片
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第24张图片

回归只能刻画发展趋势,但是具体的细节无法描述
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第25张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第26张图片

T是长期趋势,S是季节变动,C是循环变动,I是随机扰动
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第27张图片

平稳时间序列很常见
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第28张图片
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第29张图片

常用于描述随机扰动
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第30张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第31张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第32张图片

在这里插入图片描述

ARIMA模型可以通过差分使不平稳的随机序列,变成平稳的随机序列
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第33张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第34张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第35张图片

截尾就是到了某个点,基本等于0了
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第36张图片

ARIMA过程

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第37张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第38张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第39张图片

arima是一个交互式过程
nlag=12 表示要计算12个相关系数的值
在这里插入图片描述

自相关系数
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第40张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第41张图片

q=2,2阶的滑动平均模型
p=1,1阶的自回归模型
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第42张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第43张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第44张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第45张图片

因为AR1的SBC比较小,所以选择AR1模型进行预测
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第46张图片

第二个例子
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第47张图片

上证综指是非平稳的
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第48张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第49张图片

可以看到没有明显的减小情况,是非平稳的
SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第50张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第51张图片

SAS学习14、15(生存分析、lifereg过程、lifetest过程、phreg过程、时间序列、arima过程)_第52张图片

你可能感兴趣的:(数据分析,sas)