library(BayesianMonoErrorModels)
a1=read_csv(file='test/ERTEB.csv')
parms=parse_file(a1,logMIC=TRUE)
xobs=parms$xobs; yobs=parms$yobs; xcens=parms$xcens; ycens=parms$ycens
run_spline_model(xobs,yobs,xcens,ycens)
### Initialize
nobs=length(xobs)
xgrid=seq(min(xobs)-1,max(xobs)+1,length=1000)
parms=initialize_parms_spline(xobs,yobs,xcens,ycens,xgrid)
xtrue=parms$xtrue
coefs=parms$coefs
ytrue=as.numeric(parms$ytrue)
lowept=parms$lowept
upperept=parms$upperept
knotseq=parms$knotseq
bases=parms$bases
designMatrixGrid=parms$designMatrixGrid
### Run
parms=bayesian_mon_errors_spline_diag(xobs,yobs,xcens,ycens,coefs,xtrue,ytrue,xgrid,lowept,upperept,knotseq,bases,designMatrixGrid)
MICDens=parms$MICDens; fitMat=parms$fitMat; acceptCoef=parms$acceptCoef; smoothAccept=parms$smoothAccept
coefMat=parms$coefMat; xtrue_sav=parms$xtrue_sav; smoothParam_sav=parms$smoothParam_sav
spline_diagnostic(xgrid,xobs,yobs,xcens,ycens,MICDens,fitMat,acceptCoef,xtrue_sav,coefMat,smoothAccept,smoothParam_sav)
parms=bayesian_mon_errors_spline(xobs,yobs,xcens,ycens,coefs,xtrue,ytrue,xgrid,lowept,upperept,knotseq,bases,designMatrixGrid)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.