这次我们用到car这个包
install.packages("car")
library(car)
读取数据
data<-read.csv("/Users/baiyunfan/desktop/data.csv")
整理成这个形式
pdf("/Users/baiyunfan/desktop/data.pdf")
dataEllipse(data[,1],data[,2],levels=0.95,col="black",xlim=c(-1,1),ylim=c(-4,2),xlab="DNMT3A",ylab="PGLYRP2",lwd=1,
grid=F,pch=20,lty=2)
dev.off()
圆圈内表示落入置信区间95%中。