Skip to content
GitLab
Projects Groups Topics Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • mvnchol mvnchol
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributor statistics
    • Graph
    • Compare revisions
  • Issues 0
    • Issues 0
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 0
    • Merge requests 0
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Artifacts
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Container Registry
    • Model experiments
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar

Admin message

Until further notice: Mailroom does not work. Replies to git-generated emails will be silently ignored!

  • Thorsten Simon
  • mvncholmvnchol
  • Issues
  • #3
Closed
Open
Issue created Nov 15, 2021 by Thorsten Simon@c4031039Owner

Enable weights in startfun (dist_modchol)

Trying to check modchol disttree family:

devtools::load_all()
  
data(simdata)
y <- simdata$y
d <- simdata$d

df <- dist_mvn_modchol(3)
fit_dist <- df$startfun(y)

library(bamlss)
bfam <- mvnchol_bamlss(3, "modified")

# check changes in mu

eta_to_par <- function(eta, y) {
    par <- list()
    for (i in seq_along(eta)) par[[i]] <- rep(eta[[i]], nrow(y))
    names(par) <- names(fitted(b))
    for (i in 4:6) par[[i]] <- rep(exp(eta[[i]]), nrow(y))
    par
}

par <- eta_to_par(fit_dist, y)

xpar <- par

x <- runif(nrow(y), -1, 1)
xpar[["mu2"]] <- ifelse(x > 0, xpar[["mu2"]], xpar[["mu2"]] + .5)

yy <- bfam$r(xpar)

library(disttree)
dt <- disttree(yy ~ x, family = df)

But startfun() needs weights here: https://git.uibk.ac.at/c4031039/mvnchol/-/blob/master/R/dist_modchol.R#L134

Edited Nov 15, 2021 by Thorsten Simon
Assignee
Assign to
Time tracking