test/test_spline_model.R

In gdepalma/BayesianMonoErrorModels:

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)