I ran the following code to create and filter my "data" object:
data <- squareCounts(input, hg19.param, width=bin.size, filter=1)
keep <- aveLogCPM(asDGEList(data)) > 0
data <- data[keep,]
When I tried normalizing the filtered data using the following command I get the error below
data <- normOffsets(data, type="loess", se.out=TRUE)
Error in loessFit (log (mat [, x] + cont.cor.scaled [x]), ab, …) : unused argument (type = “loess” )
I tried running it on the initial data object and get the same error. Any ideas about how I can troubleshoot this error?
On quick glance, it looks like the documentation might be misleading. while the documentation indeed shows
When I examine the code the extra arguments appear to be passed to the
loessFit
function notcalcNormFactors
. So see?loessFit
for acceptable arguments to pass for...
I think if you change
type = "loess"
tomethod="loess"
it should work. It might also be worth reporting to the package maintainers to update their documentation.Thank you! changing type = "loess" to method="loess" worked out really well.