Description Usage Arguments Value Author(s) References Examples

This is the sencond version of the globaltest, the non-standardized globaltest

1 2 | ```
## a powerful variant of globaltest
gt2 (y, X, hyps, alpha = 0.05)
``` |

`y` |
The response vector (numeric vector). |

`X` |
The full design matrix, whose columns are named by the covariates. |

`hyps` |
The name vector of the covariates in the pathway of interest (character vector). |

`alpha` |
The type I error rate allowed. The default is 0.05. |

Returns the p-value, the observed and expected test statistics and the number of covariates.

Ningning Xu

Maintainer: Ningning Xu <n.xu@lumc.nl; xu15263142750@gmail.com>

Ningning Xu, Aldo solari, Jelle Goeman, Clsoed testing with global test, with applications on metabolomics data, arXiv:2001.01541, https://arxiv.org/abs/2001.01541

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | ```
#Generate the design matrix and response vector for logistic regression models
n= 100
m = 5
X = matrix(data = 0, nrow = n, ncol = m,byrow = TRUE )
for ( i in 1:n){
set.seed(1234+i)
X[i,] = as.vector(arima.sim(model = list(order = c(1, 0, 0), ar = 0.2), n = m) )
}
y = rbinom(n,1,0.6)
X[which(y==1),1:3] = X[which(y==1),1:3] + 0.8
xs = paste("x",seq(1,m,1),sep="")
colnames(X) = xs
hyps=xs[1]
#The raw p-values of globaltest
gt2(y = y, X = X, hyps = hyps, alpha = 0.05)
#p-value Statistic Expected #Cov
#7.64e-03 2.30e+02 1.24e+02 1.00e+00
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.