give.c: Work Function for 'smoothSurvReg'

View source: R/convertCDA.R

give.cR Documentation

Work Function for 'smoothSurvReg'

Description

Give a vector of all 'c' G-spline coefficients satisfying the three constrains (see below) if only (g-3) c's are given.

Usage

give.c(knots, sdspline, last.three, c.rest)

Arguments

knots

A vector of G-spline knots mu.

sdspline

Standard deviation sigma0 of the basis G-spline.

last.three

Indeces of the three 'c' G-spline coefficients which are not included in c.rest. This must be a vector of length 3 with three different numbers from 1:length(knots).

c.rest

A vector of length(knots) - 3 'c' spline coefficients corresponding to knots knots[-last.three]. They must all lie between zero and one and their sum must be lower or equal to one.

Details

The functions computes remaining three 'c' G-spline coefficients (if it is possible) such that the resulting G-spline satisfies the following:

c[1] + ... c[g] = 1,

c[1]mu[1] + ... c[g]mu[g] = 0,

c[1]mu[1]^2 + ... c[g]mu[g]^2 = 1 - sigma0^2.

Value

A vector of all length(knots) 'c' G-spline coefficients.

Note

WARNING: This function is primarily used inside smoothSurvReg. Consequently, it has very few error checks on its input arguments.

Author(s)

Arnošt Komárek arnost.komarek@mff.cuni.cz

Examples

knots <- seq(-4, 4, 0.5)
sd0 <- 0.3
ccoef <- find.c(knots, sd0, dist = "dstlogis")

last.three <- c(3, 7, 10)
c.rest <- ccoef[-last.three]
ccoef2 <- give.c(knots, sd0, last.three, c.rest)

print(ccoef)
print(ccoef2)    ## Almost no change

smoothSurv documentation built on Oct. 11, 2022, 1:05 a.m.