Dear all,
I have been trying to apply the Cy0 algorithm of the package qpcR by
creating an object "obj" with the normalized fluorescent data from a
384 plate whose characteristics were TaqMan chemistry and 45 cycles.
The import of the object was successful but when I implemented the
pcrfit model (indicating in column 1 the number of cycles and in the
successive columns the actual data) I obtained an error in
model.frame.default. Yet the dataframe is composed by 384 columns (+ 1
for the cycles) and each by 45 rows (+ 1 for the column titles).
Would you have some tips on how to debug this problem?
Many thanks,
Best regards,
Luigi
>>>
Here is a sample of the code I have written (the result database is
attached for further reference)
> obj<-pcrimport2(
+ file="cq.data.txt",
+ sep="\t",
+ dec=".",
+ header=TRUE,
+ colClasses="numeric",
+ quote=""
+ )
# j<-2:385 to be implemented with a loop cycle, to read obj columns
# i<-1:384 to be implemented with a loop cycle, to write the results
# Cy is an object with 384 values> model<-pcrfit(obj, cyc=1, j, model=l4, do.optim=TRUE, robust=FALSE)
> Cy[i]<-Cy0(model, plot=FALSE)
[1] Error in model.frame.default(formula = ~Fluo + Cycles, data DATA, weights =
WEIGHTS, :
variable lengths differ (found for '(do.optim)')