"p"=c(p,1-sum(p)), "beta"=beta, "b"=b, "weights"=weights)
save("mr", "mr_params", file=paste("res_",n,"_",d,"_",link,"_",strw,".RData",sep=""))
"p"=c(p,1-sum(p)), "beta"=beta, "b"=b, "weights"=weights)
save("mr", "mr_params", file=paste("res_",n,"_",d,"_",link,"_",strw,".RData",sep=""))