Nothing
#' @export
lambdak_gmonl <- function(times, subj, X, y, d, tau, kn, degree, lambda, gam,omega,mv){
dim = length(subj)
X = matrix(X, nrow = dim)
px = ncol(X)
qvc2_indiv = qrvcp_gmonl(times, subj, y, X, tau, kn = kn,
degree = degree, lambda = 0, d,omega,rep(1,px),mv)$hat_bt
range = NULL
for (k in 1:px) {
range[k] =(max(qvc2_indiv[seq((k -
1) * dim + 1, k * dim)]) - min(qvc2_indiv[seq((k -
1) * dim + 1, k * dim)]))^(-(gam))
}
lambdasic = Lamb_gmonl(times, subj, X, y, d, tau, kn, degree,
lambda,omega,range,mv)$lambdasic
out = list(lambdasic = lambdasic, range = range)
return(out)
}
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.