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)') ______________________________________________ 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.
