### Create a test coordinate file using ct2coord ###
ct=makeCt("((((...(((((((.........)))))))...((((.........))))...))))",
"AAAAAAAACCCCCCCCAAAGGGGGGGAUUACCCCUCCUUUAAAAGGGUUUUCCCCCCC"
)
coord=ct2coord(ct)
### add an id ###
coord$id=1
f = tempfile()
### write out test file ###
write.table( coord[,c('id','x','y','seq','num','bound')],
col.names=FALSE,row.names=FALSE,sep=",",file=f
)
### Read in the coordinate file ##
input=loadCoords(f)
### Plot the file using RNAPlot ##
RNAPlot(input)
Run the code above in your browser using DataLab