作图求助 - 数学建模社区-数学中国

x <- c(2,2,2,3,5,7,7,7,7,8,8,8,9,9,10,10,10,10,10,10,10,10,10,10,20,20,20,20,21,22,25,26,29,30,31,32,34,35,36,37,38,39,39,40,40,42,42,44,44,48,49,49,49,50,50,52,52,53,53,53,57,57,58,59,60,61,61,61,62,62,64,65,66,67,66,68,72,75,75,74,75,78,80,80,81,81,82,83,84,84,86,87,88,88,90,94,92,93,95,96,97,103,106,109);
7 ?/ {; n& E' t! z' U" f1 C
y <- c(5,20,24,20,26,25,25.5,29,31,24,26,28,29,29,20,24,27,29,30,31,33,34,38,40,41,51,55,62,40,63,49,62,62.5,74,50,61,52,50,60,70,72,58,65,64,65,68,70,69,72,74,78,80,76,74,40,65,66,66,68,68,78,76,79,80,79,92,75,74,78,79,80,78,82,85,83,79,84,86,73,90,88,87,89,89,91,90,91,90,88,94,60,88,70,90,94,93,110,66,98,90,100,106,70,90);
2 t2 ]1 z1 t$ j5 z* M
plot(x,y,pch=20,axes=T,xaxt="n",yaxt="n",ylab="methylation",xlab="Age(yrs)");
title(col.sub="red",sub="Fetal Age(wks) ",main="cg04640913");
* u% m/ G) d! Z% h2 N, p2 Z
text(98,15,"CDH22",font=3,cex=2);
, B7 Q* o1 b. U, V' \* A" f5 I
x1 <- c(x[2],x[22]);
' A" t: h$ j* ]; t
y1 <- c(y[2],y[22]);) i. y' e+ y5 M' a* \1 Z4 z! T
lines(x1,y1,lty=1,lwd=5,col=8);
! D; ~9 A% l! o
x2 <- c(x[26],x[40]);% }1 M) U9 I9 a* D# r
y2 <- c(y[26],y[40]);9 r4 l* S" H( l9 W/ G& I. F
lines(x2,y2,lty=1,lwd=5,col=8);3 N7 ?5 z5 G2 Z( k6 c1 |
x3 <- c(x[45],x[96]);
# M. z6 E2 P0 e0 Q+ N$ G5 k3 y: p
y3 <- c(y[45],y[96]);1 J1 ]5 s$ Z2 W7 Q& C
lines(x3,y3,lty=1,lwd=5,col=8);
" B4 s% K5 O0 g/ i% \
axis(side=1,at=seq(10,100,by=10),labels=c(" ","0"," ","20"," ","40"," ","60"," ","80"));( H' A5 ^8 r5 p8 T6 b" ^% A
axis(side=1,col="red",col.axis="red",at=seq(4,10,by=1),labels=c("14"," "," "," "," "," "," "));$ o/ ^/ d4 p! f5 O; X# g
axis(side=1,col="red",col.axis="red",at=c(10),labels=c("20"));$ {( {+ n3 f/ }5 G/ Y
axis(side=2,at=seq(0,100,by=25),labels=c(" ","0.4","0.5","0.6","0.7"));9 S1 f2 k! q* |

复制代码