可下载资源

本文分析的数据、代码、报告分享至会员群

作者

Qing Li
✉ 联系我们

暴露数据以15岁至65岁之间的5岁年龄区间收集。

在这里，我提供了一个简短的摘要来介绍概念和定义。

暴露-滞后-反应关联

DLNM的建模类用于描述关联，在该关联中，暴露和结果之间的依赖关系会在时间上滞后。可以使用两个不同且互补的观点来描述此过程。我们可以说，在时间t处的暴露事件确定了在时间t +l处的未来风险。使用后向视角，时间t的风险由过去在时间t-l经历的一系列风险确定。这里的l是滞后，表示暴露和测得的结果之间的滞后。

第一是与特定暴露值相关联的滞后反应曲线，定义为预测变量特定性关联。这被解释为与时间t风险相关的时间t +l的风险贡献序列。

第二是与特定滞后值相关联的暴露-反应曲线，该特定滞后值定义为滞后特定关联。这被解释为与在时间t处发生的暴露值相关联的在时间t +l处的暴露-反应关系。

第三个也是最重要的是与在考虑的滞后期内经历的整个暴露历史相关的暴露反应曲线，定义为总体累积关联。使用正向视角，这被解释为表示时间t发生的给定暴露期间[t，t+L]期间经历的净风险的暴露反应关系。

第一个未命名的参数x指定要转换的向量，而参数度设置多项式的度。定义分层函数是通过strata（）指定的。

strata(1:5,breaks=c(2,4))[,]
1 2
[1,] 0 0
[2,] 1 0
[3,] 1 0
[4,] 0 1
[5,] 0 1

结果是带有附加类别“层”的基础矩阵。转换是定义对比的虚拟参数化。参数break定义了层的右开放区间的下边界。

阈值函数通过thr（）指定。一个例子：

thr(1:5,thr.value=3,side="d")[,]
1 2
[1,] 2 0
[2,] 1 0
[3,] 0 0
[4,] 0 1
[5,] 0 2

结果是具有附加类别“ thr”的基础矩阵。参数thr.value定义一个带有一个或两个阈值的向量，而side用于指定高（“ h”，默认值），低（“ l”）或双精度（“ d”）阈值参数化。

基本转换

此函数代表以dlnm为单位进行基本转换的主要函数，适用于指定暴露-反应和滞后-反应关系。它的作用是应用选定的转换并以适用于其他函数（例如crossbasis（）和crosspred（））的格式生成基本矩阵。以下示例复制了该部分中显示的多项式变换：

onebasis(1:5,fun="poly",degree=3)
b1 b2 b3
[1,] 0.2 0.04 0.008
[2,] 0.4 0.16 0.064
[3,] 0.6 0.36 0.216
[4,] 0.8 0.64 0.512
[5,] 1.0 1.00 1.000
attr(,"fun")
[1] "poly"
attr(,"degree")
[1] 3
attr(,"scale")
[1] 5
attr(,"intercept")
[1] FALSE
attr(,"class")
[1] "onebasis" "matrix"
attr(,"range")
[1] 1 5

交叉基

这是dlnm软件包中的主要函数。它在内部调用onebasis（）来生成暴露-反应和滞后-反应关系的基矩阵，并通过特殊的张量积将它们组合起来，以创建交叉基，该交叉基在模型中同时指定了暴露-滞后-反应关联性。它的第一个参数x的类定义如何解释数据。可以使用第二个变量lag修改滞后期。

作为一个简单的示例，我模拟了2-5个滞后期内3个对象的暴露历史矩阵：它们中的每一个都将传递给onebasis（）来分别构建暴露-反应和滞后-反应关系的矩阵。仅用于时间序列数据的附加参数组定义了被视为单独无关序列的观察组，例如在季节性分析中可能有用。作为一个简单的示例，我模拟了2-5个滞后期内3个对象的暴露历史矩阵：它们中的每一个都将传递给onebasis（）来分别构建暴露-反应和滞后-反应关系的矩阵。作为一个简单的示例，我模拟了2-5个滞后期内3个对象的暴露历史矩阵：

 > hist
lag2 lag3 lag4 lag5
sub1 1 3 5 6
sub2 7 8 9 4
sub3 10 2 11 12

