lclGWAS (version 1.0.1)

betaEst: Estimate the Effect Parameter for Discrete-Time Multivariate Frailty Model for Grouped Survival Data

Description

A method to estimate the effect parameter for Multivariate Frailty Model

Usage

betaEst(fam_size, alpha, dtime, delta, g, var, m, lower, upper)

Arguments

fam_size

A vector contains the number of samples for each family.

alpha

Baseline hazard vector for each time intervals.

dtime

A vector of observed survival time for each sample.

delta

A event indicator vector: 1 means event and 0 means otherwise.

g

A vector for numerical gene type of each sample.

var

A scalar for frailty variance.

m

A integer for the number of families.

lower

A scalar for the lower bound for the effect parameter estimation.

upper

A scalar for the upper bound for the effect parameter estimation.

Value

A list with one object

Examples

# Generate dummy data	
fam_size <- rep(3,3)
alpha <- c(0.7500000, 0.6666667, 0.5000000, 0.0000000)
dtime <- c(1, 3, 3, 2, 1, 1, 2, 3, 1)
delta <- c(1, 0, 1, 1, 1, 0, 1, 0, 1)
g     <- c(0, 1, 1, 1, 2, 2, 0, 0, 0)
var   <- 0.2
m     <- 3
res   <- betaEst(fam_size, alpha, dtime, delta, g, var, m, lower = 0, upper = 2) 
res