coraallensavietta / RankingMethods

A Weighted Loss-Based Bayesian Ranking Method
0 stars 0 forks source link

running simulation efficiently #10

Open coraallensavietta opened 6 years ago

coraallensavietta commented 6 years ago

could i make this more efficient by running PostSamples only once per simulation? This seems to be the most time-consuming step.

coraallensavietta commented 6 years ago

yes, decided to do this. Need to rewrite small parts of RunSimulation to do this.

Main Questions for this Experiment

increase gaps in parameters until they're trivial. decrease until broken/impossible.

think of this as an experiment. what experimental conditions do we need to run to get a sense of behavior?

think about orders of magnitude to start (want low, medium, high)

Simulations TODO

Ranking Function Parameters:

Data Parameters

for below, think about how big average n should be to get p in CI n min: c(50) explore this to see . what differences in p1 and p2 coudl you find with this n? do a power calculation basically n max: se = sqrt(pq/n). think about how small this would have to be to think about how big n needs to be se = sqrt((.1.9)/20) #graph this out to see. If the ranking is always right, then n is too big. If ranks are uniform/terrible, then n too small. n assignment method

coraallensavietta commented 6 years ago

When I run: ranking_experiment_run_server.r Error in RunSimulation(N = n, a_p = 1, b_p = 1, n_min = n_min, n_max = n_max, : (list) object cannot be coerced to type 'double' Calls: rbind -> RunSimulation In addition: Warning message: 'rBind' is deprecated. Since R version 3.2.0, base's rbind() should work fine with S4 objects Execution halted

Trying ranking_experiment_run_server_n100.r now. MIght be more up to date. (really need to reorganize this mess!). Update: Same error message as above.

coraallensavietta commented 6 years ago

Ron suggests, "For the purposes of the simulation study, we probably don’t need to run full Bayesian inference in Stan. An empirical Bayes approach using lme4 should run much more quickly." I'll try to implement an alternate PostSamples function that implements this. started with d56774d