R/G.theo.multcomp.R In RVAideMemoire: Testing and Plotting Procedures for Biostatistics

Documented in G.theo.multcomp

```G.theo.multcomp <-
function(x,p=rep(1/length(x),length(x)),p.method="fdr") {
if (sum(p)!=1) {stop("sum of probabilities must be 1")}
theo <- integer(length(x))
G <- integer(length(x))
pval <- integer(length(x))
for (i in 1:length(x)) {
test <- G.test(c(x[i],sum(x)-x[i]),p=c(p[i],1-p[i]))
theo[i] <- as.numeric(test\$expected[1])
G[i] <- as.numeric(test\$statistic)
pval[i] <- as.numeric(test\$p.value)
}
p.adj <- p.adjust(pval,method=p.method)
comp <- data.frame("observed"=x,"expected"=theo,"G"=G,"Pr(>G)"=p.adj," "=.psignif(p.adj),
stringsAsFactors=FALSE,check.names=FALSE)
call <- match.call()
dname.x <- if(length(call\$x)==1) {call\$x} else {paste(call\$x[1],"(",paste(call\$x[-1],collapse=","),")",sep="")}
dname.p <- if(length(call\$p)==1) {call\$p} else {paste(call\$p[1],"(",paste(call\$p[-1],collapse=","),")",sep="")}
dname <- paste(dname.x," and ",dname.p,sep="")
result <- list(method="G-tests",data.name=dname,observed=x,expected=theo,p.adjust.method=p.method,statistic=G,p.value2=p.adj,p.value=comp)
class(result) <- "RV.multcomp"
return(result)
}
```

Try the RVAideMemoire package in your browser

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

RVAideMemoire documentation built on May 14, 2018, 5:07 p.m.