Running examples in ‘sn-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: selm
> ### Title: Fitting linear models with skew-elliptical error term
> ### Aliases: selm
> ### Keywords: regression univar multivariate
> 
> ### ** Examples
> 
> data(ais)
> m1 <- selm(log(Fe) ~ BMI + LBM, family="SN", data=ais)
> print(m1)
Object class: selm 
Call: selm(formula = log(Fe) ~ BMI + LBM, family = "SN", data = ais)
Number of observations: 202 
Number of covariates: 3 (counting constant term)
Number of parameters: 5 
Family: SN 
Estimation method: MLE
Log-likelihood: -175.9094 
> summary(m1)
Call: selm(formula = log(Fe) ~ BMI + LBM, family = "SN", data = ais)
Number of observations: 202 
Family: SN 
Estimation method: MLE
Log-likelihood: -175.9094 
Parameter type: CP 

CP residuals:
     Min       1Q   Median       3Q      Max 
-2.01862 -0.39918  0.03332  0.40744  1.33000 

Regression coefficients
               estimate  std.err  z-ratio Pr{>|z|}
(Intercept.CP) 2.736823 0.337376 8.112094    0.000
BMI            0.035270 0.020760 1.698934    0.089
LBM            0.009465 0.004467 2.118964    0.034

Parameters of the SEC random component
       estimate std.err
s.d.     0.5796   0.029
gamma1  -0.1742   0.180
> s<- summary(m1, "DP", cov=TRUE, cor=TRUE)
> plot(m1)
> plot(m1, param.type="DP")
> logLik(m1)
'log Lik.' -175.9094 (df=5)
> coef(m1)
(Intercept.CP)            BMI            LBM           s.d.         gamma1 
    2.73682314     0.03526978     0.00946493     0.57960516    -0.17419973 
> coef(m1, "DP")
(Intercept.DP)            BMI            LBM          omega          alpha 
    3.16595780     0.03526978     0.00946493     0.72117869    -1.11946607 
> var <- vcov(m1)
> #
> m1a <- selm(log(Fe) ~ BMI + LBM, family="SN", method="MPLE", data=ais)
> m1b <- selm(log(Fe) ~ BMI + LBM, family="ST", fixed.param=list(nu=8), data=ais)
> #
> data(barolo)
> attach(barolo)
> A75 <- (reseller=="A" & volume==75)
> logPrice <- log(price[A75],10) 
> m <- selm(logPrice ~ 1, family="ST")
Warning in st.infoUv(dp = fit$dp, NULL, x, yInfo, w, fixed.nu, symmetr,  :
  'dp' does not seem to be at MLE; score not quite 0
score(dp):  4.006562 2.255972 -1.030932 -1.301749 
norm(score)^2: 4.882399 
> summary(m)
CP does not exist. Consider param.type='DP' or 'pseudo-CP'
> plot(m, which=2, col=4, main="Barolo log10(price)")
Requested param.type='cp' evaluates to NULL.
CP of ST family requires nu>4
Error in .local(x, ...) : 
  Consider another choice of param.type (DP or pseudo-CP)
Calls: plot -> plot -> .local
Execution halted