Description Usage Arguments Details Value Note References See Also Examples

This function, when used with `da.mix`

or
`dabipf.mix`

, can be
used to create proper multiple imputations of missing data under
the general location model with or without restrictions.

1 | ```
imp.mix(s, theta, x)
``` |

`s` |
summary list of an incomplete data matrix |

`theta` |
value of the parameter under which the missing data are to be
randomly imputed. This is a parameter list such as one created
by |

`x` |
the original data matrix used to create the summary list |

This function is essentially the I-step of data augmentation.

a matrix of the same form as `x`

, but with all missing values filled in
with simulated values drawn from their predictive distribution given
the observed data and the specified parameter.

The random number generator seed must be set at least once by the
function `rngseed`

before this function can be used.

Schafer, J. L. (1996) *Analysis of Incomplete Multivariate Data.*
Chapman \& Hall, Chapter 9.

`prelim.mix`

, `da.mix`

,
`dabipf.mix`

, `rngseed`

1 2 3 4 5 6 | ```
data(stlouis)
s <- prelim.mix(stlouis,3) # do preliminary manipulations
thetahat <- em.mix(s) # ML estimate for unrestricted model
rngseed(1234567) # set random number generator seed
newtheta <- da.mix(s,thetahat,steps=100) # data augmentation
ximp <- imp.mix(s, newtheta, stlouis) # impute under newtheta
``` |

```
Steps of EM:
1...2...3...4...5...6...7...8...9...10...11...12...13...14...15...16...17...18...19...20...21...22...23...24...25...26...27...28...29...30...31...32...33...34...35...36...37...38...39...40...41...42...43...44...45...46...47...48...49...50...51...52...53...54...55...56...57...58...59...60...61...62...63...64...65...66...67...68...69...70...71...72...73...74...75...76...77...78...79...80...81...82...83...84...85...86...87...88...89...90...91...92...93...94...95...96...97...98...99...100...101...102...103...104...105...106...107...108...109...110...111...112...113...114...115...116...117...118...119...120...121...122...123...124...125...126...127...128...129...130...131...132...133...134...135...136...137...138...139...140...141...142...143...144...145...146...147...148...149...150...151...152...153...154...155...156...157...158...159...160...161...162...163...164...165...166...167...168...169...170...171...172...173...174...175...176...177...178...179...180...181...
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.