Hi all,
I am trying to run the following command to predict yhat and I am only getting missing values for yhat. Please help
the code I am running is

probit c $N $Z , r
foreach z of global Z {
gen copy_`z'=`z'
qui sum `z'
replace `z' = r(mean)
}
predict yhat,xb


the data looks like this after usign the mean for the global Z :
input byte c float(age mh)
0 25 68.35207
0 25 68.35207
0 25 68.35207
0 24 68.35207
0 22 68.35207
0 24 68.35207
0 23 68.35207
0 22 68.35207
0 25 68.35207
0 22 68.35207
0 24 68.35207
0 26 68.35207
0 22 68.35207
0 25 68.35207
0 26 68.35207
0 24 68.35207
0 26 68.35207
0 25 68.35207
0 24 68.35207
0 24 68.35207
0 25 68.35207
0 27 68.35207
0 22 68.35207
0 24 68.35207
0 26 68.35207
0 26 68.35207
0 24 68.35207
0 26 68.35207
0 22 68.35207
0 24 68.35207
0 24 68.35207
0 25 68.35207
0 24 68.35207
0 24 68.35207
0 22 68.35207
0 25 68.35207
0 25 68.35207
0 26 68.35207
0 26 68.35207
0 26 68.35207
0 24 68.35207
0 24 68.35207
0 24 68.35207
0 24 68.35207
0 26 68.35207
0 26 68.35207
0 22 68.35207
0 25 68.35207
0 24 68.35207
0 26 68.35207
0 26 68.35207
0 23 68.35207
0 22 68.35207
0 22 68.35207
0 24 68.35207
0 22 68.35207
0 26 68.35207
0 25 68.35207
0 22 68.35207
0 25 68.35207
0 25 68.35207
0 24 68.35207
0 26 68.35207
0 26 68.35207
0 25 68.35207
0 24 68.35207
0 26 68.35207
0 22 68.35207
0 26 68.35207
0 24 68.35207
0 23 68.35207
0 22 68.35207
0 24 68.35207
0 23 68.35207
0 26 68.35207
0 26 68.35207
0 25 68.35207
0 23 68.35207
0 26 68.35207
0 23 68.35207
0 23 68.35207
0 23 68.35207
0 26 68.35207
0 26 68.35207
0 24 68.35207
0 25 68.35207
0 26 68.35207
0 26 68.35207
0 26 68.35207
0 23 68.35207
0 26 68.35207
0 22 68.35207
0 27 68.35207
0 26 68.35207
0 25 68.35207
0 24 68.35207
0 24 68.35207
0 22 68.35207
0 24 68.35207
0 22 68.35207
end
Any help is appreciated. Thanks
I think it is because of overdispersion in c but i may be wrong. Thanks.

Kind Regards,
Aarushi