我正在为我的数据制作卡图,并且在自定义图时遇到问题:例如,我无法更改输出 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)