R语言画置信椭圆

这次我们用到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%中。