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);/ X9 v7 Y o1 _* ~. [* s
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); 9 ^' d- b0 h8 V _
plot(x,y,pch=20,axes=T,xaxt="n",yaxt="n",ylab="methylation",xlab="Age(yrs)");! K, v' W9 h5 x
title(col.sub="red",sub="Fetal Age(wks) ",main="cg04640913"); - d; z. G2 C& w0 ~4 I8 M: G
text(98,15,"CDH22",font=3,cex=2); 4 p4 U4 S% i% N1 ~) k I
x1 <- c(x[2],x[22]);9 c- `7 j3 F! U$ } j; X
y1 <- c(y[2],y[22]); 5 T$ w5 n7 {. \1 q+ X
lines(x1,y1,lty=1,lwd=5,col=8);9 K8 I8 U& v" ? j1 C