Hello

I tried use grofit package in our data set. We provide a subset of our
data with X iso, and 4 doses, and insect died was count each day for
long 5 days. We started with Y insects per dishes. When one is dead, it
was counted and removed. Died insect is cumulative in the next days.
i.e. day 1 died 1. day 2 no died, so, day 2 is assigned 1 died (from day
1).

Here is the script:

library(lattice)
library(grofit)
library(repmis)

url.csv <- 
https://dl.dropboxusercontent.com/u/34009642/cepajabo07_wide_acumulado.csv

data02 <- read.table(url.csv, header=TRUE, sep="\t", dec=",")

head(data02)

timepoints <- 1:5 # 5 days
time <- t(matrix(rep(timepoints, 120), c(5, 120))) # 5 days and 120 experiments
                                                   # (6 iso * 4 doses
                                                   # * 5 rep)
time

TestRun1$drFit
TestRun2$drFit

colData <- c("black", "cyan", "magenta", "blue")

plot(TestRun1$gcFit, opt = "s", colData = colData, colSpline = 1, 
     pch = 1:4, cex = 1)

plot(TestRun2$gcFit, opt = "s", colData = colData, colSpline = 1, 
     pch = 1:4, cex = 1)

plot(TestRun1$drFit$drFittedSplines[[1]], colData = colData, 
     pch = 1:4, cex = 1)

plot(TestRun2$drFit$drFittedSplines[[1]], colData = colData, 
     pch = 1:4, cex = 1)

The problem: grofit didn't deal with replicates and do a curve for each
ones.

Is it a way to get response curve with the replicates?

We are interested in LD50, and dose response curve, and graphs.

Any suggestion is very welcome!

Thank you!

-- 
Marcelo

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to