sn.2logL.profile(X=matrix(rep(1, n)), y,
param.range=c(sqrt(var(y)) * c(2/3, 3/2), -0.95, 0.95),
use.cp=TRUE, npts=31 %/% d, plot.it=TRUE, ...)
NA
s) are not allowed.col(X)
equal to
length(y). Missing values (`NA
s) are not allowed.
If X
is missing, a one-comun matrix of 1's is created.d
is set to 1, and
a 1-dimensional profile is computed and plotted, for the shape
or skewness parameter (depending on the parametrization adopted;
see belouse.cp=TRUE
(default value), the centred parametrization is used,
otherwise the direct parametrization is adopted.T
.sn.em
.d=2
, the second vector contains NA
s.param1
and param2
.2logL
.plot.it=TRUE
, a plot of the profile twice relative loglikeliood is
produced on a graphical device.sn.em
.See the reference below for explanation of the two possible parametrizations.
sn.em
, sn.mle
data(ais, package="sn")
attach(ais)
a <- sn.2logL.profile(y=bmi)
a <- sn.2logL.profile(y=bmi, use.cp=FALSE, param.range=c(3,6,1,5))
a <- sn.2logL.profile(X=cbind(1,lbm), y=bmi, param.range=c(0,0.9),
npts=50)
#
data(frontier, package="sn")
a <- sn.2logL.profile(y=frontier, param.range=c(0.8,1.6,10,30),
use.cp=FALSE, npts=11)
Run the code above in your browser using DataLab