The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

SLEEV: Semiparametric Likelihood Estimation with Errors in Variables

License: file /inst/LICENSE_sleev.txt

The complete R package logreg2ph and code for the simulation settings included in this paper.

Part of the R package TwoPhaseReg as described in this paper.

Below is a simple example for running a linear regression with sleev. The complete vignette detailing all uses of example in sleev can be found as an article.

Install

To install the package for GitHub, run the following in your R console:

  devtools::install_github("dragontaoran/sleev")

or install from CRAN:

install.packages("sleev")

Load package:

  library(sleev)

mock.vccc data

This data is a simulated two-phase sampling dataset based on VCCC data. To load this dataset, run:

  data(mock.vccc)

Preprocessing

Because of skewness, we often transform both CD4 and VL. In our analysis, CD4 was divided by 10 and square-root transformed and VL was log_10 transformed:

  mock.vccc$CD4_val_sq10 <- sqrt(mock.vccc$CD4_val/10)
  mock.vccc$CD4_unval_sq10 <- sqrt(mock.vccc$CD4_unval/10)
  mock.vccc$VL_val_l10 <- log10(mock.vccc$VL_val)
  mock.vccc$VL_unval_l10 <- log10(mock.vccc$VL_unval)

Load splines package

  library(splines)

Construct B-spline Basis

sn=20
data.linear <- spline2ph(x = "VL_unval_l10", data = mock.vccc, size = sn,
                         degree = 3,  group = "Sex")

Model fitting

The SMLEs can be obtained by running

res_linear <- linear2ph(y_unval = "CD4_unval_sq10", y = "CD4_val_sq10", 
                        x_unval = "VL_unval_l10",x = "VL_val_l10",z = "Sex", 
                        b_spline = paste0("bs", 1:sn), data = data.linear,
                        hn_scale = 1, se = TRUE, tol = 1e-04, 
                        max_iter = 1000, verbose = FALSE)

Check convergence:

c(res_linear$converge, res_linear$converge_cov)

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.