ljr01: Perform test of 0 vs 1 joinpoints.

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/ljr01.R

Description

This function tests the null hypothesis of 0 joinpoints versus the alternative of one joinpoint based on the likelihood ratio test statistic. The p-value is determined by a Monte Carlo method.

Usage

1
ljr01(y,n,tm,X,ofst,R=1000,alpha=.05)

Arguments

y

the vector of Binomial responses.

n

the vector of sizes for the Binomial random variables.

tm

the vector of ordered observation times.

X

a design matrix containing other covariates.

ofst

a vector of known offsets for the logit of the response.

R

number of Monte Carlo simulations.

alpha

significance level of the test.

Details

The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.

Value

pval

The estimate of the p-value via simulation.

Coef

A table of coefficient estimates.

Joinpoint

The estimates of the joinpoint, if it is significant.

wlik

The maximum value of the re-weighted log-likelihood.

Author(s)

The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].

References

Czajkowski, M., Gill, R. and Rempala, G. (2008). Model selection in logistic joinpoint regression with applications to analyzing cohort mortality patterns. Statistics in Medicine 27, 1508-1526.

See Also

ljr0,ljr1

Examples

1
2
3
4
 data(kcm)
 attach(kcm)
 set.seed(12345)
## Not run: ljr01(Count,Population,Year+.5,R=20)

ljr documentation built on May 29, 2017, 10:22 a.m.