Dear all,
I am trying to implement a method to read qPCR runs using the qpcR package developed by Andrej-Nikolai Spiess. I am almost there, but when I try to apply the cy0() function to my actual data I obtain the error:
Error in eff(object, ...) : could not find function "FCT"
I can create a vector with the results but the error comes when I try to make the plot of the same data. I used the function modlist() to create the model, using the argument remove = "none"; when using the argument remove ="fit" the error becomes:
Error in plot(ml[[1]]) :
error in evaluating the argument 'x' in selecting a method for function 'plot': Error in ml[[1]] : subscript out of bounds
Nonetheless when I try to extract the data out of the context and apply it to another dataframe, the things go smoothly and I can't provide a smaller dataframe to give an example.
I therefore have attached the original file itself for further reference.
the code I have written is:
# create objects
M
Hey Luigi,
you need to make declare ml, otherwise it gets written over with each loop.
wells = 2:M
ml = lapply(wells,function(j){
modlist(obj, cyc = 1, fluo=j, model = l4, norm = TRUE, remove = "none")
})
names(ml) = 2:M
##to plot
Cy0(ml[[1]][[1]],plot=TRUE)
Cy0(ml[[2]][[1]],plot=TRUE)
this should work.
Dear Gen Lin,
thank you for your reply and sorry for the delay in answering. I actually changed approach, since modlist creates a loop on its own. I therefore simplified the script with
mdl
- Badges
- Science topic