An NA error in my Gibbs sampler for mixture model

Question

I am working on a Gibbs sampler and my code is

 library(foreign)
 cognitive = read.dta("http://www.stat.columbia.edu/~gelman/arm/examples/child.iq/kidiq.