Rel 10.6.0 - Added test for Cholesky factor cov

Examples/Bernoulli/bernoulli2.jl

@@ -1,6 +1,6 @@
 ######### Stan program example  ###########
 
-using StanSample
+using StanSample, DataFrames
 
 function duplicate_tmpdir_for_windows_only(sm)
 

Project.toml

@@ -3,7 +3,7 @@ uuid = "682df890-35be-576f-97d0-3d8c8b33a550"
 license = "MIT"
 desc = "Illustrate StanJulia's interface options to cmdstan"
 authors = ["Rob J Goedman <[email protected]>"]
-version = "10.5.0"
+version = "10.6.0"
 
 [deps]
 AxisKeys = "94b1ba4f-4ee9-5380-92f1-94cde586c3c5"

test/Examples-Test-Cases/Cholesky-factor-cov/test_cholesky_factor_cov.jl

@@ -0,0 +1,35 @@
+using StanSample
+
+stan_data = (N=3, nu=13, L_Psi=[1.0 0.0 0.0; 2.0 3.0 0.0; 4.0 5.0 6.0]);
+
+model_code = "
+data {
+    int<lower=1> N;
+    real<lower=N-1> nu;
+    cholesky_factor_cov[N] L_Psi;
+}
+parameters {
+    cholesky_factor_cov[N] L_X;
+}
+model {
+    L_X ~ inv_wishart_cholesky(nu, L_Psi);
+}
+";
+
+
+sm = SampleModel("test", model_code);
+rc = stan_sample(sm; data=stan_data);
+
+if success(rc)
+    df = read_samples(sm, :dataframe)
+    ndf = read_samples(sm, :nesteddataframe)
+    display(ndf.L_X[1])
+    println()
+end
+
+for j in 1:5
+    run(pipeline(StanSample.par(sm.cmds), stdout=sm.log_file[1]));
+    ndf = read_samples(sm, :nesteddataframe)
+    display(ndf.L_X[1])
+    println()
+end    

test/runtests.jl

@@ -69,6 +69,7 @@ examples = [
   "Examples-Test-Cases/MatrixInput/input_dict_symbol.jl",
   "Examples-Test-Cases/MatrixInput/input_dict_string.jl",
   "Examples-Test-Cases/MatrixInput/input_nt.jl",
+  "Examples-Test-Cases/Cholesky-factor-cov/test_cholesky_factor_cov.jl",
 
   "Examples-Stan-Methods/Generate_Quantities/generate_quantities.jl",
   "Examples-Stan-Methods/Parse_and_Interpolate/parse.jl",