meffil.qc.summary Error and Warning about meffil HOT 4 OPEN

Sorry to hear about this error. Could you run the following code and send me the output generated?

    outlier.sd <- 3
    dat <- data.frame(
        sample.name = sapply(qc.objects, function(x) x$sample.name),
        methylated = sapply(qc.objects, function(x) x$median.m.signal),
        unmethylated = sapply(qc.objects, function(x) x$median.u.signal),
        stringsAsFactors=F
    )

    fit <- lm(methylated ~ unmethylated, dat)
    dat$resids <- residuals(fit)
    dat$methylated.lm <- predict(fit)
    interval.size <- outlier.sd*sd(dat$resids)
  
    dat$upper.lm <- dat$methylated.lm + interval.size
    dat$lower.lm <- dat$methylated.lm - interval.size

    dat$outliers <- (dat$resids > mean(dat$resids) + interval.size
                     | dat$resids < mean(dat$resids) - interval.size)
   
   na.idx <- which(is.na(dat$outliers))
   cat(nrow(dat), length(na.idx), "\n")
  if (na.idx > 3) na.idx <- na.idx[1:3]
   print(dat[na.idx,])

from meffil.

Hi,
thanks for your help.

Here is the output

`cat(nrow(dat), length(na.idx), "\n")
1 1

if (na.idx > 3) na.idx <- na.idx[1:3]
print(dat[na.idx,])
sample.name methylated unmethylated resids
206891110004_R01C01 206891110004_R01C01 2621.36 2521.409 0
methylated.lm upper.lm lower.lm outliers
206891110004_R01C01 2621.36 NA NA NA`

Thanks in advances,

Rossella

from meffil.

Ah I see. The qc.summary() is failing when it attempts to identify outliers because only one sample is being analysed. Although some quality information can be generated for a single sample, variation between samples is necessary to evaluate technical variation. Do you have more samples? Will you regularly need a QC procedure to apply to individual samples?

from meffil.

Thank you so much for your support. That was the issue.
To replay to your question no I will run QC onbatches of samples, and the tested with multiples samples worked perfectly. Thanks again

from meffil.