powerForQual = function(SampSize, ThemePrevalence, NumInstances) {
pnbinom(SampSize-NumInstances, size = NumInstances, prob = ThemePrevalence)
}
sampSizeForQual = function(Power, ThemePrevalence, NumInstances) {
qnbinom(Power, size = NumInstances, prob = ThemePrevalence) + NumInstances
}
powerForQual(30, 0.1, 2)
sampSizeForQual(0.8, 0.1, 2)
cG93ZXJGb3JRdWFsID0gZnVuY3Rpb24oU2FtcFNpemUsIFRoZW1lUHJldmFsZW5jZSwgTnVtSW5zdGFuY2VzKSB7CiAgcG5iaW5vbShTYW1wU2l6ZS1OdW1JbnN0YW5jZXMsIHNpemUgPSBOdW1JbnN0YW5jZXMsIHByb2IgPSBUaGVtZVByZXZhbGVuY2UpCn0KCnNhbXBTaXplRm9yUXVhbCA9IGZ1bmN0aW9uKFBvd2VyLCBUaGVtZVByZXZhbGVuY2UsIE51bUluc3RhbmNlcykgewogIHFuYmlub20oUG93ZXIsIHNpemUgPSBOdW1JbnN0YW5jZXMsIHByb2IgPSBUaGVtZVByZXZhbGVuY2UpICsgTnVtSW5zdGFuY2VzCn0KCnBvd2VyRm9yUXVhbCgzMCwgMC4xLCAyKQoKc2FtcFNpemVGb3JRdWFsKDAuOCwgMC4xLCAyKQo=