
normal1(lmean = "identity", lsd = "loge",
emean = list(), esd = list(), imethod = 1, zero = -2)
Links
for more choices.
Being a positive quantity, a log link is the default for the
standard deviation.earg
in Links
for general information.CommonVGAMffArguments
for more information.
If lmean = loge
then try imethod = 2
."vglmff"
(see vglmff-class
).
The object is used by modelling functions such as vglm
,
and vgam
.gaussianff
,
posnormal1
,
mix2normal1
,
Qvar
,
tobit
,
cennormal1
,
fnormal1
,
skewnormal1
,
dcennormal1
,
huber
,
studentt
,
binormal
,
dnorm
.ndata <- data.frame(x2 = rnorm(nn <- 200))
ndata <- transform(ndata, y = rnorm(nn, mean = 1-3*x2, sd = exp(1+0.2*x2)))
fit <- vglm(y ~ x2, normal1(zero = NULL), ndata, trace = TRUE)
coef(fit, matrix = TRUE)
# Generate data from N(mu = theta = 10, sigma = theta) and estimate theta.
theta <- 10
ndata <- data.frame(y = rnorm(100, m = theta, sd = theta))
fit <- vglm(y ~ 1, normal1(lsd = "identity"), ndata,
constraints = list("(Intercept)" = rbind(1, 1)))
coef(fit, matrix = TRUE)
Run the code above in your browser using DataLab