Nothing
library(robustloggamma)
do.test <- TRUE
if (!do.test) {
n <- 200
p <- 4
suppressWarnings(RNGversion("3.5.0"))
set.seed(1234)
X <- matrix(rnorm(p*n, sd=2), ncol=p)
IX <- cbind(1,X)
mu <- 2
beta <- rep(2,p)
sigma <- 2
lambda <- 1
linear <- mu + X%*%beta
y <- rloggamma(n=n, mu=linear, sigma=sigma, lambda=lambda)
cens <- rloggamma(n=n, mu=linear+2.5, sigma=sigma, lambda=lambda)
delta <- as.numeric(y <= cens)
y[delta==0] <- cens[delta==0]
set.seed(2345)
SNP <- robustloggamma:::loggammacenslmrob.S(IX,y,delta,control=robustloggamma:::loggammacenslmrob.S.control())
dput(SNP, file='Snpcensreg-savedvalues.R')
} else {
n <- 200
p <- 4
suppressWarnings(RNGversion("3.5.0"))
set.seed(1234)
X <- matrix(rnorm(p*n, sd=2), ncol=p)
IX <- cbind(1,X)
mu <- 2
beta <- rep(2,p)
sigma <- 2
lambda <- 1
linear <- mu + X%*%beta
y <- rloggamma(n=n, mu=linear, sigma=sigma, lambda=lambda)
cens <- rloggamma(n=n, mu=linear+2.5, sigma=sigma, lambda=lambda)
delta <- as.numeric(y <= cens)
y[delta==0] <- cens[delta==0]
w <- rep(1,n)
set.seed(2345)
SNPTest <- robustloggamma:::loggammacenslmrob.S(IX,y,delta,control=robustloggamma:::loggammacenslmrob.S.control())
SNP <- dget(file='Snpcensreg-savedvalues.R')
stopifnot(
all.equal(SNPTest, SNP, tol = 2e-5)
)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.