The error is near lines 496 508, in the Drinking water supply section. From the DHS, the code is:

compute h2oires=0.

if

variable labels h2oires "Public Piped into dwelling".

compute h2oyrd=0.

if

variable labels h2oyrd "Public Piped into yard/plot".

compute h2opires=0.

if

variable labels h2opires "Private Piped into dwelling".

compute h2opyrd=0.

if

variable labels h2opyrd "Private Piped into yard/plot".

You'll notice that when QH121 = 11 and 12, it's transformed twice into different dummy variables. The code should instead be:

compute h2oires=0.

if

variable labels h2oires "Public Piped into dwelling".

compute h2oyrd=0.

if

variable labels h2oyrd "Public Piped into yard/plot".

compute h2opires=0.

if

variable labels h2opires "Private Piped into dwelling".

compute h2opyrd=0.

if

variable labels h2opyrd "Private Piped into yard/plot".

This corresponds with the DHS survey question where 'Private Piped into dwelling' is coded as 13 and 'Private piped into yard/plot' is coded as 14.

Additionally, this error seems to have been incorporated into the factor analysis itself. The mean, standard deviation and component score coefficient for h2oyrd is identical to those for h2opyrd in the provided excel sheet. The same applies to h2oires and h2opires. This, of course, affects the values for the entire model.]]>

]]>