R/pvalFmla.R

Defines functions pvalFmla

Documented in pvalFmla

pvalFmla <-
function(data, i, covariable=NULL, family=binomial) {
		if(length(unique(data[,i])) == 1) { 
			pval <- 1
      	} else {
			form=CreateFormula(data,i,covariable)
			if(length(covariable)==0) { 
				fit <- glm(formula=form, data=data, family=family)
				pval <- anova(fit, test="LRT")[2,5]
				return(pval)
			}
			fit <- glm(formula=form, data=data.frame(data, covariable), family=family)
		  	pval <- anova(fit, test="LRT")[2,5]
		}
		return(pval)
}

Try the globalGSA package in your browser

Any scripts or data that you put into this service are public.

globalGSA documentation built on May 1, 2019, 9:19 p.m.