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%中。

©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容