lag=c(2,5),argvar=list(fun="poly",degree=2),
arglag=list(fun="strata",breaks=4))[,]
v1.l1 v1.l2 v2.l1 v2.l2
sub1 1.250000 0.9166667 0.4930556 0.4236111
sub2 2.333333 1.0833333 1.4583333 0.6736111
sub3 2.916667 1.9166667 2.5625000 1.8402778

该函数返回“ crossbasis”类的矩阵对象。它首先使用argvar和arglag列表中的参数调用onebasis（），以建立暴露反应空间和滞后反应空间的矩阵基础。在另一个示例中，我将crossbasis（）应用于数据集中的变量temp，该数据集表示1987-2000年期间日平均温度序列：

例如，我使用创建的交叉基矩阵cb，使用数据集时间序列数据来研究温度与心血管疾病死亡率之间的关联。首先，我将一个简单的线性模型与模型公式中包含的交叉基矩阵拟合。然后，我通过使用cross-basis和回归模型对象作为前两个参数调用crosspred（）来获得预测：

crosspred(cb,model,at=-20:30) 

结果是“ crosspred”类的列表对象，其中的存储预测和有关模型的其他信息，例如系数和与交叉基参数有关的关联（协）方差矩阵的一部分。可以为特定的预测器-滞后组合选择预测的网格。例如，我提取温度为-10°C且滞后5的预测和置信区间，然后提取25°C的整体累积预测：

 > pred$allfit["25"]
25
1.108262

第一个结果表明，在给定的一天中，-20°C的温度会在五天后导致0.95例心血管死亡的增加，或者在给定的一天中，温度为-6摄氏度时，心血管死亡的数目增加0.95。其他类型的预测可以通过crosspred（）获得。特别是，如果模型链接等于log或logit，则将自动返回取幂的预测。如果参数cum设置为TRUE，则是累积预测的矩阵cum。

crosspred（）的另一种用法是预测特定的暴露历史记录集的影响。这可以通过输入暴露历史矩阵作为参数来实现。例如，我们可以从拟合模型中预测出，在过去10天暴露于30°C和在滞后期的其余时间暴露于22°C之后，心血管死亡的总体累积增加：如果参数cum设置为TRUE，则包括增量累积预测的矩阵cum，并将其存储在组件cum中。crosspred（）的另一种用法是预测特定的暴露历史记录集的影响。这可以通过输入暴露历史矩阵作为参数来实现。例如，我们可以从拟合模型中预测出，在过去10天暴露于30°C和在滞后期的其余时间暴露于22°C之后，心血管死亡的总体累积增加：

 > crosspred(cb,model,at=histpred)$allfit
1
5.934992

dlnm软件包的主要优点之一是，用户可以使用标准回归函数执行DLNM，只需在模型公式中包括交叉基矩阵即可。函数crosspred（）自动处理来自回归函数lm（）和glm（），gam（）（软件包mgcv），coxph（）的模型。

降维

DLNM的拟合度可以降低到预测变量或滞后的一个维度，仅以诸如总累积暴露反应表达。该计算通过函数crossreduce（）进行，该函数的工作原理与crosspred（）非常相似。前两个自变量base和model指定交叉基矩阵和需要对其执行计算的模型对象。减少的类型由类型定义，带有选项“ overall”-“ lag”-“ var”，用于汇总总体累积暴露反应，滞后特异性暴露反应或预测变量特异性滞后反应。

绘图

一维或二维关联的解释通过图形表示来辅助。通过方法函数plot（），lines（）和points（）为类“ crosspred”和“ crossreduce”提供高级和低级绘图功能。例如，我使用对象pred中的预测。plot（）方法可以通过参数ptype为“ crosspred”对象生成不同类型的图。具体来说，它会生成整个二维暴露-滞后-反应关联的图形。二维关联可以绘制为3-D或等高线图，例如：

> plot(pred,ptype="3d",main="3D plot"

​​

可以通过选择不同的ptype获得定义的关联的摘要。

> plot(pred,"overall" 

在这种情况下，方法函数plot（）在内部调用函数plot.default（），如上面的示例所示，可以将其他特定参数添加到函数调用中。通过设置ptype =“ slices”，可以将滞后特异性和预测因子特异性关联分别绘制为暴露-反应和滞后-反应关系，因为它们是在3-D曲面中沿特定维度切割的切片。例如：

> plot(pred,"slices",lag=5 

​​​

这两个图分别代表了滞后5的暴露反应和特定于25°C温度的滞后反应。参数lag和var指定必须分别绘制lag和特定于预测变量的关联的值。