Nothing
## ----install_package, eval=FALSE-----------------------------------------
# install.packages("LW1949")
## ----load_package--------------------------------------------------------
library(LW1949)
## ----dose-effect_data----------------------------------------------------
conc <- c(0.0625, 0.125, 0.25, 0.5, 1, 2, 3)
numtested <- rep(8, 7)
numaffected <- c(1, 4, 4, 7, 8, 8, 8)
mydat <- dataprep(dose=conc, ntot=numtested, nfx=numaffected)
## ----dataprep()_output---------------------------------------------------
mydat
## ----fit_LW--------------------------------------------------------------
intslope <- fitLWauto(mydat)
fLW <- LWestimate(intslope, mydat)
## ----fitLWauto_output----------------------------------------------------
intslope
## ----LWestimate_output---------------------------------------------------
fLW
## ----estimate_LW---------------------------------------------------------
pctaffected <- c(25, 50, 99.9)
predlinear(pctaffected, fLW)
## ----plot_fits, fig.width=5, fig.height=5--------------------------------
plotDELP(mydat)
predLinesLP(fLW)
plotDE(mydat)
predLines(fLW)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.