- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 23 个
- 通用积分
- 3.7941
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 1151 点
- 帖子
- 63
- 精华
- 0
- 在线时间
- 85 小时
- 注册时间
- 2016-1-3
- 最后登录
- 2021-6-29
|
藤椅
苏谷壳
发表于 2017-8-18 17:07:26
- par(mfrow=c(2,2))
- #1
- bland_altman_plot <- function(the_1st_peak1,the_1st_peak2,xlab="Average testresult", ylab="Deviation of experimental test")
- {
- d1 <- ((the_1st_peak1 + the_1st_peak2)/2)
- diff1 <- the_1st_peak1 - the_1st_peak2
- plot(diff1~d1,pch=16,ylim=c(-20,160),xlab=xlab,ylab=ylab,col=cols_L)
- abline(h=mean(diff1)-c(-1.96,1.96)*sd(diff1),lty=2)
- abline(h=mean(diff1))
- }
- bland_altman_plot(the_1st_peak1,the_1st_peak2,xlab="Mean of the_1st_peak1 and the_1st_peak2",ylab="Difference the_1st_peak1 and the_1st_peak2")
- #2
- bland_altman_plot <- function(PVIF1,PVIF2,xlab="Average testresult", ylab="Deviation of experimental test")
- {
- d2 <- ((PVIF1 + PVIF2)/2)
- diff2 <- PVIF1 - PVIF2
- plot(diff2~d2,pch=16,ylim=c(-10,100),xlab=xlab,ylab=ylab,col=cols_L)
- abline(h=mean(diff2)-c(-1.96,1.96)*sd(diff2),lty=2)
- abline(h=mean(diff2))
- legend(x=60,y=100,c("1","2","3","4"),col=cols,pch=c(16),bty="n",ncol=1,cex=0.9,pt.cex=0.9,xpd=FALSE,text.width=2,x.intersp=0.6)
- }
- bland_altman_plot(PVIF1,PVIF2,xlab="Mean of PVIF1 and PVIF2",ylab="Difference PVIF1 and PVIF2")
- #3
- bland_altman_plot <- function(the_1st_peak1,PVIF1,xlab="Average testresult", ylab="Deviation of experimental test")
- {
- d3 <- ((the_1st_peak1 + PVIF1)/2)
- diff3 <- the_1st_peak1 - PVIF1
- plot(diff3~d3,pch=16,ylim=c(-100,140),xlab=xlab,ylab=ylab,col=cols_L)
- abline(h=mean(diff3)-c(-1.96,1.96)*sd(diff3),lty=2)
- abline(h=mean(diff3))
- }
- bland_altman_plot(the_1st_peak1,PVIF1,xlab="Mean of the_1st_peak1 and PVIF1",ylab="Difference the_1st_peak1 and PVIF1")
- #4
- bland_altman_plot <- function(the_1st_peak2,PVIF2,xlab="Average testresult", ylab="Deviation of experimental test")
- {
- d4 <- ((the_1st_peak2 + PVIF2)/2)
- diff4 <- the_1st_peak2 - PVIF2
- plot(diff4~d4,pch=16,ylim=c(-8,5),xlab=xlab,ylab=ylab,col=cols_L)
- abline(h=mean(diff4)-c(-1.96,1.96)*sd(diff4),lty=2)
- abline(h=mean(diff4))
- }
- bland_altman_plot(the_1st_peak2,PVIF2,xlab="Mean of the_1st_peak2 and PVIF2",ylab="Difference the_1st_peak2 and PVIF2")
复制代码
|
|