P for trend/ per 1 SD 是怎么计算得到的？

详情请点击下方：

欢迎报名！孟德尔随机化方法高级班8.3-4，药靶多组学，中介，肠道菌群。。。

---

医学研究趋势性分析，很多都知道，我们也写过论文。

临床研究趋势性分析，P trend，它是如何计算的？

其中有个P trend就是线性趋势分析。P值的算法有很多种。原始值法、等级法、等级中位数。

最近有个朋友说P trend / per 1 sd 怎么算

的确，也有一些文献应用了该方法开展趋势性分析.

如何实现per 1 sd 

Per 1 sd的实现，其实就是把原始数据进行标准化，另存为一个新的变量X，新变量X因为是被标准化后的数据，因此其均数和标准差为0和1。然后让x进入模型进行分析。请问大家此时x每增加1个单位，效应量增加的风险为HR。因为标准差为1，此时x增加1个单位，就是Per 1 sd。1=Per 1 sd。就是自变量每增加1个标准差。

R语言实现

#导入数据
exp <- foreign::read.spss(xx.sav", 
                             to.data.frame = T,
                             use.value.labels = F,
                             reencode  = "utf-8")
## re-encoding from utf-8


# 新建一列weight
df16_2$weight <- rnorm(54, 70,11)




# 进行标准化
exp$weight.scaled <- scale(exp$weight)


# 进行逻辑回归
f <- glm(y ~ weight.scaled, data = exp,family=binomial)
broom::tidy(f,conf.int=T,exponentiate=T)

如此一来，P值即为P trend per 1 sd 

风暴统计zstats平台快速计算P trend

郑老师统计平台，风暴统计zstats平台，现在同样可以实现多种方法的线性趋势分析。

• 定量暴露直接计算

• 定量暴露转为等级再计算

• 等级中位数法计算

• per 1 sd 法计算

线性趋势分析具体网址：

https://www.zstats.cn/software2/trend1/

或者百度、必应Bing搜索“风暴统计”

本平台上线的所有工具都是免费的

1.进入风暴统计平台

首先，浏览器搜索风暴统计，依次点击"风暴智能统计"——"线性趋势分析与RCS曲线"！进入分析界面后，根据提示，完成数据的导入与整理。

2.数据导入后对定量暴露因素进行转换

有多种转换：按照中位数、按照四分位数、自定义、还有数据标准化，P trend per 1 sd 就是选择“数据标准化”。

会产生新的变量，即通过标准正态分布分布Z分布产生的新变量。

3. 接下来开展回归分析, 进行线性趋势分析

我们采用logistic回归，对二分类变量“low”建模，此时自变量即刚才才生的“lwt_scale"了

就可以得到以下的分析结果

我们还可以同时开展多模型进行数据分析

4.下载结果

分析完成后，可以选择小数位数，默认情况下，P值为3位小数，其他统计量为2位小数。指定小数位数后，P值与统计量的小数位数将会统一。调整完成后，下载最终的三线表结果！

我们最后下载的word版结果也十分的清晰规范！省去了绘制表格，填写数据的繁琐。

---

详情请点击下方：

欢迎报名，孟德尔随机化1对1答疑指导学习班！