rm(list=ls())
require(VAM)
# Partie simulation
formSim <-formula( ~ (ARAInf(.4) | Weibull(.001,2.5)) & (ARA1(.7) + ARA1(.7) | Periodic(10,prob=c(.5,.5))))
simCpp <- sim.vam(formSim)
n=rep(30,1)
simDf=simulate(simCpp,n)
# Estimation un système
formMle <- update(formSim,Time & Type ~ .)
mleCpp <- mle.vam(formMle,data=simDf)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.