| maxPosterProb {glmmAK} | R Documentation |
For given G-spline basis and given data, it determines for each data point the G-spline component for which the value of the basis density is maximal.
maxPosterProb(data, intercept, std.dev, K, delta, sigma)
data |
numeric vector or matrix with data. If given as a matrix then rows correspond to observations and columns to margins. |
intercept |
numeric vector of length 1 or ncol(data) with
intercepts for each margin.
If given as a number, it is recycled. |
std.dev |
numeric vector of length 1 or ncol(data) with
standard deviations for each margin.
If given as a number, it is recycled. |
K |
numeric vector of length 1 or ncol(data) which
specifies, for each marginal G-spline, then number of knots
on each side of the zero knot. That is, the
i-th marginal G-spline has 2K[i]+1 knots.
If given as a number, it is recycled. |
delta |
numeric vector of length 1 or ncol(data)
which specifies the distance between two consecutive
knots for each marginal G-spline. That is, the i-th
marginal G-spline has the following knots
mu[i,j] = j*delta[i], j=-K[i],...,K[i].
|
sigma |
numeric vector of length 1 or ncol(data)
with basis standard deviations for marginal G-splines.
If given as a number, it is recycled. |
Matrix which specifies determined components (indeces are on scale -K[i],...,K[i]).
Arnošt Komárek arnost.komarek[AT]mff.cuni.cz
cumlogitRE, logpoissonRE.
N <- 100
intcpt <- c(0, 5, 15)
std.dev <- c(1, 0.5, 3)
data <- data.frame(b1=rnorm(N, intcpt[1], std.dev[1]),
b2=rnorm(N, intcpt[2], std.dev[2]),
b3=rnorm(N, intcpt[3], std.dev[3]))
alloc <- maxPosterProb(data=data, intercept=intcpt, std.dev=std.dev,
K=15, delta=0.3, sigma=0.2)
par(mfrow=c(1, 3), bty="n")
for (i in 1:3) hist(alloc[,i], prob=TRUE, col="seagreen3",
xlab="Allocation", breaks=(-15):15,
main=paste("Margin ", i, sep=""))