R: Provide summary of parameter estimates and standard error...
summaryParam
R Documentation
Provide summary of parameter estimates and standard error across replications
Description
This function will provide averages of parameter estimates, standard deviations of parameter estimates, averages of standard errors, and power of rejection with a priori alpha level for the null hypothesis of parameters equal 0.
The alpha level used to find the statistical power of each parameter estimate
std
If TRUE, (a) standardized coefficients and their standard errors substitute unstandardized coefficients, (b) standardized parameter values substitute parameter values, (c) confidence intervals of standardized coefficients are calculated using Wald confidence interval, and (d) all results (e.g., biases or coverage) are calculated based on standardized coefficients.
detail
If TRUE, more details about each parameter estimate are provided, such as relative bias, standardized bias, or relative standard error bias.
improper
If TRUE, include the replications that provided improper solutions
digits
The number of digits rounded in the result. If NULL, the results will not be rounded.
matchParam
If TRUE, only parameter estimates that have the same names as the parameter values will be reported. This argument is recommended when users know that the data-generation model and analysis model are the same. Then the comparison between the parameter estimates and parameter value will be valid.
Value
A data frame that provides the statistics described above from all parameters.
For using with linkS4class{SimResult}, each column means
Estimate.Average: Average of parameter estimates across all replications
Estimate.SD: Standard Deviation of parameter estimates across all replications
Average.SE: Average of standard errors across all replications
Power (Not equal 0): Proportion of significant replications when testing whether the parameters are different from zero. The alpha level can be set by the alpha argument of this function.
Average.Param: Parameter values or average values of parameters if random parameters are specified
SD.Param: Standard Deviations of parameters. Show only when random parameters are specified.
Average.Bias: The difference between parameter estimates and parameter underlying data
SD.Bias: Standard Deviations of bias across all replications. Show only when random parameters are specified.
This value is the expected value of average standard error when random parameter are specified.
Coverage: The percentage of (1-alpha)% confidence interval covers parameters underlying the data.
Rel.Bias: Relative Bias, which is (Estimate.Average - Average.Param)/Average.Param.
Hoogland and Boomsma (1998) proposed that the cutoff of .05 may be used for acceptable relative bias.
This option will be available when detail=TRUE. This value will not be available when parameter values are very close to 0.
Std.Bias: Standardized Bias, which is (Estimate.Average - Average.Param)/Estimate.SD
for fixed parameters and (Estimate.Average - Average.Param)/SD.Bias
for random parameters. Collins, Schafer, and Kam (2001) recommended that biases will be
only noticeable when standardized bias is greater than 0.4 in magnitude.
This option will be available when detail=TRUE
Rel.SE.Bias: Relative Bias in standard error, which is (Average.SE - Estimate.SD)/Estimate.SD
for fixed parameters and (Average.SE - SD.Bias)/SD.Bias
for random parameters. Hoogland and Boomsma (1998) proposed that 0.10 is the acceptable level.
This option will be available when detail=TRUE
Not Cover Below: The percentage of (1-alpha)% confidence interval does not cover the parameter and the parameter is below the confidence interval.
Not Cover Above: The percentage of (1-alpha)% confidence interval does not cover the parameter and the parameter is above the confidence interval.
Average CI Width: The average of (1-alpha)% confidence interval width across replications.
SD CI Width: The standard deviation of (1-alpha)% confidence interval width across replications.
Collins, L. M., Schafer, J. L., & Kam, C. M. (2001). A comparison of inclusive and restrictive strategies in modern missing data procedures. Psychological Methods, 6, 330-351.
Hoogland, J. J., & Boomsma, A. (1998). Robustness studies in covariance structure modeling. Sociological Methods & Research, 26, 329-367.
See Also
SimResult for the object input
Examples
showClass("SimResult")
loading <- matrix(0, 6, 1)
loading[1:6, 1] <- NA
LY <- bind(loading, 0.7)
RPS <- binds(diag(1))
RTE <- binds(diag(6))
CFA.Model <- model(LY = LY, RPS = RPS, RTE = RTE, modelType="CFA")
# We make the examples running only 5 replications to save time.
# In reality, more replications are needed.
Output <- sim(5, n=500, CFA.Model)
# Summary of the parameter estimates
summaryParam(Output)
# Summary of the parameter estimates with additional details
summaryParam(Output, detail=TRUE)
providing 
.