Neighbors have three children. Given three independent observations of boy, what is the probability that they have 3 boys? I would really appreciate some help on this one. I'm completely lost. I have no idea why my method doesnt work.

"Your new neighbors have three children. If you are told about three
  independent observations of a boy, what is the probability that they
  have three boys?"

I solve it by saying that our sample space is $S = \{bbb,ggg,bgb,bbg,ggb,gbg,gbb,bgg\}$.
If we then introduce the event $A_i = \{$observation i of boy$\}$ and call the event $B=\{bbb\}$
then $P(A_i)=7/8 $ for $i=1,2,3$ and the events $A_i$ are independent according to the problem text. Furthermore, clearly $P(B) = 1/8$.
We now want to calculate the probability
$P(B|A_1 \cap A_2 \cap A_3)$ which is (because of total law of probability/Bayes formula) equal to
$\frac{P(B)P(A_1 \cap A_2 \cap A_3 | B)}{P(A_1 \cap A_2 \cap A_3)}$
Now, $P(A_1 \cap A_2 \cap A_3) = P(A_1)P(A_2)P(A_3)$ because of their independence, and I argue that $P(A_1 \cap A_2 \cap A_3 | B) = 1$ because the probability of observing a boy three times given that there only is three boys should be one hundred percent.
This simplifies to
$\frac{8^2}{7^3}$
which is clearly wrong. The answer should be 1/2, but i have no idea how to come up with that. 
What am I doing wrong?
Thanks!

We can restate the question as follows.  Suppose there are $3$ balls in an urn, and each ball is either black or white with equal probability.  We draw one ball from the urn at random, observe the color, and replace it in the urn.  Given that after $3$ draws we observed a black ball each time, what is the probability that the urn contains no white balls?
Let the number of white balls in the urn be the random variable $X$.  Then before any data is observed, $X$ is presumed to follow a binomial distribution with parameters $n = 3$ and $p = 1/2$.  Let $Y$ represent the number of white balls drawn out of $3$ independent trials with replacement.  Then we have $$\Pr[X = 0 \mid Y = 0] = \frac{\Pr[Y = 0 \mid X = 0]\Pr[X = 0]}{\Pr[Y = 0]}$$ by Bayes theorem.  Since the prior distribution of $X$ is binomial, we have $$\Pr[X = 0] = \binom{3}{0}(1/2)^0 (1 - 1/2)^{3-0} = \frac{1}{8}.$$  Given that there are no white balls, $\Pr[Y = 0 \mid X = 0] = 1$--we cannot draw a white ball if there is none in the urn.  The only remaining quantity to be determined is the unconditional probability of $Y = 0$, which by the law of total probability is simply $$\Pr[Y = 0] = \sum_{x=0}^3 \Pr[Y = 0 \mid X = x]\Pr[X = x].$$  We can already compute the first term of this sum, since it is just the numerator in the desired probability.  But the other conditional probabilities when $X \in \{1, 2\}$ require more thought.  Specifically, the conditional distribution of $Y \mid X = x$ is binomial with parameters $n = 3$ and $p = x/3$.  That $n = 3$ is obvious; $Y$ comprises the sum of the number of white balls obtained from three independent draws with replacement from the urn.  The probability of observing a white ball in any single draw is determined by the number of white balls in the urn, which is $x$:  so if $x = 0$, all balls are black and $p = 0$, reflecting the impossibility of drawing a white ball.  If $x = 1$, then $p = 1/3$ reflects a $1/3$ chance of drawing a white ball.  Thus $$\Pr[Y = y \mid X = x]\Pr[X = x] = \binom{3}{y}(x/3)^y (1 - x/3)^{3-y} \cdot \binom{3}{x} (1/2)^x (1 - 1/2)^{3-x} = \frac{1}{216} \binom{3}{x}\binom{3}{y} x^y (3-x)^{3-y}.$$
With this in mind, we compute $$\Pr[Y = 0] = \frac{1}{216} \sum_{x=0}^3 \binom{3}{0} \binom{3}{x} x^0 (3-x)^3 = \frac{54}{216} = \frac{1}{4}.$$  It follows that $$\Pr[X = 0 \mid Y = 0] = \frac{1/8}{1/4} = \frac{1}{2}.$$
From here, the sophisticated student can easily generalize to the case where we are given $n$ balls in the urn, each with a prior probability $p$ of being white, and $m$ iid draws from the urn, in which $y$ are observed white.  Then the posterior distribution of the number of white balls is $$\Pr[X = x \mid Y = y] = \frac{\Pr[Y = y \mid X = x]\Pr[X = x]}{\Pr[Y = y]} = \frac{\binom{m}{y} (x/n)^y (1 - x/n)^{m-y} \binom{n}{x} p^x (1-p)^{n-x}}{\sum_{k=0}^n \binom{m}{y} (k/n)^y (1 - k/n)^{m-y} \binom{n}{k} p^k (1-p)^{n-k}}.$$
