现在基因沿x轴分部,y轴值不变代表染色体宽度,想要基因显示成矩形,带颜色不要边框,求助求助
这是自己写的,R刚入门的小白
setwd("E:\\keyan")
b<-read.table("chuliaa.txt",header=FALSE,sep="\t")
plot(1,xlim=c(0,12000),ylim=c(-100,0),type="n")
x<-c(b[,1],b[,1],b[,2],b[,2])
y<-c(b[,3],b[,4],b[,4],b[,3])
polygon(x,y,col =seq(2,5))
照理说想要出来3个矩形,出来的却是个大红沙漏,还说读到最后一横不完全是什么意思·····
求大神指点,最好能有详细的步骤,急求,拜谢T T
小样
前两个是基因的起点和终点。3,4是y轴的范围,求不带边框的矩形