0

我正在为我的数据制作卡图,并且在自定义图时遇到问题:例如,我无法更改输出 chi 和 lambda 限制或更改置信区间。例如,当使用普通绘图函数时,您可以编写:

plots(x,y,xlim=c(a,b)),其中 a 和 b 是您想要的输出限制。我附上了一个用伪随机数据生成的卡图示例。基本上,我希望能够将 chi 和 lambda 限制从默认值 (-1,1,-1,1) 减少到其他值。

这是我使用的代码:

 library('MASS') 
 library("asbio")
 samples = 1000
 r = 0.1
 data = mvrnorm(n=samples, mu=c(0, 0), Sigma=matrix(c(1, r, r, 1), nrow=2), empirical=TRUE)
 X = data[, 1]  # standard normal (mu=0, sd=1)
 Y = data[, 2]  # standard normal (mu=0, sd=1)
 windows()
 opar <- par(mfrow = c(2, 2),mar = rep(1.5, 4),oma = rep(4.1, 4.1))
 plot(X,Y)
 abline(lm(X ~ Y),col="darkred")
 abline(v =median(X),col ="black")
 abline(h =median(Y),col = "black")
 points(median(X),median(Y),col = "red",pch =19)
 chi.plot(X,Y)
 text(paste("r =", round(r, 2)), x = -0.8, y = 0.8)
 text(paste("N =",samples), x = 0.8, y = 0.8)

请帮助chiplots imgae

4

0 回答 0