qrcm

Running examples in ‘qrcm-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: iqr
> ### Title: Quantile Regression Coefficients Modeling
> ### Aliases: iqr
> ### Keywords: models regression
> 
> ### ** Examples
> 
> 
>   ##### Using simulated data in all examples
> 
> 
>   ##### Example 1
>   
>   n <- 1000
>   x <- runif(n)
>   y <- rnorm(n, 1 + x, 1 + x)
>   # true quantile function: Q(p | x) = beta0(p) + beta1(p)*x, with 
>     # beta0(p) = beta1(p) = 1 + qnorm(p)
>                               
>   # fit the true model: b(p) = (1 , qnorm(p))
>   m1 <- iqr(y ~ x, formula.p = ~ I(qnorm(p)))
>   # the fitted quantile regression coefficient functions are
>     # beta0(p) = m1$coef[1,1] + m1$coef[1,2]*qnorm(p)
>     # beta1(p) = m1$coef[2,1] + m1$coef[2,2]*qnorm(p)
> 
> 
> 
>   ##### Example 2 ### excluding coefficients
>   
>   n <- 1000
>   x <- runif(n)
>   qy <- function(p,x){(1 + qnorm(p)) + (1 + log(p))*x}
>   # t
...[truncated]...
 runif(n,0,5)
> 	
>   u <- runif(n)
>   beta0 <- -log(1 - u)
>   beta1 <- 0.2*log(1 - u)
>   t <- beta0 + beta1*x  # time variable
>   c <- rexp(n,2)        # censoring variable
>   y <- pmin(t,c)        # observed events
>   d <- (t <= c)         # 1 = event, 0 = censored
>   
>   iqr(Surv(y,d) ~ x, formula.p = ~ I(log(1 - p)))

Call:
iqr(formula = Surv(y, d) ~ x, formula.p = ~I(log(1 - p)))


Number of free parameters: 4 
Number of observations: 1000 
N. of events: 539

Coefficients:
             (Intercept)  I(log(1 - p))
(Intercept)   3.907e-04   -1.064e+00   
x            -4.387e-05    2.128e-01   

>   
>   ##### Example 8 (cont.) ### censored and truncated data
> 
>   z <- rexp(n,10)   # truncation variable
>   w <- which(y > z) # only observe z,y,d,x when y > z
>   z <- z[w]; y <- y[w]; d <- d[w]; x <- x[w]
> 
>   iqr(Surv(z,y,d) ~ x, formula.p = ~ I(log(1 - p)))
Error in hatF$time : $ operator is invalid for atomic vectors
Calls: iqr -> ctiqr.internal -> trans
Execution halted