df <- data.frame(x = c(seq(1, 10, by = 1), seq(12, 30, by = 2)),
                 y = c(0,0.1,0.19,0.3,0.4,0.6,0.68,
                       0.78,0.8,0.815,0.88,0.9,0.91,
                       0.92,0.91,0.93,0.93,0.92,0.915,0.93))
plot(df$x,df$y,pch = 16,cex=0,
     xlab="Soft Power Threshold",
     ylab="Scale Free Topology Model Fit,signed R^2",type="n",
     main = paste("Scale independence"),yaxt='n')
abline(h=0.8,col="black",lty = 2)
abline(v=9,col="black",lty = 2)
points(x = 9, y = 0.8,pch = 16,cex=3)
text(df$x,df$y,
     labels=df$x,cex=1,col="black")
text(x = 9, y = 0.8,labels=9,cex=1,col="white")
axis(side = 2,at = c(0,0.2,0.4,0.6,0.8,1),labels =  c(0,0.2,0.4,0.6,0.8,1))

df <- data.frame(x = c(seq(1, 10, by = 1), seq(12, 30, by = 2)),
                 y = c(0,0.1,0.19,0.3,0.4,0.6,0.68,
                       0.78,0.8,0.815,0.88,0.9,0.91,
                       0.92,0.91,0.93,0.93,0.92,0.915,0.93))
plot(df$x,df$y,pch = 16,cex=0,
     xlab="Soft Power Threshold",
     ylab="Scale Free Topology Model Fit,signed R^2",type="n",
     main = paste("Scale independence"),yaxt='n')
abline(h=0.8,col="black",lty = 2,lwd = 2)
abline(v=9,col="black",lty = 2,lwd = 2)
axis(side = 2,lwd = 2,at = c(-1,3))
axis(side = 1,lwd = 2,at = c(-1,300))
axis(side = 4,lwd = 2,at = c(-1,3))
axis(side = 3,lwd = 2,at = c(-1,300))

points(x = 9, y = 0.8,pch = 16,cex=3)
text(df$x,df$y,
     labels=df$x,cex=1,col="black")
text(x = 9, y = 0.8,labels=9,cex=1,col="white")
axis(side = 2,at = c(0,0.2,0.4,0.6,0.8,1),
     labels =  c(0,0.2,0.4,0.6,0.8,1))

df <- data.frame(x = c(seq(1, 10, by = 1), seq(12, 30, by = 2)),
                 y = c(0,0.1,0.19,0.3,0.4,0.6,0.68,
                       0.78,0.8,0.815,0.88,0.9,0.91,
                       0.92,0.91,0.93,0.93,0.92,0.915,0.93))
plot(df$x,df$y,pch = 16,cex=0,
     xlab="Soft Power Threshold",
     ylab="Scale Free Topology Model Fit,signed R^2",
     type="n",
     main = paste("Scale independence"),
     yaxt='n',
     ylim = c(0,1))
abline(h=0.8,col="black",lty = 2,lwd = 2)
abline(v=9,col="black",lty = 2,lwd = 2)
axis(side = 2,lwd = 2,at = c(-1,3))
axis(side = 1,lwd = 2,at = c(-1,300))
axis(side = 4,lwd = 2,at = c(-1,3))
axis(side = 3,lwd = 2,at = c(-1,300))

points(x = 9, y = 0.8,pch = 16,cex=3,col = 'black')
text(df$x,df$y,
     labels=df$x,cex=1,col="black")
text(x = 9, y = 0.8,labels=9,cex=1,col="white")
axis(side = 2,at = c(0,0.2,0.4,0.6,0.8,1),
     labels =  c(0,0.2,0.4,0.6,0.8,1))