gamma_fd {SIMICO}

d/d_gamma_l

Description

Calculates the gradient term for U_g for the score statistic.

Usage

gamma_fd(l, HL_array, HR_array, tpos_all, obs_all,
   temp_beta, A_i, no_l_all, gMat, a1, a2, d)

Arguments

l

Index of first outcome of interest.

HL_array

n x k matrix containing all the hazard values for the left times.

HR_array

n x k matrix containing all the hazard values for the right times.

tpos_all

n x k matrix containing a indictor for whether that time is left-censored or not.

obs_all

n x k matrix containing a indictor for whether that time is right-censored or not.

temp_beta

Vector of fitted coefficients.

A_i

Product of apply_diffs across all outcomes k summed over all quadrature nodes d.

no_l_all

n x (K - 1) matrix containing the product of apply_diffs across all outcomes K excluding the current outcome l.

gMat

n x q matrix of genetic information.

a1

First shape parameter of beta parameter.

a2

Second shape parameter of beta parameter.

d

Number of quadrature nodes.

Value

The output is a vector containing the first derivative with respect to gamma.