找回密码
 注册
查看: 3143|回复: 0

R语言求解医学参考值

[复制链接]
发表于 2012-1-3 14:52:31 | 显示全部楼层 |阅读模式
医学参考值(reference value)是指包括绝大多数正常人的人体形态、功能和代谢产物等各种生理及生化指标常数,也称正常值。由于存在个体差异,生物医学书籍并非常数而是在一定范围波动,故采用医学参考值范围(medical reference range)作为判定正常和异常的参考标准。

医学参考值反问涉及用单侧和双侧界值的问题。由于许多生物医学书籍服从或者近似服从正态分布,因此通常采用正态分布发制定医学参考值范围。
95% 医学参考值范围公式:x±1.96 · S 。其中S为标准差。
99% 医学参考值范围公式:x±2.58 · S .

R语言方法如下
#首先产生一组正太分布的数据。
> data<-rnorm(100, mean = 5.8, sd = 0.8)

#计算这组数的平均值,并保留两位有效数字。
>m<-round(mean(data,na.rm =T),digits=2);

#计算这组数据的标准差并乘以1.96。保留两位有效数字。
>std<-round(1.96*sd(data,na.rm =T),digits=2);

#95%双侧医学参考值范围

>rv=paste(m,std,sep="±")
> rv
[1] "5.69±1.41"
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2024-11-25 00:17 , Processed in 0.027031 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表