title: Unable to run ComBat script from R's sva library
body: <p>I am trying to run ComBat script on a dataset with 2 batches, but I am getting errors and I do not know how to inspect code since I am an R newbie.</p>
  
  <p>I am running ComBat method in this way:</p>
  
  <pre><code># Load sva
  library(sva)
  # Read expression values
  dat = read.table('dataset.xls', header=TRUE, sep='\t')
  # Read sample information file about batches
  sif = read.delim('sif.tsv', header=TRUE, sep='\t')
  # Call ComBat
  ComBat(dat=dat,batch=sif$Batch, mod=NULL)
  </code></pre>
  
  <p>Anyway my output is:</p>
  
  <pre><code>Found 2 batches
  Found 0  categorical covariate(s)
  Found 54675 Missing Data Values
  Standardizing Data across genes
  Error in solve(t(des) %*% des) %*% t(des) %*% y1 : 
    requires numeric/complex matrix/vector arguments
  </code></pre>
  
  <p>Data format for dat is:</p>
  
  <pre><code>probe set    &lt;sample1&gt;   ...    &lt;sampleN&gt;
  &lt;gene_name&gt;  &lt;value1&gt;    ...    &lt;valueN&gt;
  ...
  </code></pre>
  
  <p>Data format for sif is:</p>
  
  <pre><code>Array name   Sample name   Batch
  &lt;Array1&gt;     &lt;Sample1&gt;     &lt;Batch1&gt;
  ...
  </code></pre>
  
  <p>Any hint is appreciated. I'll provide more info if needed.</p>
  
  <p>Thanks</p>