# NOT RUN {
#genome-like graph structure
sim1 <- simgeno(alpha = 0.01, beta = 0.02)
plot(sim1)
#genome-like graph structure with more edges between variables in a same or different groups
sim2 <- simgeno(adjacent = 3, alpha = 0.02 , beta = 0.03)
plot(sim2)
# }
# NOT RUN {
#simulate data
D <- simgeno(p=50, n=100, g=5, k= 3, adjacent = 3, alpha = 0.06 , beta = 0.08)
plot(D)
#Reconstructing intra- and inter-chromosomal conditional interactions (LD) network
out <- netsnp(data = D$data, n.rho= 4, ncores= 1)
plot(out)
#Select an optimal graph
sel <- selectnet(out)
plot(sel, vis= "complete" )
# }
Run the code above in your browser using DataCamp Workspace