Description Usage Arguments Examples
View source: R/CreateVarCovarMix99.R
Read csv-file of the variance component estimates and construct a matrix with this estimates by parts. This matrix is checked if it is positive definite or not. If this matrix is not positive definite, 2 functions may be called. If you put the paramter psOptionRatio = FALSE, makePD2() is called . For the parameter psOptionRatio = TRUE the function make_pd_rat_ev() is called. The difference is, that with make_pd_rat_ev() you have the option to whish a maximum ratio between largest and smallest eigenvalue.
The software mix99 need a paramter-File with the variance and covariance. So this function is building this paramter-file.
1 2 3 4 5 6 7 8 9 10 11  | 
ps_input_file | 
 input csv-file  | 
ps_output_file | 
 output txt-file for the bended (co-)variances  | 
ps_raw_output_file | 
 output txt-file for the raw (co-)variances  | 
pn_ratio | 
 maximum ratio between largest and smallest eigenvalue, determines whether ratio method is used or not  | 
pn_eps | 
 lower limit of smallest eigenvalue  | 
pmat_weight | 
 weight matrix for weighted bending  | 
pn_digits | 
 number of digits to be rounded to  | 
pb_log | 
 indicator whether logs should be produced  | 
plogger | 
 log4r logger object  | 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15  | ## Not run: 
# use test input csv-file
sInputFile <- system.file("extdata","VCE_results.csv", package = "rvcetools")
# bend matrices using Schaeffer method
parameter_varCovar_mix99(ps_input_file = sInputFile, ps_output_file = 'par_varCovar_mix99_schaeffer.txt')
# bend matrices using ratio method
parameter_varCovar_mix99(ps_input_file  = sInputFile, 
                         ps_output_file = 'par_varCovar_mix99_ratio.txt',
                         pn_ratio       = 100)
# unweighted bending with minimum eigenvalue
parameter_varCovar_mix99(ps_input_file  = sInputFile, 
                         ps_output_file = 'par_varCovar_mix99_unweighted.txt',
                         pn_eps         = 1e-4)
## End(Not run)
 | 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.