inst/snippet/smoke-manual.R

row.sum <- apply(smokeTab,1,sum)
col.sum <- apply(smokeTab,2,sum)
grandTotal <- sum(smokeTab)
e <- outer(row.sum,col.sum)/grandTotal; e
o <- smokeTab
stat <- sum ( (e-o)^2/e); stat
pval <- 1 - pchisq(stat,df=2); pval

Try the fastR package in your browser

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

fastR documentation built on May 1, 2019, 11:18 p.m.