我使用“mcr”包对两个变量进行加权Deming回归。我试图为每个点添加错误条(水平和垂直。我想知道程序中是否有一个函数或一个可以与mcr包一起使用的程序,以便为x和y显示SE。我找到了一个使用箭头的临时解决方案,如下所示;但我不相信它是显示我的数据的最佳方式。
到目前为止,我发现ggplot2不适用于mcr和se.lines似乎也没有工作。任何意见都将不胜感激。
此外,如果有人能解释我在哪里可以找到X值的标准误差值,我将不胜感激。
谢谢大家的时间。
xo<- c(18.05,15.78,23.44,15.93,25.71,16.06,21.36,18.11,14.78,11.74,17.56,17.04,14.20,10.68,17.96)
yo<-c(3.07,2.32,2.63,2.50,3.02,2.24,3.13,2.46,2.23,2.31,2.55,2.69,2.32,1.77,2.97)
library(mcr)
w.dem.reg.o<-mcreg(xo,yo,method.reg="WDeming",method.ci="jackknife")
w.dem.reg.o@para #getting info on the line
MCResult.plot(w.dem.reg.o,x.lab="parameter1",y.lab="parameter2",sub="",main="Parameter 1 vs Parameter 2: Deming Regression",points.col="gold",points.cex =1,points.pch=19,ci.area = TRUE,reg=TRUE,identity=TRUE,ci.border=TRUE,add.legend = FALSE)
legend(x="topleft",c("Y = 1.07 + 0.09X", "SE of Intercept = 0.4724","SE of Slope = 0.0276")) #aesthetic.only
arrows(xo,yo+0.15692,
xo,yo-0.15692,length=0,angle=90,code=3,col="purple")
arrows(xo+0.4724,yo,
xo-0.4724,yo,length=0,angle=90,code=3,col="purple")
points(xo,yo)
找到