I have a binary outcome (success \ failure). I have 20 subjects, half of them gave 2 samples, the others just one, so I have 30 data points.

All data points without any exceptions were success (1). So my point estimate of the success rate is clearly 100%. Now I want to calculate a CI for the success rate, mainly to see what the lower limit is, I want to be able to say that with 95% confidence the success rate is higher than....(80%, 85%, whatever comes up).

The problem is, as you can see, is the clustering, I can't use n = 30 because I then ignore the correlation, I feel it is a waste to use n = 20. Is there a way to calculate the variance of a single sample proportion while taking into account the cluster ?

Thanks !