@@ -904,11 +904,11 @@ function mle_fs!{T1<:VarianceComponentModel, T2<:TwoVarCompVariateRotate}(
904904 Bcov = zeros (T, nmean, nmean)
905905 Bcov = pinv (fisher_B (vcmodel, vcdatarot, dd. vcaux))
906906 Bse = similar (vcmodel. B)
907- copy! (Bse, sqrt (diag (Bcov)))
907+ copy! (Bse, sqrt . (diag (Bcov)))
908908 Σcov = pinv (fisher_Σ (vcmodel, vcdatarot))
909909 Σse = (zeros (T, d, d), zeros (T, d, d))
910- copy! (Σse[1 ], sqrt (diag (view (Σcov, 1 : d^ 2 , 1 : d^ 2 ))))
911- copy! (Σse[2 ], sqrt (diag (view (Σcov, d^ 2 + 1 : 2 d^ 2 , d^ 2 + 1 : 2 d^ 2 ))))
910+ copy! (Σse[1 ], sqrt . (diag (view (Σcov, 1 : d^ 2 , 1 : d^ 2 ))))
911+ copy! (Σse[2 ], sqrt . (diag (view (Σcov, d^ 2 + 1 : 2 d^ 2 , d^ 2 + 1 : 2 d^ 2 ))))
912912
913913 # output
914914 maxlogl, vcmodel, Σse, Σcov, Bse, Bcov
@@ -1105,11 +1105,11 @@ function mle_mm!{T1 <: VarianceComponentModel, T2 <: TwoVarCompVariateRotate}(
11051105 # standard errors
11061106 Bcov = pinv (fisher_B (vcm, vcdatarot, vcaux))
11071107 Bse = similar (vcm. B)
1108- copy! (Bse, sqrt (diag (Bcov)))
1108+ copy! (Bse, sqrt . (diag (Bcov)))
11091109 Σcov = pinv (fisher_Σ (vcm, vcdatarot))
11101110 Σse = (zeros (T, d, d), zeros (T, d, d))
1111- copy! (Σse[1 ], sqrt (diag (view (Σcov, 1 : d^ 2 , 1 : d^ 2 ))))
1112- copy! (Σse[2 ], sqrt (diag (view (Σcov, d^ 2 + 1 : 2 d^ 2 , d^ 2 + 1 : 2 d^ 2 ))))
1111+ copy! (Σse[1 ], sqrt . (diag (view (Σcov, 1 : d^ 2 , 1 : d^ 2 ))))
1112+ copy! (Σse[2 ], sqrt . (diag (view (Σcov, d^ 2 + 1 : 2 d^ 2 , d^ 2 + 1 : 2 d^ 2 ))))
11131113
11141114 # output
11151115 logl, vcm, Σse, Σcov, Bse, Bcov
@@ -1240,7 +1240,7 @@ function fit_reml!{
12401240 # standard errors and covariance of mean parameters
12411241 Bcov = inv (fisher_B (vcmodel, vcdatarot, vcaux))
12421242 Bse = similar (vcmodel. B)
1243- copy! (Bse, sqrt (diag (Bcov)))
1243+ copy! (Bse, sqrt . (diag (Bcov)))
12441244
12451245 # output
12461246 logpdf (vcmodel, vcdatarot), vcmodel, Σse, Σcov, Bse, Bcov
0 commit comments