Nothing
logLikCR <- function(data, bw) {
p1a <- nrow(data$Sig) - 1
b <- as.vector(bw[1:p1a])
w <- bw[p1a+1]
sigma <- exp(data$W%*%data$tau + w)
sum((data$Y - data$X%*%data$beta - data$Z%*%b)^2/(2*sigma) + 0.5*log(sigma)) +
data$CH01*exp(data$X2%*%data$gamma1 + data$alpha1%*%b + data$nu1*w) +
data$CH02*exp(data$X2%*%data$gamma2 + data$alpha2%*%b + data$nu2*w) +
t(bw)%*%solve(data$Sig)%*%bw/2
}
logLikCR.learn <- function(data, bw) {
p1a <- nrow(data$Sig) - 1
b <- as.vector(bw[1:p1a])
w <- bw[p1a+1]
sigma <- exp(data$W%*%data$tau + w)
total <- sum((data$Y - data$X%*%data$beta - data$Z%*%b)^2/(2*sigma) + 0.5*log(sigma)) +
data$CH01*exp(data$X2%*%data$gamma1 + data$alpha1%*%b + data$nu1*w) +
data$CH02*exp(data$X2%*%data$gamma2 + data$alpha2%*%b + data$nu2*w) +
t(bw)%*%solve(data$Sig)%*%bw/2
if (data$D == 1) {
total <- total - log(data$HAZ01) - (data$X2%*%data$gamma1 + data$alpha1%*%b + data$nu1*w)
total
} else if (data$D == 2) {
total <- total - log(data$HAZ02) - (data$X2%*%data$gamma2 + data$alpha2%*%b + data$nu2*w)
total
} else {
total
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.