# Example 1:
require(graphics)
Dilution <- 2 + log(ToothGrowth["dose"]) / log(2)
names(Dilution) <- "Dilution"
Replicate <- rep(1:10, 6)
Data <- cbind(ToothGrowth, Replicate, Dilution)
dimnames(Data)[[2]] <- c("Response", "Sample", "Dose", "Replicate", "Dilution")
Design <- "crd"
Frame <- data2assayFrame(Data, dr = 2, design = Design)
Fits <- pla.fit(Frame, sampleLabels = c("VC", "OJ"), dr = 2,
design = Design, main = "ToothGrowth", show = TRUE)
pla.plots(Frame, design = Design,
sampleLabels = c("VC", "OJ"), main = "ToothGrowth")Run the code above in your browser using DataLab