R-based script used to simulate patterns of co-occurrence among two hypothetical species under patterns of avoidance, aggregation, and independence when the closure assumption is met or not, analyze simulated patterns with two-species occupancy models, and summarize the results.