odds=function(d) d[d %% 2==1]
x=odds(as.integer(runif(500,1,100)))
nodes=rep(0,length(x)*2)
nodes[seq(1,length(x)*2,2)]=x
nodes[seq(2,length(x)*2,2)]=sample(x+1)
g=make_bipartite_graph(rep(c(0,1),length(x)),nodes, directed = T)
g=simplify(induced_subgraph(g, V(g)[degree(g)>0]))
png("g.png",2000,2000)
plot(g, layout=matrix(c(rep(c(1,2),50),rep(1:50,2)),ncol=2), vertex.size=5)
dev.off()
v=page_rank(g)$vector
d=degree(g,mode="in")
cor(v[v>0.002],d[d>0])