Title: | Logistic Joinpoint Regression |
---|---|
Description: | Fits and tests logistic joinpoint models. |
Authors: | Michal Czajkowski, Ryan Gill, Greg Rempala |
Maintainer: | Ryan Gill <[email protected]> |
License: | GPL (>= 2) |
Version: | 1.4-0 |
Built: | 2024-12-08 06:58:09 UTC |
Source: | CRAN |
This table gives the yearly mortality counts due to neoplasms (ICD 10 codes C00-D48) and population sizes for Kentucky from 1999-2005. For more information, see http://wonder.cdc.gov/wonder/help/cmf.html.
data(kcm)
data(kcm)
A 7 by 3 data frame.
Centers for Disease Control and Prevention, National Center for Health Statistics. Compressed Mortality File 1999-2005. CDC WONDER On-line Database, compiled from Compressed Mortality File 1999-2005 Series 20 No. 2K, 2008. Accessed at http://wonder.cdc.gov/cmf-icd10.html on May 5, 2008.
Determines the maximum likelihood estimate of model coefficients in the logistic joinpoint regression model with no joinpoints.
ljr0(y,n,tm,X,ofst)
ljr0(y,n,tm,X,ofst)
y |
the vector of Binomial responses. |
n |
the vector of sizes for the Binomial random variables. |
tm |
the vector of observation times. |
X |
a design matrix containing other covariates. |
ofst |
a vector of known offsets for the logit of the response. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
Coef |
A table of coefficient estimates. |
wlik |
The maximum value of the re-weighted log-likelihood. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) ljr0(Count,Population,Year+.5)
data(kcm) attach(kcm) ljr0(Count,Population,Year+.5)
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.
ljr01(y,n,tm,X,ofst,R=1000,alpha=.05)
ljr01(y,n,tm,X,ofst,R=1000,alpha=.05)
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
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. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljr01(Count,Population,Year+.5,R=20)
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljr01(Count,Population,Year+.5,R=20)
Determines the maximum likelihood estimates of model coefficients in the logistic joinpoint regression model with one joinpoint.
ljr1(y,n,tm,X,ofst,summ=TRUE)
ljr1(y,n,tm,X,ofst,summ=TRUE)
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. |
summ |
a boolean indicator of whether summary tables should be returned. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
Coef |
A table of coefficient estimates. |
Joinpoint |
The estimate of the joinpoint. |
wlik |
The maximum value of the re-weighted log-likelihood. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) ljr1(Count,Population,Year+.5)
data(kcm) attach(kcm) ljr1(Count,Population,Year+.5)
This function performs the likelihood ratio tests to find p-values in testing the significance of each of the coefficients as well as the intercept and ordered observation times. The p-values are determined by a Monte Carlo method.
ljr11(y,n,tm,X,ofst,R=1000)
ljr11(y,n,tm,X,ofst,R=1000)
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
pvals |
The estimates of the p-values via simulation. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljr11(Count,Population,Year+.5,R=20)
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljr11(Count,Population,Year+.5,R=20)
This function performs the backward joinpoint selection algorithm with K maximum possible number of joinpoints based on the likelihood ratio test statistic. The p-value is determined by a Monte Carlo method.
ljrb(K,y,n,tm,X,ofst,R=1000,alpha=.05)
ljrb(K,y,n,tm,X,ofst,R=1000,alpha=.05)
K |
the pre-specified maximum possible number of joinpoints |
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
pvals |
The estimates of the p-values via simulation. |
Coef |
A table of coefficient estimates. |
Joinpoints |
The estimates of the joinpoint, if it is significant. |
wlik |
The maximum value of the re-weighted log-likelihood. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrb(1,Count,Population,Year+.5,R=20)
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrb(1,Count,Population,Year+.5,R=20)
This function performs the full forward joinpoint selection algorithm based on the likelihood ratio test statistic. The p-value is determined by a Monte Carlo method.
ljrf(y,n,tm,X,ofst,R=1000,alpha=.05)
ljrf(y,n,tm,X,ofst,R=1000,alpha=.05)
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
pvals |
The estimates of the p-values via simulation. |
Coef |
A table of coefficient estimates. |
Joinpoints |
The estimates of the joinpoint, if it is significant. |
wlik |
The maximum value of the re-weighted log-likelihood. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrf(Count,Population,Year+.5,R=20)
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrf(Count,Population,Year+.5,R=20)
This function tests the null hypothesis of j joinpoint(s) versus the alternative of k joinpoint(s) based on the likelihood ratio test statistic. The p-value is determined by a Monte Carlo method.
ljrjk(j,k,y,n,tm,X,ofst,R=1000,alpha=.05)
ljrjk(j,k,y,n,tm,X,ofst,R=1000,alpha=.05)
j , k
|
pre-specified number of joinpoints in the null and alternative hpyotheses (the smaller is used for the null). |
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
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. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrjk(0,1,Count,Population,Year+.5,R=20)
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrjk(0,1,Count,Population,Year+.5,R=20)
Determines the maximum likelihood estimates of model coefficients in the logistic joinpoint regression model with k joinpoints.
ljrk(k,y,n,tm,X,ofst)
ljrk(k,y,n,tm,X,ofst)
k |
the pre-specified number of joinpoints (with unknown locations). |
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
Coef |
A table of coefficient estimates. |
Joinpoints |
The estimates of the joinpoints. |
wlik |
The maximum value of the re-weighted log-likelihood. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) ljrk(1,Count,Population,Year+.5)
data(kcm) attach(kcm) ljrk(1,Count,Population,Year+.5)
This function performs the likelihood ratio tests to find p-values in testing the significance of each of the coefficients as well as the intercept and ordered observation times. The p-values are determined by a Monte Carlo method.
ljrkk(k,y,n,tm,X,ofst,R=1000)
ljrkk(k,y,n,tm,X,ofst,R=1000)
k |
the pre-specified number of joinpoints (with unknown locations). |
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. |
The re-weighted log-likelihood is the log-likelihood divided by the largest component of n.
pvals |
The estimates of the p-values via simulation. |
The authors are Michal Czajkowski, Ryan Gill, and Greg Rempala. The software is maintained by Ryan Gill [email protected].
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.
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrkk(1,Count,Population,Year+.5,R=20)
data(kcm) attach(kcm) set.seed(12345) ## Not run: ljrkk(1,Count,Population,Year+.5,R=20)